To determine whether the sampling was sufficient, I plotted rarefaction curves for four sites using rarecurve() function in R vegan package. Seems like the curves have not yet "converged" (i.e. reached asymptote), and each site has different sample size. I would like to add fittings on the curves, so that I can comment on how similar the sites are in terms of their need for further sampling. I assume then I would also be able to comment on which site is subset of which.
One way to do this I guess was explained in this paper (http://www.wsl.ch/info/mitarbeitende/wohlgemu/pdf/koellner_etal2004.pdf), but I could not quite understand the approach here.
I appreciate any suggestions.