[R] predict.gls choking on levels of factor

Ben Bolker bbolker at gmail.com
Sat May 21 17:39:10 CEST 2011


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



More information about the R-help mailing list