使用带有cph（）对象的幸存图的风险表编号

Question

这对于rms包非常具体。当使用带有cph对象的幸存图时，根据协变量，风险表中的数字不会分层。使用npsurv（）可以正确执行此操作。

library(survival);library(rms)
data(lung)
fit <- cph(Surv(time, status==2) ~ sex,data = lung,surv = TRUE,x=TRUE, y=TRUE)
survplot(fit, n.risk=TRUE)
# compared to:
 survplot(npsurv(fit$sformula, data = lung), conf= "none",n.risk=TRUE)

对于通常的用例我想用npsurv重新运行模型会做，但对于其他人，比如想要使用fit.mult.impute（我的具体用例），你想要使用创建的“原始”拟合通过cph（）。 有解决方法吗？

Answer 1

所以我应该使用

strat(sex)

作为术语。