data(hockey)
## fit a regularization path for 20 significant players
## team effects are unpenalized
fit <- gamlr(covars=cbind(Xteam,Xplayer), response=Ygoals,
free=1:30, verb=1, family="binomial", scale=FALSE,
stoparg=list(rule="maxvar", val=20))
plot(fit)
summary(fit)
Run the code above in your browser using DataLab