#!/usr/bin/env python # vim: tw=0 from microarray import DataSet inputs = [ 'pickles/control.1/000.pkl', 'pickles/control.1/030.pkl', 'pickles/control.1/060.pkl', 'pickles/control.1/180.pkl' ] experiments = [ DataSet.restore(input) for input in inputs ] header = '{0.path}' template = '{0.id}\t{0.normed_ratio}' DataSet.tabulate(header, template, *experiments)
#!/usr/bin/env python from __future__ import division import math from microarray import DataSet pickles = ["pickles/A+D.000.pkl", "pickles/A+D.030.pkl", "pickles/A+D.060.pkl", "pickles/A+D.180.pkl"] def log_ratio(feature): return feature.log_ratio def too_extreme(feature): return abs(feature.log_ratio) > 15 header = "{0.path} (R/G = {0.intensity_ratio})" feature = "{0.id:<15} {0.log_ratio}" timepoints = [DataSet.restore(path) for path in pickles] for timepoint in timepoints: timepoint.prune(too_extreme) timepoint.sort(log_ratio, reverse=True) timepoint.truncate(50) DataSet.tabulate(header, feature, *timepoints)