R dpnorm 穩定評估正常 c.d.f.s 之間的差異

說明

評估兩個 N(0,1) 累積分布函數之間的差異，避免取消錯誤。

用法

dpnorm(x0,x1)

參數

x0	用於評估標準正態分布函數的較低值向量。
x1	用於評估標準正態分布函數的上限值向量。

細節

相當於 pnorm(x1)-pnorm(x0) ，但當 x0 和 x1 值非常接近或位於標準法線的上尾時穩定。

例子

require(mgcv)
x <- seq(-10,10,length=10000)
eps <- 1e-10
y0 <- pnorm(x+eps)-pnorm(x) ## cancellation prone
y1 <- dpnorm(x,x+eps)       ## stable
## illustrate stable computation in black, and
## cancellation prone in red...
par(mfrow=c(1,2),mar=c(4,4,1,1))
plot(log(y1),log(y0),type="l")
lines(log(y1[x>0]),log(y0[x>0]),col=2)
plot(x,log(y1),type="l")
lines(x,log(y0),col=2)

作者

Simon N. Wood simon.wood@r-project.org