二项式和泊松随机变量的总和


10

如果我们有两个独立的随机变量X 2P ø 小号λ ,什么是的概率质量函数X 1 + X 2X1Binom(n,p)X2Pois(λ)X1+X2

注意:这不是我的作业。


我猜你试过回旋了吗?en.wikipedia.org/wiki / ...您在哪里被卡住?我认为是没有封闭的形式,否则该解决方案将可能是在这里:en.wikipedia.org/wiki/...
斯蒂芬Kolassa

3
是的,这就是我尝试过的方法,但是也许我在这里找到了答案:mathstatica.com/SumBinomialPoisson Kummer融合超几何函数..休
Matteo Fasiolo 2012年

1
我已经根据该站点上的用途读取了作业标签。干杯。:-)
主教

2
新颖是指新的(以前未知或未出版)。我也不同意使用已知的方法解决新问题会使其成为家庭作业-对于大多数在发行版上发布结果的期刊文章也可以这样说。
wolfies

2
与统计中的其他许多情况一样,超几何函数与积分参数一起出现,如果愿意,您可以将其理解为卷积中隐式(有限)和的简写形式。这种表达式的优点在于,有无数种方法可以将其操纵为更简单的形式,并且常常可以在不实际进行求和的情况下对其进行评估。
whuber

Answers:


7

pX1+X2(k)0k<nkni=0npX1(i)zkj=0pX2(j)zjpX1+X2(k)zk


1

P(X1+X2=k)=x1=0min(n,k)(nx1)px1(1p)nx1eλλkx1(kx1)!=(1p)neλλkΓ(k+1)2F0(k,n; ;p(p1)λ)


0

Dilip Sarwate 7年前曾表示不可能简化,尽管在评论中对此提出了挑战。但是,我认为有必要指出,即使不做任何简化,任何电子表格或编程语言中的计算也非常简单。

这是R中的实现:

# example parameters
n <- 10
p <- .3
lambda <- 5

# probability for just a single value
x <- 10  # example value
sum(dbinom(0:x, n, p) * dpois(x:0, lambda))

# probability function for all values
x0  <- 0:30   # 0 to the maximum value of interest
x   <- outer(x0, x0, "+")
db  <- dbinom(x0, n, p)
dp  <- dpois(x0, lambda)
dbp <- outer(db, dp)
aggregate(as.vector(dbp), by=list(as.vector(x)), sum)[1:(max(x0)+1),]

1
迪利普(Dilip)没有表明不可能简化总和:他说了这样的断言(而且断言似乎是不正确的)。如果您遵循OP提供的链接,则会提供有关Kummer融合超几何功能的解决方案。
Wolfies '19

@wolfies-在这个老问题的新答案中,这将是非常有趣的一点。可能比我的有趣。
Pere

1
对于二项式中较大的n和较大的lambda而言,一种可能更快的方法将涉及快速傅立叶变换(或类似的变换)。我已经成功地将其用于许多实际问题,这些问题卷积在代数上并不方便,但是数值答案就足够了,并且添加了多个独立变量。
Glen_b-恢复莫妮卡

1
nλdpoisxxx<-qpois(0:1+c(1,-1)*1e-6, lambda)dpoisxzapsmalln

确实。我对自己的应用程序执行了类似的操作-进行足够远的操作可以准确地提供所需的分位数。
Glen_b-恢复莫妮卡
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.