### R code from vignette source 'effects.Rnw'

###################################################
### code chunk number 1: effects.Rnw:5-6
###################################################
options(show.signif.stars=FALSE,digits=4,width=70)


###################################################
### code chunk number 2: effects.Rnw:10-16
###################################################
load(url("http://www.stat.umn.edu/~sandy/courses/8053/Data/Achiv.Rda"))
meanses <- with(Achiv, tapply(ses, school, mean))
head(meanses)
length(meanses)
Achiv$meanses <- meanses[as.character(Achiv$school)]
Achiv$indivses <- with(Achiv, ses - meanses)


###################################################
### code chunk number 3: effects.Rnw:19-24
###################################################
library(effects)
library(lme4)
m5 <- lmer(mathach ~ indivses*sector + indivses*meanses + (indivses | school),
        data=Achiv)
print(m5, corr=FALSE)


###################################################
### code chunk number 4: eff1
###################################################
print(plot(effect("sector", m5) ))


###################################################
### code chunk number 5: eff2
###################################################
print(plot(effect("indivses*sector", m5)))


###################################################
### code chunk number 6: eff3
###################################################
print(plot(effect("indivses:meanses", m5)))


###################################################
### code chunk number 7: eff4
###################################################
print(plot(effect("indivses:sector:meanses", m5), factor.names=FALSE))