Я скопировал конкурирующую модель риска в Книга Терри стр. 53-4 и получаю сообщение об ошибке: идентификатор объекта не найден . Как я могу решить эту проблему? Вот код репликации
library(survival)
crdata <- mgus2
crdata$etime <- pmin(crdata$ptime, crdata$futime)
crdata$event <- ifelse(crdata$pstat==1, 1, 2*crdata$death)
crdata$event <- factor(crdata$event, 0:2, c("censor", "PCM", "death"))`
cfit <- coxph(Surv(etime, event) ~ I(age/10) + sex + mspike, id = id, crdata)