data(dative)
# analysis with CART tree
library(rpart)
# ---- initial tree
dative.rp = rpart(RealizationOfRecipient ~ .,
data = dative[ ,-c(1, 3)]) # exclude the columns with subjects, verbs
plot(dative.rp, compress = TRUE, branch = 1, margin = 0.1)
text(dative.rp, use.n = TRUE, pretty = 0)
# ---- pruning the initial tree
plotcp(dative.rp)
dative.rp1 = prune(dative.rp, cp = 0.041)
plot(dative.rp1, compress = TRUE, branch = 1, margin = 0.1)
text(dative.rp1, use.n = TRUE, pretty = 0)
# analysis with logistic regression
# ---- logistic regression with the Design package
library(Design)
dative.dd = datadist(dative)
options(datadist = 'dative.dd')
dative.lrm = lrm(RealizationOfRecipient ~
AccessOfTheme + AccessOfRec + LengthOfRecipient + AnimacyOfRec +
AnimacyOfTheme + PronomOfTheme + DefinOfTheme + LengthOfTheme+
SemanticClass + Modality, data = dative)
anova(dative.lrm)
par(mfrow = c(4,3))
plot(dative.lrm)
par(mfrow = c(1,1))
# ---- mixed-effects logistic regression with the lme4 package
library(lme4)
dative.lmer = lmer(RealizationOfRecipient ~ AccessOfTheme +
AccessOfRec + LengthOfRecipient + AnimacyOfRec +
AnimacyOfTheme + PronomOfTheme + DefinOfTheme + LengthOfTheme +
SemanticClass + Modality + (1|Verb),
data = dative, family = "binomial")
print(dative.lmer, corr=FALSE)
Run the code above in your browser using DataLab