与R中的pnorm()相比,atan()是否具有计算优势?

出于好奇想看它

首先定义函数

PNORM <- function(x) { 1/(exp(-358/23*x + 111*atan(37*x/294)) + 1) }

x <- seq(-4, 4, .01)
plot(x, pnorm(x)-PNORM(x), type="l", lwd=3, ylab="Difference")

这导致了这张图

因此,差异很小,但在某些应用中可能不足以忽略。YMMV。如果我们看一下计算时间,那么它们大致等于近似值,似乎略快一些

> microbenchmark::microbenchmark(pnorm(x), PNORM(x))
Unit: microseconds
     expr    min      lq     mean  median      uq    max neval cld
 pnorm(x) 34.703 34.8785 36.54254 35.1820 38.3150 47.786   100   b
 PNORM(x) 24.293 24.4625 27.07660 24.8875 28.9035 59.216   100  a