Jessica Myers
2011-Apr-27 16:16 UTC
[R] treatment of factors and errors in ridge() function with coxph
I am trying to fit a large Cox model with many predictors. Because there are many predictors, I would like to use the ridge() function to get penalized ml estimates for all coefficients. The problems are that: 1. When I include a factor (like race) in the ridge() function, dummy variables are not created. The resulting model has a single coefficient for the race variable, and I have no idea what is the interpretation of this coefficient. Obviously, I can create the dummy variables myself, but it would be better not to. 2. When I include a large number of predictors in ridge() [> ~ 35], coxph does not complete and I get the error: Error in model.matrix.default(Terms, m) : model frame and formula mismatch in model.matrix() Using a smaller subset of the same covariates, I don't get this error and it doesn't seem to matter which subset of covariates I use. Any help with these problems would be greatly appreciated! Thanks, Jessica Myers Instructor in Medicine Division of Pharmacoepidemiology Brigham and Women's Hospital The information in this e-mail is intended only for the ...{{dropped:7}}