#5.1 Editando alguns parâmetros gráficos
esaligna <- read.table("esaligna.csv", header=T, sep=",", as.is=F)
str(esaligna)
plot(esaligna$dap~esaligna$ht, xlab="Altura (m)", ylab="Dap (cm)", tcl=0.5, bty="l", main="Relacao Dap~h", font.axis=4, font.main=2, font.lab=3)

#5.2 Dois gráficos juntos
esaligna <- read.table("esaligna.csv", header=T, sep=",", as.is=F)
str(esaligna)
par(mfrow=c(1,2))
boxplot(esaligna$dap~esaligna$talhao, xlab="Talhao", ylab="Dap(cm)")
text(1,22, "a", cex=1.8) 
mean.plot<- aggregate(x=esaligna$dap, by=list(esaligna$talhao),FUN=mean)
mean.plot
sd.plot <- aggregate(x=esaligna$dap, by=list(esaligna$talhao),FUN=sd)
sd.plot
plot(mean.plot[,2], xlab="Talhao", ylab="Media Dap(cm)", ylim=c(0,25), pch=19, xaxt="n")
y0=mean.plot[,2]-sd.plot[,2] 
y1=mean.plot[,2]+sd.plot[,2]
x0=c(1:6)
segments(x0, y0, x0, y1)
axis(side=1, at=c(1,2,3,4,5,6), labels=c(16,17,18,22,23,32))
text(5.5,25, "b", cex=1.8)

#5.3 Adivinhando o código
exerc5.3 <- read.table("exercicio5.3.txt", header=T, sep=",", as.is=F)
str(exerc5.3)
par(mfrow=c(1,2))
plot(exerc5.3$x1~exerc5.3$y1, xlab="Log(Patch size)(ha)", ylab="Euclidean Distances", xlim=c(0.5,2.5), ylim=c(0,3), pch=17, tcl=0.5, bty="l", cex.axis=1.5, family="serif", cex.lab=1.5, family="serif")
exerc5.3.lm <- lm(exerc5.3$x1~exerc5.3$y1)
abline(exerc5.3.lm, lwd=2)
text(2.4,3, "a", cex=1.8)
exerc5.3$y2 <- as.factor(exerc5.3$y2)
exerc5.3$y2
boxplot(exerc5.3$x2~exerc5.3$y2, xlab="", ylab="", par(tcl=0.5), ylim=c(0,3), par(bty="l"), xaxt="n", family="serif", cex.axis=1.2)
axis(side=1, at=c(1,2,3,4,5,6), labels=c("Small", "Medium Edge", "Medium Interior", "Large Edge", "Large Interior", "Control"), cex.axis=1.2, family="serif")
text(c(1:6,6.5),3, c("*","*","**","*","***","b"))