def _merge(ks, outpath): k = Kyub.merge(*ks) def _mkTC(hu): p = hu.split('=') h = p[0] u = p[1] if len(p) > 1 else None return Column(h, iterable=k.get_treatment_column(hu), units=u) def _mkDC(h): return Column(h, iterable=k.get_column(h)) tc = [_mkTC(f) for f in k.factors] dc = [_mkDC(f) for f in k.readouts] datapflex.write_datapflex(outpath, tc, dc)
k = Kyub.merge(*ks) def _mkTC(hu): p = hu.split('=') h = p[0] u = p[1] if len(p) > 1 else None return Column(h, iterable=k.get_treatment_column(hu), units=u) def _mkDC(h): return Column(h, iterable=k.get_column(h)) # def _mkDC(h): # m = k.get_column(h) # s = k.get_column(h + '=stdev') # return MSColumn(h, iterable=zip(m, s)) tc = [_mkTC(f) for f in k.factors] dc = [_mkDC(f) for f in k.readouts] #dc = [_mkDC(f) for f in k.readouts if not f.endswith('=stdev')] datapflex.write_datapflex(outpath, tc, dc) # for p in sorted(glob('scans/linkfarm/*')): # assay = op.basename(p) # for z in ('GF', 'CK'): # f0, f1 = fs = [op.join(d, '%s_%s.csv' % (assay, z)) for d in dpdirs] # if not op.exists(f1): # continue # diffs, dkyub = diff(f0, f1) # out = []