Почему оценка кумулятивной базовой опасности после Кокса 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)

Stata базовая опасность

R базовая опасность

Я прочитал много онлайн-постов о базовой опасности, но не могу понять, что я сделал неправильно.

0 ответов

Другие вопросы по тегам