###5.1 Editando alguns parâmetros gráficos
##5.Crie um gráfico de dispersão entre “dap”1) e “ht”2) com: 

# Legendas dos eixos com nomes das variáveis e suas unidades
# Marcações do eixos (ticks) para dentro da área do gráfico 
# Apenas dois eixos (formato “L”)
# Título informativo
# Tamanho das fontes maiores que o padrão

e.saligna<-read.csv("esaligna.csv")
str(e.saligna)
head(e.saligna)
summary(e.saligna)
names(e.saligna)

#nomes
dap<-"Diâmetro à altura do peito(cm)"
tit.saligna<-"Relação DAP x Altura"

par(cex=1.2)
plot(e.saligna[,4]~e.saligna[,5], xlab="Altura(m)",ylab=dap,
          tcl=0.4,bty="l",main=tit.saligna)

##########################
###5.2 Dois gráficos juntos

# Use as variáveis “dap” e “talhao” para construir dois gráficos, 
#colocando-os lado a lado. O primeiro deve ser um gráfico de desenho de 
#caixa (boxplot) da variável “dap” em função do fator “talhão”. 
#O segundo deve ter apenas a média e uma barra de desvio-padrão do dap, 
#para cada talhão.
# Insira também uma letra para dizer qual é o gráfico “a” e qual é o “b”
# (tanto faz, quem é um e quem é outro).

par(mfrow=c(1,2))
names(e.saligna)
as.factor(e.saligna[,3])

#primeiro
boxplot(e.saligna[,4]~e.saligna[,3],ylab="DAP(cm)",xlab="Talhão")
text(6.5,23,"a")

#segundo
media.dap<-tapply(e.saligna[,4],INDEX=e.saligna[,3],FUN=mean)
desvio.dap<-tapply(e.saligna[,4],INDEX=e.saligna[,3],FUN=sd)

(dap.talhao<-matrix(c(media.dap,desvio.dap),nrow=6,ncol=2))
rownames(dap.talhao)<-c(16,17,18,22,23,32)

desv.pos<-dap.talhao[,1]+dap.talhao[,2]
desv.neg<-dap.talhao[,1]-dap.talhao[,2]

plot(dap.talhao[,1],xaxt="n",ylab="DAP(cm)",ylim=c(0,22),xlab="Talhão")

axis(side=1,labels=c(16,17,18,22,23,32), at=1:6)
text(6,22,"b")

arrows(seq(1:6),dap.talhao[,1],seq(1:6),desv.pos,angle=90, code=2,length=0.05)
arrows(seq(1:6),dap.talhao[,1],seq(1:6),desv.neg,angle=90,code=2,length=0.05)

#####################

#### 5.3Adivinhando o código
##Figura 1

ex.3<-read.csv("exercicio3.csv")
str(ex.3)
head(ex.3)
summary(ex.3)
names(ex.3)

par(mfrow=c(1,2),mar=c(5,4,2,1))
ED<-"Euclidean distances"

par(bty="L",tcl=0.3,pch=17,cex=1,cex.axis=0.8)
plot(ex.3[,1],ex.3[,2],xlab="Log(Patch size)(ha)",ylab=ED,ylim=c(0,3.0),
       xlim=c(0.5,2.0),cex.lab=1.1)

#model.ex.3 <-  lm(ex.3[,1]~ex.3[,2])
#abline(model.ex.3)#não é identica# intersepto e inclinação diferentes

abline(2.5,-0.68,lwd=2.5)

text(2.05,3,"a")

###Figura 2

ex.3$y2<-as.factor(ex.3[,4])
class(ex.3[,4])

plot(ex.3[,4],ex.3[,3])
identify(ex.3[,4],ex.3[,3])

#27- linha que esta o ponto,excluir essa linha
#não deu certo excluir linha pq boxplot muda a forma

plot(ex.3[,4],ex.3[,3],ylim=c(0,3.0),,xaxt="n",tcl=0.3)

plot.box<-c("","Medium\nEdge","Medium\nInterior",
            "Large\nEdge","Large\nInterior","")

mtext("Small",side=1,at=1,line=0.5,cex=0.9)
axis(side=1, at=1:6,label=plot.box)
mtext("Control",side=1,at=6,line=0.5,cex=0.9)

#inserir figura pra ficar acima do ponto circular (ponto da linha 27)
rect(4.7,1.6,5.2,1.8,col="white",border=NA)

text(1:5,3,c("*","*","**","*","***"))
text(6.5,3,"b")



##################         ##############################