Tony,
gls() fits the model by ML or REML (the default), and the parameters in
the correlation structure (e.g., lambda) are estimated if needed or
specified.
Cheers,
Emmanuel
Anthony R Ives wrote on 21/04/2011 17:45:
Emmanuel,
I don't know gls, so I'd have to look at this. Technically,
Scott,
I would love it if somebody would. If I were doing it now, I'd
simultaneously use an OU or maybe lambda transform. Several people
have code that does, for example, regression while assuming residual
variation has some non-Brownian phylogenetic structure (e.g., Lavins
et al.