我是Stata用户,现在正在教我自己R.我一直试图概括使用R在Stata中完成的分析但是已经卡住了。 Stata和R中的多变量Cox回归给出相同的结果
stcox ag1 ag2 sz nd grade screen , offset(hc1) nohr
mva <- coxph(fit ~ ag1 + ag2 + sz + nd + grade + screen +
offset(hc1), data = er1)
然而,当我得出累积基线危险函数时,它们是不同的
predict h0, basechazard
h0 <- basehaz(mva, centered = F)
R基线危害
我已经阅读了大量关于基线危害的在线帖子,但无法弄清楚我做错了什么。