library(lme4)
load("/home/bates/afs/src/Rlibs/AppLong/data/fev1.rda")
fev1$loght <- log(fev1$ht)
fev1$logbht <- log(fev1$baseht)
(fm1 <- lmer(logfev1 ~ age + loght + baseage + logbht + (age|id), fev1))
(fm1a <- lmer(logfev1 ~ age + loght + baseage + logbht + (age|id), fev1, subset = id != 197))
deviance(fm1a)
(fm2 <- lmer(logfev1 ~ age + loght + baseage + logbht + (loght|id), fev1))
(fm2a <- lmer(logfev1 ~ age + loght + baseage + logbht + (loght|id), fev1, subset = id != 197))
deviance(fm2a)
(fm3 <- lmer(logfev1 ~ age + loght + baseage + logbht + (age + loght|id), fev1, verb = 1))
(fm3a <- lmer(logfev1 ~ age + loght + baseage + logbht + (age + loght|id), fev1, verb = 1, subset = id != 197))
deviance(fm3a)