########## OLIDAN POCIUS __ EXERCICIO 5
 ## Editando alguns parâmetros gráficos
esaligna <- read.csv("C:/Users/OPocius/Desktop/RAULAS/Aula5/esaligna.csv")
View(esaligna)
plot(esaligna$ht~esaligna$dap)
plot(esaligna$ht~esaligna$dap,xlab="DAP(cm)",ylab="Altura(m)")
plot(esaligna$ht~esaligna$dap,xlab="DAP(cm)",ylab="Altura(m)",tck=0.03)
plot(esaligna$ht~esaligna$dap,xlab="DAP(cm)",ylab="Altura(m)",tck=0.03,bty="l")
title(main="DAP vs Altura")
plot(esaligna$ht~esaligna$dap,xlab="DAP(cm)",ylab="Altura(m)",tck=0.03,bty="l",cex=2, cex.main=2, cex.lab=2, cex.axis=2)

## Dois gráficos juntos
media.dap<-tapply(esaligna$dap,esaligna$talhao,mean)
media.dap
desviopadrao.dap <- tapply(esaligna$dap,esaligna$talhao,sd)
desviopadrao.dap

boxplot(esaligna$dap~esaligna$talhao,xlab="Talhao",ylab="DAP")
table(esaligna$talhao)
boxplot(esaligna$dap~esaligna$talhao,xlab="Talhao",ylab="DAP",ylim=c(5,25))
title(main="A")
media.dap
plot(media.dap,main="B", xlab="Talhao",ylab="Mean DAP")
arrows(x0=1:6, y0=media.dap[1:6]-desviopadrao.dap[1:6], x1=1:6, y1=media.dap[1:6]+desviopadrao.dap[1:6],angle=90,code=3)

par(mfrow=c(1,2))
boxplot(esaligna$dap~esaligna$classe,xlab="Talhao",ylab="DAP",pch="*",ylim=c(5,25))
title(main="A")
plot(media.dap,main="B", xlab="Talhao",ylab="Mean DAP",tck=0.03, ylim=c(5,25), xlim=c(0,7))
arrows(x0=1:6, y0=media.dap[1:6]-desviopadrao.dap[1:6], x1=1:6, y1=media.dap[1:6]+desviopadrao.dap[1:6],angle=90,code=3)

## Adivinhando o Código