Estimadores de riqueza no R

Estão disponíveis diversos pacotes e estimadores de riqueza no R, cada um com sua vantagem e desvantagem. Mas, para facilitar, desenvolvi uma função para calcular, de uma única vez, os estimadores Jackknife 1 e 2, ICE, ACE, Bootstrap, e Chao 1 e 2, além de plotar um gráfico para facilitar a visualização e comparação dos resultados.

A função criada se chama “bio_estimate” (basta copiar e colar em um script no R):

Para utilizar, basta executar a função com o título (entre aspas) que deseja para o gráfico, dessa forma:

Como resultado, será gerado o gráfico:

O gráfico também mostra a riqueza amostrada, além dos demais estimadores citados. Para que a função seja executada corretamente são necessários os pacotes ggplot2 e vegan.

Quem sabe essa função não vire um pacote um dia?

Dúvidas? Visite nosso site ou mande um e-mail para viniciusbrbio@gmail.com.