Я хочу вычислить предсказанную функцию выживания для модели пропорциональных рисков Кокса с ненадежными терминами [используя пакет выживания]. Похоже, что когда в модели присутствуют слабые члены, предсказанная функция выживания не может быть вычислена.
## Example
require(survival)
data(rats)
## Create fake weight
set.seed(90989)
rats$weight<-runif(nrow(rats),0.2,0.9)
## Cox model with gamma frailty on litter
fit <- coxph(Surv(time, status) ~ rx+weight+frailty(litter,dist="gamma"),
data = rats)
## Compute survival curve from the cox model for rx=0 and weight=0.5 kg
plot(survfit(fit, newdata=data.frame(rx=0,weight=0.5)),xlab = "time",
ylab="Survival")
## Running this line, I get following error message:
Error in survfit.coxph(fit, newdata = data.frame(rx = 0, weight = 0.5)) :
Newdata cannot be used when a model has sparse frailty terms
Я пытался использовать как разреженные, так и не разреженные методы вычислений с использованием
sparse=TRUE, Sparse =FALSE, sparse =0, sparse=5
опций. Однако никто не работал.
Как рассчитать прогнозную кривую выживаемости на основе моей модели слабости?