library(languageR)


4.1
pt(-2,10)
2*pt(-2,10)
2*(1-pt(2,10))
2*(1-pt(abs(-2),10))
2*(1-pt(abs(2),10))

4.1.1
plot(density(ver$Frequency))
ver$Frequency=log(ver$Frequency)
plot(density(ver$Frequency))

qqnorm(rnorm(length(ver$Frequency),4,3))
abline(v=qnorm(0.025),col="grey")
abline(h=qnorm(0.025,4,3),col="grey")

qqnorm(ver$Frequency)

shapiro.test(ver$Frequency)

ks.test(ver$Frequency,"pnorm",mean(ver$Frequency),sd(ver$Frequency))

ks.test(jitter(ver$Frequency),"pnorm",mean(ver$Frequency),sd(ver$Frequency))

ver$Frequency[1:5]
jitter(ver$Frequency[1:5])
ks.test(jitter(ver$Frequency),"pnorm",mean(ver$Frequency),sd(ver$Frequency))

intro=c(75,68,45,40,39,39,38,33,24,24)
names(intro)=c("the","to","of","you","is","a","and","in","that","data")
intro

chisq.test(intro)

x=c(37,21,26,30,23,26,41,26,37,33)
chisq.test(x)


4.1.2
meanLengthN=mean(durationsOnt$DurationPrefixNasal)
meanLengthN

t.test(durationsOnt$DurationPrefixNasal,mu=0.053)

2*(1-pt(abs(-1.5038),101))

mean(durationsOnt$DurationPrefixPlosive)

shapiro.test(durationsOnt$DurationPrefixPlosive)

wilcox.test(durationsOnt$DurationPrefixPlosive,mu=0.044)

4.2
4.2.1
ver$Frequency=log(ver$Frequency)
ver.transp=ver[ver$SemanticClass=="transparent",]$Frequency
ver.opaque=ver[ver$SemanticClass=="opaque",]$Frequency

ver.transp.d = density(ver.transp)
ver.opaque.d = density(ver.opaque)
xlimit = range(ver.transp.d$x, ver.opaque.d$x)
ylimit = range(ver.transp.d$y, ver.opaque.d$y)
plot(ver.transp.d, lty = 1, col = "black", xlab = "frequency", ylab = "density", xlim = xlimit,
ylim = ylimit, main = "")
lines(ver.opaque.d, col = "darkgrey")

ks.test(jitter(ver.transp), jitter(ver.opaque))

4.2.2
library(lattice)
bwplot(Frequency ~ Class | Complex, data = ratings)

simplex = ratings[ratings$Complex == "simplex",]
freqAnimals = simplex[simplex$Class == "animal",]$Frequency
freqPlants = simplex[simplex$Class == "plant",]$Frequency
t.test(freqAnimals, freqPlants)

t.test(simplex[simplex$Class=="animal",]$Frequency,simplex[simplex$Class=="plant",]$Frequency,conf.level=0.99)

wilcox.test(ver.opaque,ver.transp)

tapply(verbs$LengthOfTheme,verbs$AnimacyOfRec,mean)

t.test(LengthOfTheme ~ AnimacyOfRec, data = verbs)

4.2.3
x = rnorm(50, mean = 0, sd = 2)
y = rnorm(30, mean = 1, sd = 1)

var.test(x, y)

var(x)/var(y)

2 * (1 - pf(var(x)/var(y), 49, 29))