##5.1 Editando alguns parametros graficos
dir()
esaligna<-read.table("esaligna.csv", header=T,sep=",")
head(esaligna)
str(esaligna)
par(cex=0.8)
plot(esaligna$dap~esaligna$ht,xlab="Altura(m)",ylab="DAP(cm)",main="Relação entre DAP e Altura",tcl=0.3,bty="L")

##5.2 Dois graficos juntos
head(esaligna)
par(mfrow=c(2,1),mar=c(6,4,2.5,2))
boxplot(esaligna$dap~esaligna$talhao,xlab="Tamanho talhão(cm)",ylab="Dap(cm)")
mtext("a",side=2,line=1,at=30,cex=0.8,adj=0,padj=1,las=1)
summary(esaligna)
sd(esaligna$dap)
par(mar=c(6,3,2,2))
plot(x=NULL,y=NULL,xlim=c(15,33),ylim=c(0,22),bty="L",xlab="Tamanho do talhao(cm)",ylab="")
media.por.talhao<-tapply(esaligna$dap, INDEX=esaligna$talhao,mean)
talhao<-c(16,17,18,22,23,32)
points(x=talhao,y=media.por.talhao)
sd.por.talhao<-tapply(esaligna$dap,INDEX=esaligna$talhao,sd)
segments(x0=talhao,y0=media.por.talhao-sd.por.talhao,x1=talhao,y1=media.por.talhao)
segments(x0=talhao,y0=media.por.talhao,x1=talhao,y1=sd.por.talhao+media.por.talhao)
sd.por.talhao
media.por.talhao
mtext("b",side=2,line=1,at=30,cex=0.8,adj=0,padj=1,las=1)

##5.3 Adivinhando o codigo
dir()
exercicio5.3<-read.table("exercicio3.csv", header=T,sep=",")
head(exercicio5.3)
exercicio5.3
par(mfrow=c(1,2))
modelo<-lm(exercicio5.3$x1~exercicio5.3$y1)
par(mar=c(5,5,4,2))
plot(exercicio5.3$x1,exercicio5.3$y1,ylab="Euclidiean distances",xlab="Log(Patch size)(ha)",pch=17,xlim=c(0.5,2.25),ylim=c(0,3),bty="L",tcl=0.3,cex.lab=1)
abline(modelo,lwd=2)
text(2.2,3,"a")
eixo<-c(rep("Small",6),rep("Medium\nEdge",4),rep("Medium\nInterior",8),rep("Large\nEdge",4),rep("Large\nInterior",5),rep("Control",4))
exercicio5.3$eixo.x<-factor(eixo,levels=c("Small","Medium\nEdge","Medium\nInterior","Large\nEdge","Large\nInterior","Control"))
levels(exercicio5.3$eixo.x)
head(exercicio5.3)
par(bty="L",tcl=0.3,cex.axis=1,mar=c(3,1,4,0.5))
boxplot(exercicio5.3$x2~exercicio5.3$eixo.x,ylim=c(0,3),outline=F)
text(1,3, "*",cex=1.5)
text(2,3,"*",cex=1.5)
text(2.9,3,"*",cex=1.5)
text(3,3,"*",cex=1.5)
text(4,3,"*",cex=1.5)
text(4.9,3,"*",cex=1.5)
text(5,3,"*",cex=1.5)
text(5.1,3,"*",cex=1.5)
text(6.3,3,"b")