Как отменить объект выживания (вернуть фрейм данных)
После генерации модели с использованием cph()
от rms
звонит model$y
вернет объект выживания. Есть ли функция, которая будет "отменять" объект выживания и возвращать фрейм данных?
Я хотел бы иметь возможность использовать объект выживания в качестве аргумента для функции, которую я пишу, но мне также нужны данные ответа. Я пытаюсь избежать использования данных в качестве аргумента и создания модели внутри функции.
Минимальный рабочий пример представлен ниже:
library(rms)
# generate data
time <- c(82, 73, 89, 79, 72, 87, 103, 83, 100, 79)
event <- c(0, 0, 1, 0, 1, 0, 0, 0, 1, 1)
covar <- c(15, 11, 11, 20, 12, 13, 10, 11, 10, 14)
df <- data.frame(time, event, covar)
# Cox model
dd <- datadist(df)
options(datadist = 'dd')
model <- cph(Surv(time, event) ~ covar, x=TRUE, y=TRUE, surv=TRUE, data=df)
# returns a survival object
model$y
# what I want is a data frame
want <- data.frame(time, event)
want
1 ответ
Я думаю, что вы можете получить то, что вы хотите, используя as.matrix
:
pl <- as.matrix(model$y)
as.data.frame(pl)