####5.1######

dir()
getwd()
setwd("C:/Users/Mauro Brum-Jr/Documents/Arquivos_R/Exercicio5")
dir()
esaligna<-read.table("esaligna.csv", header=T, sep=",")
esaligna
plot(dap~ht, data=esaligna)
#colocando nomems das variáveis e suas unidades
plot(dap~ht, data=esaligna, xlab="Altura total (m)", ylab="DAP (cm)")
#colocando os ticks para dentro da área do gráfico
plot(dap~ht, data=esaligna, xlab="Altura total (m)", ylab="DAP (cm)", tcl=0.3)
#colocando apenas dois eixos no formato L
plot(dap~ht, data=esaligna, xlab="Altura total (m)", ylab="DAP (cm)", tcl=0.3, bty="l")
#colocando título informativo
plot(dap~ht, data=esaligna, xlab="Altura total (m)", ylab="DAP (cm)", tcl=0.3, bty="l", main="Eucalyptus saligna")

#alterando o tamanho da fonte maior que o padrão
#COMANDO FINAL DO EXERCICIO
par(las=1)
plot(dap~ht, data=esaligna, xlab="Altura total (m)", ylab="DAP (cm)", tcl=0.3, bty="l", main="Eucalyptus saligna", cex.lab=1.5, cex.main=1.5)


###########5.2#############
##### 5.2 #############
getwd()
dir()
esaligna<-read.table("esaligna.csv", header=T, sep="," )
str(esaligna)
talhao<-factor(esaligna$talhao)
esaligna$talhao<-talhao
str(esaligna)
#boxplot
par(mfrow=c(1,2))
par(las=1, tcl=0.2, cex.axis=0.6)
boxplot(esaligna$dap~esaligna$talhao, xlab="Talhão", ylab="DAP (cm)")
text(1,22,labels="a")
#segundo gráfico aplicar os desvios padrões em valores de média
# em gráficos. O uso da funçao segment()
head(esaligna)
str(esaligna)
grafico.media<-aggregate(esaligna$dap, by=list(esaligna$talhao), FUN=mean)
grafico.desvio<-aggregate(esaligna$dap, by=list(esaligna$talhao), FUN=sd)
grafico.media$sd<-grafico.desvio$x
colnames(grafico.media)<-c("talhao","media","desvpad")
grafico.media
str(grafico.media)
plot(media~talhao, data=grafico.media, ylim=c(0,25), xlab="Talhão", ylab="DAP (cm)")
#segments()
##a funçao segments adiciona uma linha de um segmento em um 
#par de pontos. Os primeiros argumentos é a origem da reta 
#a ser desenhada. 
segments(x0=1:6, y0=grafico.media$media-grafico.media$desvpad, x1=1:6, y1=grafico.media$media+grafico.media$desvpad)
#colocando indices A e B nos gráficos
text(1,24,labels="b")

#########################5.3###################################
######5.3######
getwd()
setwd("C:/Users/Mauro Brum-Jr/Documents/Arquivos_R/Exercicio5")
misterio<-read.table("exercicio3.csv", header=T, sep="," )
str(misterio)
head(misterio)
tail(misterio)
misterio
par(mfrow=c(1,2), family="serif", mar=c(5,3.5,3,1))
#plotando o gráfico
plot(y1~x1, data=misterio, 
     xlab="Log(Patch size)(ha)", ylab="Euclidean distances", 
     family="serif", pch=17, bty="l", ylim=c(0,3), xlim=c(0.5,2.5),
     tcl=0.3, cex.lab=1.5, mgp=c(2.3,1,0))

#adicionando a linha do modelo linear
modelo<-lm(misterio$y1~misterio$x1)
abline(modelo) 
#adicionando o indice "a" na figura
text(2.5,3,"a", cex=1.2)

#plotando o boxplot
par(bty="l", tcl=0.3)
boxplot(x2~y2, data=misterio, bty="l", ylim=c(0,3),outline=F, names=F, mar=c(5,0,3,1))
mtext(text=c("Small", "Medium\nEdge", "Medium\nInferior",
             "Large\nEdge","Large\nInferior", "Control"),
      side=1, at=1:6, cex=1.2, mgp=c(3,1,3), line=1.5)
asterisco<-c("*","*","**","*","***","b")
text(1:6,3, labels=asterisco)