Почему оценка кумулятивной базовой опасности после Кокса PH в Стате отличается от R
Я - пользователь Stata, который сейчас изучает себя R. Я пытался подвести итоги анализов, выполненных в Stata, используя R, но застрял. Многовариантная регрессия Кокса в Stata и R дает идентичные результаты
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 базовая опасность
Я прочитал много онлайн-постов о базовой опасности, но не могу понять, что я сделал неправильно.