#5.1 #leyendo esaligna esaligna<-read.table("C:\\Users\\USER\\Desktop\\curso R\\Exe3\\3.6\\esaligna.csv",header=T,as.is=T,sep=",") esaligna str(esaligna) #creando un grafico de dispersao plot(esaligna$dap,esaligna$ht,xlab="Diametro altura peito(cm)",ylab="altura(m)") par(cex.axis=1.5) #aumenta os carateres par(las=1,tcl=0.3,bty="L")#marcadores dentro del grafico plot(esaligna$dap,esaligna$ht,main="Diagrama de dispersion",xlab="Diametro altura peito(cm)",ylab="altura(m)") #################################################### #5.2 #use as variaves dap e talhao para fazer 2grafi par(mfrow=c(1,2)) boxplot(esaligna$dap~esaligna$talhao) text(16,20,"a") men<-aggregate(esaligna$dap,by=list(talhao=esaligna$talhao),mean) men colnames(men)[2]<-"dap medio" men sd<-aggregate(esaligna$dap,by=list(talhao=esaligna$talhao),sd) colnames(sd)[2]<-"Desvio padro" sd #junt<-matrix(c(mean,sd),ncol=2,nrow=6) junt barplot(datos$men,datos$sd) abline(sd) ###################################### #5.3 dad<-read.table("exercicio3.csv",header=T,sep=",",as.is=T) dad names(dad) is.na(dad) is.na(dad$x1) is.na(dad$y1) is.na(dad$x2) is.na(dad$x2) #guardando em temp1 los na temp1<-dad[is.na(dad[,1])|is.na(dad[,2]),] dad$x1[is.na(dad$x1)]<-0 dad$y1[is.na(dad$y1)]<-0 dad# puso 0 donde era NA #calculando distancia euclidiana dist<-sqrt((dad$x1-dad$x2)^2+(dad$y1-dad$y2)^2) dist #tamano<-(abs(dad$x1-dad$x2)*abs(dad$y1-dad$y2))*10000 #tamano #log<-log(tamano) #log datos<-data.frame(dad,dist) range(datos$x1) par(mfrow=c(1,2),bty="l",family="serif",cex.lab=1.5,cex.axis=1.3,las=1,tcl=0.3) plot(datos$x1,datos$dist,xlab="Log(Patch size)(ha)",ylab="Euclideam distance",cex=1.5,pch=17,xaxp=c(0.5,2,3),ylim=c(0,3)) model<-lm(datos$x1~dist) abline(model,lwd=2) text(2.5,3,"a",cex=1.8) boxplot(datos$x1,datos$dist) text(3,3,"a",cex=1.8)