#O exercicio foi postado primeiro no wiki!!! #5.1 setwd("C:/Users/Camila Hohlenwerger/Desktop/aula4.exercicio_aula5") esaligna<- read.table(file="esaligna.csv", header=TRUE, sep=",", dec=".", as.is=TRUE) esaligna plot(esaligna$ht~esaligna$dap, xlab="Dap (cm)", ylab="Altura (m)") dev.off() dev.off() ?mtext plot(esaligna$ht~esaligna$dap, xlab="Dap (cm)", ylab="Altura (m)", , bty="l", tcl=0.3, main="Relação entre Dap e Altura", cex.lab=1.5, cex.axis=1.5) savePlot #5.2 setwd("C:/Users/Camila Hohlenwerger/Desktop/aula4.exercicio_aula5") esaligna<- read.table(file="esaligna.csv", header=TRUE, sep=",", dec=".", as.is=TRUE) esaligna par(mfrow=c(1,2)) boxplot(esaligna$dap~esaligna$talhao) mtext("a", side=4, line=0.9, at=20,cex=2, family="serif") esaligna$talhao<- as.factor(esaligna$talhao) #mudou pra fator dentro da tabela. mediatalhao<-tapply(X=esaligna$dap, INDEX=esaligna$talhao, FUN=mean) sd<- tapply(X=esaligna$dap, INDEX=esaligna$talhao, FUN=sd) sd matrix<- matrix(c(mediatalhao, sd), ncol=2) matrix mediatalhao plot(matrix[,1], xlab="talhão", ylab="dap", ylim=c(5,25), xlim=c(1,6)) segments(x0=c(1:6) , y0=mediatalhao-sd, x1=c(1:6), y1=mediatalhao+sd) mtext("b", side=4, line=0.9, at=20,cex=2, family="serif") #5.3 codigo<- read.table(file="exercicio3.csv", header=TRUE, sep=",", dec=".", as.is=TRUE) head(codigo) codigo3<-codigo[,3:4] codigo3 #x e y do bloxplot codigo2<-na.omit(codigo [,1:2]) codigo2 #x e y do plot par(mfrow=c(1,2)) plot(codigo2$y1~codigo2$x1, xlab="Log(Patch size)(ha)",bty="l", ylab="Euclidean distances", ylim=c(0, 3.0), xlim=c(0.5, 2.3), tcl=0.3, pch=17, family="serif") modelo<- lm(codigo2$y1~codigo2$x1) summary(modelo) #abline(modelo, lwd=2.0) segments(x0=0.5, y0= -0.5173*0.5+2.4594 , x1=2.5, y1=-0.5173*2.5+2.4594, lwd=2.0 ) text(x=2.5, y=3.0, "a", bty="l", cex=1, family="serif") #x2= y, x1= x par(bty="l") par(tcl=0.3) boxplot(codigo3$x2~codigo3$y2,ylim=c(0, 3.0), outline=FALSE, names= c("Small", "Medium Edge", "Medium Interior", "Large Edge", "Large Interior", "Control")) text(x=6, y=3.0, "b", family="serif") text(x=1, y=3.0, "*", family="serif") text(x=2, y=3.0, "*", family="serif") text(x=3, y=3.0, "**", family="serif") text(x=4, y=3.0, "*", family="serif") text(x=5, y=3.0, "***", family="serif")