setwd("C:/Users/Fabio Caltabellotta/Documents/Disciplina R")
e.saligna<- read.csv("arvore.csv", header=TRUE, sep=",")
e.saligna
names(e.saligna)
str(e.saligna)
summary(e.saligna)

### Editando alguns parâmetros gráficos ###

plot(dap~ht, data=e.saligna, xlab="altura (m)", ylab="diâmetros à altura do peito (cm)",
 tck=0.020, bty="l", main="Dispersão entre dap~ht", cex.lab=1.25)

### Dois gráficos juntos ###

boxplot(dap~talhao, data=e.saligna, xlab="talhão", ylab="altura (m)", ylim=c(5,25)) ## gráfico 1 em uma só janela

media<- tapply(e.saligna$dap,factor(e.saligna$talhao),FUN=mean)### cálculo das médias de dap por talhão (OBS: talhão em fator)
desv.pad<- tapply(e.saligna$dap,factor(e.saligna$talhao),FUN=sd) ### cálculo dos desvios padrão de dap por talhão (OBS: talhão em fator)
media ### médias
desv.pad ### desvios padrão
matriz<- matrix(c(media, desv.pad),nrow=6,nco=2, dimnames=list(c("16","17","18","22","23","32"), 
 c("media","desv.pad"))) ### criando uma matriz com a média e o desvio padrão calculados por talhão
matriz

par(mfrow=c(1,2)) ## dois gráficos na mesma janela
boxplot(dap~talhao, data=e.saligna, main="a", xlab="talhão", ylab="altura (m)", ylim=c(5,25),xlim=c(1,6), xaxt="n",
 pars = list(boxwex = 0.3, staplewex = 0.4, outwex = 0.3)) ### gráfico 1
axis(1, 1:6, c("16","17","18","22","23","32")) 
plot(matriz[,"media"], xlab= "talhão", ylab="altura (m)", main="b", xlim=c(1,6), ylim=c(5,25), xaxt="n", pch=20) ### plot das médias com o desvio padrão por talhão  
axis(1, 1:6, c("16","17","18","22","23","32"))  
arrows(1:6,matriz[,"media"]-matriz[,"desv.pad"], 1:6,matriz[,"media"]+matriz[,"desv.pad"],code=3,angle=90,length=0.04)

### Adivinhando o código ###

code <- read.csv("adivinhe.csv", header=TRUE, sep=",", dec = ".", as.is=TRUE, na.strings="NA")
code
head(code)
tail(code)
str(code)
summary(code)

par(mfrow=c(1,2), pch=17, tcl=0.15, family="serif", bty="l", mar=c(5.2,5.2,5.2,2.2)) ### mfrow: divide as janelas
### pch: estilo do ponto ### tcl: tamanho do ponto ### family: tipo de letra ###  bty: contornos ### mar: tamanho da margem 
plot(code$y1~code$x1,  ylim=c(0,3), xlim=c(0.5,2.5), xlab="", ylab="") ### plota os dados y1~x1 ### ylim e xlim: limites ### xlab e ylab:legendas
mtext("Log(Patch size)(ha)",line=3,side=1,cex=1.3, family="serif") ## essa função insere "Log(Patch size)(ha)" e os args controlam posição e tamanho
##
mtext("Euclidean distance",side=2,line=3,cex=1.3, family="serif") ### essa função insere "Euclidean distance" e os args controlam posição e tamanho
code.lm<- lm(code$y1~code$x1) ### regressão linear
abline(code.lm,lwd=3)  ## a função abline coloca a linha da regressão e lwd muda a grossura da linha
text(2.3,3.0,"a", cex=1.3) ## coloca um texto na area do grafico definindo o tamanho e a sua posição

code$y2<- factor(code$y2)
str(code)
boxplot(code$x2~code$y2,outline=FALSE, ylim=c(0,3), names=c("","","","","","")) ### gráfico de bigodes ## outline=FALSE -> tira os valores extremos
nomes=c("Small","Medium","Medium","Large","Large","Control")  
nome=c("","Edge","Interior","Edge","Interior","")
mtext(nomes,side=1,line=1,at=1:6, cex=1.0)  ### defini qual é a posição de cada nome
mtext(nome,side=1,line=2,at=1:6, cex=1.0) ### defini qual é a posição de cada nome
text(6.5,3.0,"b", cex=1.8) ## identificação da figura (b)
text(x=1,y=3,"*",cex=1.5) ## coloca o *
text(2,3,"*",cex=1.5) ## coloca o *
text(3,3,"**",cex=1.5) ## coloca o *
text(4,3,"*",cex=1.5) ## coloca o *
text(5,3,"***",cex=1.5) ## coloca o *