class TestDisplayResults(TestCaseSuppressOutput):
samples = list(get_samples(controls, sample_groups, qadata, matrix))
def test01(self):
display_results(self.samples, sys.stdout)
def test02(self):
display_results(self.samples, sys.stdout, style='file')
def test03(self):
display_results(self.samples, sys.stdout, grouped=False)
def action(args):
controls, sample_groups = get_input(args.infile, split_description = args.split_desc)
qadata = qa_from_csv(qafile)
matrix = read_matrix(matrix_file)
style = 'screen' if args.outfile == sys.stdout else 'file'
outfile = get_outfile(args, label = 'results', ext = 'csv')
if args.calculate_ion_ratios:
qadata = add_ion_ratios(qadata, controls)
samples = get_samples(controls, sample_groups, qadata, matrix)
display_results(samples, outfile, style = style, grouped = args.grouped)
if args.outfile is None:
outfile.close()
def action(args):
matrix = read_matrix(matrix_file)
qadata = qa_from_csv(qafile)
compound_ids, compound_codes = map(list, zip(*COMPOUND_CODES))
writer = csv.DictWriter(
args.outfile if hasattr(args.outfile, 'write') else open(args.outfile, 'w'),
fieldnames = ['infile','label'] + compound_ids,
extrasaction = 'ignore')
# create headers for the first row
d = dict(label = 'label', infile = 'infile')
d.update(dict(COMPOUND_CODES))
writer.writerow(d)
counter = count(1)
seen = defaultdict(list) # keep only the first instance of each accession
for infile in args.infiles:
d = dict(infile = infile)
controls, sample_groups = get_input(infile, split_description = args.split_desc)
if args.calculate_ion_ratios:
qadata = add_ion_ratios(qadata.copy(), controls)
patient_samples = get_samples(controls, sample_groups, qadata, matrix, quantitative = True)
for samples in patient_samples:
samples = list(samples)
# exclude any malformed samples
if any(any(c.malformed for c in sample.compounds.values()) for sample in samples):
log.warning('malformed samples in %s' % infile)
continue
# keep only the first instance of each sample
if samples[0].row_label() not in seen:
d['label'] = samples[0].row_label() if args.keep_phi else counter.next()
d.update(dict((s.COMPOUND_id, s.result) for s in samples))
writer.writerow(sanitize(d))
seen[samples[0].row_label()].append(infile)
def test01(self):
_qadata = add_ion_ratios(qadata, controls)
rows = get_samples(controls, sample_groups, _qadata, matrix)
rows.next()
def test01(self):
_qadata = add_ion_ratios(qadata, controls)
rows = get_samples(controls, sample_groups, _qadata, matrix)
rows.next()
