[miso-users] Error with MISO

[Erin Wissink]

Hello!

I've been trying to run MISO on my data, but I keep getting error messages
when I try to find differential events. I'm using Python 2.7 and samtools
1.3.1. The messages I get when I run compare_miso are:

Retrieving MISO files in sample directories...

Computing sample comparison between R2_C_0hr_SE_single and
R2_C_1hr_SE_single...

  - No. of events in R2_C_0hr_SE_single: 3928

  - No. of events in R2_C_1hr_SE_single: 3391

Creating comparisons parent directory:
comparisons/R2_C_0hr_SE_single_vs_R2_C_1hr_SE_single

Traceback (most recent call last):

  File "/u/home/e/ewissink/.local/bin/compare_miso", line 9, in <module>

    load_entry_point('misopy==0.5.3', 'console_scripts', 'compare_miso')()

  File
"/u/home/e/ewissink/.local/lib/python2.7/site-packages/misopy/compare_miso.py",
line 94, in main

    use_compressed=use_compressed)

  File
"/u/home/e/ewissink/.local/lib/python2.7/site-packages/misopy/hypothesis_test.py",
line 289, in output_samples_comparison

    confidence_level=alpha)

  File
"/u/home/e/ewissink/.local/lib/python2.7/site-packages/misopy/credible_intervals.py",
line 22, in format_credible_intervals

    cred_interval = compute_credible_intervals(samples)

  File
"/u/home/e/ewissink/.local/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
Do you know what's causing this error and how I can fix it?

Thanks!
Erin
-- 
Erin Wissink, PhD
Postdoctoral Fellow, UCLA
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.mit.edu/pipermail/miso-users/attachments/20170313/7f8a68c6/attachment.html