dados.esaligna<-read.csv("esaligna.csv",header=TRUE,sep=",")
str(dados.esaligna)
head(dados.esaligna)
summary(dados.esaligna)
novacoluna1 <- dados.esaligna$tronco + dados.esaligna$folha
dados.esaligna2 <- data.frame(dados.esaligna,novacoluna1)
novacoluna2 <- pi*((dados.esaligna2$dap/2)^2)
dados.esaligna.fin <- data.frame(dados.esaligna2,novacoluna2)
dados.esaligna.10cm <- dados.esaligna.fin[dados.esaligna$dap>10,]
area.basal <- by(dados.esaligna.10cm$novacoluna2,dados.esaligna.10cm$talhao,FUN=sum)
media.area.basal <- by(dados.esaligna.10cm$novacoluna2,dados.esaligna.10cm$talhao,FUN=mean)
write.table (c(dados.esaligna.10cm, area.basal, media.area.basal), file="novoarquivo.txt",row.names=T, sep="\t")