Rarefaction curves showing combined samples

Is it possible to add a feature in future releases where rarefaction.single command takes the design file and outputs an average value for certain samples (i.e. combining treatment and control samples)?

You can kind of do this with summary.single(shared=, subsample=) and then use good old fashion statistics from there. Is that what you’re looking for?