taim=rep(c("a","b"),c(6,6))
kiirus=c(23,24,23,26,27,29,34,30,29,32,34,29)
röövik=data.frame(taim,kiirus)
#algandmed valmis
t.test(kiirus~taim,data=röövik,var.equal=T)
t.test(kiirus~taim,data=röövik,var.equal=F)
#kas dispersioonid on võrdsed? fligneri test on mittenormaalsuse suhtes robustne
fligner.test(kiirus,as.factor(taim))

enne=c(23,26,25,27,23,24,30,40)
pärast=c(34,27,32,29,25,26,37,45)
vahe=pärast-enne
t.test(vahe)
#t.test(pärast,enne,paired=T)

taim=rep(c("a","b","c"),c(6,6,6))
kiirus=c(23,24,23,26,27,29,34,30,29,32,34,29,23,19,15,20,24,25)
röövik=data.frame(taim,kiirus)
fligner.test(kiirus,as.factor(taim))
#teeb ainult mudeli
mudel=lm(kiirus~as.factor(taim))
#ise tuleb välja trükkida
anova(mudel)
summary(mudel)
#kruskal wallise test (SASis Wilcoxoni test)
kruskal.test(kiirus,as.factor(taim))