Hey dear mothur community,
I hope somebody can help me realize where the error is. Hope enough detail is provided.
I am using the last version to analyze V4 18S MiSeq samples, based on MiSeq SOP. I have 15 samples, group in 5 sets of three samples each. I want to compare OW/SI and ACC/BS/WASW, but after merging my design file is not recognized correctly and mothur crashes.
After getting my occurrence table (create.database) I subsampled the shared file and merged groups according to my design file to run AMOVA and Unifrac.
Here is how my commands look like:
make.shared(list=18Sfcm.pick.pick.filter.opti_mcc.unique_list.0.02.pick.list, count=18S5fcm.pick.pick.pick.count_table, label=0.02)
dist.shared(shared=18Sfcm.pick.pick.filter.opti_mcc.unique_list.0.02.pick.0.02.subsample.merge.shared, calc=braycurtis, output=lt)
When I try to run:
unifrac.unweighted(tree=18Sfcm.pick.pick.filter.opti_mcc.unique_list.0.02.pick.0.02.subsample.merge.braycurtis.0.02.lt.tre, random=F, groups=ACC-BS-WASW, group=15x.design)
It send an error claiming that my design file does not include ‘ACC’, etc.
Please, if anybody can give me guidance on what can be going wrong I will greatly appreciate it.