function sp.averages to create CT results table

with mean, median and mode of estimated species richness for each site