我已经审查了其他类似的问题,我不相信我找到了答案。
我正在尝试使用两个协变量进行 Cox 回归 - 性别和疾病状态。
原始数据框看起来有点像这样:
Patient ID: 1001, 1002
Age: 56, 60
Sex: Male, Female
Mortality event: 1 0
Follow up years: 6,7
我已经打电话给cxmod <- coxph(Surv(Mortality event, time) ~ Disease_status + Sex, data = original data)
我已经根据协变量的这个包的说明设置了一个dummy_df
作为网格:
Sex Male Disease_status 0,
Sex Female Disease_status 0,
Sex Male Disease status 1,
Sex Female Disease status 1
我已将行名重命名为字母,因为我知道这是需要的。
但是,当我打电话时:
cxsf <- survfit(cxmod, data= orginal_data_frame, newdata = dummy_df, conf.type = "none")
我收到以下错误消息:
Warning message:
'newdata' had 4 rows but variables found have 500000 rows
此外,如果我调用surv_summary(cxsf)
来帮助可视化情节 - R 会话会因遇到致命错误而终止。
谁能就出了什么问题提供建议?
您不需要提供原始数据集来survfit,因为模型已经使用它进行了训练,并且现在参数存在于模型本身(cxmod(中。
survfit(cxmod , newdata = dummy_df, conf.type = "none")