I am new to this topic. Please don’t hesitate to tell me what’s what.
I am more-or-less following the MiSeq SOP (http://www.mothur.org/wiki/MiSeq_SOP) for 16S (Illumina V3 chemistry). I have generated my phylogenetic tree with clearcut, and I would like to compare each of my groups with each other group in the sequencing run without treatments labelling those groups. The Mothur wiki page on design files does not seem to cover how to generate a design file with this property (http://www.mothur.org/wiki/Design_File), and perhaps it is not possible.
What I have tried
One idea I had was formatting the design file with a unique treatment label for each group (as below), but I am not sure sure I trust the output of the weighted unifrac given that all the wscores=1. Although the scores were all equal, I noticed that the significance scores were not. I am not sure whether this is valid, or whether it is plausible for the weighted fraction of branch lengths to be entirely unshared.
What I would like to learn
- How to use UniFrac functions to compare each of my sequencing groups
- Whether a weighted score of 1 in every group comparison indicates that I did something wrong given that I want to compare each of my sequencing groups to each other without grouping into treatments.