Displaying 1 result from an estimated 1 matches for "lme2a".
Did you mean:
lme2
2010 Sep 16
1
Help for an absolutely r-noob
...brary(lattice)
xyplot(resp ~ amp|id, dat, type=c("g","p","smooth"))
##### (3) MIXED-EFFECTS REGRESSION MODELS #####
library(nlme)
lme1 <- lme(resp ~ amp, dat, ~1|id, method="ML")
lme2 <- lme(resp ~ amp + I(amp^2), dat, ~1|id, method="ML")
lme2a <- lme(resp ~ amp + I(amp^2), dat, ~amp|id, method="ML")
lme3 <- lme(resp ~ amp + I(amp^2) + I(amp^3), dat, ~1|id, method="ML")
lme5 <- lme(resp ~ factor(amp), dat, ~1|id, method="ML")
anova(lme1, lme2)
anova(lme2, lme3, lme5)
summary(lme2)
##### (4) obs...