Ejemplo n.º 1
0
#!/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)

Ejemplo n.º 2
0
#!/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)