2
$\begingroup$

I can not put my company data online, but I can provide a reproducible example here.

We're modelling Insurance's frequency using Poisson distribution with exposure as offset.

Here in this example, we want to model the number of claim Claims ($y_i$) with exposure Holders ($e_i$)

In the traditional GLM model, we can dirrectly model $y_i$ and put $e_i$ in the offset term. This option is not available in xgboost. So the alternative is to model the rate $\frac{y_i}{e_i}$, and put $e_i$ as a the weight term (reference)

When I do that with a lot of iteractions, the results are coherent ($\sum y_i = \sum \hat{y_i}$). But it is not the case when nrounds = 5. I think that the equation $\sum y_i = \sum \hat{y_i}$ must be satisfied after the very first iteration.

The following code is an extreme example for the sake of reproducibility. In my real case, I performed a CV on the training set (optimizing MAE), I obtained nrounds = 1200, training MAE = testing MAE. Then I re-run a xgboost on the whole data set with 1200 iteration, I see that $\sum y_i \ne \sum \hat{y_i}$ by a large distance, this doesn't make sense, or am I missing something?

So my questions are:

  1. Am I correctly specify parameters for Poisson regression with offset in xgboost?
  2. Why such a high bias at the first iterations?
  3. Why after tuning nrounds using xgb.cv, I still have high bias?

Here is the graphics plotting the ratio $\frac{\sum \hat{y_i}}{\sum y_i}$ by nrounds

enter image description here

Code edited after the comment of @JonnyLomond

library(MASS)
library(caret)
library(xgboost)
library(dplyr)

#-------- load data --------#
data(Insurance)

#-------- data preparation --------#

#small adjustments
Insurance$rate = with(Insurance, Claims/Holders)
temp<-dplyr::select(Insurance,District, Group, Age, rate)
temp2= dummyVars(rate ~ ., data = temp, fullRank = TRUE) %>% predict(temp)

#create xgb matrix
xgbMatrix <- xgb.DMatrix(as.matrix(temp2), 
                         label = Insurance$Claims)
setinfo(xgbMatrix, "base_margin",log(Insurance$Holders))


#-------------------------------------------#
#      First model with small nround
#-------------------------------------------#
bst.1 = xgboost(data    = xgbMatrix, 
                objective ='count:poisson', 
                nrounds   = 5)


pred.1 = predict(bst.1, xgbMatrix)
sum(Insurance$Claims) #3151
sum(pred.1) #12650.8 fails



#-------------------------------------------#
#      Second model with more iteractions
#-------------------------------------------#

bst.2 = xgboost(data    = xgbMatrix, 
                objective = 'count:poisson', 
                nrounds   = 100)


pred.2 = predict(bst.2, xgbMatrix)
sum(Insurance$Claims) #3151
sum(pred.2) #same
$\endgroup$
  • 3
    $\begingroup$ Why would you expect any algorithm to converge at first iteration..? $\endgroup$ – Tim Jan 25 '18 at 11:09
  • $\begingroup$ @Tim : I expect the nrounds = 1 to be equivalent to fitting a tree on the data, so sum of observations should be equal to sum of predictions. Am I missing something? $\endgroup$ – Metariat Jan 25 '18 at 11:10
  • $\begingroup$ perhaps "converge" is not the right word, but I hope the question is clear. $\endgroup$ – Metariat Jan 26 '18 at 8:54
  • 2
    $\begingroup$ If fitting single tree worked out-of-the-box, we wouldn't use random forests and boosting at all, don't you think? $\endgroup$ – Tim Jan 26 '18 at 8:58
  • 1
    $\begingroup$ Why do 2 people vote to close this question? Can you give me the reason? $\endgroup$ – Metariat Jan 26 '18 at 17:16
4
$\begingroup$

First a few technical things:

  1. You can use an offset in xgboost for Poisson regression, by setting the base_margin value in the xgb.DMatrix object.

  2. You will not get the same results with your above code as if you use the base_margin term. (You get the same results in a GLM, but this is not a GLM. I think the weight term in xgboost means something different.)

For your question:

Sum of predictions will not equal sum of observations after a small number of rounds for several reasons:

  1. Xgboost is regularizing predictions in the nodes (shrinking them toward 0). This will happen by default

  2. Xgboost is scaling predictions from each tree (by a positive number less than 1). This will happen by default

  3. Xgboost is randomly sampling rows every round, not fitting on the whole data set. I think this will not happen by default.

  4. Some other behaviors.

Basically: It is not true that fitting one round in xgboost is the same as fitting a basic decision tree in the usual way.

$\endgroup$
  • $\begingroup$ The first point 2: I have the same doubt as well, can you give more details? For the rest: I completely agree with you, but those arguments only explain the fact that the ratio $\frac{\sum \hat{y_i}}{\sum y_i} \ne 1$ but should be close to 1. I have this ratio equal to 1.5 after 1200 iterations. $\endgroup$ – Metariat Jan 26 '18 at 15:39
  • 2
    $\begingroup$ I dont think the result is even true for a single unregularized tree. The op is generalizing results about glms with canonical links to domains where they do not apply. The trees in gradient boosting are not even fit to the response. $\endgroup$ – Matthew Drury Jan 26 '18 at 15:40
  • $\begingroup$ @MatthewDrury thanks for you suggestions. But I still don't get it why with nrounds = 1200 obtained by minimizing CV MAE, this ratio is still 1.5? $\endgroup$ – Metariat Jan 26 '18 at 16:17
  • $\begingroup$ Are you adjusting your predictions for exposure? If not, what is your average exposure? $\endgroup$ – Matthew Drury Jan 26 '18 at 16:20
  • $\begingroup$ Yes, I adjust my predictions for exposure, when I try to increase nrounds to 2000, the equation is satisfied, but not with nrounds = 1200 found by CV. My average exposure is about 0.4 $\endgroup$ – Metariat Jan 29 '18 at 8:59

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.