Uninformative crash with low read depth samples
Opened this issue · 0 comments
From Aatzli Coria:
Here is the command I entered:
python /data/data2/Programs/RingMapper-master/pairmapper.py --modified_parsed Pipeline_Modified_XXXX_parsed.mut --untreated_parsed Pipeline_Untreated_XXXX_parsed.mut --profile Pipeline_XXXX_profile.txt --out XXXX
Here is the error I got:
python /data/data2/Programs/RingMapper-master/pairmapper.py --modified_parsed Pipeline_Modified_XXXX_parsed.mut --untreated_parsed Pipeline_Untreated_XXXX_parsed.mut --profile Pipeline_XXXX_profile.txt --out XXXX
Traceback (most recent call last):
File "/data/data2/Programs/RingMapper-master/pairmapper.py", line 569, in
profile.normalize(DMS=True)
File "/data/data2/Programs/RingMapper-master/ReactivityProfile.py", line 521, in normalize
ac_nfac, ac_nerr = self.norm90( nprof[acmask] )
File "/data/data2/Programs/RingMapper-master/ReactivityProfile.py", line 704, in norm90
bnds = np.percentile(finitedata, [90., 99.]) # 99
File "/usr/lib/python2.7/dist-packages/numpy/lib/function_base.py", line 4274, in percentile
interpolation=interpolation)
File "/usr/lib/python2.7/dist-packages/numpy/lib/function_base.py", line 4016, in _ureduce
r = func(a, **kwargs)
File "/usr/lib/python2.7/dist-packages/numpy/lib/function_base.py", line 4391, in _percentile
x1 = take(ap, indices_below, axis=axis) * weights_below
File "/usr/lib/python2.7/dist-packages/numpy/core/fromnumeric.py", line 134, in take
return _wrapfunc(a, 'take', indices, axis=axis, out=out, mode=mode)
File "/usr/lib/python2.7/dist-packages/numpy/core/fromnumeric.py", line 57, in _wrapfunc
return getattr(obj, method)(*args, **kwds)
IndexError: cannot do a non-empty take from an empty axes.