library(languageR)
affixes.pr=prcomp(affixProductivity[, 1:(ncol(affixProductivity)-3)])
summary(affixes.pr)
affixes.pr$rotation
biplot(affixes.pr)

ilm=read.table("http://www.tlu.ee/~jaagup/andmed/ilm/harkutund.txt", header=TRUE, sep=",")
ilmtervik=ilm[complete.cases(ilm), ]
ilmskaleeritud=scale(ilmtervik[, 4:11])
k=prcomp(ilmskaleeritud)
summary(k)
k$rotation
barplot(k$sdev**2)
plot(k)
biplot(k)

ilmskaleeritud=scale(ilmtervik[(ilmtervik$Kuu==7) & (ilmtervik$Paev==3), 5:11])
k=prcomp(ilmskaleeritud)
summary(k)
k$rotation
biplot(k)


affixes.fac=factanal(affixProductivity[, 1:27], factors=3)
affixes.fac
loadings(affixes.fac)
plot(loadings(affixes.fac))

affixes.fac=factanal(affixProductivity[, 1:27], factors=4, rotation="promax")
affixes.fac
plot(loadings(affixes.fac))


ilm=read.table("http://www.tlu.ee/~jaagup/andmed/ilm/harkutund.txt", header=TRUE, sep=",")
ilmtervik=ilm[complete.cases(ilm), ]
ilmskaleeritud=scale(ilmtervik[, 4:11])
t1=factanal(ilmskaleeritud, factors=3, rotation="promax")
t1


oldFrench
oldFrenchMeta
oldFrench.ca=corres.fnc(oldFrench)
summary(oldFrench.ca)
plot(oldFrench.ca)
plot(oldFrench.ca, rlabels=oldFrenchMeta$Genre)


variationLijk
chisq.test(variationLijk)
k=corres.fnc(variationLijk)
plot(k)


ilm=read.table("http://www.tlu.ee/~jaagup/andmed/ilm/harkutund.txt", header=TRUE, sep=",")
ilmtervik=ilm[complete.cases(ilm), ]

k=corres.fnc(ilmtervik[, 4:11])
plot(k, rlabels = ilmtervik$Kuu, rcex = 0.5)

k=corres.fnc(ilmtervik[ilmtervik$Kuu==7, 4:11])
plot(k, rlabels = substr(ilmtervik$Kell, 1, 2), rcex = 0.3)