dap<-runif(100,10,80) e <- exp(1) especies<-rep(c("especie1","especie2","especie3","especie4","especie5","especie6","especie7","especie8","especie9","especie10"),times=1,each=10) biomassa<-(e^-1.7953)*(dap^2.2974) log(biomassa)-((-2.6464)+(1.9960*log(dap)))/0.7558->x altura<-e^x factor(especies) prefixo<-"ind" sufixo<-1:100 individuos<-paste(prefixo,sufixo,sep="") dados<-data.frame(especies,dap,biomassa,altura,row.names=individuos) dados.decr<-dados[order(-dap),] t1<-dados$dap<=50 dap50<-table(t1)["TRUE"] t2<-dados$dap>=50&dados$dap<=60 dap50.60<-table(t2)["TRUE"] t3<-dados$dap<=30|dados$dap>=70 dap30.70<-table(t3)["TRUE"]