##exercicio 5
##exercicio 5.1
getwd()
setwd("C:/Users/Osmarina Marinho/Documents/teste")
esaligna<-read.table("esaligna.csv",header=T,sep=",", as.is=T)
##5.1 Editando alguns parâmetros gráficos
plot(esaligna$dap~esaligna$h)
plot(esaligna$dap~esaligna$ht,xlab="ht(m)",ylab="Dap(cm)")
par(tcl=0.3)
par(las=1)
par(cex.axis=1.5)
par(cex.lab=1.5)

plot(esaligna$dap~esaligna$ht,xlab="ht(m)",ylab="Dap(cm)",main="DAPxAltura", xlim=c(0,30),ylim=c(0,30), bty="n", xaxt="n", yaxt="n")
axis(side=1,at=seq(0,30,by=10),labels=c("",10,20,30),pos=0)
axis(side=2,at=seq(0,30,by=10),labels=seq(0,30,by=10),pos=0)

##5.2 Dois gráficos juntos
par(mfrow=c(1,2))
par(las=1)
boxplot(esaligna$dap~esaligna$talhao,ylab="DAP(cm)")
mtext("Talhão",side=1,cex=1.5,line = 2.5)
#####locator(1) usar locator
text(x= 0.956,y=21.8147,"a",cex=1)

esaligna$talhao=as.factor(esaligna$talhao)
media.dap<-aggregate(esaligna$dap, by= list(esaligna$talhao),FUN=mean)

sd.dap<-aggregate(esaligna$dap, by= list(esaligna$talhao),FUN=sd)

plot(media.dap$x,xaxt="n", ylim=c(0,30),ylab = "Média DAP(cm)",xlab="")
axis(side=1,at=1:6,labels = c(16,17,18,22,23,32))
segments(x0=1:6,x1=1:6,y0=(media.dap$x-sd.dap$x),y1=(media.dap$x+sd.dap$x))
mtext("Talhão",side=1,cex=1.5,line = 2.5)
###locator(1)
text(x=1.696,y=27.9205,"b",cex=1)


##5.3Adivinhando o código
###pch escolhe o tipo de ponto, tcl controla o tamanho e direção dos 
# risquinhos do eixo, family cntrola o tipo de letra, bty a quantidade de contornos, mar o tamanho da margem
##Plot x1 e y1
par(mfrow=c(1,2))
arquivo<-read.table("arquivo.csv", header=T,sep=",", as.is=T)
par()$mar
par(las=1)
plot(arquivo$y1~arquivo$x1,xlab="Log(Patch size(ha)",ylab="Euclidian distances", xlim=c(0.5,2.3),ylim=c(0,3), bty="l", xaxt="n", yaxt="n",pch=17)

axis(side=1,at= seq(0,2,by=0.5),tcl=0.3)
axis(side=2,at=seq(0,3,by=0.5),tcl=0.3, cex=0.3)

model<-lm(arquivo$y1~arquivo$x1)
abline(model,lwd=2)
segments(x0=0.49,y0=2.204,x1=2.3084,y1=1.2804)
text(x=2.32,y=2.900,"a", cex=1)


##boxplot 
arquivo$y2<-as.factor(arquivo$y2)
class(arquivo$y2)
###outline=false remove o outline
par(tcl=0.3)
par(bty="l")
par(las=1)
par(cex.axis=0.8,cex.lab=0.8)
boxplot(arquivo$x2~arquivo$y2, xaxt="n", outline=FALSE)
nomex<-c("","Small", "Medium\n edge", "Medium\n interior", "large\n edge", "Large\n interior", "Control","")
axis(side=1,at=0:7,labels = nomex,tcl=0.3,cex=1)

##colocabdo texto dentro do grafico
text(x=1.01,y=2.907, "*",cex=1)
text(x=1.99,y=2.907, "*",cex=1)
text(x=2.99,y=2.907, "**",cex=1)
text(x=3.98,y=2.907,"*",cex=1)
text(x=5.00,y=2.907,"***",cex=1)
text(x=5.96,y=2.707,"b", cex=1)