-1

我用“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")
4

1 回答 1

1

阅读?predict.gbm帮助页面,尤其是参数type。默认情况下,预测是在链接规模上的。

于 2018-09-02T16:57:20.090 回答