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

###Section 10.2


##Main text

install.packages(c("ggplot2", "corrgram") #if you haven't installed these packages yet

library(Rling); library(ggplot2); library(corrgram)

data(ditr)

head(ditr)

a <- ditr$Freq_VC

b <- ditr$Freq_V - a

c <- 667 - a

d <- 104162 - (a + b + c)

head(c)

aExp <- (a + b)*(a + c)/(a + b + c + d)

head(aExp)

Attr <- 100*a/(a + c)

Rel <- 100*a/(a + b)

names(Attr) <- rownames(ditr)

head(Attr)

sort(Attr, decreasing = TRUE)[1:5]

names(Rel) <- rownames(ditr)

sort(Rel, decreasing = TRUE)[1:5]

plot(Attr, Rel, type = "n", main = "Attraction and Reliance scores of verbs in Russian Ditransitive Cx")

text(Attr, Rel, rownames(ditr), cex = 0.7)

plot(sort(Attr, decreasing = TRUE), type = "l", main = "Attraction", ylab = "Attraction, in %")

plot(sort(Rel, decreasing = TRUE), type = "l", main = "Reliance", ylab = "Reliance, in %")

dP.cueVerb <- a/(a + b) - c/(c + d)

dP.cueCx <- a/(a + c) - b/(b + d)

pvF <- pv.Fisher.collostr(a, b, c, d)

round(head(pvF), 3)

logpvF <- ifelse(a < aExp, log10(pvF), -log10(pvF))

round(head(logpvF), 3)

-log10(0.05)

LL <- LL.collostr(a, b, c, d)

LL1 <- ifelse(a < aExp, -LL, LL)

random <- sample(0:100, length(a), replace = TRUE)

head(random)

assoc <- cbind(Attr, Rel, dP.cueCx, dP.cueVerb, logpvF, LL1, random)
assoc.cor <- cor(assoc)

round(assoc.cor, 3)

corrgram(assoc, order = TRUE, lower.panel = panel.shade, upper.panel = panel.pie)


corrgram(assoc, order = TRUE, lower.panel = panel.ellipse, upper.panel = panel.pts)

assoc.cor1 <- cor(assoc, method = "kendall")

round(assoc.cor1, 3)


##Boxes with additional information

ggplot(data.frame(Attr, Rel), aes(x = Attr, y = Rel)) + geom_text(label = names(Attr), size = 3)