growth=read.table("http://user.mendelu.cz/drapela/Forest_Biometry/Data/growth.txt",header=TRUE)
growth
attach(growth)

id=as.factor(id)
species=as.factor(species)
time.point=as.factor(time.point)

library(psych)
library(ggpubr)
library(ez)

describeBy(height,species)
describeBy(height,time.point)
bxp <- ggboxplot(growth, x = "time.point", y = "height", color="species",palette="jco")
bxp
ggqqplot(growth, "height", facet.by = "time.point")
ggdensity(growth, "height", rug=T,color="species", facet.by = "time.point")

repANOVA = ezANOVA(growth,dv = .(height), wid = .(id), within = .(time.point),between = .(species), detailed = TRUE)
repANOVA

detach(growth)