मुझे इस वितरण का नाम नहीं पता, लेकिन आप इसे कुल संभावना के कानून से निकाल सकते हैं। मान लीजिए कि X,Y प्रत्येक में क्रमशः मापदंडों (r1,p1) और साथ नकारात्मक द्विपद वितरण (r2,p2)है। मैं पैरामीटर का उपयोग कर रहा हूं जहां क्रमशः आर 1 'वें, और आर 2 ' वें विफलताओं X,Yसे पहले की सफलताओं की संख्या का प्रतिनिधित्व करते हैं। फिर,r1r2
P(X−Y=k)=EY(P(X−Y=k))=EY(P(X=k+Y))=∑y=0∞P(Y=y)P(X=k+y)
हम जानते है
P(X=k+y)=(k+y+r1−1k+y)(1−p1)r1pk+y1
तथा
P(Y=y)=(y+r2−1y)(1−p2)r2py2
इसलिए
P(X−Y=k)=∑y=0∞(y+r2−1y)(1−p2)r2py2⋅(k+y+r1−1k+y)(1−p1)r1pk+y1
यह सुंदर नहीं है (yikes!)। एकमात्र सरलीकरण जो मैं देख रहा हूं वह सही है
pk1(1−p1)r1(1−p2)r2∑y=0∞(p1p2)y(y+r2−1y)(k+y+r1−1k+y)
जो अभी भी बहुत बदसूरत है। मुझे यकीन नहीं है कि यह मददगार है, लेकिन यह भी फिर से लिखा जा सकता है
pk1(1−p1)r1(1−p2)r2(r1−1)!(r2−1)!∑y=0∞(p1p2)y(y+r2−1)!(k+y+r1−1)!y!(k+y)!
p
मैंने सिमुलेशन के साथ सत्यापित किया कि उपरोक्त गणना सही है। इस सामूहिक समारोह की गणना करने और कुछ सिमुलेशन करने के लिए यहां एक कच्चा आर फ़ंक्शन है
f = function(k,r1,r2,p1,p2,UB)
{
S=0
const = (p1^k) * ((1-p1)^r1) * ((1-p2)^r2)
const = const/( factorial(r1-1) * factorial(r2-1) )
for(y in 0:UB)
{
iy = ((p1*p2)^y) * factorial(y+r2-1)*factorial(k+y+r1-1)
iy = iy/( factorial(y)*factorial(y+k) )
S = S + iy
}
return(S*const)
}
### Sims
r1 = 6; r2 = 4;
p1 = .7; p2 = .53;
X = rnbinom(1e5,r1,p1)
Y = rnbinom(1e5,r2,p2)
mean( (X-Y) == 2 )
[1] 0.08508
f(2,r1,r2,1-p1,1-p2,20)
[1] 0.08509068
mean( (X-Y) == 1 )
[1] 0.11581
f(1,r1,r2,1-p1,1-p2,20)
[1] 0.1162279
mean( (X-Y) == 0 )
[1] 0.13888
f(0,r1,r2,1-p1,1-p2,20)
[1] 0.1363209
rp1,p21−p1,1−p2