Plotting a rarefaction curve using mean species richness in R

I have bird community matrices (presence-absence) for 10 sites (with 72 samples each) and would like to plot a single averaged sample-based rarefaction curve for all these sites. At the moment I've merely averaged individual estimates (for each site) generated by the specaccum function, but am unsure whether this is a sensible approach.