Я хотел бы использовать функцию pregnetSurvProb из пакета pe c, чтобы предсказать вероятности выживания для каждого наблюдения моего фрейма данных. Я применил модель Кокса к своим данным и применяю это:
predictSurvProb(object = model, newdata = new_data, times = c(200))
Приведенный выше код рассчитывает вероятности выживания через 200 дней для каждой строки new_data , используя результаты примерная модель .
Теперь проблема заключается в следующем. Я использую 200 дней для вычисления вероятностей, но эти 200 применяются ко всем наблюдениям. Но что произойдет, если я захочу использовать разные номера для каждой строки? Я попробовал это с для l oop, но у меня так много наблюдений, что это занимает много времени.
surv_probs <- data.frame()
times_for_each_observations = #vector with times for each observations in new_data I want to use
for(i in 1:nrow(new_data))
{
surv_probs_aux <- predictSurvProb(object = model, newdata = new_data[i,], times = c(times_for_each_observations[i]))
colnames(surv_probs_aux) <- c('prob')
surv_probs <- rbind(surv_probs, surv_probs_aux)
}
Может кто-нибудь помочь мне вычислить это в простой Кстати?
Я хотел бы сделать это, потому что все мои наблюдения имеют очень разные кривые выживания, и я хочу изменить значение времени для каждого из них.
Большое спасибо.