#1. Criando e Editando um gráfico

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

    *
      Legenda em ambos os eixos
    *
      Marcações (ticks) para dentro
    *
      Com apenas dois eixos (formato L)
    *
      Título do gráfico
    *
      Tamanho de fonte


	esaligna = read.csv("esaligna.csv")
	str(esaligna)
	head(esaligna)

#Plot de DAP e H

	plot(esaligna$dap~esaligna$ht)

#Inserindo legenda em ambos os eixos

	plot(esaligna$dap~esaligna$ht, xlab="Altura (m)", ylab="Diametro Altura do Peito (cm)")

#Com os ticks para dentro

	par(tcl=0.3)
	plot(esaligna$dap~esaligna$ht, xlab="Altura (m)", ylab="Diametro Altura do Peito (cm)",)

#Com apenas dois eixos

	plot(esaligna$dap~esaligna$ht, xlab="Altura (m)", ylab="Diametro Altura do Peito (cm)", bty="L")

#Titulo do grafico

	plot(esaligna$dap~esaligna$ht, xlab="Altura (m)", ylab="Diametro Altura do Peito (cm)", bty="L", main="Grafico de DAP x Altura")

#Tamanho da fonte

	par(c(cex.lab=1.5,cex.main=2,cex.axis=1.2))
	par(cex.lab=1)
	par(cex.main=1)
	par(cex.axis=1)

	plot(esaligna$dap~esaligna$ht, xlab="Altura (m)", ylab="Diametro Altura do Peito (cm)", bty="L", main="Grafico de DAP x Altura", cex=0.5)


#2. - Construindo dois gráficos

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

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).

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


	esaligna = read.csv("esaligna.csv")
	str(esaligna)
	head(esaligna)
	
#Criando a Matriz de Media e Desvio Padrao para cada talhao

	esaligna$talhao = as.factor(esaligna$talhao)

	talhao = sort(unique(esaligna$talhao))
	talhao = as.data.frame(talhao)
	media = tapply(esaligna$dap,esaligna$talhao,mean)
	var = tapply(esaligna$dap,esaligna$talhao,var)
	talhao$media = as.vector(media)
	talhao$var = as.vector(var)
	talhao

$Inserindo o grafico

	par(mfrow=c(1,2))
	boxplot(esaligna$dap~esaligna$talhao, main="Boxplot DAP X Talhões")
	text(16,20,"a",cex=1.8)
	plot(talhao$media~talhao$var, ylab="Media dos Talhões", xlab="Variança dos Talhões", 	main="Plot Var X Talhoes", bty="L")
	text(15,15,"b",cex=1.8)
	par(mfrow=c(1,1))

#3. - Adivinhando o Código

Leia o arquivo exercicio3 e usando os dados “x1 e y1” e “x2 e y2” tente reproduzir a figura.
	points = read.csv("exercicio3.csv")
	str(points)
	head(points)

	par(mfrow=c(1,2))
	par(bty="L")	
	par(cex.lab=1)
	par(cex.axis=1)
	par(mar=c(5,4,5,1))

# Plot x1 e y1
	plot(points$y1~points$x1, ylab="Euclidian Distances", xlab="Log(Patch Size)(ha)", ylim=c(0,3), xlim=c(0,2), pch=17)
	model <-  lm(points$y1~points$x1)
	abline(model, lwd=2)
	text(2,3,"a",cex=1.8)

#Boxplot de x2 e y2

	data = as.vector(c("Small","Medium \n Edge","Medium \n Interior","Large \n Edge","Large \n Interior","Control"))
	boxplot(points$x2~points$y2,names=data, 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)