How to run core_diversity_py after using multiple_split_libraries.py (What mapping file to use?)

16 views
Skip to first unread message

Mgall00

unread,
Jul 29, 2016, 5:10:30 AM7/29/16
to Qiime 1 Forum
Hello,

Because of my already demultiplexed data, I split up my mapping file up and made one for each fastq file. multiple_split_libraries.py worked, and I proceeded (as in the illumina tutorial) with pick_open_reference_otus.py

However, I am not sure how to get the rarefaction curves or the plots from the tutorial, because core_diversity_analyses.py requires a (single?) mapping file. As mentioned, I have one for each fastq. I tried giving the script my master-mapping file, the one containing info for all before I split it up, but that failed. Any idea of what I could do?

Here is the error raised, by the way, but it seems pretty clear that it's just not happy with the mapping file.

qiime@qiime-190-virtual-box:~/Desktop/Shared_Folder/16S_Roskilde_2016/joined_reads/LABPrimer$ core_diversity_analyses.py -o cdout/ -i otus/otu_table_mc2_w_tax_no_pynast_failures.biom -m MappingOverview.tsv  -t otus/rep_set.tre -e 23751
/usr/local/lib/python2.7/dist-packages/skbio/stats/ordination/_principal_coordinate_analysis.py:107: RuntimeWarning: The result contains negative eigenvalues. Please compare their magnitude with the magnitude of some of the largest positive eigenvalues. If the negative ones are smaller, it's probably safe to ignore them, but if they are large in magnitude, the results won't be useful. See the Notes section for more details. The smallest eigenvalue is -0.0100882321604 and the largest is 1.02269957114.
  RuntimeWarning
Traceback (most recent call last):
  File "/usr/local/bin/core_diversity_analyses.py", line 202, in <module>
    main()
  File "/usr/local/bin/core_diversity_analyses.py", line 199, in main
    status_update_callback=status_update_callback)
  File "/usr/local/lib/python2.7/dist-packages/qiime/workflow/core_diversity_analyses.py", line 252, in run_core_diversity_analyses
    status_update_callback=status_update_callback)
  File "/usr/local/lib/python2.7/dist-packages/qiime/workflow/downstream.py", line 183, in run_beta_diversity_through_plots
    close_logger_on_success=close_logger_on_success)
  File "/usr/local/lib/python2.7/dist-packages/qiime/workflow/util.py", line 122, in call_commands_serially
    raise WorkflowError(msg)
qiime.workflow.util.WorkflowError: 

*** ERROR RAISED DURING STEP: Make emperor plots, weighted_unifrac)
Command run was:
 make_emperor.py -i cdout//bdiv_even23751//weighted_unifrac_pc.txt -o cdout//bdiv_even23751//weighted_unifrac_emperor_pcoa_plot/ -m MappingOverview.tsv 
Command returned exit status: 2
Stdout:

Stderr
Error in make_emperor.py: None of your sample identifiers match between the mapping file and the coordinates file. Verify you are using a coordinates file and a mapping file that belong to the same dataset.

Reply all
Reply to author
Forward
0 new messages