LUDMILA MARIA RATTIS TEIXEIRA

Exercício 5:

Usando o conjunto de dados Conjunto de Dados: Dados de Biomassa de Árvores de Eucalyptus Saligna, construa os seguintes gráficos: 

##Criando e editando um gráfico

Um gráfico de dispersão entre “dap” e “ht” editando:

1-)Legenda em ambos os eixos

2-)Marcações (ticks) para dentro 

3-)Com apenas dois eixos (formato L)

4-)Título do gráfico

R:
  plot(esaligna$dap~esaligna$ht, xlab="altura", ylab="DAP", bty="l", tcl=0.3,
  title(main="Altura de E.saligna em função do DAP"))


5-)Tamanho de fonte

R:
  par(cex.axis=1.5)
  par(cex.lab=2)


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

##Construindo dois gráficos

Dois gráficos lado a lado entre “dap” e “talhao”:

1-)Construir um gráfico de desenho de caixa (boxplot) entre as variáveis
“dap” e “classe” e ao lado um plot com apenas a média e desvio-
padrão. E insira também uma letra para dizer qual é o gráfico “a” e
qual é o “b” (tanto faz, quem é um e quem é outro). 

##Problema no enunciado

esaligna$classe = as.factor(esaligna$classe)
classe = sort(unique(esaligna$classe))
classe = as.data.frame(classe)
media = tapply(esaligna$dap,esaligna$classe,mean)
sd = tapply(esaligna$dap,esaligna$classe,sd)
classe$media = as.vector(media)
classe$sd = as.vector(sd)
classe
par(mfrow=c(1,2))
boxplot(esaligna$dap~esaligna$classe, main="Boxplot DAP X Classes")
text(15,15,"a",cex=1.8)
plot(classe$media~classe$sd, ylab="Media das Classes", xlab="Desvio das Classes", main="Plot Desvio X Classes", bty="L")
text(15,15,"b",cex=1.8)
par(mfrow=c(1,1))

Dica: vocês terão que calcular a média e os devio padrão do dap das
árvores em cada classe. Depois crie uma matriz com estes valores e
crie o plot destes valores. 
###################################################################

##Adivinhando o código

1-)Leia o arquivo exercicio3 e usando os dados “x1 e y1” e “x2 e y2”
tente reproduzir esta figura: 

R:

##Gráfico 1

desc=read.table("exercicio3.csv",as.is=T,header=T,sep=",")
par(mfrow=c(1,2))
str(desc)
plot(desc$y1~desc$x1,xlab="Log(PatchSize)ha",ylab="Euclidean distances",xaxp=c(0.5,2.5,5),yaxp=c(0,3,8),bty="l",pch=17)
model=lm(desc$y1~desc$x1)
abline(model,lwd=1.5)
text(3,3,"a",cex=3.8)

##Gráfico 2

legenda.x = as.vector(c("Small","Medium \n Edge","Medium \n Interior","Large \n Edge","Large \n Interior","Control"))
par(bty="l")
par(cex.lab=1)
par(cex.axis=1)
par(mar=c(5,4,5,1)) 
boxplot(desc$x2~desc$y2,names=legenda.x, cex.axis=0.7, ylim=c(0,3))
text(7,3,"b",cex=1,5)
text(1,3,"*",cex=2)
text(2,3,"*",cex=2)
text(3,3,"**",cex=2)
text(4,3,"*",cex=2)
text(5,3,"***",cex=2)
par(mfrow=c(1,1)