################################################### ### Logistic regression: BresDative data ########## ################################################### # Open data file. dat <- read.table("BresDative.txt", header=T) head(dat) SwitchDat <- subset(dat,mod=="switchboard") ####################################################### ### Logistic regression model of dative alternation ## ###################################################### dat.glm <- glm(real~class+accrec+accth+prorec+proth+defrec+defth+ animrec+ldiff, family=binomial, data=SwitchDat) summary(dat.glm) #################################################################### ### Mixed effects logistic regression model of dative alternation ## #################################################################### library(lmerTest) dat.glmer01 <- glmer(real~ class+accrec+accth+prorec+proth+defrec+defth+ animrec+ldiff +(1|vsense), family=binomial, data=SwitchDat) # model failed to converge # use a different optimizer dat.glmer02 <- glmer(real~ class+accrec+accth+prorec+proth+defrec+defth+ animrec+ldiff +(1|vsense), family=binomial, data=SwitchDat, control=glmerControl(optimizer="bobyqa", optCtrl=list(maxfun=100000))) summary(dat.glmer02)