使用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
答案 0 :(得分:1)
我认为您可以使用as.matrix
:
pl <- as.matrix(model$y)
as.data.frame(pl)