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)
sub.sample(shared=18Sfcm.pick.pick.filter.opti_mcc.unique_list.0.02.pick.shared, size=89233)
merge.groups(shared=18Sfcm.pick.pick.filter.opti_mcc.unique_list.0.02.pick.0.02.subsample.shared, design=15x.design)
count.groups(shared=18Sfcm.pick.pick.filter.opti_mcc.unique_list.0.02.pick.0.02.subsample.merge.shared)
dist.shared(shared=18Sfcm.pick.pick.filter.opti_mcc.unique_list.0.02.pick.0.02.subsample.merge.shared, calc=braycurtis, output=lt)
clearcut(phylip=18Sfcm.pick.pick.filter.opti_mcc.unique_list.0.02.pick.0.02.subsample.merge.braycurtis.0.02.lt.dist, verbose=t)
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)
amova(phylip=18Sfcm.pick.pick.filter.opti_mcc.unique_list.0.02.pick.0.02.subsample.merge.braycurtis.0.02.lt.dist,group=15x.design)
It send an error claiming that my design file does not include ‘ACC’, etc.
Design file:
group categorie
12_0207.15 SI
12_0409.15 SI
12_0509.15 SI
12_0202.15 OW
12_0505.15 OW
12_0701.15 OW
13_0502.15 ACC
13_0701.15 ACC
13_0902.15 ACC
13_0201.15 WASW
13_1107.15 WASW
13_1909.15 WASW
13_0808.15 BS
13_0909.15 BS
13_1314.15 BS
Please, if anybody can give me guidance on what can be going wrong I will greatly appreciate it.
Again, thanks,
Carla