内容来源
保险风险与破产(原书第二版)科学出版社
定义
如果一个计数分布的分布律满足
p n = ( a + b n ) p n − 1 , n = 1 , 2 , ⋯ p_n=\left(a+\frac{b}{n}\right)p_{n-1},n=1,2,\cdots pn=(a+nb)pn−1,n=1,2,⋯
其中 a , b a,b a,b 均为常数, p 0 > 0 p_0>0 p0>0
则其属于 ( a , b , 0 ) (a,b,0) (a,b,0) 类,其中 0 0 0 表示递推计算的起点
在 ( a , b , 0 ) (a,b,0) (a,b,0) 类中有三种非平凡分布,分别是泊松分布、二项分布和负二项分布
通过 a , b a,b a,b 的取值判别分布
不难看出 a + b ⩾ 0 a+b\geqslant0 a+b⩾0,否则 p n p_n pn 可能取负值
a + b = 0 a+b=0 a+b=0
此时 p n = 0 p_n=0 pn=0,又 ∑ i = 0 ∞ p i = 1 \sum^\infty_{i=0} p_i=1 ∑i=0∞pi=1,所以 p 0 = 1 p_0=1 p0=1
从而分布在 0 0 0 处退化
a = 0 a=0 a=0
此时 p n = ( b / n ) p n − 1 p_n=(b/n)p_{n-1} pn=(b/n)pn−1
因为 ∑ i = 0 ∞ p i = 1 \sum^\infty_{i=0} p_i=1 ∑i=0∞pi=1,而且
∑ i = 0 ∞ p i = p 0 ∑ i = 0 ∞ b n n ! = p 0 e b \sum^\infty_{i=0} p_i=p_0\sum^\infty_{i=0}\frac{b^n}{n!}=p_0e^b i=0∑∞pi=p0i=0∑∞n!bn=p0eb
所以 p 0 = e − b p_0=e^{-b} p0=e−b
所以,当 a = 0 a=0 a=0 时,所得分布是参数为 b b b 的泊松分布
a > 0 , a ≠ − b a>0,a\ne-b a>0,a=−b
通过重复利用递推公式得
p n = ( a + b n ) ⋯ ( a + b 1 ) p 0 = ( a n + b ) ⋯ ( a ⋅ 1 + b ) n ! p 0 = a n n ! ( n + b a ) ⋯ ( 1 + b a ) p 0 \begin{align*} p_n&=\left(a+\frac{b}{n}\right)\cdots\left(a+\frac{b}{1}\right)p_0\\ &=\frac{(an+b)\cdots(a\cdot1+b)}{n!}p_0\\ &=\frac{a^n}{n!}\left(n+\frac{b}{a}\right)\cdots \left(1+\frac{b}{a}\right)p_0\\ \end{align*} pn=(a+nb)⋯(a+1b)p0=n!(an+b)⋯(a⋅1+b)p0=n!an(n+ab)⋯(1+ab)p0
将 1 + b / a 1+b/a 1+b/a 记为 α \alpha α,则
p n = a n n ! ( n − 1 + α ) ⋯ ( 1 + α ) α p 0 = ( α + n − 1 n ) a n p 0 \begin{align*} p_n&=\frac{a^n}{n!}\left(n-1+\alpha\right)\cdots \left(1+\alpha\right)\alpha p_0\\ &=\left(\begin{matrix}\alpha+n-1\\n\end{matrix}\right)a^np_0\\ \end{align*} pn=n!an(n−1+α)⋯(1+α)αp0=(α+n−1n)anp0
因为级数 ∑ i = 0 ∞ p i = 1 \sum^\infty_{i=0} p_i=1 ∑i=0∞pi=1 绝对收敛
由比式判别法知
lim n → ∞ ∣ p n p n − 1 ∣ = lim n → ∞ ∣ a 2 n + a − a n ∣ < 1 \lim_{n\rightarrow\infty}\left|\frac{p_n}{p_{n-1}}\right|= \lim_{n\rightarrow\infty}\left|\frac{a^2}{n}+a-\frac{a}{n}\right|<1 n→∞lim pn−1pn =n→∞lim na2+a−na <1
所以 ∣ a ∣ < 1 |a|<1 ∣a∣<1,结合 a > 0 a>0 a>0,条件退化为 0 < a < 1 0<a<1 0<a<1
又因为
p 0 + p 0 ∑ n = 1 ∞ ( α + n − 1 n ) a n = 1 p_0+p_0\sum^\infty_{n=1}\left(\begin{matrix}\alpha+n-1\\n\end{matrix}\right)a^n=1 p0+p0n=1∑∞(α+n−1n)an=1
对比负二项分布的密度函数,得 p 0 = ( 1 − a ) α p_0=(1-a)^\alpha p0=(1−a)α
此时分布为 N B ( α , 1 − a ) NB(\alpha,1-a) NB(α,1−a)
a + b > 0 , a < 0 a+b>0,a<0 a+b>0,a<0
一定存在某一正整数 c c c 满足
a + b c ⩾ 0 , a + b c + 1 < 0 a+\frac{b}{c}\geqslant0,a+\frac{b}{c+1}<0 a+cb⩾0,a+c+1b<0
取 d = − ( 1 + b / a ) d=-(1+b/a) d=−(1+b/a),结合
p n = a n n ! ( n + b a ) ⋯ ( 1 + b a ) p 0 p_n=\frac{a^n}{n!}\left(n+\frac{b}{a}\right)\cdots \left(1+\frac{b}{a}\right)p_0\\ pn=n!an(n+ab)⋯(1+ab)p0
得
p n = a n n ! ( − d + n − 1 ) ⋯ ( − d ) p 0 = ( − 1 ) n a n n ! ( d − n + 1 ) ⋯ d p 0 = ( − 1 ) n a n n ! d ! ( d − n ) ! p 0 = ( − a ) n ( d n ) p 0 \begin{align*} p_n&=\frac{a^n}{n!}(-d+n-1)\cdots(-d)p_0\\ &=(-1)^n\frac{a^n}{n!}(d-n+1)\cdots dp_0\\ &=(-1)^n\frac{a^n}{n!}\frac{d!}{(d-n)!}p_0\\ &=(-a)^n\left(\begin{matrix}d\\n\end{matrix}\right)p_0 \end{align*} pn=n!an(−d+n−1)⋯(−d)p0=(−1)nn!an(d−n+1)⋯dp0=(−1)nn!an(d−n)!d!p0=(−a)n(dn)p0
令 A = − a > 0 A=-a>0 A=−a>0,加上 ∑ i = 0 ∞ p i = 1 \sum^\infty_{i=0} p_i=1 ∑i=0∞pi=1,得
p 0 ∑ i = 0 d ( d n ) A n = 1 p_0\sum^d_{i=0}\left(\begin{matrix}d\\n\end{matrix}\right)A^n=1 p0i=0∑d(dn)An=1
记 A = p / 1 − p A=p/1-p A=p/1−p
则 p = A / ( 1 + A ) = a / ( a − 1 ) p=A/(1+A)=a/(a-1) p=A/(1+A)=a/(a−1), 0 < p < 1 0<p<1 0<p<1
上式变为
p 0 ∑ i = 0 d ( d n ) p n ( 1 − p ) − n = 1 p_0\sum^d_{i=0}\left(\begin{matrix}d\\n\end{matrix}\right) p^n(1-p)^{-n}=1 p0i=0∑d(dn)pn(1−p)−n=1
对比二项分布的密度函数,可得 p 0 = ( 1 − p ) d p_0=(1-p)^d p0=(1−p)d
此时分布为 B ( d , p ) B(d,p) B(d,p)