Repository navigation
Expand file tree
/
Copy pathmatplotlib_venn_example.py
More file actions
16 lines (10 loc) · 15.5 KB
/
Copy pathmatplotlib_venn_example.py
File metadata and controls
16 lines (10 loc) · 15.5 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
from matplotlib import pyplot as plt
import numpy
import scipy
from matplotlib_venn import venn2, venn2_circles
Expected_Scans = set([ 5, 6, 82, 83, 85, 86, 88, 89, 91, 92, 94, 95, 97, 98, 100, 101, 103, 104, 106, 107, 109, 110, 112, 113, 115, 116, 118, 119, 121, 122, 124, 125, 127, 128, 130, 131, 133, 134, 136, 137, 139, 140, 142, 143, 145, 146, 148, 149, 151, 152, 154, 155, 157, 158, 160, 161, 163, 164, 166, 167, 169, 170, 172, 173, 175, 176, 178, 179, 181, 182, 184, 185, 187, 188, 190, 191, 193, 194, 196, 197, 199, 200, 202, 203, 205, 206, 208, 209, 211, 212, 214, 215, 217, 218, 220, 221, 223, 224, 226, 227, 229, 230, 232, 233, 235, 236, 238, 239, 241, 242, 244, 245, 247, 248, 250, 251, 253, 254, 256, 257, 259, 260, 262, 263, 265, 266, 268, 269, 271, 272, 274, 275, 277, 280, 282, 283, 285, 286, 288, 289, 291, 292, 294, 295, 297, 298, 300, 301, 303, 304, 306, 307, 309, 310, 312, 313, 315, 316, 318, 319, 321, 322, 324, 325, 327, 329, 330, 332, 333, 335, 336, 338, 339, 341, 342, 344, 345, 347, 348, 350, 351, 353, 354, 356, 357, 359, 360, 362, 363, 365, 366, 368, 369, 371, 372, 374, 375, 377, 378, 380, 381, 383, 384, 386, 387, 389, 390, 392, 393, 395, 396, 398, 399, 401, 402, 404, 405, 407, 408, 410, 411, 413, 414, 416, 417, 419, 420, 422, 423, 425, 426, 428, 429, 431, 432, 435, 436, 438, 439, 441, 443, 444, 446, 447, 449, 450, 452, 453, 455, 456, 458, 459, 461, 462, 464, 465, 467, 468, 470, 471, 473, 474, 476, 478, 479, 481, 482, 484, 486, 487, 489, 491, 492, 494, 498, 499, 502, 504, 505, 507, 508, 510, 511, 513, 514, 516, 517, 519, 520, 522, 523, 525, 526, 528, 530, 533, 535, 537, 539, 542, 545, 547, 551, 558, 560, 561, 563, 565, 566, 569, 571, 573, 576, 578, 579, 581, 582, 584, 585, 587, 589, 592, 594, 595, 599, 604, 607, 608, 610, 611, 613, 615, 617, 618, 620, 621, 624, 625, 629, 634, 638, 639, 641, 642, 644, 647, 649, 650, 652, 653, 655, 656, 658, 659, 661, 664, 665, 667, 668, 673, 675, 677, 679, 680, 682, 683, 685, 687, 689, 690, 692, 694, 695, 697, 699, 700, 702, 704, 706, 707, 709, 711, 714, 717, 719, 721, 723, 724, 726, 727, 729, 731, 732, 734, 736, 737, 739, 741, 743, 744, 747, 748, 750, 753, 754, 757, 758, 760, 762, 763, 769, 771, 772, 774, 778, 779, 781, 782, 784, 785, 787, 788, 790, 792, 793, 795, 796, 798, 802, 803, 805, 807, 809, 810, 812, 814, 816, 817, 820, 821, 823, 824, 827, 832, 834, 835, 837, 838, 843, 844, 845, 847, 848, 850, 851, 853, 855, 858, 862, 863, 865, 867, 868, 870, 872, 875, 876, 878, 879, 884, 886, 888, 889, 891, 893, 894, 896, 898, 899, 901, 903, 904, 906, 907, 909, 910, 912, 913, 915, 916, 918, 919, 921, 923, 924, 926, 929, 930, 932, 933, 935, 936, 938, 940, 941, 943, 944, 946, 947, 949, 950, 952, 953, 955, 956, 958, 959, 961, 962, 964, 965, 967, 968, 970, 972, 973, 975, 976, 978, 979, 981, 982, 984, 985, 987, 988, 990, 991, 995, 1001, 1003, 1004, 1006, 1008, 1011, 1013, 1015, 1017, 1019, 1021, 1026, 1028, 1029, 1031, 1032, 1034, 1035, 1038, 1039, 1041, 1042, 1044, 1045, 1047, 1048, 1050, 1052, 1054, 1056, 1057, 1059, 1061, 1063, 1065, 1067, 1068, 1070, 1071, 1073, 1074, 1076, 1077, 1079, 1080, 1082, 1083, 1085, 1086, 1088, 1089, 1091, 1092, 1094, 1095, 1097, 1098, 1100, 1101, 1103, 1104, 1106, 1107, 1109, 1110, 1112, 1114, 1116, 1117, 1119, 1120, 1122, 1125, 1126, 1128, 1130, 1136, 1137, 1139, 1140, 1142, 1143, 1145, 1146, 1148, 1149, 1151, 1152, 1154, 1155, 1157, 1158, 1160, 1161, 1163, 1164, 1166, 1167, 1169, 1170, 1172, 1173, 1175, 1176, 1178, 1179, 1181, 1182, 1184, 1185, 1187, 1188, 1190, 1191, 1193, 1195, 1196, 1198, 1202, 1204, 1207, 1210, 1212, 1213, 1215, 1217, 1219, 1221, 1226, 1228, 1231, 1234, 1238, 1241, 1243, 1245, 1247, 1249, 1251, 1253, 1255, 1257, 1258, 1261, 1263, 1266, 1267, 1269, 1270, 1272, 1273, 1275, 1276, 1278, 1279, 1283, 1295, 1296, 1298, 1301, 1303, 1305, 1307, 1311, 1320, 1327, 1329, 1331, 1332, 1334, 1335, 1337, 1338, 1340, 1341, 1343, 1344, 1346, 1347, 1350, 1352, 1353, 1355, 1361, 1363, 1364, 1366, 1367, 1369, 1370, 1372, 1374, 1376, 1377, 1379, 1381, 1383, 1385, 1386, 1388, 1389, 1391, 1392, 1394, 1395, 1397, 1398, 1400, 1401, 1403, 1404, 1406, 1408, 1410, 1411, 1413, 1414, 1416, 1418, 1420, 1422, 1423, 1425, 1426, 1428, 1429, 1431, 1432, 1434, 1435, 1438, 1440, 1443, 1449, 1453, 1454, 1456, 1458, 1460, 1462, 1463, 1465, 1466, 1468, 1469, 1471, 1472, 1474, 1476, 1481, 1483, 1485, 1487, 1488, 1490, 1491, 1493, 1495, 1496, 1499, 1504, 1505, 1507, 1508, 1510, 1511, 1513, 1514, 1516, 1517, 1519, 1520, 1522, 1523, 1525, 1526, 1528, 1529, 1531, 1532, 1534, 1535, 1538, 1539, 1542, 1547, 1548, 1550, 1552, 1553, 1555, 1557, 1558, 1560, 1565, 1567, 1568, 1572, 1578, 1580, 1581, 1585, 1586, 1588, 1593, 1595, 1599, 1601, 1602, 1604, 1605, 1607, 1608, 1610, 1611, 1613, 1615, 1616, 1618, 1619, 1621, 1622, 1624, 1625, 1627, 1628, 1630, 1631, 1633, 1634, 1636, 1637, 1639, 1645, 1650, 1651, 1660, 1664, 1667, 1673, 1676, 1678, 1687, 1696, 1701, 1706, 1709, 1710, 1713, 1714])
X_found = set([116,101,104,100,103,711,101,104,711,909,711,907,910,907,910,924,923,329,332,935,1017,478,484,625,83,83,137,137,599,1017,599,251,251,692,695,1056,1056,1076,1056,624,736,561,571,1275,1279,1337,1275,1279,1092,1092,94,97,1181,1148,1151,1154,1148,1151,1154,1038,1041,795,798,86,655,658,238,241,244,655,658,378,795,798,378,769,459,1103,1340,1103,467,461,404,1103,679,679,679,679,426,432,426,432,697,1676,1709,1301,1303,578,519,528,107,607,1195,607,372,374,1267,1267,199,955,1267,1267,1272,510,516,510,673,675,573,823,594,731,262,262,1651,98,89,92,98,98,98,450,1471,89,92,1212,89,92,1142,1139,1283,1283,1283,1283,1142,1142,1418,1182,1184,1397,1400,1520,1522,1523,1526,1520,1522,1523,1526,1305,1397,1400,514,517,1261,1261,1334,436,1266,1363,1367,1608,1519,1398,1401,581,587,581,587,306,300,304,1045,1048,1065,1068,300,1247,1249,1504,1504,462,462,1191,1191,1091,1091,1604,1466,1468,176,1460,1462,1463,1465,1466,1468,1460,1462,1463,1465,1633,1631,1631,1714,1548,1217,566,579,1098,1100,1420,1607,1621,1624,1622,151,155,1599,116,478,484,1525,1507,1510,1525,1525,478,484,625,1017,254,599,599,625,629,748,692,748,889,889,695,1056,889,1667,1664,1667,1664,1664,1667,1667,1097,1097,380,721,1092,1092,1021,1021,1181,1148,1151,1154,1181,795,798,769,679,1103,1340,464,1272,1031,268,656,1355,1212,396,560,393,289,292,889,795,798,286,285,286,1103,1103,1709,1466,1468,1460,1462,1463,1465,1466,1468,1714,1714,1714,1608,1181,1195,1143,482,924,923,907,910,936,1077,1077,157,1565,1713,322,1532,1714,1283,375,1182,1184,1340,107,1516,1516,1516,1519,1257,510,1337,723,1337,262,1181,697,1241,719,124,124,1453,1531,1531,1531,482,384,384,396,1516,130,128,5,1337,106,109,116,101,104,1355,1352,1352,1355,106,109,125,125,110,110,1095,1110,1471,748,1103,1103,1625,1628,719,1625,1628,1625,1628,912,915,1471,172,375,1182,1184,1182,1184,1182,1184,1182,1184,1182,1184,143,130,127,134,128,128,468,128,1295,1602,1601,206,345,345,310,309,1251,257,769,823,1026,1029,1031,1028,1295,1520,1522,1523,1526,118,115,112,112,1092,729,726,947,950,656,656,119,113,112,113,316,313,309,316,288,286,285,446,1067,1070,330,333,316,319,309,1608,286,285,288,919,307,306,301,282,335,341,285,330,333,303,203,330,333,307,306,1678,1567,1120,1119,301,282,335,341,203,1567,285,307,306,285,303,301,1713,291,294,203,301,1678,275,1676,1709,1713,307,306,1120,1119,1567,303,291,294,1391,304,304,1713,300,306,303,304,1329,304,300,306,304,300,1079,1074,1079,1079,1713,300,304,304,300,300,471,304,300,322,100,103,107,297,322,295,297,938,940,929,932,297,295,280,295,297,148,184,188,193,200,197,190,1347,1338,1344,1714,295,297,297,295,1052,1054,1044,1047,1032,1034,1032,1034,101,104,683,1026,1029,1031,1028,387,390,629,106,109,112,110,930,1056,1157,1160,1169,1172,1157,1160,1169,1172,447,1157,1160,1340,1469,1469,525,197,611,145,145,148,145,1122,1125,677,510,1337,1550,1552,1637,1604,1613,1615,1611,1593,1636,1637,1613,1615,736,695,1619,1613,1615,1611,1619,1613,1615,1593,1627,1613,1615,1627,1627,226,230,1627,1630,1634,1611,1613,1615,1619,1627,1627,1631,1633,930,1631,1630,1634,1633,1627,94,97,1631,1630,1634,1633,1631,1633,1631,1631,1633,1633,1395,1406,413,417,1660,1667,732,1528,1531,459,455,444,446,446,429,1170,901,899,190,193,187,200,197,145,145,145,148,149,1152,1155,1181,970,1181,455,461,459,452,582,410,414,110,1196,1198,1196,1226,1228,462,470,461,459,455,479,1166,417,413,410,414,404,1094,1101,1103,1125,1080,1083,479,467,464,471,470,468,1558,1560,154,154,152,151,155,158,154,1645,1639,128,152,154,151,155,154,1013,1031,1028,151,155,452,1063,82,329,332,321,324,325,327,329,332,298,1519,904,906,904,906,891,893,913,916,904,906,149,285,289,292,298,140,146,137,1388,1398,1401,1391,1519,285,289,292,83,659,655,658,661,131,137,145,133,136,139,140,146,149,976,981,982,985,988,145,157,145,145,143,130,133,136,139,140,146,131,148,137,149,1059,1057,1065,1068,128,127,131,1418,1414,1416,131,148,699,680,682,702,695,692,148,152,847,850,805,785,787,779,782,832,844,778,781,1238,1243,1278,1272,1276,1267,1272,1276,824,790,820,1508,970,967,929,932,912,915,926,933,921,116,83,82,384,386,398,407,389,392,987,990,384,393,386,389,392,384,245,232,236,239,238,241,244,242,251,254,85,88,91,86,83,1495,1035,1039,248,904,906,938,940,913,916,608,615,85,88,91,86,618,621,248,416,419,422,431,461,423,479,470,436,438,439,441,359,435,420,428,458,425,416,419,422,449,537,547,563,423,431,248,251,1678,1676,85,88,91,86,401,405,404,359,435,416,419,422,420,401,405,426,432,426,432,426,432,1376,1379,1381,1383,1413,1410,1376,1379,1381,1383,1413,85,88,91,450,453,456,467,464,107,113,426,432,107,106,109,1136,1140,107,224,214,217,209,196,209,199,263,259,245,253,256,257,254,251,247,250,263,489,492,478,484,487,491,510,516,507,511,494,513,502,505,508,814,851,452,834,837,814,851,1178,1185,1176,1178,814,851,507,511,502,505,508,504,522,1145,1149,1173,1175,1173,1175,254,1555,1557,1555,1557,344,1210,1219,968,1618,823,238,241,244,248,1397,1400,1403,1152,1155,85,88,91,181,185,184,188,178,182,191,194,181,185,85,88,91,1434,1422,1426,1420,1422,1443,1411,1434,1397,1400,1403,1411,1385,1374,1377,1370,1372,1420,1422,731,734,1420,1422,1188,1443,1440,1438,1454,351,357,338,347,353,348,342,329,332,89,92,95,98,94,97,202,200,89,92,94,97,98,95,94,97,410,414,89,92,95,94,97,98,94,97,107,89,92,95,98,206,89,92,95,98,1425,1423,1429,1432,1429,368,369,371,380,372,374,1429,1432,94,97,89,92,1604,1608,94,97,89,92,128,127,143,142,149,130,137,128,128,507,511,504,522,1483,1485,176,184,188,178,182,179,191,194,1431,1428,191,194,179,184,188,178,182,285,283,354,377,360,351,357,362,365,378,359,435,356,363,347,353,354,377,384,360,329,332,116,122,121,121,116,122,121,122,1074,464,467,471,470,462,465,271,724,717,727,731,734,736,747,744,757,739,743,731,734,271,122,218,220,212,223,229,217,214,206,215,221,1307,1311,1298,1305,1301,1303,1307,266,289,292,263,280,280,274,277,265,269,280,274,277,298,205,206,200,1261,265,262,266,257,263,1305,1301,1303,827,802,809,823,805,779,785,787,782,803,812,1019,1017,1021,1021,1017,1021,1019,426,432,1019,1021,1017,408,1021,1019,1021,1017,1529,1017,1021,1019,411,408,1019,260,253,256,247,250,251,257,254,232,236,239,247,250,233,235,253,256,247,250,338,330,333,336,339,329,332,1513,1517,921,919,903,918,342,360,351,357,356,363,347,353,350,344,354,377,366,350,336,1513,1517,721,709,723,736,699,702,697,719,706,714,709,704,721,732,704,673,675,679,697,680,315,318,312,272,265,271,268,187,175,176,744,1534,1532,747,760,744,750,128,1532,1534,375,372,374,380,847,850,858,862,865,867,462,452,1535,468,476,462,474,468,452,468,476,473,481,487,491,474,486,468,449,453,456,450,455,452,443,449,452,1095,1110,1082,1086,1528,1311,1320,351,357,347,353,336,339,342,344,1346,1343,903,918,868,878,888,898,891,893,921,924,876,863,879,884,868,878,494,513,566,579,516,561,571,494,513,530,533,479,470,478,484,498,510,482,467,462,389,392,498,396,566,579,542,576,393,561,571,530,533,479,470,478,484,510,482,573,516,535,545,396,387,390,558,408,395,494,513,464,467,498,566,579,401,405,447,393,566,579,561,571,464,467,470,530,533,478,484,516,510,401,405,396,395,408,494,513,462,479,447,498,482,396,396,380,387,390,447,408,494,513,479,467,470,478,484,498,482,395,464,396,380,387,390,489,492,610,617,624,638,647,641,625,629,396,625,629,624,523,523,504,522,561,571,592,592,625,629,641,641,169,170,173,566,579,604,625,581,587,629,1050,594,611,566,579,561,571,585,589,578,625,599,573,537,547,563,565,560,525,523,519,528,504,522,507,511,585,604,589,629,599,560,561,571,578,507,511,525,519,528,561,571,523,504,522,578,561,571,584,604,589,578,599,585,584,607,629,823,816,604,1604,232,236,239,233,235,1056,1038,1041,1044,1047,1056,1045,1048,1251,1221,1234,1076,1077,1065,1068,1181,1076,1056,1182,1184,1038,1041,1044,1047,1163,1167,1056,1056,233,235,1003,1011,1006,685,664,685,667,1221,1215,1234,1258,1255,1261,1247,1249,1221,1234,1221,1213,1261,401,405,399,402,398,407,401,405,410,414,1161,1158,446,450,218,220,205,206,208,211,215,221,217,214,1337,1338,1340,1341,1340,208,211,215,221,205,206,1616,1263,224,223,229,224,232,236,239,227,223,229,215,221,226,230,1152,1155,1148,1151,1154,196,202,200,197,193,191,194,1148,1151,1154,348,607,1332,1335,224,1107,1109,196,827,841,834,837,838,848,1119,823,848,760,737,771,762,760,747,753,737,771,762,741,754,757,758,848,870,1097,1106,1095,1110,1091,1104,1098,1100,1112,1089,1114,1122,1120,1071,1073,1091,1089,1095,1110,1106,1114,1112,1098,1100,1104,661,649,653,655,658,652,661,655,658,1091,1106,1095,1110,1104,1098,1100,1114,639,685,665,665,649,653,668,655,658,679,652,1095,1110,1092,1085,1088,1103,561,571,520,526,530,533,514,517,510,502,505,508,516,450,453,456,1471,446,1092,1085,1088,1092,1085,1088,926,933,736,748,711,689,889,907,910,711,748,729,689,719,736,697,889,907,910,692,736,689,711,692,695,1706,697,719,729,924,923,909,907,910,912,915,889,695,706,714,692,729,711,889,909,907,910,924,923,689,1709,1466,1468,1460,1462,1463,1465,1471,706,714,729,1706,909,924,923,907,910,889,689,1456,1460,1462,1463,1465,1466,1468,1458,706,714,924,923,907,910,909,907,910,924,923,889,1460,1462,1463,1465,1466,1468,1474,1471,1474,909,638,647,650,641,924,923,558,641,930,936,935,936,935,936,1469,699,702,694,700,1469,1472,1622,1621,1624,1472,175,172,179,1622,1621,1624,160,166,161,164,163,167,157,158,175,172,175,158,160,166,151,155,163,157,161,164,263,257,160,166,161,164,163,387,390,380,387,390,378,381,383,384,1217,1191,1193,1202,1179,1185,1210,1219,1182,1184,1198,1196,1257,1251,1202,449,991,984,972,978,982,985,988,979,975,309,976,981,972,978,979,975,995,982,985,988,984,991,972,978,170,173,1610,1604,1607,1607,1604,1605,1139,1142,1158,1143,1139,1142,1143,1142,1139,1139,1142,1143,820,824,941,949,938,940,962,1015,952,944,943,946,967,955,903,918,919,938,940,896,1152,1155,788,779,784,792,816,785,787,782,807,1599,1602,1601,1052,1054,1602,1601,1599,785,787,1146,1145,1149,1148,1151,1154,953,956,958,965,961,964,959,973,970,1182,1184,1295,1181,1212,1298,1301,1303,1195,1701,1251,1350,1353,1257,1696,1231,1334,1195,1065,1068,1035,1039,1212,1195,814,851,613,620,617,611,638,647,624,610,608,615,618,621,639,607,617,613,620,584,595,814,851,812,868,878,868,878,855,827,823,772,863,795,798,769,835,817,737,771,778,781,875,821,832,844,853,888,898,903,918,758,772,736,769,762,737,771,781,778,795,798,741,754,769,741,754,762,737,771,772,737,771,769,1269,1278,1278,1283,1270,1273,1269,1266,1267,1275,1279,1278,1238,1243,1266,1269,1267,1247,1249,1245,1253,1270,1273,1296,1550,1552,1493,1496,1374,1377,1366,1369,1364,1367,1363,1370,1372,638,647,649,653,642,644,1332,1385,1386,1529,1528,1389,1520,1522,1523,1526,1525,1532,1520,1522,1523,1526,1519,1507,1510,1514,1525,1519,1520,1522,1523,1526,1525,1514,1516,1519,1507,1510,1525,1520,1522,1523,1526,1520,1522,1523,1526,1519,1516,1507,1510,1488,1490,1520,1522,1523,1526,1520,1522,1523,1526,1163,1167,1166,1608,1334,1329,1331,1163,1167,1212,1181,1173,1175,1187,1164,1215,1204,1166,1163,1167,1394,1392,1408,1404,1397,1400,1391,1398,1401,1350,1353,1664,1097,941,949,827,446])
inter = Expected_Scans.symmetric_difference (X_found)
print (inter)
venn2([Expected_Scans,X_found],["Expected Scans","XTRACT Targets"])
plt.show()