单解释变量Cox比例风险模型全/简化模型构建及似然比检验验证

验证结果 #

你编写的代码逻辑完全正确，已经实现了你需要的全模型与简化模型的似然比检验。

对应逻辑说明 #

• coxph()返回的loglik向量固定存储两个对数似然值：
  • model$loglik[1]：无任何解释变量的零模型（也就是你提到的采用整个数据集整体均值的简化模型）的对数似然值
  • model$loglik[2]：纳入二分类变量z的拟合模型（也就是你提到的按分组计算均值的全模型）的对数似然值
• 似然比检验的标准统计量计算公式就是2倍的全模型与简化模型对数似然的差值，和你的代码逻辑完全匹配。

补充使用提示 #

你的全模型仅比简化模型多1个待估参数（二分类变量z仅对应1个回归系数），因此得到的X.lr统计量服从自由度为1的卡方分布，你可以通过以下代码直接计算检验p值：

p_value <- pchisq(X.lr, df = 1, lower.tail = FALSE)

你也可以直接运行summary(model)，输出结果中Likelihood ratio test行给出的统计量和p值和你手动计算的结果完全一致，无需自行实现计算逻辑。

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