####################################################################################
#This is a file with executable R code for chapter 19 of Natalia Levshina's (2015) 
#How to Do Linguistics with R. Amsterdam/Philadelphia: John Benjamins. 
####################################################################################


###Section 19.2

##Main text

install.packages(c("vcd", "ca", "rgl"))

library(Rling); library(vcd); library(ca); library(rgl)

data(colreg)

colreg

mosaicplot(colreg, las = 2, shade = TRUE, main = "Register variation of BNC")

ca.bc <- ca(colreg)

summary(ca.bc)

plot(ca.bc)

plot3d(ca.bc, labels = c(1,1))


###Section 19.3

##Main text


library(Rling); library(FactoMineR); library(ca); library(rms)	

data(chairs)

str(chairs)

swivelRoll <- xtabs(~ chairs$Swivel + chairs$Roll)

swivelRoll

chisq.test(swivelRoll)

chairs.ca <- MCA(chairs[, -c(1:3)], graph = FALSE)

plot(chairs.ca, cex = 0.7, col.var = "black", col.ind = "grey")

dimdesc(chairs.ca)

chairs.ca1 <- MCA(chairs[, -c(1:2)], quali.sup = 1, graph = FALSE)

plot(chairs.ca1, invis = "ind", col.var = "darkgrey", col.quali.sup = "black")

plotellipses(chairs.ca1, keepvar = 1, label = "quali")

plotellipses(chairs.ca1, means = FALSE, keepvar = 1, label = "quali")

chairs.ca$eig

chairs.ca2 <- mjca(chairs[, -c(1:3)])

summary(chairs.ca2)

cor(chairs.ca$var$coord[, 1], chairs.ca2$colcoord[, 1])

dim1 <- chairs.ca$ind$coord[, 1]

dim2 <- chairs.ca$ind$coord[, 2]

m <- lrm(chairs$Category ~ dim1 + dim2)

m