From 411c86d7a494e0b6244ba5dc2502a902b9883731 Mon Sep 17 00:00:00 2001
From: gpratt <gpratt@gmail.com>
Date: Wed, 27 Jan 2016 12:08:54 -0800
Subject: [PATCH] first commit of parse clip analysis

---
 clipper/src/parse_clip_analysis.py | 28 ++++++++++++++++++++++++++++
 1 file changed, 28 insertions(+)
 create mode 100644 clipper/src/parse_clip_analysis.py

diff --git a/clipper/src/parse_clip_analysis.py b/clipper/src/parse_clip_analysis.py
new file mode 100644
index 0000000..954f745
--- /dev/null
+++ b/clipper/src/parse_clip_analysis.py
@@ -0,0 +1,28 @@
+__author__ = 'gpratt'
+
+import pandas as pd
+import cPickle as pickle
+import os
+import numpy as np
+import argparse
+
+
+def output_region_density(fn):
+    with open(fn) as f_handle:
+        result = pickle.load(f_handle)
+    name = os.path.basename(fn).split(".")[0]
+    regions_dict = {}
+    for region in result['distributions']:
+        region_type_dict = {}
+        for region_type in ['total', 'individual']:
+            count, bins = np.histogram(result['distributions'][region][region_type], bins=100, normed=True)
+            region_type_dict[region_type] = {key: value for key, value in zip(np.arange(0,1, .01), count)}
+        regions_dict[region] = pd.DataFrame(region_type_dict).T
+    pd.concat(regions_dict).to_csv(name + ".distributions.csv")
+
+if __name__ == '__main__':
+    parser = argparse.ArgumentParser(description="Takes a pickle file from the output of clip_analysis and outputs the region distributions as csv in histogram form")
+    parser.add_argument("-i", "--input", help="pickle file from clip analysis to parse")
+    args = parser.parse_args()
+
+    output_region_density(args.input)