<div dir="ltr">Hello!<div><br></div><div>I&#39;ve been trying to run MISO on my data, but I keep getting error messages when I try to find differential events. I&#39;m using Python 2.7 and samtools 1.3.1. The messages I get when I run compare_miso are:</div><div><br><div>Retrieving MISO files in sample directories...<br></div><div>
<p class="gmail-p1"><span class="gmail-s1">Computing sample comparison between R2_C_0hr_SE_single and R2_C_1hr_SE_single...</span></p>
<p class="gmail-p1"><span class="gmail-s1">  - No. of events in R2_C_0hr_SE_single: 3928</span></p>
<p class="gmail-p1"><span class="gmail-s1">  - No. of events in R2_C_1hr_SE_single: 3391</span></p>
<p class="gmail-p1"><span class="gmail-s1">Creating comparisons parent directory: comparisons/R2_C_0hr_SE_single_vs_R2_C_1hr_SE_single</span></p>
<p class="gmail-p1"><span class="gmail-s1">Traceback (most recent call last):</span></p>
<p class="gmail-p1"><span class="gmail-s1">  File &quot;/u/home/e/ewissink/.local/bin/compare_miso&quot;, line 9, in &lt;module&gt;</span></p>
<p class="gmail-p1"><span class="gmail-s1">    load_entry_point(&#39;misopy==0.5.3&#39;, &#39;console_scripts&#39;, &#39;compare_miso&#39;)()</span></p>
<p class="gmail-p1"><span class="gmail-s1">  File &quot;/u/home/e/ewissink/.local/lib/python2.7/site-packages/misopy/compare_miso.py&quot;, line 94, in main</span></p>
<p class="gmail-p1"><span class="gmail-s1">    use_compressed=use_compressed)</span></p>
<p class="gmail-p1"><span class="gmail-s1">  File &quot;/u/home/e/ewissink/.local/lib/python2.7/site-packages/misopy/hypothesis_test.py&quot;, line 289, in output_samples_comparison</span></p>
<p class="gmail-p1"><span class="gmail-s1">    confidence_level=alpha)</span></p>
<p class="gmail-p1"><span class="gmail-s1">  File &quot;/u/home/e/ewissink/.local/lib/python2.7/site-packages/misopy/credible_intervals.py&quot;, line 22, in format_credible_intervals</span></p>
<p class="gmail-p1"><span class="gmail-s1">    cred_interval = compute_credible_intervals(samples)</span></p>
<p class="gmail-p1"><span class="gmail-s1">  File &quot;/u/home/e/ewissink/.local/lib/python2.7/site-packages/misopy/credible_intervals.py&quot;, line 54, in compute_credible_intervals</span></p>
<p class="gmail-p1"><span class="gmail-s1">    cred_interval = [samples[lower_bound_indx], samples[upper_bound_indx]]</span></p>
<p class="gmail-p1"><span class="gmail-s1">IndexError: only integers, slices (`:`), ellipsis (`...`), numpy.newaxis (`None`) and integer or boolean arrays are valid indices</span></p><div>Do you know what&#39;s causing this error and how I can fix it?</div><div><br></div><div>Thanks!</div><div>Erin</div>-- <br><div class="gmail_signature"><div dir="ltr">Erin Wissink, PhD<div>Postdoctoral Fellow, UCLA</div></div></div>
</div></div></div>