I need to make sure that I get a posterior of region-wide richness, and that I don't average over the iterations prior to that.
This is somewhat cumbersome if trying to get site-specific and region-wide richness, because the iteration-averaging always occurs at the last step!
E.g., for some Stan code I had previously taken each species' average (across posterior samples) chance of being present in the region, then summed across those chances to get richness. But I really should have summed across species within each iteration so that I could get a distribution of region-wide richness.
I need to make sure that I get a posterior of region-wide richness, and that I don't average over the iterations prior to that.
This is somewhat cumbersome if trying to get site-specific and region-wide richness, because the iteration-averaging always occurs at the last step!
E.g., for some Stan code I had previously taken each species' average (across posterior samples) chance of being present in the region, then summed across those chances to get richness. But I really should have summed across species within each iteration so that I could get a distribution of region-wide richness.