Possible interfaces: - Add a `.groups` argument to `mutate_sample_data()` etc that functions like that in `dplyr::with_groups()` - Add a `with_groups()` method that works on phyloseq objects