Just discovered another bug… I think… It seems that get.oturep has problems dealing with missing labels, unlike most other functions.
Here’s what I did to make get.oturep crash:
mothur > get.oturep(fasta=IVSCS_QC_R.fasta, list=IVSCS_QC_R.nr.aligned.fn.list, label=unique-0.03-0.05-0.15-0.30-0.40)
unique 1148
0.03 402
0.05 275
0.15 84
0.30 15
Segmentation fault
get.sharedotu doesn’t have this problem:
mothur > get.sharedotu(list=IVSCS_QC_R.nr.aligned.fn.list, group=IVSCS_QC_R.groups, fasta=IVSCS_QC_R.fasta, label=unique-0.03-0.05-0.15-0.30-0.40)
unique
0.03
0.05
0.15
0.30
Your file does not include the label 0.40. I will use 0.38.
0.38
I can fix get.oturep by specifying a label that does exist:
mothur > get.oturep(fasta=IVSCS_QC_R.fasta, list=IVSCS_QC_R.nr.aligned.fn.list, label=unique-0.03-0.05-0.15-0.30-0.38)
unique 1148
0.03 402
0.05 275
0.15 84
0.30 15
0.38 3
Should be a relatively easy fix since get.sharedotu already deals with this correctly. Thanks!