Re: [R] R^2 from lme function

2007-05-14 Thread Douglas Bates
On 5/14/07, Martin Henry H. Stevens <[EMAIL PROTECTED]> wrote: > Hi Cleber, > By "full" I simply meant "not REML." the function assumes that the > fixed effects were estimated using REML criteria, and using update() > simply changes that to ML. If the model was fit originally with ML, > it shouldn

Re: [R] R^2 from lme function

2007-05-14 Thread Martin Henry H. Stevens
Hi Cleber, By "full" I simply meant "not REML." the function assumes that the fixed effects were estimated using REML criteria, and using update() simply changes that to ML. If the model was fit originally with ML, it shouldn't make any difference. I am reasonably sure that it should not mat

Re: [R] R^2 from lme function

2007-05-14 Thread Cleber Borges
Hi Martin, many thanks for your tip! but,{ :-( } what it 'full MLE' ? how to calculate? it is a saturated model??? and it is valid for 'no-intercept model? Many thanks again... Cleber > Hi Cleber, > I have been using this function I wrote for lme

Re: [R] R^2 from lme function

2007-05-14 Thread Martin Henry H. Stevens
Hi Cleber, I have been using this function I wrote for lmer output. It should be easy to convert to lme. As with everything, buyer beware. Note that it requires (full) maximum likelihood estimates. Rsq <- function(reml.mod) { ## Based on ## N. J. D. Nagelkerke. A note on a general defini

[R] R^2 from lme function

2007-05-14 Thread Cleber Borges
Hello allR How to access R^2 from lme object? or how to calculate it? ( one detail: my model do not have a intercept ) thanks in advanced Cleber ___ Experimente já e veja as novidades.