lifelines.CoxPHFitter p值计算所用检验方法及适用规则咨询

• 首先要区分两类不同场景下的p值，二者对应检验方法不同，不存在混用情况：
  • 协变量系数对应的p值（也就是结果表中每个变量单独一行展示的p值），全部采用Wald检验计算：计算逻辑为估计得到的协变量系数除以对应标准误，构造z统计量后计算双侧检验p值。该计算逻辑在单变量Cox回归、多变量Cox回归场景下完全一致，不会因为纳入协变量的数量切换检验方法。
  • 模型整体显著性对应的p值，采用似然比检验计算：对比当前拟合模型与无协变量的空模型的对数似然差值，基于卡方分布得到p值，该值属于模型整体拟合优度指标，不会和单变量系数的p值混在同一列展示。
• 关于验证依据：

不需要依赖额外的公开说明文档，直接对照拟合输出的字段和手动复现即可确认：拟合完成后调用summary属性查看结果表，单系数对应的统计量列明确标注为z，对应Wald检验的z统计量，而非似然比检验使用的卡方统计量；手动按照Wald检验、似然比检验的公式复现计算，结果和库输出值完全匹配。

• 补充说明：从lifelines v0.2x至今的稳定版本中，该检验逻辑从未发生过变动，不存在版本差异导致的方法不一致问题。

内容的提问来源于stack exchange，提问作者Mischa