Hello -
I’ve run into a seemingly not-uncommon error when attempting to run the beta-group-significance command:
qiime diversity beta-group-significance
–i-distance-matrix Iceland-2018-1003-core-metrics/Iceland-2018-1003-unweighted-unifrac-distance-matrix.qza
–m-metadata-file Iceland-2018-1003-no-mock.tsv
–m-metadata-column group
–o-visualization no-mock-unweighted-unifrac-sector-significance.qzv
–p-pairwise
Plugin error from diversity:
All values in the grouping vector are unique. This method cannot operate on a grouping vector with only unique values (e.g., there are no ‘within’ distances because each group of objects contains only a single object).
Debug info has been saved to /tmp/qiime2-q2cli-err-c84pvu6j.log
After researching, it appears that my problem is that I have two samples in my data set that each belong to one group. Here are my samples and their groups:
sample-id | . barcode-sequence | . group |
---|---|---|
categorical. | categorical. | categorical |
Solo-spot2 | . TAGATGCCGTCC | . two |
Solo-spot3 | . CGCTCCTTCCTT | . two |
Solo-spot4 | . CGAGCGCCGTCC | . two |
Solo-spot5 | . TCGAGAGTCCGG | . three |
Solo-spot6 | . GCTCGGCCGGCC | . three |
Solo-spot8 | . AGCTCAGTCCGG | . three |
Solo-spot9 | . ACTGCCGCGTCC | . three |
Solo-spot10 | . TTACTTCTCCGG | . three |
Solo-spot12 | . TGACGCGCGGCC | . four |
Solo-spot16 | . AATGATCTCCGG | . one |
Solo-spot17 | . TTATCAGTCCTT | . one |
Solo-spot18 | . AAGCGTACGTCC | . one |
Mock1 | . TAAGAGAAGTCC | . mock1 |
As you can see, the Mock belongs to its own group (mock1), and Spot12 belongs to its own group (four). All well and good. I have run the following command to remove them from my data:
qiime feature-table filter-samples
–i-table Iceland-2018-1003-table-dada2.qza
–m-metadata-file 2018-1003-Solo-Metadata.tsv
–p-where “[group]=‘three’ AND [group]=‘four’ AND [group]=’one’”
–o-filtered-table Iceland-2018-1003-filtered-table.qza
My question now is, how do I incorporate that into beta diversity? As far as I can tell, beta diversity analysis doesn’t use the table.qza as input, and still needs the metadata file, which was causing the error in the first place.
I am running Qiime2 2020.2 on a cluster installed with conda.
Thanks much!