## Exercício 5

#Editando alguns parâmetros gráficos
e.salina=file.choose()
esa=read.table(e.salina, header=TRUE, sep=",",dec=".",as.is=TRUE)
head(esa)
str(esa)
plot(esa$h~esa$dap)

plot(esa$h~esa$dap, main="E.salina", xlab="dap (cm)", ylab="altura (m)")# grafico espalhagrama !

# explicação de Tick está dentro de plot-par
#colocando tick para dentro - sinal positivo

par(bty="l")# tick
par(tcl=0.3)# L 
plot(esa$h~esa$dap)# ver alterações, os nomes sumiram !

# colocando os nomes novamente ! 
plot(esa$h~esa$dap,main="E. salina", xlab="dap (cm)", ylab="altura (m)") 

# Fontes

par(cex=1.5)# aumentando a fonte - mais do que 1.5 não fica bom !
plot(esa$h~esa$dap, main="E. salina", xlab="dap (cm)", ylab="altura (m)") 

# Dois gráficos juntos

#(boxplot) da variável “dap” em função do fator “talhão
head(esa)
boxplot(esa$dap~esa$talhao)

#média e uma barra de desvio-padrão do dap, para cada talhão
summary(esa)
table(esa$talhao)
# tapply para devolver todos os tipos de talhao

mean.tl= tapply(esa$dap, esa$talhao, FUN= mean)#calculo da media por tal
sd.tl=tapply(esa$dap, esa$talhao, FUN= sd)#desvio por tal
mean.tl
sd.tl

par(mfrow=c(1,2))


tal=c(mean.tl,sd.tl) # juntando em um objeto
tal# ok
help(matrix)
mat=matrix(tal,6,2)#fazendo uma matrix com os dados!
mat

help(boxplot)


#Colocando a e b !
boxplot(mat,names=c("a","b"),)


#Adivinhando o código

adivinhando=file.choose()
adv=read.table(adivinhando,header=TRUE, sep=",", dec=".", as.is=TRUE)
adv

# o primeiro grafico é um espalhagrama das distancias em função do log do tamanho de algo !

#x2 são as dictancias ! comparando os valores com o grafico 
#x1 são os tamanhos do patch

str(adv)
class(adv$x1)
# logo é o grafico q é o plot de x2~(log)x1 
plot(adv$x2~adv$x1))


#Adivinhando o código

adivinhando=file.choose()
adv=read.table(adivinhando,header=TRUE, sep=",", dec=".", as.is=TRUE)
adv

# o primeiro gráfico é um espalhagrama das distancias em função do log do tamanho de algo !

#x2 são as dictancias ! comparando os valores com o grafico 
#x1 são os tamanhos do patch

str(adv)
class(adv$x1)
# logo é o grafico q é o plot de x2~(log)x1 

par(mfrow=c(1,2), bty="l", tck=0.05) # mfrow - dividir área no meio; tck - ticks do grafico, bty - forma do gráf (l)

#xlim, ylim - limites do gráfico, , pch - forma do ponto (triangula fechado)

plot(adv$x2~adv$x1,xlab="Log(Patch size)",ylab="Euclidean distances", ylim=c(0,3.5),xlim=c(0,2.8),pch=17)
help(abline)
abline(lm(adv$x2~adv$x1))
help(points)
args(text)
help(text)
text(2.5,3.5,"a",cex=1.5)# texto dentro 
# no gráfico do exercício um dos pontos estava cortado por causa do limite configurado para o eixo x
# aumentei esse limite para esse ponto aparecer !


# gráfico 2 - boxplot de distancia (x2) e categorias (y2)
#boxplot de 
par(cex = 0.8,bty="l")
boxplot(adv$x2~adv$y2, xlim=c(0.5,6.5),ylim=c(0,4), names=c("small","medium ed","medium in", "Large ed", "Large in", "control") )
text(1,3,"*",cex=1.5)
text(2,3,"*",cex=1.5)
text(3,3,"**",cex=1.5)
text(4,3,"*",cex=1.5)
text(5,3,"***",cex=1.5)
mtext("b",side=3,at=6)# texto nas margens