X,Y(r1,p1)(r2,p2)X,Yr1r2
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
এটি সুন্দর নয় (হাই!)। আমি এখনই দেখছি কেবল সরলকরণ
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
আমি সিমুলেশন দিয়ে যাচাই করেছি যে উপরের গণনাটি সঠিক is এই ভর ফাংশনটি গণনা করতে এবং কয়েকটি অনুকরণ চালানোর জন্য এখানে একটি অশোধিত আর ফাংশন
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