From: Jon Olav Vik <jonovik_at_start.no>

Date: Sat 06 May 2006 - 22:21:09 EST

intervals(lme1) # confidence for random effect anova(lm1,lme1)

# models with response "NULL" removed because

*# response differs from model 1 in: anova.lmlist(object, ...)
*

anova(lme1,lm1)

logLik(lme1,REML=TRUE)

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.html Received on Sat May 06 22:26:25 2006

Date: Sat 06 May 2006 - 22:21:09 EST

library(nlme) summary(lm1 <- lm(travel~1,data=Rail)) # no random effect summary(lme1 <- lme(fixed=travel~1,random=~1|Rail,data=Rail)) # randomeffect

intervals(lme1) # confidence for random effect anova(lm1,lme1)

## Outputs warning message:

# models with response "NULL" removed because

anova(lme1,lm1)

## Output: Can I trust this?

# Model df AIC BIC logLik Test L.Ratio p-value

# lme1 1 3 128.1770 130.6766 -61.08850# lm1 2 2 162.6815 164.3479 -79.34075 1 vs 2 36.50451 <.0001

## Various log likelihoods: logLik(lm1,REML=FALSE) logLik(lm1,REML=TRUE) logLik(lme1,REML=FALSE)

logLik(lme1,REML=TRUE)

Any help is highly appreciated.

Best regards,

Jon Olav Vik

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.html Received on Sat May 06 22:26:25 2006

Archive maintained by Robert King, hosted by
the discipline of
statistics at the
University of Newcastle,
Australia.

Archive generated by hypermail 2.1.8, at Mon 08 May 2006 - 08:10:02 EST.

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