INNER CODE UNIT · Python
main
getzlab/rnaseqc · test_data/approx_diff.py:7
def main(args):
if args.mode == 'tables':
df = pd.read_csv(args.input1, index_col=0, header=2, sep='\t').join(
pd.read_csv(args.input2, index_col=0, header=2, sep='\t'),
how='outer',
rsuffix='_'
)
elif args.mode == 'metrics' or args.mode == 'fragments':
df = pd.read_csv(args.input1, sep='\t', index_col=0).join(
pd.read_csv(args.input2, sep='\t', index_col=0),
how='outer',
rsuffix='_',
)
assert not (df[args.columns[0]].isna() ^ df[args.columns[1]].isna()).any(), df[df.isna().any(1)]
assert len(df[np.abs(df[args.columns[0]] - df[args.columns[1]]) > args.tolerance]) == 0, df[np.abs(df[args.columns[0]] - df[args.columns[1]]) > args.tolerance].head()