diff --git a/output/analyze/analyze_sum_matrix/1M/threshold-0.1/output_sum_matrix_1M_1CDX1_v1.txt b/output/analyze/analyze_sum_matrix/1M/threshold-0.1/output_sum_matrix_1M_1CDX1_v1.txt new file mode 100644 index 0000000..ceb6530 --- /dev/null +++ b/output/analyze/analyze_sum_matrix/1M/threshold-0.1/output_sum_matrix_1M_1CDX1_v1.txt @@ -0,0 +1,234 @@ +bin1 bin2 count +3 1542 51 +4 1542 62 +5 1542 52 +6 1542 60 +7 1542 37 +8 1542 35 +9 1542 50 +10 1542 54 +11 1542 34 +12 1542 40 +13 1542 36 +14 1542 29 +15 1542 40 +16 1542 40 +18 1542 33 +19 1542 31 +20 1542 29 +22 1542 31 +23 1542 28 +24 1542 29 +25 1542 31 +28 1542 33 +34 1542 28 +35 1542 28 +36 1542 28 +90 201 30 +90 1542 70 +111 964 38 +136 1199 30 +201 1542 61 +202 1542 40 +203 1542 30 +204 1542 40 +205 1542 32 +206 1542 37 +207 1542 34 +208 1542 44 +209 1542 50 +210 1542 45 +211 1542 38 +212 1542 43 +213 1542 28 +214 1542 36 +215 1542 32 +217 1542 32 +221 1542 30 +222 1542 28 +254 1359 36 +267 833 32 +383 1542 36 +384 1542 29 +385 1542 31 +386 1542 28 +387 1542 28 +388 1542 29 +389 1542 33 +394 1542 28 +518 570 31 +543 852 30 +543 1155 30 +543 1542 33 +545 1542 31 +546 1542 33 +547 1542 32 +548 1542 37 +549 1542 34 +550 1542 35 +551 1542 35 +554 1542 32 +555 1542 38 +557 1542 33 +559 1542 29 +677 738 30 +695 2364 34 +699 1542 31 +701 1542 32 +729 1542 82 +771 1525 40 +786 1366 33 +826 2257 31 +843 2272 47 +852 1542 75 +852 2272 105 +853 1542 57 +854 1542 50 +855 1542 45 +856 1542 30 +857 1542 42 +859 1542 30 +860 1542 28 +861 1542 35 +862 1542 36 +863 1542 37 +864 1542 34 +865 1542 28 +866 1542 35 +867 1542 30 +868 1542 29 +870 1542 35 +871 1542 31 +872 1542 28 +875 1542 31 +886 1688 45 +915 1585 31 +919 2447 39 +1002 1542 80 +1003 1542 56 +1004 1542 36 +1005 1542 49 +1011 1542 36 +1012 1542 45 +1013 1542 28 +1015 1542 35 +1016 1542 33 +1018 1542 30 +1019 1542 35 +1024 1542 28 +1025 1542 30 +1029 1542 33 +1134 2465 28 +1155 1542 93 +1155 2272 89 +1156 1542 68 +1157 1542 56 +1158 1542 54 +1159 1542 45 +1160 1542 51 +1161 1542 57 +1162 1542 44 +1163 1542 41 +1164 1542 37 +1165 1542 46 +1166 1542 40 +1167 1542 29 +1168 1542 32 +1169 1542 34 +1170 1542 44 +1175 1542 33 +1177 1542 35 +1178 1542 36 +1179 1542 31 +1180 1542 32 +1256 1713 32 +1260 1466 102 +1277 1665 32 +1283 2328 28 +1287 1542 52 +1288 1542 53 +1289 1542 42 +1290 1542 33 +1291 1542 52 +1292 1542 35 +1293 1542 30 +1294 1542 28 +1295 1542 44 +1296 1542 38 +1297 1542 35 +1298 1542 28 +1302 1542 32 +1304 1542 28 +1306 1542 32 +1402 2183 32 +1412 1542 29 +1413 1542 30 +1414 1542 29 +1415 1542 39 +1416 1542 34 +1417 1542 30 +1421 1542 36 +1422 1542 31 +1542 1786 48 +1542 1787 42 +1542 1788 45 +1542 1789 31 +1542 1790 39 +1542 1791 38 +1542 1792 38 +1542 1793 30 +1542 1794 41 +1542 1795 34 +1542 1800 30 +1542 1803 29 +1542 1806 28 +1542 1841 33 +1542 1913 32 +1542 1914 34 +1542 1916 31 +1542 1917 28 +1542 1919 37 +1542 1921 37 +1542 1922 33 +1542 1923 54 +1542 1924 30 +1542 2033 51 +1542 2034 33 +1542 2035 32 +1542 2036 43 +1542 2037 36 +1542 2038 47 +1542 2039 37 +1542 2040 41 +1542 2041 36 +1542 2042 38 +1542 2043 29 +1542 2044 34 +1542 2045 28 +1542 2046 38 +1542 2047 30 +1542 2048 32 +1542 2049 31 +1542 2050 28 +1542 2052 28 +1542 2053 28 +1542 2129 254 +1542 2137 38 +1542 2138 46 +1542 2139 32 +1542 2140 35 +1542 2141 28 +1542 2142 33 +1542 2146 29 +1542 2236 73 +1542 2237 53 +1542 2238 32 +1542 2242 39 +1542 2243 36 +1542 2244 31 +1542 2245 33 +1542 2248 28 +1542 2332 52 +1647 2272 68 +1660 2481 28 +1841 2129 254 diff --git a/output/analyze/analyze_sum_matrix/1M/threshold-0.1/output_sum_matrix_1M_1CDX2_v1.txt b/output/analyze/analyze_sum_matrix/1M/threshold-0.1/output_sum_matrix_1M_1CDX2_v1.txt new file mode 100644 index 0000000..71eec46 --- /dev/null +++ b/output/analyze/analyze_sum_matrix/1M/threshold-0.1/output_sum_matrix_1M_1CDX2_v1.txt @@ -0,0 +1,190 @@ +bin1 bin2 count +3 1542 47 +4 1542 61 +5 1542 43 +6 1542 41 +7 1542 40 +9 1542 51 +10 1542 41 +11 1542 33 +13 1542 39 +34 1542 37 +37 1542 31 +55 2462 36 +90 201 48 +90 1542 67 +111 964 32 +136 1199 42 +201 1542 58 +202 1542 62 +203 1542 47 +207 1542 36 +208 1542 38 +209 1542 41 +210 1542 33 +211 1542 37 +215 1542 37 +216 1542 37 +217 1542 34 +218 1542 33 +222 1542 36 +254 1359 45 +267 833 46 +271 645 37 +319 1163 33 +352 1280 39 +355 1496 32 +379 2272 38 +383 1542 42 +385 1542 33 +387 1542 36 +388 1542 44 +389 1542 38 +479 2063 35 +518 570 41 +543 852 33 +543 1155 31 +543 1542 44 +544 1542 33 +545 1542 34 +546 1542 33 +549 1542 31 +550 1542 35 +551 1542 31 +574 1542 31 +677 738 37 +695 2364 65 +699 1227 33 +699 1542 37 +704 1542 31 +716 1542 31 +729 1542 102 +771 1525 55 +786 1366 34 +826 2257 48 +843 2272 80 +852 1542 68 +852 2272 152 +853 1542 47 +854 1542 36 +856 1542 36 +857 1542 32 +858 1542 31 +860 1542 32 +861 1542 37 +864 1542 36 +866 1542 33 +870 1542 33 +871 1542 33 +877 1542 31 +878 1542 31 +879 1542 32 +886 1688 51 +916 2419 37 +919 2447 36 +925 2182 32 +940 1103 32 +1002 1542 100 +1003 1542 60 +1004 1542 32 +1005 1542 45 +1006 2648 34 +1018 1542 32 +1095 2327 35 +1134 2465 37 +1155 1542 113 +1155 2272 116 +1156 1542 57 +1157 1542 31 +1158 1542 44 +1159 1542 44 +1160 1542 42 +1161 1542 44 +1162 1542 40 +1163 1542 43 +1164 1542 33 +1165 1542 32 +1166 1542 31 +1170 1542 36 +1175 1542 42 +1178 1542 31 +1239 1339 36 +1256 1713 35 +1260 1466 151 +1277 1665 49 +1283 2328 43 +1287 1542 60 +1288 1542 49 +1289 1542 42 +1290 1542 40 +1291 1542 43 +1292 1542 31 +1293 1542 34 +1294 1542 32 +1296 1542 32 +1297 1542 31 +1299 1542 34 +1301 1542 31 +1304 1542 40 +1306 1542 32 +1415 1542 54 +1416 1542 39 +1417 1542 35 +1430 1786 34 +1542 1786 50 +1542 1787 43 +1542 1789 35 +1542 1790 36 +1542 1791 31 +1542 1792 38 +1542 1793 35 +1542 1795 37 +1542 1797 31 +1542 1841 45 +1542 1913 33 +1542 1916 38 +1542 1917 35 +1542 1918 35 +1542 1919 32 +1542 1920 33 +1542 1921 32 +1542 1922 38 +1542 1923 57 +1542 1924 55 +1542 1925 33 +1542 1926 33 +1542 2033 48 +1542 2034 36 +1542 2035 32 +1542 2036 39 +1542 2037 38 +1542 2038 39 +1542 2041 34 +1542 2042 40 +1542 2044 33 +1542 2048 35 +1542 2049 39 +1542 2057 39 +1542 2129 274 +1542 2137 40 +1542 2138 45 +1542 2236 83 +1542 2237 47 +1542 2238 44 +1542 2242 32 +1542 2243 42 +1542 2245 40 +1542 2246 34 +1542 2248 32 +1542 2250 33 +1542 2332 56 +1542 2489 33 +1543 2166 34 +1571 2057 31 +1571 2444 35 +1616 2055 55 +1647 2272 116 +1653 2455 42 +1660 2481 34 +1771 1888 42 +1841 2129 287 diff --git a/output/analyze/analyze_sum_matrix/1M/threshold-0.1/output_sum_matrix_1M_1CDX3_v1.txt b/output/analyze/analyze_sum_matrix/1M/threshold-0.1/output_sum_matrix_1M_1CDX3_v1.txt new file mode 100644 index 0000000..e8b6149 --- /dev/null +++ b/output/analyze/analyze_sum_matrix/1M/threshold-0.1/output_sum_matrix_1M_1CDX3_v1.txt @@ -0,0 +1,159 @@ +bin1 bin2 count +3 1542 34 +4 1542 50 +4 2502 37 +5 1542 44 +6 1542 37 +7 1542 30 +10 1542 27 +12 1542 29 +16 1542 28 +17 1542 27 +55 2462 41 +90 201 28 +90 1542 61 +111 964 35 +136 1199 44 +201 1542 60 +202 1542 40 +203 1248 30 +203 1542 40 +204 1542 39 +207 1542 30 +208 1542 31 +209 1542 35 +211 1542 29 +212 1542 27 +216 1542 34 +218 1542 28 +254 1359 45 +267 833 59 +271 645 35 +326 635 31 +329 2001 34 +351 1955 30 +352 1280 29 +383 1542 28 +388 1542 31 +389 1542 29 +399 694 27 +415 1916 35 +464 1013 36 +479 2063 41 +518 570 30 +543 852 33 +543 1155 60 +543 1542 36 +544 1542 31 +546 1542 30 +547 1542 28 +548 1542 35 +551 1542 33 +611 2429 28 +677 738 33 +684 2156 37 +694 2006 29 +695 2364 58 +699 1542 34 +729 1542 71 +771 1525 65 +826 1938 33 +826 2257 42 +843 2272 91 +852 1542 48 +852 2137 30 +852 2272 127 +853 1542 34 +854 1542 33 +855 1542 31 +856 1542 28 +857 1542 31 +886 1688 36 +902 2478 27 +905 1831 28 +915 1585 30 +916 2419 50 +919 2447 69 +940 1103 28 +957 2187 33 +996 1496 40 +1002 1542 77 +1003 1542 43 +1005 1542 45 +1012 1542 35 +1023 1196 30 +1095 2327 28 +1120 1462 34 +1134 2465 37 +1155 1542 86 +1155 2272 110 +1156 1542 40 +1157 1923 33 +1160 1542 35 +1161 1542 49 +1162 1542 37 +1163 1542 35 +1165 1542 27 +1166 1542 35 +1170 1542 28 +1175 1542 31 +1239 1339 43 +1260 1466 114 +1283 2328 59 +1287 1542 47 +1288 1542 33 +1289 1542 30 +1291 1542 32 +1292 1542 35 +1293 1542 28 +1295 1542 29 +1296 1542 30 +1297 1542 35 +1299 1542 28 +1301 1542 27 +1304 1542 30 +1402 2183 42 +1415 1542 47 +1416 1542 31 +1421 1542 28 +1422 1542 29 +1475 1788 55 +1542 1786 54 +1542 1787 33 +1542 1791 30 +1542 1792 32 +1542 1794 29 +1542 1795 29 +1542 1806 27 +1542 1841 36 +1542 1916 29 +1542 1918 34 +1542 1923 74 +1542 2033 35 +1542 2034 31 +1542 2036 31 +1542 2038 39 +1542 2041 34 +1542 2042 41 +1542 2129 235 +1542 2137 32 +1542 2138 27 +1542 2236 61 +1542 2237 33 +1542 2247 27 +1542 2332 42 +1543 2166 27 +1571 2057 27 +1571 2444 37 +1571 2445 47 +1616 2055 43 +1622 1838 27 +1647 2272 118 +1653 2455 27 +1660 2480 35 +1660 2481 48 +1669 2283 31 +1771 1888 35 +1841 2129 242 +2137 2332 27 +2267 2454 30 diff --git a/output/analyze/analyze_sum_matrix/1M/threshold-0.1/output_sum_matrix_1M_1CDX4_v1.txt b/output/analyze/analyze_sum_matrix/1M/threshold-0.1/output_sum_matrix_1M_1CDX4_v1.txt new file mode 100644 index 0000000..21cc43b --- /dev/null +++ b/output/analyze/analyze_sum_matrix/1M/threshold-0.1/output_sum_matrix_1M_1CDX4_v1.txt @@ -0,0 +1,214 @@ +bin1 bin2 count +3 1542 68 +4 1542 66 +5 1542 54 +6 1542 61 +7 1542 48 +8 1542 38 +9 1542 52 +10 1542 56 +11 1542 40 +12 1542 38 +13 1542 46 +14 1542 35 +17 1476 36 +35 1542 38 +55 2462 38 +66 1352 50 +90 201 44 +90 202 35 +90 543 37 +90 1542 82 +111 964 62 +136 1199 62 +201 1542 71 +202 1542 64 +203 1248 33 +203 1542 75 +204 1542 50 +205 1542 37 +207 1542 33 +208 1542 34 +209 1542 40 +210 1542 37 +215 1542 33 +216 1542 44 +244 2010 35 +254 1359 63 +267 833 59 +271 645 36 +297 1482 41 +297 2575 36 +329 2001 37 +337 2007 42 +351 1955 53 +352 1280 45 +355 1496 33 +360 602 35 +383 1542 46 +388 1542 49 +389 1542 47 +468 1237 39 +479 2063 49 +516 1940 33 +518 570 44 +543 852 64 +543 1155 77 +543 1542 61 +543 2033 37 +544 1542 45 +545 1542 40 +546 1542 37 +547 1542 42 +548 1542 54 +549 1542 40 +550 1542 35 +551 1542 51 +555 1542 39 +677 738 62 +695 2364 59 +699 1227 33 +699 1542 46 +701 1542 34 +728 1542 39 +729 1542 126 +771 1525 72 +786 1366 55 +826 1938 36 +826 2257 60 +843 2272 88 +852 1542 75 +852 2137 36 +852 2272 156 +853 1542 57 +854 1542 57 +855 1542 46 +857 1542 42 +858 1542 41 +859 1542 35 +866 1542 41 +867 1542 33 +875 1182 44 +877 1542 34 +886 1688 65 +915 1585 55 +916 2419 57 +919 2447 62 +925 2182 55 +955 2521 36 +996 1496 48 +1002 1542 110 +1002 2129 35 +1003 1542 81 +1004 1542 34 +1005 1542 67 +1016 1542 36 +1023 1196 40 +1092 1977 34 +1095 2327 46 +1102 2203 33 +1108 1331 33 +1120 1462 33 +1134 2465 36 +1155 1542 148 +1155 2272 127 +1156 1542 90 +1157 1542 63 +1158 1542 54 +1159 1542 45 +1160 1542 80 +1161 1542 66 +1162 1542 60 +1163 1542 47 +1164 1542 41 +1165 1542 45 +1166 1542 38 +1167 1542 33 +1170 1542 33 +1175 1542 35 +1177 1542 34 +1239 1339 35 +1256 1713 42 +1260 1466 173 +1277 1665 35 +1283 2328 62 +1287 1542 78 +1288 1542 46 +1289 1542 45 +1290 1542 35 +1291 1542 47 +1292 1542 33 +1293 1542 37 +1294 1542 40 +1296 1542 36 +1297 1542 39 +1298 1542 48 +1299 1542 38 +1304 1542 39 +1305 1542 42 +1366 1822 39 +1402 2183 59 +1414 1542 43 +1415 1542 66 +1416 1542 39 +1418 1542 33 +1421 1542 38 +1427 1542 33 +1439 2407 39 +1475 1788 63 +1542 1786 71 +1542 1787 55 +1542 1788 46 +1542 1789 35 +1542 1790 38 +1542 1791 41 +1542 1792 44 +1542 1794 33 +1542 1797 36 +1542 1841 60 +1542 1916 41 +1542 1917 39 +1542 1918 41 +1542 1919 37 +1542 1920 40 +1542 1921 47 +1542 1922 52 +1542 1923 93 +1542 1924 53 +1542 2033 49 +1542 2034 34 +1542 2035 35 +1542 2036 34 +1542 2037 40 +1542 2038 42 +1542 2041 33 +1542 2042 38 +1542 2044 39 +1542 2045 37 +1542 2129 303 +1542 2137 44 +1542 2138 38 +1542 2139 34 +1542 2140 40 +1542 2141 36 +1542 2144 37 +1542 2236 68 +1542 2237 51 +1542 2238 46 +1542 2248 34 +1542 2332 89 +1542 2493 33 +1571 2444 45 +1571 2445 37 +1616 2055 79 +1622 1838 34 +1647 2272 144 +1653 2455 45 +1660 2480 50 +1660 2481 43 +1669 2283 37 +1771 1888 56 +1841 2129 307 +2137 2332 57 +2137 2423 49 +2267 2454 45 diff --git a/output/analyze/analyze_sum_matrix/500k/threshold-0.1/output_sum_matrix_500k_1CDX1_v1.txt b/output/analyze/analyze_sum_matrix/500k/threshold-0.1/output_sum_matrix_500k_1CDX1_v1.txt new file mode 100644 index 0000000..1a2f5b5 --- /dev/null +++ b/output/analyze/analyze_sum_matrix/500k/threshold-0.1/output_sum_matrix_500k_1CDX1_v1.txt @@ -0,0 +1,72 @@ +bin1 bin2 count +6 3082 33 +7 3082 30 +8 3082 32 +9 3082 37 +10 3082 33 +12 3082 33 +13 3082 35 +19 3082 34 +20 3082 28 +21 3082 28 +24 3082 28 +180 3082 65 +222 1928 36 +401 3082 45 +408 3082 29 +417 3082 32 +508 2717 30 +1085 2309 28 +1389 4721 31 +1457 3082 75 +1541 3048 36 +1572 2731 28 +1685 4538 47 +1703 3082 37 +1703 4538 105 +1704 3082 48 +1705 3082 30 +1706 3082 32 +1713 3082 31 +1771 3375 41 +1829 3168 28 +1837 4887 32 +2003 3082 66 +2005 3082 29 +2006 3082 29 +2009 3082 28 +2309 3082 69 +2310 3082 45 +2310 4538 87 +2311 3082 47 +2312 3082 35 +2313 3082 30 +2314 3082 33 +2315 3082 35 +2320 3082 33 +2321 3082 35 +2322 3082 32 +2323 3082 28 +2330 3082 28 +2340 3082 28 +2511 3424 28 +2520 2931 85 +2554 3329 29 +2573 3082 33 +2575 3082 30 +2581 3082 30 +2582 3082 29 +2592 3082 28 +2803 4361 30 +3082 3570 31 +3082 3679 30 +3082 3842 29 +3082 4061 34 +3082 4254 254 +3082 4465 51 +3082 4466 32 +3082 4467 33 +3082 4656 46 +3293 4538 66 +3319 4954 28 +3679 4254 254 diff --git a/output/analyze/analyze_sum_matrix/500k/threshold-0.1/output_sum_matrix_500k_1CDX2_v1.txt b/output/analyze/analyze_sum_matrix/500k/threshold-0.1/output_sum_matrix_500k_1CDX2_v1.txt new file mode 100644 index 0000000..ec09359 --- /dev/null +++ b/output/analyze/analyze_sum_matrix/500k/threshold-0.1/output_sum_matrix_500k_1CDX2_v1.txt @@ -0,0 +1,61 @@ +bin1 bin2 count +8 3082 33 +9 3082 35 +14 3082 33 +19 3082 36 +111 4916 32 +180 401 38 +180 3082 65 +273 2397 41 +401 3082 44 +403 3082 36 +404 3082 36 +406 3082 31 +508 2717 38 +534 1666 42 +542 1290 36 +704 2560 35 +758 4538 37 +957 4122 32 +1035 1140 36 +1389 4721 60 +1398 2454 31 +1457 3082 95 +1541 3048 47 +1652 4507 31 +1685 4538 80 +1703 3082 40 +1703 4538 152 +1704 3082 33 +1771 3375 46 +1837 4887 34 +1879 2206 32 +2003 3082 86 +2005 3082 35 +2011 5287 34 +2189 4647 31 +2309 3082 86 +2310 3082 45 +2310 4538 116 +2311 3082 38 +2520 2930 48 +2520 2931 126 +2554 3329 48 +2573 3082 38 +2575 3082 32 +2828 3082 35 +2859 3569 33 +3082 3570 35 +3082 3679 40 +3082 3842 34 +3082 4254 274 +3082 4465 56 +3082 4466 37 +3082 4468 31 +3082 4656 50 +3231 4106 43 +3293 4538 116 +3305 4902 38 +3319 4954 33 +3541 3774 32 +3679 4254 287 diff --git a/output/analyze/analyze_sum_matrix/500k/threshold-0.1/output_sum_matrix_500k_1CDX3_v1.txt b/output/analyze/analyze_sum_matrix/500k/threshold-0.1/output_sum_matrix_500k_1CDX3_v1.txt new file mode 100644 index 0000000..709137b --- /dev/null +++ b/output/analyze/analyze_sum_matrix/500k/threshold-0.1/output_sum_matrix_500k_1CDX3_v1.txt @@ -0,0 +1,69 @@ +bin1 bin2 count +9 3082 38 +9 4996 30 +111 4916 39 +180 3082 57 +222 1928 31 +273 2397 38 +401 3082 40 +402 3082 28 +406 2495 27 +508 2717 38 +534 1666 45 +542 1290 28 +651 1269 29 +657 3999 34 +797 1388 27 +830 3828 29 +957 4122 38 +1035 1140 27 +1085 1703 27 +1085 2309 50 +1086 3082 27 +1367 4306 28 +1388 4008 27 +1389 4721 53 +1457 3082 65 +1541 3048 57 +1652 3873 30 +1652 4507 34 +1685 4538 91 +1703 4538 127 +1771 3375 30 +1810 3660 27 +1831 4831 42 +1837 4887 58 +1914 4369 30 +1991 2990 36 +2003 3082 61 +2009 3082 29 +2239 2922 33 +2268 4922 33 +2309 3082 62 +2310 3082 40 +2310 4538 110 +2321 3082 36 +2478 2677 35 +2520 2930 31 +2520 2931 100 +2566 4649 42 +2573 3082 28 +2803 4361 40 +2829 3082 30 +2948 3574 52 +3082 3570 42 +3082 3679 35 +3082 3842 45 +3082 3843 41 +3082 4080 30 +3082 4254 234 +3082 4465 42 +3082 4656 41 +3141 4881 30 +3141 4882 40 +3231 4106 36 +3293 4538 118 +3319 4953 28 +3319 4954 46 +3541 3774 27 +3679 4254 241 diff --git a/output/analyze/analyze_sum_matrix/500k/threshold-0.1/output_sum_matrix_500k_1CDX4_v1.txt b/output/analyze/analyze_sum_matrix/500k/threshold-0.1/output_sum_matrix_500k_1CDX4_v1.txt new file mode 100644 index 0000000..a3812e7 --- /dev/null +++ b/output/analyze/analyze_sum_matrix/500k/threshold-0.1/output_sum_matrix_500k_1CDX4_v1.txt @@ -0,0 +1,110 @@ +bin1 bin2 count +6 3082 42 +7 3082 38 +8 3082 39 +9 3082 44 +11 3082 35 +12 3082 35 +14 3082 33 +19 3082 36 +20 3082 36 +26 3082 35 +111 4916 37 +132 2703 40 +180 3082 79 +222 1928 56 +273 2397 57 +401 3082 52 +402 3082 35 +403 3082 39 +405 3082 42 +406 3082 46 +408 3082 34 +487 4017 34 +508 2717 51 +534 1666 49 +593 5141 34 +594 2962 33 +674 4011 36 +702 3907 51 +704 2560 40 +720 1203 33 +935 2474 36 +957 4122 38 +1035 1140 39 +1085 1703 54 +1085 2309 60 +1085 4061 33 +1086 3082 41 +1096 3082 39 +1101 3082 33 +1353 1476 33 +1389 4721 55 +1457 3082 113 +1541 3048 63 +1572 2731 45 +1652 3873 33 +1652 4507 42 +1685 4538 88 +1703 3082 41 +1703 4538 156 +1704 3082 42 +1771 3375 56 +1829 3168 46 +1831 4831 40 +1837 4887 48 +1849 4358 44 +1991 2990 39 +2003 3082 87 +2004 3082 42 +2005 3082 53 +2006 3082 44 +2009 3082 33 +2010 3082 44 +2189 4647 39 +2309 3082 107 +2310 3082 85 +2310 4538 127 +2311 3082 64 +2312 3082 42 +2313 3082 41 +2314 3082 33 +2316 3082 33 +2317 3082 34 +2320 3082 57 +2321 3082 44 +2322 3082 37 +2323 3082 39 +2511 3424 36 +2520 2930 52 +2520 2931 149 +2565 4649 34 +2573 3082 46 +2574 3082 40 +2587 3082 33 +2732 3641 38 +2803 4361 53 +2828 3082 38 +2829 3082 35 +2948 3574 54 +3082 3570 55 +3082 3571 38 +3082 3574 36 +3082 3679 56 +3082 3842 59 +3082 3843 46 +3082 3844 39 +3082 4254 303 +3082 4465 50 +3082 4656 85 +3141 4881 35 +3231 4106 62 +3293 4538 144 +3305 4902 39 +3319 4953 36 +3319 4954 43 +3541 3774 42 +3679 4254 307 +4268 4656 41 +4268 4838 39 +4527 4900 39 diff --git a/output/cytoscape-fig/SumMatrix_500k_1M.cys b/output/cytoscape-fig/SumMatrix_500k_1M.cys new file mode 100644 index 0000000..d0acff7 Binary files /dev/null and b/output/cytoscape-fig/SumMatrix_500k_1M.cys differ diff --git a/src/analyze/AnalyzeSumMatrix.py b/src/analyze/AnalyzeSumMatrix.py index 580ee1c..c7f0cc8 100644 --- a/src/analyze/AnalyzeSumMatrix.py +++ b/src/analyze/AnalyzeSumMatrix.py @@ -61,8 +61,11 @@ def main(): matrix_file = "analyze-{}/SumMatrix/sum_matrix_{}_{}.txt".format(bin_size, bin_size, cell_type) chrom_bin_range = "metadata/chrom_bins_range_{}.txt".format(bin_size) - output_file = "output/analyze/analyze_sum_matrix/{}/output_sum_matrix_{}_{}.txt".format(bin_size, bin_size, - cell_type) + # output_file = "output/analyze/analyze_sum_matrix/{}/output_sum_matrix_{}_{}_v1.txt".format(bin_size, bin_size, + # cell_type) + output_file = "output/analyze/analyze_sum_matrix/{}/threshold-0.1/output_sum_matrix_{}_{}_v1.txt".format(bin_size, + bin_size, + cell_type) bin_range = {} with open(chrom_bin_range) as f: @@ -121,22 +124,31 @@ def main(): # unique, counts = np.unique(data, return_counts=True) M = ((count_N - count_X) ** 2 - sum_reduction) / 2 + # print(M) + # sys.exit(0) + count1 = 0 p_max = calculate_pmax(M, sum_max) # my_list = [] data = np.triu(data, k=-1) out = open(output_file, "w") - out.write("{} {} {} {}\n".format("bin1", "bin2", "count", "p_value")) + out.write("{} {} {}\n".format("bin1", "bin2", "count", "p_value")) + # out.write("{} {} {} {}\n".format("bin1", "bin2", "count", "p_value")) for i in range(0, rows): for j in range(0, cols): if data[i][j] != 0: value = int(data[i][j]) p_value = calculate_f(number_of_cells, value, p_max) - if p_value <= threshold(M): + # if p_value <= threshold(M): + # count1 += 1 + # out.write("{} {} {} {}\n".format(i, j, value, p_value)) + # my_list.append(value) + + if value >= number_of_cells / 10: count1 += 1 - out.write("{} {} {} {}\n".format(i, j, value, p_value)) - # my_list.append(value) + out.write("{} {} {}\n".format(i, j, value, p_value)) + out.close() # my_arr = np.array(my_list)