我试图计算正态分布密度的积分,期望值为200,标准差为20.从-Inf到Inf,这应该是1.
我得到以下内容:
> integrate(dnorm, mean=200, sd=20,-Inf, Inf)$value
[1] 1.429508e-08
对于低于169的预期值,我得到正确的值,1。 如何获得更大预期值的正确值?
答案 0 :(得分:4)
可选地
21474839
见here。
要了解发生了什么,请注意以下细分数量:
integrate(dnorm, mean=200, sd=20, lower= -Inf, upper= Inf, abs.tol = 0)$value
[1] 1
答案 1 :(得分:3)
设置间隔有限似乎有帮助
integrate(dnorm, mean=200, sd=20, -1e4, 1e4)
# 1 with absolute error < 2e-07