Skip to content

Uninformative crash with low read depth samples #3

Description

@anthonymustoe

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.

Activity

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

      Milestone

      No milestone

      Relationships

      None yet

      Development

      No branches or pull requests

      Issue actions