richness function

Estimate posterior distributions of site richness