Я использую функцию cph в пакете RMS, чтобы оценить связь ВЕСА с результатом события. Я использую сплайн cubi c с 3 узлами.
Как я могу использовать модель для прогнозирования степени опасности для весов, которые не были включены в исходные данные? т.е. как я могу использовать модель для интерполяции или экстраполяции?
Вот пример:
set.seed(123)
library(rms)
WT <- rnorm(10, 30, 10)
EVENT <- sample(c(0,1), replace=TRUE, size=10)
TIME <- c(seq(1,10,1))
df <- as.data.frame(cbind(TIME,EVENT,WT))
fit <- cph(Surv(TIME, EVENT==1) ~ rcs(WT, 3), data = df)
fit
d = datadist(df)
options(datadist = 'd')
Predict(fit, WT) # this predicts hazard for WT included in the data
Как я могу использовать модель fit для прогнозирования хазрада на WT==70 для пример?
Я использую R studio.