[R] mboost: Proportional odds boosting model - how to specify the offset?
Nussbaum Madlene
madlene.nussbaum at env.ethz.ch
Fri Mar 13 14:17:58 CET 2015
Dear R team
The package mboost allows for boosting of proportional odds models.
However, I would like to include an offset for every observation. This
produces an error - no matter how I put the offset (as response
probabilities or as response link).
Fitting gamboost-models with offset works satisfactory with family =
Gaussian() or Multinomial().
Questions:
1) How do I need to specify the offset with family = PropOdds()?
2) Where in the mboost-object do I find the Theta's (response category
dependent intercept)?
# --- minimal example with iris data ---
library(MASS)
library(mboost)
data(iris)
iris$Species <- factor(iris$Species, ordered = T)
p.iris <- polr(Species ~ Sepal.Length, data = iris)
mlp <- gamboost(Species ~ bols(Sepal.Length) + bols(Sepal.Width),
data = iris, family = PropOdds(),
offset = fitted(p.iris) )
Error in tmp[[i]] : subscript out of bounds
Thank you
M. Nussbaum
--
ETH Zürich
Madlene Nussbaum
Institut für Terrestrische Ökosysteme
Boden- und Terrestrische Umweltphysik
CHN E 37.2
Universitätstrasse 16
CH-8092 Zürich
Telefon + 44 632 73 21
Mobile + 79 761 34 66
madlene.nussbaum at env.ethz.ch
www.step.ethz.ch
More information about the R-help
mailing list