diff repmatch_gff3.py @ 0:d33030c8e2cc draft

Uploaded

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/repmatch_gff3.py	Tue Nov 17 14:26:08 2015 -0500
@@ -0,0 +1,51 @@
+# repmatch.py
+#
+# Replicate matching - matches paired peaks from two or more replicates
+#
+# Input: one or more gff files (simple output from cwpair2, each a list of paired peaks from a replicate
+#
+# Output: list of matched groups and list of unmatched orphans
+# Files: key.tabular (file to replicate IDsummary.tabular, detail.tabular, orphans.tabular
+
+import argparse
+import repmatch_gff3_util
+
+if __name__ == '__main__':
+    parser = argparse.ArgumentParser()
+    parser.add_argument('--input', dest='inputs', action='append', nargs=2, help="Input datasets")
+    parser.add_argument('--method', dest='method', default='closest', help='Method of finding match')
+    parser.add_argument('--distance', dest='distance', type=int, default=50, help='Maximum distance between peaks in different replicates to allow merging')
+    parser.add_argument('--step', dest='step', type=int, default=0, help='Step size of distance for each iteration')
+    parser.add_argument('--replicates', dest='replicates', type=int, default=2, help='Minimum number of replicates that must be matched for merging to occur')
+    parser.add_argument('--low_limit', dest='low_limit', type=int, default=-1000, help='Lower limit for c-w distance filter')
+    parser.add_argument('--up_limit', dest='up_limit', type=int, default=1000, help='Upper limit for c-w distance filter')
+    parser.add_argument('--output_files', dest='output_files', default='simple', help='Restrict output dataset collections.')
+    parser.add_argument('--plot_format', dest='plot_format', default=None, help='Output format for graph')
+    parser.add_argument('--output_summary', dest='output_summary', help='Matched groups in gff format')
+    parser.add_argument('--output_orphan', dest='output_orphan', default=None, help='Orphans in tabular format')
+    parser.add_argument('--output_detail', dest='output_detail', default=None, help='Details in tabular format')
+    parser.add_argument('--output_key', dest='output_key', default=None, help='Keys in tabular format')
+    parser.add_argument('--output_histogram', dest='output_histogram', default=None, help='Histogram in plot_format')
+
+    args = parser.parse_args()
+
+    dataset_paths = []
+    hids = []
+    for (dataset_path, hid) in args.inputs:
+        dataset_paths.append(dataset_path)
+        hids.append(hid)
+    repmatch_gff3_util.process_files(dataset_paths,
+                                     hids,
+                                     args.method,
+                                     args.distance,
+                                     args.step,
+                                     args.replicates,
+                                     args.up_limit,
+                                     args.low_limit,
+                                     args.output_files,
+                                     args.plot_format,
+                                     args.output_summary,
+                                     args.output_orphan,
+                                     args.output_detail,
+                                     args.output_key,
+                                     args.output_histogram)