caixeta <- read.csv("caixeta.csv", as.is=T) names(caixeta) table(caixeta$especie) hist(caixeta$cap/pi) agr=aggregate(caixeta$h, by=list(local=caixeta$local,arvores=caixeta$arvore), FUN=mean) ##Árvores com mais de um fuste são agregadas pela altura média, mantendo-se a informação do local. chauas=agr[grep("chauas", agr$local),] jureia=agr[grep("jureia", agr$local),] retiro=agr[grep("retiro", agr$local),] par(mfrow=c(3,1)) #Visualização dos três histogramas na mesma página hist(chauas$x) #Frequência de árvores por altura em chauas hist(jureia$x) #Frequência de árvores por altura em jureia hist(retiro$x) #Frequência de árvores por altura em retiro par(mfrow=c(1,1)) #Visualização de gráficos para um por página. #Há diferenças entre as estruturas (distribuição de tamanhos) dos caixetais? summary(chauas) summary(jureia) summary(retiro) par(mfrow=c(3,1)) ##Ajusta a visualização dos três histogramas na mesma página. hist(chauas$x, ylim=c(0,0.015), xlim=c(0,500), prob=T, breaks=21) ##breaks é necessário para igualar a divisão de barras dos outros gráficos. lines(density(chauas$x), col="blue") ##Ajuste de densidade probabilística. curve(expr = dnorm(x,mean=mean(chauas$x), sd=sd(chauas$x)), add=T, col="red") ##Curva normal. hist(jureia$x, ylim=c(0,0.015), xlim=c(0,500), prob=T) #Ajuste o eixo x com o eixo do gráfico anterior. lines(density(jureia$x), col="blue") ##Ajuste de densidade probabilística. curve(expr = dnorm(x,mean=mean(jureia$x), sd=sd(jureia$x)), add=T, col="red") ##Curva normal. hist(retiro$x, ylim=c(0,0.015), xlim=c(0,500), prob=T) #Ajuste o eixo x com o eixo do gráfico de chauas. lines(density(retiro$x), col="blue") ##Ajuste de densidade probabilística. curve(expr = dnorm(x,mean=mean(retiro$x), sd=sd(retiro$x)), add=T, col="red") ##Curva normal. par(mfrow=c(1,1)) ##Ajusta a visualização dos três histogramas na mesma página. ##Apesar de apresentarem médias diferentes e medianas não discrepantes em relação às primeiras, os locais avaliados não apresentaram grandes alterações na estrutura.