Re: [R] R^2 from lme function

From: Cleber Borges <klebyn_at_yahoo.com.br>
Date: Mon, 14 May 2007 12:55:17 -0300

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 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 definition
> ## of the coefficient of determination. Biometrika, 78:691–692, 1991.
> ml.mod <- update(reml.mod, method="ML")
> l.B <- logLik(ml.mod)
> l.0 <- logLik( lm(ml.mod_at_y ~ 1) )
> Rsq <- 1 - exp( - ( 2/length(ml.mod_at_y) ) * (l.B - l.0) )
> Rsq[1]
> }
>
> Hank
>
>
>
>
>> 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.



R-help_at_stat.math.ethz.ch 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. Received on Mon 14 May 2007 - 16:02:26 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Mon 14 May 2007 - 17:31:33 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.