Mercurial > repos > mheinzl > hd
comparison hd.py @ 30:46bfbec0f9e6 draft
planemo upload for repository https://github.com/monikaheinzl/duplexanalysis_galaxy/tree/master/tools/hd commit 033dd7b750f68e8aa68f327d7d72bd311ddbee4e-dirty
author | mheinzl |
---|---|
date | Wed, 07 Aug 2019 04:01:32 -0400 |
parents | 6b15b3b6405c |
children | 8beced3064e3 |
comparison
equal
deleted
inserted
replaced
29:6b15b3b6405c | 30:46bfbec0f9e6 |
---|---|
1078 minHD_tags = numpy.concatenate([item[4] for item in diff_list_a]) | 1078 minHD_tags = numpy.concatenate([item[4] for item in diff_list_a]) |
1079 minHD_tags_zeros1 = numpy.concatenate([item[7] for item in diff_list_a]) | 1079 minHD_tags_zeros1 = numpy.concatenate([item[7] for item in diff_list_a]) |
1080 minHD_tags_zeros2 = numpy.concatenate([item[7] for item in diff_list_b]) | 1080 minHD_tags_zeros2 = numpy.concatenate([item[7] for item in diff_list_b]) |
1081 | 1081 |
1082 chimera_tags1 = sum([item[10] for item in diff_list_a], []) | 1082 chimera_tags1 = sum([item[10] for item in diff_list_a], []) |
1083 chimera_tags2 = numpy.concatenate([item[10] for item in diff_list_b]) | 1083 chimera_tags2 = sum([item[10] for item in diff_list_b], []) |
1084 | 1084 |
1085 rel_Diff = [] | 1085 rel_Diff = [] |
1086 diff_zeros = [] | 1086 diff_zeros = [] |
1087 minHD_tags_zeros = [] | 1087 minHD_tags_zeros = [] |
1088 diff = [] | 1088 diff = [] |
1089 chimera_tags = [] | 1089 chimera_tags = [] |
1090 for d1, d2, rel1, rel2, zeros1, zeros2, tag1, tag2, ctag1, ctag2 in \ | 1090 for d1, d2, rel1, rel2, zeros1, zeros2, tag1, tag2, ctag1, ctag2 in \ |
1091 zip(diff1, diff2, rel_Diff1, rel_Diff2, diff_zeros1, diff_zeros2, minHD_tags_zeros1, minHD_tags_zeros2, | 1091 zip(diff1, diff2, rel_Diff1, rel_Diff2, diff_zeros1, diff_zeros2, minHD_tags_zeros1, minHD_tags_zeros2, |
1092 chimera_tags1, chimera_tags2): | 1092 chimera_tags1, chimera_tags2): |
1093 rel_Diff.append(max(rel1, rel2)) | 1093 relatives = numpy.array([rel1, rel2]) |
1094 diff.append(max(d1, d2)) | 1094 absolutes = numpy.array([d1, d2]) |
1095 max_idx = numpy.argmax(relatives) | |
1096 rel_Diff.append(relatives[max_idx]) | |
1097 diff.append(absolutes[max_idx]) | |
1095 | 1098 |
1096 if all(i is not None for i in [zeros1, zeros2]): | 1099 if all(i is not None for i in [zeros1, zeros2]): |
1097 diff_zeros.append(max(zeros1, zeros2)) | 1100 diff_zeros.append(max(zeros1, zeros2)) |
1098 minHD_tags_zeros.append(str(tag1)) | 1101 minHD_tags_zeros.append(str(tag1)) |
1099 tags = [ctag1, ctag2] | 1102 tags = [ctag1, ctag2] |