Skip to content
GitLab
Projects
Groups
Snippets
/
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
Blaise LI
bioinfo_utils
Commits
51089311
Commit
51089311
authored
Jan 09, 2020
by
Blaise Li
Browse files
More info in warning logs.
parent
88623070
Changes
1
Hide whitespace changes
Inline
Side-by-side
small_RNA-seq/small_RNA-seq.snakefile
View file @
51089311
...
@@ -2788,6 +2788,12 @@ rule merge_bigwig_reps:
...
@@ -2788,6 +2788,12 @@ rule merge_bigwig_reps:
warn(f"{msg}:\n")
warn(f"{msg}:\n")
#raise
#raise
warn(f"The bigwig files without {chrom} will be skipped.\n")
warn(f"The bigwig files without {chrom} will be skipped.\n")
except AssertionError as e:
chrom_lens = []
for bw in bws:
chrom_lens.append(bw.chroms()[chrom])
warn(f"chrom_len for {chrom}: {chrom_len}\nchrom_lens: {chrom_lens}")
raise
to_use = [bw for bw in bws if chrom in bw.chroms()]
to_use = [bw for bw in bws if chrom in bw.chroms()]
if to_use:
if to_use:
means = np.nanmean(np.vstack([bw.values(chrom, 0, chrom_len, numpy=True) for bw in to_use]), axis=0)
means = np.nanmean(np.vstack([bw.values(chrom, 0, chrom_len, numpy=True) for bw in to_use]), axis=0)
...
...
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment