mobile wallpaper 1mobile wallpaper 2mobile wallpaper 3mobile wallpaper 4
665 字
2 分钟
生存分析中四个基本函数的数学关系
2026-09-18

1. 记号与基本假设#

设随机变量 X≥0X\ge 0 表示个体从研究起点到事件发生的时间。以下推导主要针对连续型生存时间,并假设:

  • XX 的密度函数为 f(x)f(x),分布函数为 F(x)F(x);
  • 生存函数 S(x)>0S(x)>0 且可微;
  • S(0)=1S(0)=1,即事件不会在研究起点之前或恰好于起点发生;
  • 所涉及的积分有限,并满足必要的极限条件,例如 lim⁡t→∞(t−x)S(t)=0\lim_{t\to\infty}(t-x)S(t)=0。

四个基本函数分别为:

  1. 密度函数

    f(x)=ddxF(x).f(x)=\frac{d}{dx}F(x).
  2. 生存函数

    S(x)=P(X>x)=1−F(x).S(x)=P(X>x)=1-F(x).
  3. 风险函数(瞬时风险率)

    h(x)=lim⁡Δx→0+P(x≤X<x+Δx∣X≥x)Δx.h(x)=\lim_{\Delta x\to 0^+} \frac{P(x\le X<x+\Delta x\mid X\ge x)}{\Delta x}.
  4. 平均剩余寿命函数(mean residual life)

    mrl(x)=E(X−x∣X>x).mrl(x)=E(X-x\mid X>x).

为简化记号,下面将 mrl(x)mrl(x) 记为 m(x)m(x)。对于连续型 XX,条件中的 X>xX>x 与 X≥xX\ge x 没有区别。


2. 密度函数 f(x)f(x) 与生存函数 S(x)S(x)#

由生存函数的定义,

S(x)=P(X>x)=∫x∞f(t) dt.S(x)=P(X>x)=\int_x^\infty f(t)\,dt.

对 xx 求导,根据微积分基本定理,

f(x)=−S′(x).f(x)=-S'(x).

因此二者可以相互转换:

S(x)=∫x∞f(t) dt,f(x)=−dS(x)dx.S(x)=\int_x^\infty f(t)\,dt, \qquad f(x)=-\frac{dS(x)}{dx}.

由于 f(x)≥0f(x)\ge 0,所以 S′(x)≤0S'(x)\le 0,即生存函数必为单调不增函数。


3. 风险函数 h(x)h(x) 与 f(x),S(x)f(x),S(x)#

从风险函数的定义出发:

h(x)=lim⁡Δx→0+P(x≤X<x+Δx∣X≥x)Δx=lim⁡Δx→0+P(x≤X<x+Δx)Δx P(X≥x)=lim⁡Δx→0+F(x+Δx)−F(x)Δx S(x)=f(x)S(x).\begin{aligned} h(x) &=\lim_{\Delta x\to0^+} \frac{P(x\le X<x+\Delta x\mid X\ge x)}{\Delta x}\\ &=\lim_{\Delta x\to0^+} \frac{P(x\le X<x+\Delta x)}{\Delta x\,P(X\ge x)}\\ &=\lim_{\Delta x\to0^+} \frac{F(x+\Delta x)-F(x)}{\Delta x\,S(x)}\\ &=\frac{f(x)}{S(x)}. \end{aligned}

所以

h(x)=f(x)S(x).h(x)=\frac{f(x)}{S(x)}.

再利用 f(x)=−S′(x)f(x)=-S'(x),得到

h(x)=−S′(x)S(x)=−ddxln⁡S(x).h(x)=-\frac{S'(x)}{S(x)} =-\frac{d}{dx}\ln S(x).

即

h(x)=−ddxln⁡S(x).h(x)=-\frac{d}{dx}\ln S(x).

3.1 累积风险函数#

定义累积风险函数

H(x)=∫0xh(u) du.H(x)=\int_0^x h(u)\,du.

对 h(x)=−(ln⁡S(x))′h(x)=-(\ln S(x))' 积分,并使用 S(0)=1S(0)=1,可得

H(x)=−ln⁡S(x).H(x)=-\ln S(x).

所以

S(x)=e−H(x)=exp⁡[−∫0xh(u) du].S(x)=e^{-H(x) =\exp\left[-\int_0^x h(u)\,du\right]}.

继而

f(x)=h(x)S(x)=h(x)exp⁡[−∫0xh(u) du].f(x)=h(x)S(x) =h(x)\exp\left[-\int_0^x h(u)\,du\right].

这说明:只要给定风险函数 h(x)h(x),原则上就能恢复整个生存时间分布。

注意:h(x)h(x) 本身不是概率,而是单位时间内的瞬时条件事件率;h(x)Δxh(x)\Delta x 才近似表示已经生存到 xx 时,在接下来长度为 Δx\Delta x 的短时间内发生事件的条件概率。


4. 平均剩余寿命 m(x)m(x) 与 f(x),S(x)f(x),S(x)#

根据条件期望的定义,

m(x)=E(X−x∣X>x)=∫x∞(t−x)f(t) dtS(x).\begin{aligned} m(x) &=E(X-x\mid X>x)\\ &=\frac{\displaystyle\int_x^\infty (t-x)f(t)\,dt}{S(x)}. \end{aligned}

因此

m(x)=∫x∞(t−x)f(t) dtS(x).m(x)=\frac{\displaystyle\int_x^\infty (t-x)f(t)\,dt}{S(x)} .

利用 f(t)=−S′(t)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.\begin{aligned} \int_x^\infty (t-x)f(t)\,dt &=\int_x^\infty (t-x)[-S'(t)]\,dt\\ &=\left[-(t-x)S(t)\right]_x^\infty +\int_x^\infty S(t)\,dt\\ &=\int_x^\infty S(t)\,dt. \end{aligned}

所以平均剩余寿命也可写为

m(x)=∫x∞S(t) dtS(x).m(x)=\frac{\displaystyle\int_x^\infty S(t)\,dt}{S(x)} .

当 x=0x=0 时,

m(0)=E(X)=∫0∞S(t) dt,m(0)=E(X)=\int_0^\infty S(t)\,dt,

即平均寿命等于生存曲线下的总面积。

还可把上式写成条件生存概率的积分。因为

P(X>t∣X>x)=S(t)S(x),t≥x,P(X>t\mid X>x)=\frac{S(t)}{S(x)},\qquad t\ge x,

所以

m(x)=∫x∞S(t)S(x) dt.m(x)=\int_x^\infty \frac{S(t)}{S(x)}\,dt.

5. 平均剩余寿命 m(x)m(x) 与风险函数 h(x)h(x)#

令

A(x)=∫x∞S(t) dt,m(x)=A(x)S(x).A(x)=\int_x^\infty S(t)\,dt, \qquad m(x)=\frac{A(x)}{S(x)}.

由

A′(x)=−S(x),S′(x)=−h(x)S(x),A'(x)=-S(x), \qquad S'(x)=-h(x)S(x),

对 m(x)=A(x)/S(x)m(x)=A(x)/S(x) 求导:

m′(x)=A′(x)S(x)−A(x)S′(x)S2(x)=−S2(x)+A(x)h(x)S(x)S2(x)=−1+h(x)m(x).\begin{aligned} m'(x) &=\frac{A'(x)S(x)-A(x)S'(x)}{S^2(x)}\\ &=\frac{-S^2(x)+A(x)h(x)S(x)}{S^2(x)}\\ &=-1+h(x)m(x). \end{aligned}

于是得到核心微分关系

m′(x)=h(x)m(x)−1.m'(x)=h(x)m(x)-1.

当 m(x)>0m(x)>0 时,可反解出风险函数:

h(x)=1+m′(x)m(x).h(x)=\frac{1+m'(x)}{m(x)}.

该式也说明:

m′(x){>0,h(x)m(x)>1,=0,h(x)m(x)=1,<0,h(x)m(x)<1.m'(x) \begin{cases} >0, & h(x)m(x)>1,\\ =0, & h(x)m(x)=1,\\ <0, & h(x)m(x)<1. \end{cases}

因此,平均剩余寿命在某一点增加还是减少,取决于该点的风险率与平均剩余寿命之积。

5.1 由风险函数直接计算平均剩余寿命#

由

S(t)S(x)=exp⁡[−∫xth(u) du],\frac{S(t)}{S(x)} =\exp\left[-\int_x^t h(u)\,du\right],

可得

m(x)=∫x∞exp⁡[−∫xth(u) du]dt.m(x)=\int_x^\infty \exp\left[-\int_x^t h(u)\,du\right]dt .

6. 由平均剩余寿命反推 S(x)S(x) 与 f(x)f(x)#

由上一节可知

h(x)=1+m′(x)m(x)=1m(x)+m′(x)m(x).h(x)=\frac{1+m'(x)}{m(x)} =\frac{1}{m(x)}+\frac{m'(x)}{m(x)}.

又因为 h(x)=−(ln⁡S(x))′h(x)=-(\ln S(x))',所以

ddxln⁡S(x)=−1m(x)−m′(x)m(x).\frac{d}{dx}\ln S(x) =-\frac{1}{m(x)}-\frac{m'(x)}{m(x)}.

从 00 到 xx 积分,并使用 S(0)=1S(0)=1:

ln⁡S(x)=−∫0xdum(u)−∫0xm′(u)m(u) du=−∫0xdum(u)−ln⁡m(x)+ln⁡m(0).\begin{aligned} \ln S(x) &=-\int_0^x\frac{du}{m(u)} -\int_0^x\frac{m'(u)}{m(u)}\,du\\ &=-\int_0^x\frac{du}{m(u)} -\ln m(x)+\ln m(0). \end{aligned}

因此

S(x)=m(0)m(x)exp⁡[−∫0xdum(u)].S(x)=\frac{m(0)}{m(x)} \exp\left[-\int_0^x\frac{du}{m(u)}\right].

再利用 f(x)=h(x)S(x)f(x)=h(x)S(x),得到

f(x)=m(0)[1+m′(x)]m2(x)exp⁡[−∫0xdum(u)].f(x)= \frac{m(0)[1+m'(x)]}{m^2(x)} \exp\left[-\int_0^x\frac{du}{m(u)}\right] .

所以,在 m(x)>0m(x)>0、可微、m′(x)≥−1m'(x)\ge -1 且所得 S(x)S(x) 满足适当边界条件时,平均剩余寿命函数同样能够唯一确定生存分布。

分享

如果这篇文章对你有帮助,欢迎分享给更多人!

生存分析中四个基本函数的数学关系
https://umi-phi.vercel.app/posts/survival_analysis_01/
作者
Umi
发布于
2026-09-18
许可协议
Unlicensed

目录