Significado dos termos de saída no pacote gbm?

15

Estou usando o pacote gbm para classificação. Como esperado, os resultados são bons. Mas estou tentando entender a saída do classificador. Existem cinco termos na saída.

`Iter   TrainDeviance   ValidDeviance   StepSize   Improve`

Alguém poderia explicar o significado de cada termo, especialmente o significado de Melhorar .

Yoanh27
fonte

Respostas:

11

Você deve achar que estão relacionados à determinação do melhor valor para o número de funções básicas - ou seja, iterações - ou seja, número de árvores no modelo aditivo. Não consigo encontrar documentação descrevendo exatamente o que são, mas aqui está o meu melhor palpite e talvez alguém possa comentar.

Siga o seguinte no manual:

library(gbm)
# A least squares regression example
# create some data
N <- 1000
X1 <- runif(N)
X2 <- 2*runif(N)

X3 <- ordered(sample(letters[1:4],N,replace=TRUE),levels=letters[4:1])
X4 <- factor(sample(letters[1:6],N,replace=TRUE))
X5 <- factor(sample(letters[1:3],N,replace=TRUE))
X6 <- 3*runif(N)
mu <- c(-1,0,1,2)[as.numeric(X3)]
SNR <- 10 # signal-to-noise ratio
Y <- X1**1.5 + 2 * (X2**.5) + mu
sigma <- sqrt(var(Y)/SNR)
Y <- Y + rnorm(N,0,sigma)
# introduce some missing values
X1[sample(1:N,size=500)] <- NA
X4[sample(1:N,size=300)] <- NA
data <- data.frame(Y=Y,X1=X1,X2=X2,X3=X3,X4=X4,X5=X5,X6=X6)
# fit initial model
gbm1 <- gbm(Y~X1+X2+X3+X4+X5+X6, # formula
data=data, # dataset
var.monotone=c(0,0,0,0,0,0), # -1: monotone decrease,
# +1: monotone increase,
# 0: no monotone restrictions
distribution="gaussian", # bernoulli, adaboost, gaussian,
# poisson, coxph, and quantile available
n.trees=3000, # number of trees
shrinkage=0.005, # shrinkage or learning rate,
# 0.001 to 0.1 usually work
interaction.depth=3, # 1: additive model, 2: two-way interactions, etc.
bag.fraction = 0.5, # subsampling fraction, 0.5 is probably best
train.fraction = 0.5, # fraction of data for training,
# first train.fraction*N used for training
n.minobsinnode = 10, # minimum total weight needed in each node
cv.folds = 5, # do 5-fold cross-validation
keep.data=TRUE, # keep a copy of the dataset with the object
verbose=TRUE) # print out progress

O número de iterações ( Iter) é 3000, que é o número de árvores selecionadas para serem construídas (1 a 3000, embora nem todas sejam mostradas). O processo completo é repetido 5 vezes a propósito, porque selecionamos cv.folds = 5.

StepSize é a taxa de encolhimento ou aprendizado selecionada (0,005 aqui).

Acredito que Improveseja a redução no desvio (função de perda) adicionando outra árvore e é calculada usando os registros out-of-bag (OOB) (observe que não será calculado se bag.fraction não for <1).

Então, para cada iteração, TrainDeviance ValidDevianceé o valor da função de perda nos dados de treinamento e retém os dados (um único conjunto de retenções). O ValidDeviance não será calculado se train.fractionnão for <1.

Você já viu isso descrevendo os três tipos de métodos para determinar o número ideal de árvores?

B_Miner
fonte