如何从R语言coxph对象中提取se(coef)值

你通过以下代码拟合了Cox比例风险模型：

library(survival)
library(survey)

res.cox <- coxph(Surv(time, status) ~ sex, data = lung)

已知模型结果里的se(coef)值为0.1672，使用svycontrast函数计算得到的SE值0.0983是错误的——因为svycontrast是为**复杂调查设计的模型（如svycoxph）**设计的，不适用于普通的coxph对象。

正确提取方法 #

直接从coxph对象中提取标准误，有两种简单方式：

• 方式一：直接访问对象元素
  coxph拟合结果对象自带se.coef元素，直接调用即可：

# 提取系数标准误
se_coef <- res.cox$se.coef
print(se_coef)
# 输出结果：0.1672（与模型结果一致）

• 方式二：从汇总结果中提取
  如果需要从模型汇总表中获取，可以使用summary()函数生成汇总后，提取对应列的值：

# 生成模型汇总
summary_res <- summary(res.cox)
# 提取se(coef)列
se_coef_from_summary <- summary_res$coefficients[, "se(coef)"]
print(se_coef_from_summary)

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