egrandis <- read.csv("egrandis.csv", as.is=T, sep=";")
boxplot(dap~rotacao, egrandis)
boxplot(dap~regiao, egrandis)
qqnorm(egrandis$ht)
qqline(egrandis$ht,col="red",lwd=2)