forked from palday/lme4-convergence
-
Notifications
You must be signed in to change notification settings - Fork 0
/
models.R
21 lines (16 loc) · 851 Bytes
/
models.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
library(lme4)
library(lme4.0)
c. <- function(x) scale(x, center=TRUE, scale=FALSE)
reml <- FALSE
modeldata <- read.table("data.tab",header=T,
colClasses=c(mean="numeric",subj="factor",item="factor",
roi="factor",win="factor",
sdiff="numeric", dist="numeric",signdist="numeric",ambiguity="factor"))
modeldata <- subset(modeldata, roi == 'Left-Posterior')
modeldata.n400 <- subset(modeldata,win=="N400")
sdiff.new <- lme4::lmer(mean ~ ambiguity * c.(sdiff) + (1+c.(sdiff)|item) + (1+c.(sdiff)|subj), data=modeldata.n400, REML=reml)
print("New lme4")
print(summary(sdiff.new))
sdiff.old <- lme4.0::lmer(mean ~ ambiguity * c.(sdiff) + (1+c.(sdiff)|item) + (1+c.(sdiff)|subj), data=modeldata.n400, REML=reml)
print("Old lme4")
print(summary(sdiff.old))