# Exercicio 5 - Criação de gráficos
# Camila Galheigo Coelho


setwd("C:/Users/Camila/Documents/PESQUISA/1_DOUTORADO/Disciplinas/R")
getwd()

read.table("esaligna.csv", header=T, sep=",")
esaligna <- read.table("esaligna.csv", header=T, sep=",")
summary(esaligna)


# Criando e editando um gráfico de dispersão entre “dap” e “ht” #


plot(esaligna$dap~esaligna$ht)

#Q1: legenda em ambos os eixos
plot(esaligna$dap~esaligna$ht,xlab="Altura (m)", ylab="DAP (cm)")

#Q2:
?plot  #length of tick marks as a fraction of the smaller of the width or height of the plotting region. Default setting (tck = NA) is to use tcl = -0.5.
plot(esaligna$dap~esaligna$ht,xlab="Altura (m)",ylab="DAP (cm)",tck=0.03)  

#Q3:
?par  #bty: A character string which determined the type of box which is drawn about plots.
#If bty is one of "o" (the default), "l", "7", "c", "u", or "]" the resulting box resembles the corresponding upper case letter. "n" suppresses the box.
plot(esaligna$dap~esaligna$ht,xlab="Altura (m)",ylab="DAP (cm)",tck=0.03,bty="l") 

#Q4:
args(title)
title(main="Relação Altura da árvore e DAP")

#Q5: 
plot(esaligna$dap~esaligna$ht,xlab="Altura (m)",ylab="DAP (cm)",tck=0.03,bty="l",cex=2, cex.main=2, cex.lab=2, cex.axis=2)
title(main="Relação Altura da árvore e DAP")


# Construindo dois gráficos lado a lado entre “dap” e “classe” #

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="Talhão",ylab="DAP")
?boxplot
?par
table(esaligna$talhao)
boxplot(esaligna$dap~esaligna$talhao,xlab="Talhão",ylab="DAP",ylim=c(5,25))
title(main="A")
media.dap
plot(media.dap,main="B", xlab="Talhão",ylab="Média DAP")
plot(media.dap,main="B", xlab="Talhão",ylab="Média DAP",tck=0.03, ylim=c(2,22), xlim=c(0,7), xaxis=c(16,17,18,22,23,31))
?arrows
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="Talhão",ylab="DAP",pch="*",ylim=c(5,25))
title(main="A")
plot(media.dap,main="B", xlab="Talhão",ylab="Média 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)

save.image()

# Adivinhando o código #
adv.cod<-read.csv("exercicio3.csv", header=TRUE, sep=",")
adv.cod
par(mfrow=c(1,2))
plot(adv.cod$x1~adv.cod$y1,main="a", xlab="Log(Patch size)(ha)",ylab="Euclidean distances",tck=0.03, ylim=c(0,3),xlim=c(0.5,3),bty="l",pch=17)
?abline
lines(    #???
?par
boxplot(adv.cod$x2~adv.cod$y2,main="b",tck=0.03, ylim=c(0,3),bty="l")
?boxplot
?par

save.image()