生存曲线(三):统计分析方法这么多,到底选哪个?

前面花了两期内容,讨论了如何使用软件(包括GraphPad Prism、SPSS和Origin)绘制生存曲线。

生存曲线只是数据呈现的方式之一,其核心必须是统计结果。涉及统计,就不可避免地需要选择使用哪种统计方法。有时,你的数据没有统计意义,还真得回头看看统计方法是否选择正确了。

打开SPSS统计软件,录好数据后,点击生存分析,我们会发现软件中有4种可供使用的统计方法,此时到底该选择哪个呢?

本期主要以SPSS统计分析为例进行讨论,毕竟这个属于专业的统计软件,更加全面。

① 寿命表法

SPSS软件官方解释是这样描述的:用从每个区间估计的概率估计在不同时间点发生该事件的整体概率。“寿命表的基本概念是将观察区间划分为较小的时间区间。对于每个区间,使用所有观察至少该时长的人员计算该区间内发生终端事件的概率。然后使用从每个区间估计的概率估计在不同时间点发生该事件的整体概率。”

上来就整这么一段话,估计打击有点懵。

举个例子。假如一个临床实验需要随访2000人的生存时间(按月计算),此时会出现两个问题。第一,我们很难做到及时、有效地随访每一位病人,不可能每个月都去问病人是否还活着?第二,如此大量的数据,我们在最终统计时,如果挨个录入,费时费力,极易出错。

别急,寿命表法为我们提供了有力帮助。

寿命表法本质是通过人为划分时间区间的方式,以频数方式呈现数据,采用加权法进行生存率分析的方法。注意,此时的时间变量是不连续的。

例如上面的例子,我们可以人为划分时间区间为3个月,那么只要病人死亡发生在这个3个月内,都将其归入此区间,极大地减轻了随访难度。同时,以频数的方式呈现数据,无论多少个病人,只要其死亡时间在这个区间,频数就是多少,不用再细分每个死亡患者的具体死亡时间,可降低整理数据的难度。

因此,寿命表法适用于临床研究中的终点事件无法被有效、及时随访和大样本数据的生存分析。

② Kaplan-Meier法(又称乘积极限法)

SPSS官方解释又说了:Kaplan-Meier法与寿命表法的关键区别在于时间变量为连续变量,状态变量可以是分类变量或连续变量,因子和层次变量应为分类变量。

基于此,我们反推一下,可发现Kaplan-Meier法的要求需要时间变量为连续状态,就是说需要我们获取完整的、精确的生存资料。Kaplan-Meier法本质是一个生存时间函数。

再举个例子,动物实验过程中,我们非常容易控制是动物数量、品系、性别的一致性,同时也冷非常准确及时地获取动物生存时间。有时,前脚刚出动物房,管理员就来电话说动物死了,这效率,杠杠的。此外,动物实验涉及的样本量一般不会特别大,就算一组动物20只,分5组,也才100只。

因此,Kaplan-Meier法比较适合动物实验这种小样本量的、观测记录较为完整的数据类或者大样本但未划分时间段的数据。

③ Cox 回归分析

SPSS官方解释贴心指出:Cox 回归数据分析时,其时间变量应是定量变量,但状态变量可以是分类或连续变量。

Cox回归分析本质是多因素分析,与之对应的上面的单因素分析(Kaplan-Meier法)。

理想状态下,我们非常希望有且只有1个因素影响生存率,便于我们更好地评价治疗作用。但现实往往不如人愿。临床研究中,经常会出现多个因素影响患者生存率,如果简单地套用单因素分析法,非常容易得出错误的研究结论。因此,我们必须要做Cox回归分析,充分评价观测因素以外的其它因素在患者生存率上回带来多大的影响。

举个例子。评价新药和市售药治疗糖尿病时患者的生存率差异。我们都知道很多糖尿病患者都存在心梗、脑梗的风险,这些风险发生时均会导致患者死亡,必然影响我们最终的研究结论。此时即可采用Cox分析法。

重点来了。尽管Cox回归分析看起来很符合实际操作情况,但是Cox回归分析也称为Cox 比例风险模型

SPSS广泛统计中强调:观察值应是独立的,风险比应是时间恒定值;即,各个个案风险的比率不应随时间变化。后一个假设称为Cox 比例风险模型。

什么叫等比例风险?预后因素对死亡风险的作用强度随着时间的变化需保持一致。

即,假设A、B两组糖尿病患者,分别给予新药和市售药,随访5年生存率。第1年A组风险死亡率是6%,B组是3%。第2年A组风险死亡率是10%,B组需为5%。此时才刚刚满足Cox风回归分析的前提。

验证是否为等比例风险的方法大致分两类。分类变量采用K-M生存曲线法查看是否存在交叉(若不交叉,则说明等比例风险);连续型变量采用残差散点图来判定(P>0.05表示)。(图例来自网络,侵删)

(K-M生存曲线法)

(残差散点图)

有人又会说了,这要求忒严格了吧,如果不是等比例风险呢?见下方。

④ Cox 依时协变量分析

若判定后发现不是等比例风险,说明影响影响因素在随着时间发生变化,此时可能需要采用带依时协变量的 Cox 过程。

感觉越说越远了......懵了吗朋友们

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 212,222评论 6 493
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 90,455评论 3 385
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 157,720评论 0 348
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 56,568评论 1 284
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 65,696评论 6 386
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 49,879评论 1 290
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,028评论 3 409
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 37,773评论 0 268
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,220评论 1 303
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 36,550评论 2 327
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 38,697评论 1 341
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,360评论 4 332
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,002评论 3 315
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 30,782评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,010评论 1 266
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 46,433评论 2 360
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 43,587评论 2 350