From c2c1d8533ddf94e9fc0af2ac5b92482907244156 Mon Sep 17 00:00:00 2001 From: Lingfei Wang Date: Fri, 18 Aug 2023 16:25:11 -0400 Subject: [PATCH] . --- src/dictys/chromatin.py | 5 ++--- 1 file changed, 2 insertions(+), 3 deletions(-) diff --git a/src/dictys/chromatin.py b/src/dictys/chromatin.py index 7903d9a..b9c3997 100644 --- a/src/dictys/chromatin.py +++ b/src/dictys/chromatin.py @@ -183,14 +183,13 @@ def _motif_postproc(dret,fi_exp:str,fo_bed:str,fo_wellington:str,fo_homer:str)-> t1=set(namet) namem=[x.split('_') for x in namem] namem=['_'.join([','.join(list(filter(lambda y:y in t1,x[0].split(','))))]+x[1:]) for x in namem] - if len(namem)!=len(set(namem)): - raise ValueError('Found non-unique motif name suffices. Each motif name is recommended to contain a unique suffix.') - dh.columns=dw.columns=namem namem=np.array(namem) #Remove motifs having no TF in current dataset t1=[not x.startswith('_') for x in namem] dw,dh=[x[:,t1] for x in [dw,dh]] namem,=[x[t1] for x in [namem]] + if len(namem)!=len(set(namem)): + raise ValueError('Found non-unique motif name suffices. Each motif name is recommended to contain a unique suffix.') assert dw.shape==(len(namep),len(namem)) and dh.shape==(len(namep),len(namem)) assert np.isfinite(dw).all() and np.isfinite(dh).all() assert (dw>=0).all() and (dh>=0).all()