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.nussb...@env.ethz.ch
www.step.ethz.ch
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
