Commit c34d2612 authored by Venkat Malladi's avatar Venkat Malladi

Move logic ahead of redo for blacklist.

parent f5c0ad5c
Pipeline #7235 passed with stages
in 494 minutes and 49 seconds
......@@ -242,14 +242,15 @@ def main():
for experiment, df_experiment in design_peaks_df.groupby('experiment_id'):
replicated_peak, chr_peak = overlap(experiment, df_experiment)
design_diff.loc[design_diff.experiment_id == experiment, "peak"] = replicated_peak
design_anno.loc[experiment] = [experiment, chr_peak]
# Remove blacklist regions; if blacklist = True
if blacklist and os.path.exists(blacklist):
bl_peaks, bl_chr_peak = blacklist_peaks(experiment, blacklist)
design_diff.loc[design_diff.experiment_id == experiment, "peak"] = bl_peaks
design_anno.loc[design_anno.Condition == experiment, "Peaks"] = bl_chr_peak
else:
design_diff.loc[design_diff.experiment_id == experiment, "peak"] = replicated_peak
design_anno.loc[experiment] = [experiment, chr_peak]
# Write out file
design_diff.columns = ['SampleID',
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment