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

remove errors from generate_design function

parent e56be52c
Branches
Tags
1 merge request!54Resolve "Single control or experiment fails at call peaks macs"
Pipeline #4371 failed with stages
in 7 seconds
...@@ -322,22 +322,6 @@ def generate_design(paired, cutoff_ratio, design_df, cwd, no_reps, no_unique_con ...@@ -322,22 +322,6 @@ def generate_design(paired, cutoff_ratio, design_df, cwd, no_reps, no_unique_con
tmp_metadata['tag_align'] = path_to_file tmp_metadata['tag_align'] = path_to_file
design_new_df = design_new_df.append(tmp_metadata) design_new_df = design_new_df.append(tmp_metadata)
# Add in pool experiment
tmp_metadata['sample_id'] = experiment_id + '_pooled' else:
path_to_pool_control = cwd + '/' + pool_control
design_new_df['control_tag_align'] = path_to_pool_control
# Add in pseudo replicates
tmp_metadata = design_new_df.loc[0].copy()
tmp_metadata['control_tag_align'] = path_to_pool_control
for rep, pseudorep_file in pool_pseudoreplicates_dict.items():
tmp_metadata['sample_id'] = experiment_id + '_pr' + str(rep)
tmp_metadata['replicate'] = str(rep) + '_pr'
tmp_metadata['xcor'] = 'Calculate'
path_to_file = cwd + '/' + pseudorep_file
tmp_metadata['tag_align'] = path_to_file
design_new_df = design_new_df.append(tmp_metadata)
# Add in pool experiment # Add in pool experiment
tmp_metadata['sample_id'] = experiment_id + '_pooled' tmp_metadata['sample_id'] = experiment_id + '_pooled'
tmp_metadata['replicate'] = 'pooled' tmp_metadata['replicate'] = 'pooled'
...@@ -346,16 +330,6 @@ def generate_design(paired, cutoff_ratio, design_df, cwd, no_reps, no_unique_con ...@@ -346,16 +330,6 @@ def generate_design(paired, cutoff_ratio, design_df, cwd, no_reps, no_unique_con
tmp_metadata['tag_align'] = path_to_file tmp_metadata['tag_align'] = path_to_file
design_new_df = design_new_df.append(tmp_metadata) design_new_df = design_new_df.append(tmp_metadata)
# Write out new dataframe
design_new_df.to_csv(experiment_id + '_ppr.tsv',
header=True, sep='\t', index=False)
tmp_metadata['replicate'] = 'pooled'
tmp_metadata['xcor'] = 'Calculate'
path_to_file = cwd + '/' + pool_experiment_se
tmp_metadata['tag_align'] = path_to_file
design_new_df = design_new_df.append(tmp_metadata)
return design_new_df return design_new_df
......
...@@ -5,7 +5,6 @@ import pandas as pd ...@@ -5,7 +5,6 @@ import pandas as pd
from io import StringIO from io import StringIO
import os import os
import pool_and_psuedoreplicate import pool_and_psuedoreplicate
import sys
test_output_path = os.path.dirname(os.path.abspath(__file__)) + \ test_output_path = os.path.dirname(os.path.abspath(__file__)) + \
'/../output/design/' '/../output/design/'
......
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