警告：'newdata'具有x行，但是找到的变量具有> x行且带有coxph和survfit

Question

我已经审查了其他类似的问题，但我不相信我已经找到答案了。

我正在尝试使用两个协变量进行性别与疾病状态的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会话会遇到致命错误而终止。

任何人都可以对发生的问题提出建议吗？

Answer 1

您不需要提供原始数据集就可以生存，因为已经使用它训练了模型，现在参数本身就存在于模型本身（cxmod）中。

survfit(cxmod , newdata = dummy_df, conf.type = "none")