a<-c(1,2,3,4,4,5,5,5,4,4,3,2,1)
b<-c(1,1,2,3,3,4,5,4,3,3,2,2,1)
c<-c(1,3,3,4,4,5,5,4,4,3,3,2,2)

data<-data.frame("hodnota"=c(a,b,c),
                 "kategorie"=rep(c("a","b","c"),each=13),
                 "opakovani"=rep(1:13,3))

boxplot(hodnota~kategorie,data=data)

model<-aov(hodnota~kategorie,data=data)
summary(model)

model_zavisle<-aov(hodnota~kategorie+Error(opakovani/(kategorie)),data=data)

summary(model_zavisle)


hist(c)