我用“gbm”R 包分析了我的数据。我的数据基于一项队列研究。因此,我根据“coxph”结果运行了“gbm”模型。
建好一个模型后,我想看看这个模型如何能很好地预测。另一方面,就像下面的代码一样,预测的值是负的。所以,我很难理解这种现象。请让我知道如何解释这个值。
这是我的代码。
install.packages("survival")
install.packages("randomForestSRC")
install.packages("gbm")
library(survival)
library(randomForestSRC)
library(gbm)
data(pbc, package="randomForestSRC")
data <- na.omit(pbc)
exposure <- names(data[, names(data.model) !=c("days", "status")])
formula <- as.formula(paste("Surv(days, status)~", paste(exposure, collapse="+")))
set.seed(123)
ex <- gbm(Surv(days, status)~.,
data=data,
distribution="coxph",
cv.folds=5,
shrinkage=.01,
n.trees=1000)
set.seed(123)
pred <- predict(ex, n.trees=1000, type="response")