### 5.1

esaligna <- read.table("esaligna.csv", header=TRUE, sep=",")
head(esaligna)

plot(esaligna$dap~esaligna$ht, xlab="Altura (m)", ylab="DAP (cm)", tcl=0.3,
	bty="l", main="Relação DAP/altura",cex.lab=1.3, cex.axis=1.3, cex.main=1.3)



### 5.2

par(mfrow=c(1, 2))

boxplot(esaligna$dap~esaligna$talhao, xlab="Talhão", ylab="DAP (cm)")
text(x=1, y=23, "a")

esaligna$talhao <- as.factor(esaligna$talhao)
(media <- aggregate(esaligna$dap, by=list(esaligna$talhao), FUN=mean))
(desvio <- aggregate(esaligna$dap, by=list(esaligna$talhao), FUN=sd))
(matriz <- matrix(c(media$x, desvio$x), nrow=6, ncol=2))
rownames(matriz) <- c(16:18,22,23,32)
plot(matriz[,1], ylim=c(5,23), xlab="Talhão", ylab="DAP (cm)")
segments(x0=c(1:6), y0=media$x-desvio$x, x1=c(1:6), y1=media$x+desvio$x)
text(x=1.5, y=23, "b")



### 5.3

ex <- read.csv("exercicio3.csv", as.is=T, header=T)
head(ex)
str(ex)
ex$y2 <- factor(ex$y2)
ex <- ex[,-2]
ex1 <- na.omit(ex[,1:2])

par(mfrow=c(1,2), bty="l", tcl=0.3)
plot(ex1$x2~ex1$x1, tcl=0.3, pch=17, cex=1.5, xlim=c(0.5,2.3), ylim=c(0,3), 
	xlab="Log(Patch size)(ha)", ylab="Euclidean distances")
modelo <- lm(ex1$x2~ex1$x1)
abline(modelo, lwd=3)
text(2.3,3, "a")

boxplot(ex$x2~ex$y2, ylim=c(0,3), outline=FALSE, 
	names=c("Small","Medium\nEdge","Medium\nInterior","Large\nEdge",
	"Large\nInterior","Control"))

text(1,3, "*")
text(2,3, "*")
text(3,3, "**")
text(4,3, "*")
text(5,3, "***")
text(6.5,3, "b")

savePlot(filename="ex3",type="jpeg")