diff table-union.py @ 2:8c5fb0c5e560

planemo upload
author jpayne
date Fri, 26 Jan 2018 16:38:41 -0500
parents f1f2497301d3
children 4eaafbdfb8bf
line wrap: on
line diff
--- a/table-union.py	Mon Jan 08 11:33:27 2018 -0500
+++ b/table-union.py	Fri Jan 26 16:38:41 2018 -0500
@@ -8,13 +8,40 @@
 def main(files):
 	header = []
 	items = []
+	possible_identity_headers = None
 	for fi in files:
 		with open(fi, 'rU') as table:
-			rows = csv.DictReader(table, delimiter='\t', dialect='excel-tab')
-			for field in rows.fieldnames:
+			reader = csv.DictReader(table, delimiter='\t', dialect='excel-tab')
+			rows = list(reader)
+			for field in reader.fieldnames:
 				if field not in set(header):
 					header.append(field)
+				#try to find identity columns in the files, to use to join
+				if possible_identity_headers is None:
+					possible_identity_headers = set(reader.fieldnames)
+				#winnow down the shared columns in each file by whether they're present in all, and all their values are unique in each file and not null
+				#because these are the most likely to be shared keys
+				possible_identity_headers = possible_identity_headers.intersection(filter(lambda f: len(set([r[f] for r in rows])) == len(rows) and all([r[f] is not None for r in rows]), reader.fieldnames))
 			items.extend(rows)
+	
+	# if len(possible_identity_headers) > 1:
+	# 	#if there's more than one, we need to check that joining on any one of them produces the same results
+
+	# 	#finally
+	# 	possible_identity_headers = set((possible_identity_headers.pop(), ))
+
+	#if we found an identity column, then try to join rows
+	if possible_identity_headers:
+		key_column = possible_identity_headers.pop()
+		keys = set([r[key_column] for r in items])
+		merged_rows = []
+		for key in sorted(keys):
+			new_row = {}
+			for row in filter(lambda r: r[key_column] == key, items):
+				new_row.update(row)
+			merged_rows.append(new_row)
+		items = merged_rows
+
 	wr = csv.DictWriter(sys.stdout, delimiter='\t', dialect='excel-tab', fieldnames=header)
 	wr.writeheader()
 	wr.writerows(items)