yarden / MISO

MISO: Mixture of Isoforms model for RNA-Seq isoform quantitation
http://genes.mit.edu/burgelab/miso/index.html
132 stars 74 forks source link

summarize_miso error #109

Closed wxzheng001 closed 7 years ago

wxzheng001 commented 7 years ago

When I summarize the output, an IndexError occurs as below.What's the problem? ``summarize_miso --summarize-samples gene-psi-output summaries MISO (Mixture of Isoforms model) Summarize MISO output to get Psi values and confidence intervals. Use --help argument to view options.

Loading events from: /home/wxzheng/bin/myenv/lib/python2.7/site-packages/misopy/test-output/gene-psi-output Writing summary to: /home/wxzheng/bin/myenv/lib/python2.7/site-packages/misopy/test-output/summaries/summary/gene-psi-output.miso_summary Traceback (most recent call last): File "/home/wxzheng/bin/myenv/bin/summarize_miso", line 11, in load_entry_point('misopy==0.5.3', 'console_scripts', 'summarize_miso')() File "/home/wxzheng/bin/myenv/lib/python2.7/site-packages/misopy/summarize_miso.py", line 93, in main use_compressed=use_compressed) File "/home/wxzheng/bin/myenv/lib/python2.7/site-packages/misopy/samples_utils.py", line 309, in summarize_sampler_results output_fields = format_credible_intervals(event_name, samples) File "/home/wxzheng/bin/myenv/lib/python2.7/site-packages/misopy/credible_intervals.py", line 22, in format_credible_intervals cred_interval = compute_credible_intervals(samples) File "/home/wxzheng/bin/myenv/lib/python2.7/site-packages/misopy/credible_intervals.py", line 54, in compute_credible_intervals cred_interval = [samples[lower_bound_indx], samples[upper_bound_indx]] IndexError: only integers, slices (:), ellipsis (...), numpy.newaxis (None) and integer or boolean arrays are valid indices