[R] predict.gls choking on levels of factor
I've got a gls formula that's a mix of continuous and ordered variables. I wanted to use gls because I wanted to use the varIdent structure. Anyway, attempts to use predict.gls choke with the error that the levels I use are not allowed for one of them -- the first one alphabetically, so I'd guess the second would have the same problem. So I have three linked questions -- answering any of the three would set my mind at ease. (1) is it possible to make predict.gls friendlier to my model structure? (2) Is it possible to get the varIdent into lm or ols from Design? (3) is Yee's smartpred package, last updated for R 1.7, reliable for current R? Thanks! __ R-help@r-project.org mailing list 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.
Re: [R] predict.gls choking on levels of factor
Mitchell Maltenfort mmalten at gmail.com writes: I've got a gls formula that's a mix of continuous and ordered variables. I wanted to use gls because I wanted to use the varIdent structure. Anyway, attempts to use predict.gls choke with the error that the levels I use are not allowed for one of them -- the first one alphabetically, so I'd guess the second would have the same problem. Can you give us the precise error message? [snip] (1) is it possible to make predict.gls friendlier to my model structure? Need a reproducible example (see below) (2) Is it possible to get the varIdent into lm or ols from Design? No. (3) is Yee's smartpred package, last updated for R 1.7, reliable for current R? Can we have a reproducible example please? Ben Bolker __ R-help@r-project.org mailing list 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.
