egrandis <- read.table("egrandis.csv",header=TRUE,sep=";")

summary(egrandis)
head(egrandis)

names(egrandis)

par(mfrow (1,2))
boxplot (egrandis$dap~egrandis$regiao, data=egrandis)
boxplot (egrandis$dap~egrandis$rotacao, data=egrandis)
dev.off()


par(mfrow (1,1))
qqnorm(egrandis$ht)
qqline(egrandis$ht)