require(lme4) bySubj <- read.csv("bySubj.csv") byItem <- read.csv("byItem.csv") model.bySubj <- lmer(elNP2 ~ (t1+t2)*lang*cons*group*conflict + (t1+t2||ID) + (t1+t2||ID:lang:cons), bySubj, weights=1/wtsNP2) model.byItem <- lmer(elNP2 ~ (t1+t2)*lang*cons*group + (t1+t2||Item) + (t1+t2||Item:lang:cons:group), bySubj, weights=1/wtsNP2)