乐于分享
好东西不私藏

【更新】肝帝新增功能说明

【更新】肝帝新增功能说明

关联分析研究数据分析-肝帝  

新增功能说明:Table 3B 马尔科夫多状态生存分析

一、这个功能叫什么?

在关联分析的肝帝模块里,当您为用户配置了随访时间并完成分析后,系统会在原来的多因素回归表(现更名为 Table 3A)之外,自动多出一部分结果:Table 3B — 马尔科夫多状态生存分析(Markov Multi-State Survival Analysis)。

可以理解为:在「时间 + 事件」的前提下,多了一把用来描述「人随时间在各状态之间怎么流动」的统计工具,作为 Cox 等多因素生存分析结果的补充,而不是替代。

二、一句话讲清楚它在干什么?

把每个研究对象想成一开始都在「尚未发生目标结局」的状态;随着时间推移,有的人一直保持在该状态(相当于仍未发生您关心的那件「事」),有的人进入目标结局,在部分数据里还有人进入其它类型的重要事件(竞争风险)。Table 3B 用马尔科夫多状态模型来描述这些状态随时间怎么变化,并给出主暴露因素对「走向哪一类结局」的影响(用类似风险比的方式报告),同时还能画出未来各时点「还有多少人仍在风险集、多少人已发生结局、多少人走了竞争事件」的比例图。

三、原理(通俗易懂版)

1.三个「状态」是什么?

状态 0:起点——随访开始时,大家算在「尚未发生目标结局」的池子里(仍可能被观察到事件发生)。

状态 1:目标结局发生了(您最关心的那类事件)。

状态 2:竞争风险——若您的结局变量编码里包含「另一类会「终结」随访、且与目标结局互斥」的事件(例如死亡或其它 major event),模型会把它当作竞争事件;若没有这类编码,系统会自动退化成只有两个方向有意义的简化模型(相当于没有第二条「竞争」出路)。

状态 1 和 2 在模型里被设为吸收态:一旦发生,就像「已经到站」,不再回到状态 0——这符合临床上「发生过的事件通常不可逆」的常识。

2.「分段常数强度」是什么意思?

现实里,风险不一定随时间完全「一条直线」不变。模型把随访时间轴切成若干时间段,在每一小段里假设基准转移强度是常数(像阶梯一样一节一节变),这样能在一定程度上反映风险随时间变化,比简单假设「从头到尾比例风险一成不变」更灵活一些,也更适合作为方法学上的补充叙述。

3.暴露和协变量怎么进来?

在每个时间段、每条「从状态 0 走出去」的路径上,转移强度可以随主暴露 X和协变量 Z按指数形式变化(与 Cox 里「乘性风险」的直觉相近)。系统用数据做极大似然估计,给出暴露对走向状态 1(以及在有竞争风险时对走向状态 2)的作用,并报告 HR、95% 置信区间、P 值等,便于和 Table 3A 对照阅读。

4.「状态占据概率」和 stacked 面积图

在估计好模型后,可以算出:在随访第 t 时刻,人群中平均还有多大比例「仍在状态 0」、多大比例「已进入结局 1」、多大比例「已进入竞争事件 2」。把这些比例随时间画成堆叠面积图,就很直观:例如可以对比暴露组与参考组在未来若干时间点上,仍处于未发病与已发病的构成差异——这是 Table 3A 单独给 HR 时不容易一眼看出来的「动态画面」。

四、有什么用?(给临床与合作者的价值)

与 Cox / 多因素生存结果形成互补

Table 3A 仍回答「调整协变量后,暴露对事件风险的比值型关联」;Table 3B 在此基础上补充多状态结构和随时间的状态比例,便于写进方法学或讨论:我们不仅报告 HR,还展示了状态占用随时间的演变。

竞争风险场景下语义更清晰

若结局编码支持「主事件 + 竞争事件」,模型会同时考虑通向两类吸收态的路径,避免把所有非主事件简单等同为「截尾」而忽略其临床含义(具体仍依赖您的变量定义与编码规范)。

沟通成本更低

堆叠面积图对非统计背景的医生、审稿人更友好:强调的是「随时间的人群在各状态上的占比」,而不是只有一张系数表。

全自动触发,不增加配置负担

只要用户在肝帝里配置了随访时间列,完整分析跑完后 Table 3B 会自动出现;无需单独勾选「马尔科夫」开关(与产品当前实现一致)。

五、和 Table 3A 怎么称呼?

Table 3A:多因素回归(含随访时为 Cox 等)——原来的「Table 3」已更名为 3A,避免与附表里的 Supplemental Table 3 混淆。

Table 3B:马尔科夫多状态生存分析——排在 3A 之后,作为补充表格与图示。

六、使用与解读时注意什么?(简要提醒)

结局变量需按约定编码:0 = 删失,1 = 目标结局,2 = 竞争风险(若有);若用肝帝其它模块,请与团队统一编码说明。

样本量过小时,复杂模型可能不稳定,界面可能不展示 Table 3B 或结果需谨慎解读。

Table 3B 是补充性、方法学增强型输出,主结论仍建议以预设主分析(如 Table 3A)与方案为准。