From: Beatrijs Moerkerke <Beatrijs.Moerkerke_at_UGent.be>

Date: Wed 04 May 2005 - 20:03:59 EST

Date: Wed 04 May 2005 - 20:03:59 EST

The code I use for calculating the likelihood ratio test statistics is as follows:

a1<-summary(lme(y~x,random=~1|gr,method="ML"))$logLik
a2<-logLik(lm(y~x))

(-2*(a2-a1))

I don't know how I can simulate the null distribution in R using lme.

Kind regards,

Beatrijs Moerkerke

-- Beatrijs Moerkerke Department of Applied Mathematics and Computer Science Ghent University Krijgslaan 281 - S9 B-9000 GENT Tel: +32-(0)9-264.47.56 Fax: +32-(0)9-264.49.95 E-mail: Beatrijs.Moerkerke@UGent.be ______________________________________________ R-help@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.htmlReceived on Wed May 04 22:10:22 2005

*
This archive was generated by hypermail 2.1.8
: Fri 03 Mar 2006 - 03:31:35 EST
*