[R-sig-ME] using weights=varIdent in lme: how to change reference level of grouping factor?

[Emma]

Ben Bolker <bbolker at ...> writes:

> 
>   I have spent some time delving through the guts of lme()
> but haven't figured it out yet.  I'm somewhat baffled,
> as the default order of factor() should be Female, Male in
> any case, so I don't see where the order could be getting
> messed up.
> 
> 


Hi Ben,
Thanks for having a look at this.
I think lme is using the category with the largest number of observations as
the reference level.
I have no idea how to change this.
Emma