泊松分布（dpois）积分R.

Question

我正在尝试将泊松分布（dpois）整合到R中但是我得到了一个错误的答案（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累积分布函数。