Add function for determining the species richness within groups

This is a general function for doing things like:
1. Determining the richness at each site
2. Determining the richness at each site-year combination
3. Determining the richness within each genus