[转]生存分析-1

转自生存分析(1)

临床研究中,生存分析对于一项干预措施或者是危险因素的评估是一种关键方法。生存分析对应于一组统计方法,用于调查感兴趣事件发生所花费的时间。生存分析的概念是广泛的,可以是指通常意义下的生存时间,也可以是指所关心的某个现象(如疾病愈后,合格品使用后)。多见到的生存分析概念有总体生存分析(overrall survival analysis,OS),无病进展生存期(Progression-Free Survival,PFS)。

1. 肿瘤研究中,典型的问题

  • 某些临床特征对患者生存的影响是什么
  • 一个人生存3年的概率是多少?
  • 两组患者的生存率是否存在差异?

2.几个概念
2.1 总体生存分析
总体生存是指,从病人确认患有疾病开始至因任何原因引起死亡的时间。该指标常常被认为是肿瘤临床试验中最佳的疗效终点。确认病人因病或其他因素引起死亡的日期通常几乎没有困难,并且死亡的时间有其独立的因果关系。如果在生存期上有小幅度的提高,可以认为是有意义的临床受益证据。

2.2 无病进展生存期
所谓无病进展生存期(Progression-Free Survival,PFS)通常定义为病人经过治疗,随机选择某个时间直到肿瘤复发或因各种原因出现死亡,病人总的生存时间。

PFS的优点在于它能反映肿瘤的生长(这个现象可能反映了肿瘤相关疾病或死亡的因果联系),可以于生存获益证实前被评价,不会受到后续治疗的潜在的易混淆的指标或症状影响。而且PFS的结果比生存期结果出现得更早,治疗过程中,病人一旦出现了症状,肿瘤复发了,过了无病进展生存期就要采取其他积极治疗手段,从而进一步改善患者的症状,延长生存时间。

2.3 无复发生存时间
无复发生存时间( relapse-free survival time) 是指从治疗到疾病复发之间的时间。也称为无病生存时间(disease-free survival time) 和无事件生存时间 (event-free survival time) 。

3. 生存函数和风险函数
使用生存概率和危险概率这两个概念来描述生存数据。

3.1 生存函数
生存函数用S(t)表示, 是指观察对象生存时间T大于某时刻t的概率。

S(t) = Prob(T > t)

在具体问题中,该函数t时刻的取值可以用以下公式来估计

S(t) = 生存时间长于t的观察对象人数 / 观察对象总人数

显然,S(t) 是一个随时间增加而下降的函数,表示观察对象随访到t时刻的累计生存率。

3.2 风险函数
风险函数用h(t)表示,指的是是在时间t被观察的个体在该时间发生事件的概率。

风险函数随着时间的延长可呈现为递增、递减或者其他的波动形式。如果风险函数为常数时,表示死亡速率不随着时间而加速;如果风险函数随着时间上升,则表示死亡速率随时间而加速;风险函数随着时间下降,则表示死亡速率随时间减速。

4. Kaplan-Meier 生存分析
Kaplan-Meier 生存分析,即就是乘积极限法,是一种非参数方法,用于根据观察到的生存时间估算生存概率(Kaplan和Meier,1958年)。

时刻的生存概率计算如下:

S ( t i ) = S ( t i − 1 ) ( 1 − d i n i ) S(t_i)=S(t_{i-1})(1-\frac{d_i}{n_i}) S(ti)=S(ti1)(1nidi)

S ( t i − 1 ) S(t_{i-1}) S(ti1) = 在 t i − 1 t_{i-1} ti1存活的概率
n i n_i ni = t i t_i ti之前还活着的病人数量
d i d_i di是指再时刻事件的发生的书目
t 0 t_0 t0=0, S ( 0 ) = 1 S(0)=1 S(0)=1

估计概率(S(t))是仅在每个事件发生时才改变值的阶跃函数。同样也可以计算生存概率的置信区间。

KM生存曲线是KM生存概率与时间的关系图,提供了可用于估计度量值(例如中位生存时间)的数据。

5.Cox回归分析
5.1 什么是Cox回归模型
Cox模型是目前生存分析多因素预后评价中较好的统计学方法。Cox将模型以半参数方式出现,适用于许多分布未知的资料的多因素分析,可以再许多因素共存的情况下,排除混杂因子的影响,提高预后分析的质量,并能处理截尾数据。

上面提到Cox是半参数模型,什么是半参数模型

参数模型回归分析:已知生存时间服从特定的参数模型时,拟合相应的参数模型,更准确地分析确定变量之间的变化规律。
缺点:使用参数法进行生存分析需要事先知道生存时间的分布。 优点:1.
能比较两组或多组生存时间分布函数。2.能分析危险因素对生存时间的影响。3.能建立生存时间与危险因素之间的依存关系的模型。

非参数检验:检验分组变量各水平所对应的生存曲线是否一致。 对生存时间的分布没有要求。常用的方法有乘积极限法(PL法)和寿命表法(LT法)。

半参数模型回归分析:不需要对生存时间的分布作出假定,但却可以通过一个模型来分析生存时间的分布规律,以及危险因素对生存时间的影响。

Cox比例风险回归分析法是这种半参数模型的代表,在表达形式上与参数模型相似,但在对模型中的各个参数进行估计时又不依赖于特定分布假设。

  • 1
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
在 MATLAB 中进行生存分析通常涉及使用统计函数来估计生存曲线、制作 Kaplan-Meier 曲线和进行 Cox 比例风险回归等分析。以下是一些常用的 MATLAB 函数和步骤: 1. 导入数据:可以从文件中导入生存数据,其中包含每个观察的生存时间和事件状态(存活/死亡)。 2. 估计生存曲线:使用 `survfit` 函数来估计生存曲线(Kaplan-Meier 曲线)。该函数根据观察到的生存时间和事件状态计算出生存曲线的估计值。 3. 绘制生存曲线:使用 `plot` 函数将估计的生存曲线可视化。 ```matlab % 导入数据 data = readtable('survival_data.csv'); % 将数据换为适合进行生存分析的格式 time = data.Time; status = data.Status; % 估计生存曲线 sf = survfit(time, status); % 绘制生存曲线 plot(sf); ``` 4. 进行 Cox 比例风险回归:使用 `coxphfit` 函数拟合 Cox 比例风险模型,并使用 `coxphval` 函数进行预测。 ```matlab % 导入数据 data = readtable('survival_data.csv'); % 将数据换为适合进行 Cox 比例风险回归的格式 time = data.Time; status = data.Status; covariates = data.Covariates; % 拟合 Cox 比例风险模型 coxModel = coxphfit(covariates, time, status); % 进行预测 predictedSurvival = coxphval(coxModel, covariates); ``` 这些是一些基本的 MATLAB 函数和步骤,可用于进行生存分析。根据具体的需求,你可能需要进一步探索和使用其他相关函数来完成更复杂的分析。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值