Skip to content

Commit

Permalink
Fix ac.tl.nucleosome_signal (#105)
Browse files Browse the repository at this point in the history
  • Loading branch information
gtca committed Jun 7, 2023
1 parent 19f98d2 commit 0484ace
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions muon/_atac/tools.py
Original file line number Diff line number Diff line change
Expand Up @@ -1153,14 +1153,14 @@ def nucleosome_signal(

for i in tqdm(range(n), desc="Reading Fragments"):
try:
f = fr.next()
f = next(fr)
length = f.end - f.start
row_ind = d[f.name]
if length < nucleosome_free_upper_bound:
mat[row_ind, 0] += 1
elif length < mononuleosomal_upper_bound:
mat[row_ind, 1] += 1
except:
except KeyError:
pass
# if i % 1000000 == 0:
# print(f"Read {i/1000000} Mio. fragments.", end='\r')
Expand Down

0 comments on commit 0484ace

Please sign in to comment.