如果我使用lung
中的survival
数据集,我可以很容易地获得累积危害和幸存者估计值
lung.example <- survfit(Surv(time, status) ~ 1, data = lung)
# for survivor estimates
fortify(lung.example)
# for cumulative hazard estimates
fortify(lung.example, fun = "cumhaz")
是否有类似的“简便方法”提取我的危险功能?我已经搜索了很长时间,而且看起来很复杂。 我不确定为什么这种方式比累积危害要复杂。
谢谢