Skip to content
Snippets Groups Projects
Commit efc98ccc authored by Jeremy Mathews's avatar Jeremy Mathews
Browse files

redo replicate

parent 59962d33
Branches
Tags
1 merge request!54Resolve "Single control or experiment fails at call peaks macs"
Pipeline #4388 failed with stages
in 4 seconds
...@@ -229,14 +229,14 @@ def generate_design(paired, cutoff_ratio, design_df, cwd, no_reps, no_unique_con ...@@ -229,14 +229,14 @@ def generate_design(paired, cutoff_ratio, design_df, cwd, no_reps, no_unique_con
design_new_df.at[3, 'tag_align'] = design_new_df.at[0, 'tag_align'] design_new_df.at[3, 'tag_align'] = design_new_df.at[0, 'tag_align']
# Make 2 self psuedoreplicates # Make 2 self psuedoreplicates
self_pseudoreplicates_dict = {} pool_pseudoreplicates_dict = {}
for rep, tag_file in zip(design_df['replicate'], design_df['tag_align']): for rep, tag_file in zip(design_df['replicate'], design_df['tag_align']):
replicate_prefix = experiment_id + '_' + str(rep) replicate_prefix = experiment_id + '_' + str(rep)
self_pseudoreplicates_dict = \ pool_pseudoreplicates_dict = \
self_psuedoreplication(tag_file, replicate_prefix, paired) self_psuedoreplication(tag_file, replicate_prefix, paired)
# Update design to include new self pseudo replicates # Update design to include new self pseudo replicates
for rep, pseudorep_file in self_pseudoreplicates_dict.items(): for rep, pseudorep_file in pool_pseudoreplicates_dict.items():
path_to_file = cwd + '/' + pseudorep_file path_to_file = cwd + '/' + pseudorep_file
replicate = rep + 1 replicate = rep + 1
design_new_df.loc[replicate, 'tag_align'] = path_to_file design_new_df.loc[replicate, 'tag_align'] = path_to_file
...@@ -262,6 +262,7 @@ def generate_design(paired, cutoff_ratio, design_df, cwd, no_reps, no_unique_con ...@@ -262,6 +262,7 @@ def generate_design(paired, cutoff_ratio, design_df, cwd, no_reps, no_unique_con
replicate_prefix = experiment_id + '_' + str(rep) replicate_prefix = experiment_id + '_' + str(rep)
pr_dict = self_psuedoreplication(tag_file, replicate_prefix, paired) pr_dict = self_psuedoreplication(tag_file, replicate_prefix, paired)
pseudoreplicates_dict[rep] = pr_dict pseudoreplicates_dict[rep] = pr_dict
# Merge self psuedoreplication for each true replicate # Merge self psuedoreplication for each true replicate
pseudoreplicates_df = pd.DataFrame.from_dict(pseudoreplicates_dict) pseudoreplicates_df = pd.DataFrame.from_dict(pseudoreplicates_dict)
pool_pseudoreplicates_dict = {} pool_pseudoreplicates_dict = {}
......
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