泊松分布（dpois）積分R

Question

我正在嘗試在R中集成泊松分布（dpois），但得到的答案不正確（0，絕對錯誤0）和21條警告。 我不明白R如何消化我的便餐，以及為什么它會吐出21條警告。

    dpoisd1 <- function(x) {dpois(x, 0.0001)}
    dpoisd1(1:20)
    integrate(dpoisd1, lower = 1, upper = 20)

它產生的絕對誤差<0的0和一些21警告。 如果有人可以告訴我我的錯誤，我將不勝感激。

Answer 1

使用warnings來查看警告：

warnings()
#Warning messages:
#1: In dpois(x, 1e-04) : non-integer x = 10.500000
#<snip>

dpois的第一個參數必須是非負整數（請參閱help("dpois") ）。 integrate將非整數值傳遞給它。 實際上，不清楚要計算什么。 您正在嘗試集成離散密度函數。 可能需要ppois （累積分布函數）。