outfiles_list = []
for d in batch_parameters:
dirname = c.experiment_name + '_clustered_reads'
outfile = c.experiment_name + '_clustered_reads'
if 'c_thresh' in d:
dirname = dirname + '-c{}'.format(int(d['c_thresh']*100))
outfile = outfile + '-c{}'.format(int(d['c_thresh']*100))
if 'n_filter' in d:
dirname = dirname + '-n{}'.format(d['n_filter'])
outfile = outfile + '-n{}'.format(d['n_filter'])
if 'maskN' in d:
dirname = dirname + '-maskN'
outfile = outfile + '-maskN'
path = os.path.join(c.clusters_outpath, dirname)
outfiles_list.append(os.path.join(path, outfile))
for cluster_file in outfiles_list:
name = os.path.split(cluster_file)[1]
out6 = cluster_summary_plot(cluster_file, plot_hist = 0)
hist_counter(out6[1], bins=5000, range =(1,10000),label=name)
''' Created on 28 Feb 2013 @author: musselle ''' import os import sys import numpy as np from plot_utils import cluster_summary_plot, hist_counter infileL6 = '/space/musselle/datasets/gazelles-zebras/clusters/L6clustered_reads' infileL8 = '/space/musselle/datasets/gazelles-zebras/clusters/L8clustered_reads' outL6 = cluster_summary_plot(infileL6, plot_hist = 0) outL8 = cluster_summary_plot(infileL8, plot_hist = 0) hist_counter(outL6[1], bins=5000, range =(1,10000),label='Lane 6') hist_counter(outL8[1], bins=5000, range =(1,10000),label='Lane 8')