#You will have to have languageR and Design installed to run this code.

library(languageR)
library(Design)

lrm = lrm(RealizationOfRec~AnimacyOfRec+AnimacyOfTheme+LengthOfTheme,data=verbs)
lrm
glm = glm(RealizationOfRec~AnimacyOfRec+AnimacyOfTheme+LengthOfTheme,data=verbs,family="binomial")
summary(glm)
lmer = lmer(RealizationOfRec~AnimacyOfRec+AnimacyOfTheme+LengthOfTheme+(1|Verb),data=verbs,family="binomial")
lmer

#LRM Stepwise

anova(lrm)

#LRM Comparison

lrm = lrm(RealizationOfRec~AnimacyOfRec+AnimacyOfTheme+LengthOfTheme,data=verbs)
lrm.length = lrm(RealizationOfRec~LengthOfTheme,data=verbs)
lrm.animac = lrm(RealizationOfRec~AnimacyOfRec+AnimacyOfTheme,data=verbs)

dchisq(deviance(lrm.length)[2]-deviance(lrm)[2],lrm$stat[4]-lrm.length$stat[4])
dchisq(deviance(lrm.animac)[2]-deviance(lrm)[2],lrm$stat[4]-lrm.animac$stat[4])

#LMER Comparison

lmer = lmer(RealizationOfRec~AnimacyOfRec+AnimacyOfTheme+LengthOfTheme+(1|Verb),data=verbs,family="binomial")
lmer.length = lmer(RealizationOfRec~LengthOfTheme+(1|Verb),data=verbs,family="binomial")
lmer.animac = lmer(RealizationOfRec~AnimacyOfRec+AnimacyOfTheme+(1|Verb),data=verbs,family="binomial")

anova(lmer,lmer.length)
anova(lmer,lmer.animac)

d=read.csv(file="PluralComparisonMontreal.csv",header=T)

head(d)

#Your turn! Go nuts!

lmer(resp~(1|as.factor(d$subj))+(1|as.factor(d$item))+(1|as.factor(d$order))
+MatuRuys
+catMean
+probMean
,data=d,family="binomial")