1. 记号与基本假设#
设随机变量 X≥0 表示个体从研究起点到事件发生的时间。以下推导主要针对连续型生存时间,并假设:
- X 的密度函数为 f(x),分布函数为 F(x);
- 生存函数 S(x)>0 且可微;
- S(0)=1,即事件不会在研究起点之前或恰好于起点发生;
- 所涉及的积分有限,并满足必要的极限条件,例如
limt→∞(t−x)S(t)=0。
四个基本函数分别为:
-
密度函数
f(x)=dxdF(x).
-
生存函数
S(x)=P(X>x)=1−F(x).
-
风险函数(瞬时风险率)
h(x)=Δx→0+limΔxP(x≤X<x+Δx∣X≥x).
-
平均剩余寿命函数(mean residual life)
mrl(x)=E(X−x∣X>x).
为简化记号,下面将 mrl(x) 记为 m(x)。对于连续型 X,条件中的 X>x 与 X≥x 没有区别。
2. 密度函数 f(x) 与生存函数 S(x)#
由生存函数的定义,
S(x)=P(X>x)=∫x∞f(t)dt.对 x 求导,根据微积分基本定理,
f(x)=−S′(x).因此二者可以相互转换:
S(x)=∫x∞f(t)dt,f(x)=−dxdS(x).由于 f(x)≥0,所以 S′(x)≤0,即生存函数必为单调不增函数。
3. 风险函数 h(x) 与 f(x),S(x)#
从风险函数的定义出发:
h(x)=Δx→0+limΔxP(x≤X<x+Δx∣X≥x)=Δx→0+limΔxP(X≥x)P(x≤X<x+Δx)=Δx→0+limΔxS(x)F(x+Δx)−F(x)=S(x)f(x).所以
h(x)=S(x)f(x).再利用 f(x)=−S′(x),得到
h(x)=−S(x)S′(x)=−dxdlnS(x).即
h(x)=−dxdlnS(x).3.1 累积风险函数#
定义累积风险函数
H(x)=∫0xh(u)du.对 h(x)=−(lnS(x))′ 积分,并使用 S(0)=1,可得
H(x)=−lnS(x).所以
S(x)=e−H(x)=exp[−∫0xh(u)du].继而
f(x)=h(x)S(x)=h(x)exp[−∫0xh(u)du].这说明:只要给定风险函数 h(x),原则上就能恢复整个生存时间分布。
注意:h(x) 本身不是概率,而是单位时间内的瞬时条件事件率;h(x)Δx 才近似表示已经生存到 x 时,在接下来长度为 Δx 的短时间内发生事件的条件概率。
4. 平均剩余寿命 m(x) 与 f(x),S(x)#
根据条件期望的定义,
m(x)=E(X−x∣X>x)=S(x)∫x∞(t−x)f(t)dt.因此
m(x)=S(x)∫x∞(t−x)f(t)dt.利用 f(t)=−S′(t) 对分子作分部积分:
∫x∞(t−x)f(t)dt=∫x∞(t−x)[−S′(t)]dt=[−(t−x)S(t)]x∞+∫x∞S(t)dt=∫x∞S(t)dt.所以平均剩余寿命也可写为
m(x)=S(x)∫x∞S(t)dt.当 x=0 时,
m(0)=E(X)=∫0∞S(t)dt,即平均寿命等于生存曲线下的总面积。
还可把上式写成条件生存概率的积分。因为
P(X>t∣X>x)=S(x)S(t),t≥x,所以
m(x)=∫x∞S(x)S(t)dt.
5. 平均剩余寿命 m(x) 与风险函数 h(x)#
令
A(x)=∫x∞S(t)dt,m(x)=S(x)A(x).由
A′(x)=−S(x),S′(x)=−h(x)S(x),对 m(x)=A(x)/S(x) 求导:
m′(x)=S2(x)A′(x)S(x)−A(x)S′(x)=S2(x)−S2(x)+A(x)h(x)S(x)=−1+h(x)m(x).于是得到核心微分关系
m′(x)=h(x)m(x)−1.当 m(x)>0 时,可反解出风险函数:
h(x)=m(x)1+m′(x).该式也说明:
m′(x)⎩⎨⎧>0,=0,<0,h(x)m(x)>1,h(x)m(x)=1,h(x)m(x)<1.因此,平均剩余寿命在某一点增加还是减少,取决于该点的风险率与平均剩余寿命之积。
5.1 由风险函数直接计算平均剩余寿命#
由
S(x)S(t)=exp[−∫xth(u)du],可得
m(x)=∫x∞exp[−∫xth(u)du]dt.
6. 由平均剩余寿命反推 S(x) 与 f(x)#
由上一节可知
h(x)=m(x)1+m′(x)=m(x)1+m(x)m′(x).又因为 h(x)=−(lnS(x))′,所以
dxdlnS(x)=−m(x)1−m(x)m′(x).从 0 到 x 积分,并使用 S(0)=1:
lnS(x)=−∫0xm(u)du−∫0xm(u)m′(u)du=−∫0xm(u)du−lnm(x)+lnm(0).因此
S(x)=m(x)m(0)exp[−∫0xm(u)du].再利用 f(x)=h(x)S(x),得到
f(x)=m2(x)m(0)[1+m′(x)]exp[−∫0xm(u)du].所以,在 m(x)>0、可微、m′(x)≥−1 且所得 S(x) 满足适当边界条件时,平均剩余寿命函数同样能够唯一确定生存分布。