# 3. TYÖARKKI
#
# Yksisuuntainen varianssianalyysi ja qqnorm-kuvio residuaalien normaalisuuden
# tarkistamiseksi. Tarkastellaan "potaska"-aineistoa. Siinä vertaillaan viiden
# eri käsittelyn vaikutusta puuvillakuitujen lujuuteen. Käsittelyt ovat käytettävän
# lannoitteen, potaskan, määriä. Jokainen käsittely toistetaan kolme kertaa.
# Saadaan 15 havaintoa. Mallinnetaan lujuuden riippuvuus potaskan
# määrästä.
#
# sink("A:\\TMallit\\Rscripts\\TA3.txt")
# komennolla voidaan tulostus ohjata tiedostoon
#
y <- scan("A:\\TMallit\\Datat\\Potaska.txt")
y
# Annataan koesuunnitelma
x <- scan()
36 36 36
54 54 54
72 72 72
108 108 108
144 144 144
 
# x on käsittely (paunaa/aari). Pauna on n. 1/2 kiloa. y on lujuus.
# Tyhjä rivi on loppumerkki aineiston syötössä.
tapply(y,x,mean)   # y:n keskiarvo x:n eri tasoilla 
plot(x,y)
regr <- lm(y ~ x); summary(regr)
potaska <- factor(x); potaska 
plot(potaska,y)
teeny <- lm(y ~ potaska)
anova(teeny
#summary(teeny) #Huom. Regr. on teeny:n erikoistapaus. teeny:llä parempi selitysaste
names(teeny)        
coefficients(teeny) # tulkitse
par(mfrow=c(2,1))
qqnorm(resid(teeny))
qqline(resid(teeny))
hist(resid(teeny))
plot(fitted(teeny),resid(teeny))
#
# Generoidaan vielä 100 alkion otos viiden vap. asteen t-jakaumasta.
# Arvioidaan poikkemaa normaalijakaumasta.
#
y <- rt(100,5)
hist(y)
qqnorm(y);qqline(y)