setwd("c:/Users/brozoval/Desktop/6. cvičení_ANOVA/")

#install.packages("readxl")
library("readxl")

trenink <- read_excel("trenink.xlsx",sheet = "data",col_names=T)
trenink<-as.data.frame(trenink)
attach(trenink)

### overeni predpokladu
#normalita v kazde kombinaci skupin

windows()
par(mfrow=c(3,2))
for (i in 1:2) 
{for (j in 1:3)
{
  spoj<-Cas[which(Zpusob_treninku==unique(Zpusob_treninku)[j] & Suplementace==unique(Suplementace)[i])]
  hist(spoj,
       main=paste0(unique(Zpusob_treninku)[j],"+",unique(Suplementace)[i],", p=",round(shapiro.test(spoj)$p.value,3)),
       freq = F,
       xlab="Cas")
  lines(seq(min(spoj), max(spoj), by=.1), 
        dnorm(seq(min(spoj), max(spoj), by=.1),mean(spoj), sd(spoj)), col="blue")
}}

# homogenita rozptylu

library("car")
leveneTest(Cas ~ (Suplementace*Zpusob_treninku),center=mean)

# nahled na data
boxplot(Cas ~(Suplementace*Zpusob_treninku))

# ANOVA s interakci
anova_res<-aov(Cas ~ Zpusob_treninku*Suplementace)
anova_res<-aov(Cas ~ Zpusob_treninku+Suplementace+Zpusob_treninku:Suplementace) # jinak zapsano
summary(anova_res)

# Post hoc testovani
TukeyHSD(anova_res, which = "Zpusob_treninku")
TukeyHSD(anova_res, which = "Zpusob_treninku:Suplementace")