In [1]:
from pytadbit.mapping.analyze import eig_correlate_matrices, correlate_matrices
from pytadbit import load_hic_data_from_reads
from cPickle import load
from matplotlib import pyplot as plt
In [2]:
reso = 1000000
base_path = 'results/{0}/03_filtering/valid_reads12_{0}.tsv'
bias_ice_path = 'results/{1}/04_normalizing/biases_{0}_{1}.pick'
bias_dry_path = 'results/{1}/04_normalizing/biases_dryhic_{0}_{1}.tsv'
bads_path = 'results/{1}/04_normalizing/bad_columns_{0}_{1}.pick'
Write a little function to load HiCData obeject
In [3]:
def my_load_hic_data(renz, reso, which='ice'):
hic_data = load_hic_data_from_reads(base_path.format(renz), resolution=reso)
if which=='ice':
hic_data.bias = load(open(bias_ice_path.format(reso, renz)))
else:
hic_data.bias = dict([(int(l.split()[0]), float(l.split()[1]))
for l in open(bias_dry_path.format(reso, renz))])
hic_data.bads = load(open(bads_path.format(reso, renz)))
return hic_data
In [4]:
renz1 = 'HindIII'
renz2 = 'MboI'
In [5]:
hic_data1_ice = my_load_hic_data(renz1, reso, which='ice')
hic_data2_ice = my_load_hic_data(renz2, reso, which='ice')
In [6]:
hic_data1_dry = my_load_hic_data(renz1, reso, which='dry')
hic_data2_dry = my_load_hic_data(renz2, reso, which='dry')
In [8]:
%matplotlib inline
In [9]:
## this part is to "tune" the plot ##
plt.figure(figsize=(9, 6))
axe = plt.subplot()
axe.grid()
axe.set_xticks(range(0, 55, 5))
axe.set_xticklabels(['%d Mb' % int(i * 0.2) if i else '' for i in range(0, 55, 5)], rotation=-45)
#####################################
_ = correlate_matrices(hic_data1_ice, hic_data2_ice, max_dist=50, show=False, axe=axe, normalized=True)
In [10]:
## this part is to "tune" the plot ##
plt.figure(figsize=(9, 6))
axe = plt.subplot()
axe.grid()
axe.set_xticks(range(0, 55, 5))
axe.set_xticklabels(['%d Mb' % int(i * 0.2) if i else '' for i in range(0, 55, 5)], rotation=-45)
#####################################
_ = correlate_matrices(hic_data1_dry, hic_data2_dry, max_dist=50, show=False, axe=axe, normalized=True)
In [12]:
reso = 1000000
hic_data1_ice = my_load_hic_data(renz1, reso, which='ice')
hic_data2_ice = my_load_hic_data(renz2, reso, which='ice')
hic_data1_dry = my_load_hic_data(renz1, reso, which='dry')
hic_data2_dry = my_load_hic_data(renz2, reso, which='dry')
In [13]:
## this part is to "tune" the plot ##
plt.figure(figsize=(9, 6))
axe = plt.subplot()
axe.grid()
axe.set_xticks(range(0, 55, 5))
axe.set_xticklabels(['%d Mb' % int(i * 0.2) if i else '' for i in range(0, 55, 5)], rotation=-45)
#####################################
_ = correlate_matrices(hic_data1_ice, hic_data2_ice, max_dist=50, show=False, axe=axe, normalized=True)
In [14]:
## this part is to "tune" the plot ##
plt.figure(figsize=(9, 6))
axe = plt.subplot()
axe.grid()
axe.set_xticks(range(0, 55, 5))
axe.set_xticklabels(['%d Mb' % int(i * 0.2) if i else '' for i in range(0, 55, 5)], rotation=-45)
#####################################
_ = correlate_matrices(hic_data1_dry, hic_data2_dry, max_dist=50, show=False, axe=axe, normalized=True)
In [15]:
corrs = eig_correlate_matrices(hic_data1_ice, hic_data2_ice, show=True, aspect='auto', normalized=True)
for cor in corrs:
print ' '.join(['%5.3f' % (c) for c in cor]) + '\n'
0.613 0.282 0.326 0.045 0.085 0.064
0.740 0.391 0.222 0.077 0.002 0.109
0.109 0.821 0.323 0.191 0.069 0.024
0.068 0.187 0.308 0.797 0.244 0.084
0.037 0.087 0.257 0.206 0.095 0.885
0.092 0.086 0.472 0.374 0.179 0.285
In [17]:
corrs = eig_correlate_matrices(hic_data1_dry, hic_data2_dry, show=True, aspect='auto', normalized=True)
for cor in corrs:
print ' '.join(['%5.3f' % (c) for c in cor]) + '\n'
0.755 0.363 0.020 0.130 0.024 0.238
0.529 0.099 0.179 0.231 0.272 0.373
0.226 0.364 0.022 0.068 0.765 0.172
0.040 0.288 0.689 0.076 0.201 0.510
0.151 0.503 0.616 0.043 0.377 0.367
0.121 0.563 0.153 0.069 0.026 0.416
Once agreed that experiments are similar, they can be merged.
Here is a simple way to merge valid pairs. Arguably we may want to merge unfiltered data but the difference would be minimal specially with non-replicates.
In [18]:
from pytadbit.mapping import merge_2d_beds
In [24]:
! mkdir -p results/both
! mkdir -p results/both/
! mkdir -p results/both03_filtering/
! mkdir -p results/both/03_filtering/
In [25]:
renz1 = 'HindIII'
renz2 = 'MboI'
hic_data1 = 'results/{0}/03_filtering/valid_reads12_{0}.tsv'.format(renz1)
hic_data2 = 'results/{0}/03_filtering/valid_reads12_{0}.tsv'.format(renz2)
hic_data = 'results/both/03_filtering/valid_reads12.tsv'
merge_2d_beds(hic_data1, hic_data2, hic_data)
Out[25]:
9751233
In [26]:
from pytadbit.mapping.analyze import hic_map
from cPickle import dump
In [27]:
! mkdir -p results/both/04_normalizing
All in one loop to:
all at diferent resolutions, and for both time points
In [28]:
for reso in [1000000, 300000, 100000]:
print ' *', reso
# load hic_data
hic_data = load_hic_data_from_reads(
'results/both/03_filtering/valid_reads12.tsv', reso)
# filter columns
hic_data.filter_columns(draw_hist=False, min_count=10, by_mean=True)
# normalize
hic_data.normalize_hic(iterations=0)
# save biases to reconstruct normalization
out = open('results/both/04_normalizing/biases_{0}.pick'.format(reso), 'w')
dump(hic_data.bias, out)
out.close()
# save filtered out columns
out = open('results/both/04_normalizing/bad_columns_{0}.pick'.format(reso), 'w')
dump(hic_data.bads, out)
out.close()
# save data as raw matrix per chromsome
hic_map(hic_data, by_chrom='intra', normalized=False,
savedata='results/both/04_normalizing/{0}_raw'.format(reso))
# save data as normalized matrix per chromosome
hic_map(hic_data, by_chrom='intra', normalized=True,
savedata='results/both/04_normalizing/{0}_norm'.format(reso))
# if the resolution is low save the full genomic matrix
if reso > 500000:
hic_map(hic_data, by_chrom=False, normalized=False,
savefig ='results/both/04_normalizing/{0}_raw.png'.format(reso),
savedata='results/both/04_normalizing/{0}_raw.mat'.format(reso))
hic_map(hic_data, by_chrom=False, normalized=True,
savefig ='results/both/04_normalizing/{0}_norm.png'.format(reso) ,
savedata='results/both/04_normalizing/{0}_norm.mat'.format(reso))
-
WARNING: removing columns having less than 20 counts:
124 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 742 931
932 1295 1563 1571 1586 1587 1588 1589 1590 1591 1592 1593 1594 1595 1596 1597 1598 1599 1600 1601
1721 1867 1868 2085 2086 2087 2088 2089 2090 2091 2092 2093 2094 2095 2096 2097 2098 2099 2100 2101
2102 2200 2201 2202 2203 2204 2205 2206 2207 2208 2209 2210 2211 2212 2213 2214 2215 2216 2217 2307
2308 2309 2310 2311 2312 2313 2314 2315 2316 2317 2318 2319 2320 2321 2322 2323 2324 2447 2448 2449
2450 2451 2452 2453 2454 2455 2601 2602 2603 2604 2691 2692 2752 2790 2791 2792 2793 2794 2801 2837
2838 2839 2840 2841 2842 2843 2844 2845 2846 2850 2851 2888 2889 2947 2948 2949 3044 3045 3046 3047
3052 3053 3054 3063 3067 3068 3069 3070 3072 3073 3074 3075 3076 3077 3078 3079 3080 3081 3082 3083
3084 3085 3086 3087 3088 3089 3090 3091 3092 3093 3094 3095 3096 3097 3098 3099 3100 3102
WARNING: removing columns having less than 1082.318 counts:
1 121 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140
141 142 143 144 145 147 337 340 341 342 492 584 585 691 742 882 930 931 932 953
1064 1124 1294 1295 1296 1403 1440 1441 1541 1563 1564 1565 1566 1567 1571 1572 1581 1582 1585 1586
1587 1588 1589 1590 1591 1592 1593 1594 1595 1596 1597 1598 1599 1600 1601 1602 1603 1604 1605 1606
1607 1608 1609 1680 1721 1866 1867 1868 1950 1986 1987 2085 2086 2087 2088 2089 2090 2091 2092 2093
2094 2095 2096 2097 2098 2099 2100 2101 2102 2199 2200 2201 2202 2203 2204 2205 2206 2207 2208 2209
2210 2211 2212 2213 2214 2215 2216 2217 2218 2307 2308 2309 2310 2311 2312 2313 2314 2315 2316 2317
2318 2319 2320 2321 2322 2323 2324 2325 2326 2327 2329 2446 2447 2448 2449 2450 2451 2452 2453 2454
2455 2500 2524 2525 2584 2601 2602 2603 2604 2605 2665 2666 2667 2691 2692 2752 2789 2790 2791 2792
2793 2794 2795 2796 2797 2801 2802 2837 2838 2839 2840 2841 2842 2843 2844 2845 2846 2847 2850 2851
2852 2888 2889 2890 2947 2948 2949 3044 3045 3046 3047 3048 3049 3050 3051 3052 3053 3054 3055 3056
3057 3058 3059 3060 3061 3062 3063 3064 3065 3066 3067 3068 3069 3070 3071 3072 3073 3074 3075 3076
3077 3078 3079 3080 3081 3082 3083 3084 3085 3086 3087 3088 3089 3090 3091 3092 3093 3094 3095 3096
3097 3098 3099 3100 3101 3102
T60
* 1000000
Found 266 of 3102 columns with poor signal
iterative correction
- copying matrix
- computing baises
rescaling to factor 1
- getting the sum of the matrix
=> 3198.486
- rescaling biases
* 300000
WARNING: removing columns having less than 20 counts:
410 411 412 413 414 415 416 419 420 421 422 423 424 425 426 427 428 429 430 431
432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451
452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471
472 473 474 475 476 477 481 483 1133 1134 1945 1947 1948 2466 2467 2468 2469 2470 2934 3093
3094 3095 3096 3097 3098 3099 3100 3101 3167 3540 3736 3737 3738 3739 3740 4306 4307 4308 4309 4310
4311 4312 4313 5197 5198 5199 5200 5201 5202 5203 5204 5205 5206 5207 5208 5209 5210 5211 5212 5213
5214 5223 5224 5225 5226 5227 5228 5229 5230 5231 5261 5272 5273 5274 5275 5276 5277 5278 5279 5280
5281 5282 5283 5284 5285 5286 5287 5288 5289 5290 5291 5292 5293 5294 5295 5296 5297 5298 5299 5300
5301 5302 5303 5304 5305 5306 5307 5308 5309 5310 5311 5312 5313 5314 5315 5316 5317 5318 5319 5320
5321 5322 5323 5324 5325 5326 5327 5329 5333 5334 5337 5588 5721 5722 5723 5724 5725 5726 6205 6206
6207 6209 6210 6211 6212 6213 6214 6215 6931 6932 6933 6934 6935 6936 6937 6938 6939 6940 6941 6942
6943 6944 6945 6946 6947 6948 6949 6950 6951 6952 6953 6954 6955 6956 6957 6958 6959 6960 6961 6962
6963 6964 6965 6966 6967 6968 6969 6970 6971 6972 6973 6974 6975 6976 6977 6978 6979 6980 6981 6982
6983 6984 6985 6986 6987 6988 6989 6990 7313 7314 7315 7316 7317 7318 7319 7320 7321 7322 7323 7324
7325 7326 7327 7328 7329 7330 7331 7332 7333 7334 7335 7336 7337 7338 7339 7340 7341 7342 7343 7344
7345 7346 7347 7348 7349 7350 7351 7352 7353 7354 7355 7356 7357 7358 7359 7360 7361 7362 7363 7364
7365 7366 7367 7368 7369 7370 7371 7372 7670 7671 7672 7673 7674 7675 7676 7677 7678 7679 7680 7681
7682 7683 7684 7685 7686 7687 7688 7689 7690 7691 7692 7693 7694 7695 7696 7697 7698 7699 7700 7701
7702 7703 7704 7705 7706 7707 7708 7709 7710 7711 7712 7713 7714 7715 7716 7717 7718 7719 7720 7721
7722 7723 7724 7725 7731 7732 7733 7734 8132 8133 8134 8135 8136 8137 8138 8139 8140 8141 8142 8143
8144 8145 8146 8147 8148 8149 8150 8151 8152 8153 8154 8155 8156 8157 8158 8159 8160 8161 8162 8163
8311 8643 8644 8645 8646 8647 8648 8649 8650 8651 8652 8653 8654 8655 8656 8657 8940 8941 8942 8943
8944 8945 8946 8947 9142 9143 9144 9145 9146 9147 9148 9269 9270 9271 9272 9273 9274 9275 9276 9277
9278 9279 9280 9281 9282 9283 9284 9287 9288 9289 9292 9294 9296 9306 9307 9308 9309 9310 9311 9425
9426 9427 9428 9429 9430 9431 9432 9433 9434 9435 9436 9437 9438 9439 9440 9441 9442 9443 9444 9445
9446 9447 9448 9449 9450 9451 9452 9453 9454 9455 9456 9457 9458 9459 9468 9469 9470 9471 9472 9473
9474 9595 9596 9597 9598 9599 9600 9601 9602 9603 9791 9792 9793 9794 9795 9796 9797 9798 9799 9800
9801 9802 10115 10116 10117 10118 10119 10120 10121 10122 10123 10124 10125 10128 10129 10130 10133 10136 10137 10138
10139 10140 10141 10142 10143 10144 10145 10146 10147 10148 10151 10152 10156 10157 10158 10159 10160 10162 10163 10164
10165 10166 10167 10168 10169 10170 10171 10172 10173 10174 10175 10176 10177 10178 10179 10180 10183 10184 10185 10187
10188 10189 10190 10191 10192 10193 10194 10195 10196 10197 10198 10199 10200 10201 10202 10203 10205 10206 10207 10208
10209 10210 10211 10212 10213 10214 10215 10216 10217 10218 10219 10220 10221 10222 10223 10224 10225 10226 10227 10228
10229 10230 10231 10232 10233 10234 10235 10236 10237 10238 10239 10240 10241 10242 10243 10244 10245 10246 10247 10248
10249 10250 10251 10252 10253 10254 10255 10256 10257 10258 10259 10260 10261 10262 10263 10264 10265 10266 10267 10268
10269 10270 10271 10272 10273 10274 10275 10276 10277 10278 10279 10280 10281 10282 10283 10284 10285 10286 10287 10288
10289 10290 10291 10292 10293 10294 10295 10296 10297 10298 10299 10300 10301 10302 10303 10304 10306
WARNING: removing columns having less than 286.661 counts:
1 2 5 401 403 407 408 409 410 411 412 413 414 415 416 418 419 420 421 422
423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442
443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462
463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 480 481 482 483
484 485 488 496 498 499 1122 1123 1132 1133 1134 1135 1138 1145 1155 1197 1638 1944 1945 1946
1947 1948 1949 1950 2466 2467 2468 2469 2470 2471 2934 3090 3092 3093 3094 3095 3096 3097 3098 3099
3100 3101 3167 3168 3169 3170 3171 3172 3539 3540 3736 3737 3738 3739 3740 3741 4110 4305 4306 4307
4308 4309 4310 4311 4312 4313 4315 4322 4354 4360 4361 4642 4668 4669 4790 4791 4792 4793 4794 4795
5196 5197 5198 5199 5200 5201 5202 5203 5204 5205 5206 5207 5208 5209 5210 5211 5212 5213 5214 5222
5223 5224 5225 5226 5227 5228 5229 5230 5231 5257 5258 5259 5260 5261 5262 5265 5266 5269 5270 5271
5272 5273 5274 5275 5276 5277 5278 5279 5280 5281 5282 5283 5284 5285 5286 5287 5288 5289 5290 5291
5292 5293 5294 5295 5296 5297 5298 5299 5300 5301 5302 5303 5304 5305 5306 5307 5308 5309 5310 5311
5312 5313 5314 5315 5316 5317 5318 5319 5320 5321 5322 5323 5324 5325 5326 5327 5328 5329 5330 5331
5332 5333 5334 5335 5337 5338 5341 5342 5343 5344 5345 5346 5347 5348 5349 5350 5351 5352 5353 5354
5588 5721 5722 5723 5724 5725 5726 5727 5748 6035 6205 6206 6207 6208 6209 6210 6211 6212 6213 6214
6215 6216 6602 6603 6604 6605 6606 6607 6608 6609 6930 6931 6932 6933 6934 6935 6936 6937 6938 6939
6940 6941 6942 6943 6944 6945 6946 6947 6948 6949 6950 6951 6952 6953 6954 6955 6956 6957 6958 6959
6960 6961 6962 6963 6964 6965 6966 6967 6968 6969 6970 6971 6972 6973 6974 6975 6976 6977 6978 6979
6980 6981 6982 6983 6984 6985 6986 6987 6988 6989 6990 7312 7313 7314 7315 7316 7317 7318 7319 7320
7321 7322 7323 7324 7325 7326 7327 7328 7329 7330 7331 7332 7333 7334 7335 7336 7337 7338 7339 7340
7341 7342 7343 7344 7345 7346 7347 7348 7349 7350 7351 7352 7353 7354 7355 7356 7357 7358 7359 7360
7361 7362 7363 7364 7365 7366 7367 7368 7369 7370 7371 7372 7373 7374 7376 7377 7669 7670 7671 7672
7673 7674 7675 7676 7677 7678 7679 7680 7681 7682 7683 7684 7685 7686 7687 7688 7689 7690 7691 7692
7693 7694 7695 7696 7697 7698 7699 7700 7701 7702 7703 7704 7705 7706 7707 7708 7709 7710 7711 7712
7713 7714 7715 7716 7717 7718 7719 7720 7721 7722 7723 7724 7725 7727 7728 7729 7730 7731 7732 7733
7734 7735 7739 7741 7742 7743 7747 7765 7771 7777 8064 8071 8108 8119 8120 8121 8126 8131 8132 8133
8134 8135 8136 8137 8138 8139 8140 8141 8142 8143 8144 8145 8146 8147 8148 8149 8150 8151 8152 8153
8154 8155 8156 8157 8158 8159 8160 8161 8162 8163 8164 8311 8373 8375 8388 8389 8390 8391 8392 8393
8394 8396 8397 8399 8585 8643 8644 8645 8646 8647 8648 8649 8650 8651 8652 8653 8654 8655 8656 8657
8658 8858 8859 8860 8861 8862 8869 8873 8904 8940 8941 8942 8943 8944 8945 8946 8947 9053 9142 9143
9144 9145 9146 9147 9148 9154 9268 9269 9270 9271 9272 9273 9274 9275 9276 9277 9278 9279 9280 9281
9282 9283 9284 9285 9287 9288 9289 9290 9291 9292 9294 9296 9297 9306 9307 9308 9309 9310 9311 9425
9426 9427 9428 9429 9430 9431 9432 9433 9434 9435 9436 9437 9438 9439 9440 9441 9442 9443 9444 9445
9446 9447 9448 9449 9450 9451 9452 9453 9454 9455 9456 9457 9458 9459 9468 9469 9470 9471 9472 9473
9474 9475 9476 9486 9487 9594 9595 9596 9597 9598 9599 9600 9601 9602 9603 9604 9790 9791 9792 9793
9794 9795 9796 9797 9798 9799 9800 9801 9802 10109 10114 10115 10116 10117 10118 10119 10120 10121 10122 10123
10124 10125 10126 10127 10128 10129 10130 10131 10132 10133 10134 10135 10136 10137 10138 10139 10140 10141 10142 10143
10144 10145 10146 10147 10148 10149 10150 10151 10152 10154 10155 10156 10157 10158 10159 10160 10161 10162 10163 10164
10165 10166 10167 10168 10169 10170 10171 10172 10173 10174 10175 10176 10177 10178 10179 10180 10181 10182 10183 10184
10185 10186 10187 10188 10189 10190 10191 10192 10193 10194 10195 10196 10197 10198 10199 10200 10201 10202 10203 10204
10205 10206 10207 10208 10209 10210 10211 10212 10213 10214 10215 10216 10217 10218 10219 10220 10221 10222 10223 10224
10225 10226 10227 10228 10229 10230 10231 10232 10233 10234 10235 10236 10237 10238 10239 10240 10241 10242 10243 10244
10245 10246 10247 10248 10249 10250 10251 10252 10253 10254 10255 10256 10257 10258 10259 10260 10261 10262 10263 10264
10265 10266 10267 10268 10269 10270 10271 10272 10273 10274 10275 10276 10277 10278 10279 10280 10281 10282 10283 10284
10285 10286 10287 10288 10289 10290 10291 10292 10293 10294 10295 10296 10297 10298 10299 10300 10301 10302 10303 10304
10306
Found 901 of 10306 columns with poor signal
iterative correction
- copying matrix
- computing baises
rescaling to factor 1
- getting the sum of the matrix
=> 10667.025
- rescaling biases
* 100000
WARNING: removing columns having less than 20 counts:
2 4 5 1203 1218 1220 1222 1226 1227 1228 1229 1230 1231 1232 1233 1234 1235 1236 1237 1238
1239 1240 1241 1242 1243 1244 1245 1246 1247 1248 1253 1254 1255 1256 1257 1258 1259 1260 1261 1262
1263 1264 1265 1266 1267 1268 1269 1270 1271 1272 1273 1274 1275 1276 1277 1278 1279 1280 1281 1282
1283 1284 1285 1286 1287 1288 1289 1290 1291 1292 1293 1294 1295 1296 1297 1298 1299 1300 1301 1302
1303 1304 1305 1306 1307 1308 1309 1310 1311 1312 1313 1314 1315 1316 1317 1318 1319 1320 1321 1322
1323 1324 1325 1326 1327 1328 1329 1330 1331 1332 1333 1334 1335 1336 1337 1338 1339 1340 1341 1342
1343 1344 1345 1346 1347 1348 1349 1350 1351 1352 1353 1354 1355 1356 1357 1358 1359 1360 1361 1362
1363 1364 1365 1366 1367 1368 1369 1370 1371 1372 1373 1374 1375 1376 1377 1378 1379 1380 1381 1382
1383 1384 1385 1386 1387 1388 1389 1390 1391 1392 1393 1394 1395 1396 1397 1398 1399 1400 1401 1402
1403 1404 1405 1406 1407 1408 1409 1410 1411 1412 1413 1414 1415 1416 1417 1418 1419 1420 1421 1422
1423 1424 1425 1426 1427 1428 1429 1430 1431 1432 1433 1434 1435 1441 1442 1443 1444 1447 1448 1449
1454 1462 1463 1488 1496 3385 3395 3396 3397 3398 3399 3400 3401 3402 3403 3404 3434 3435 3465 3589
3590 3594 3595 3596 5829 5831 5832 5833 5834 5835 5836 5837 5838 5839 5840 5841 5842 5844 5845 5846
5848 6895 6988 7393 7394 7395 7396 7397 7398 7399 7400 7401 7402 7403 7404 7405 7406 7407 7408 7409
7411 7412 7413 8798 9264 9266 9267 9271 9272 9273 9274 9275 9276 9277 9278 9279 9280 9281 9282 9283
9284 9285 9286 9287 9288 9289 9290 9291 9292 9293 9294 9295 9296 9297 9298 9299 9495 9496 9497 9498
9499 9500 9501 9502 9503 9504 9506 9507 9508 9510 9511 10612 10613 10614 11200 11201 11202 11203 11204 11205
11206 11207 11208 11209 11210 11211 11212 11213 11214 11215 11216 12323 12906 12907 12909 12910 12911 12912 12913 12914
12915 12916 12917 12918 12919 12920 12921 12922 12923 12924 12925 12926 12927 12928 12929 12930 12931 12932 12937 12939
12944 13054 13072 13075 13076 13077 13918 13919 13990 13993 13994 13995 13998 14312 14358 14359 14360 14361 14364 14365
14366 14367 14368 14370 14371 14373 14374 14375 14376 14775 15369 15578 15579 15580 15581 15582 15583 15584 15585 15586
15587 15588 15589 15590 15591 15592 15593 15594 15595 15596 15597 15598 15599 15600 15601 15602 15603 15604 15605 15606
15607 15608 15609 15610 15611 15612 15613 15614 15615 15616 15617 15618 15619 15620 15621 15622 15623 15624 15625 15626
15627 15628 15629 15630 15631 15632 15633 15634 15656 15657 15658 15659 15660 15661 15662 15663 15664 15665 15666 15667
15668 15669 15670 15671 15672 15673 15674 15675 15676 15677 15678 15679 15680 15681 15682 15683 15684 15685 15762 15763
15766 15769 15772 15773 15774 15776 15781 15796 15804 15805 15806 15807 15808 15809 15810 15811 15812 15813 15814 15815
15816 15817 15818 15819 15820 15821 15822 15823 15824 15825 15826 15827 15828 15829 15830 15831 15832 15833 15834 15835
15836 15837 15838 15839 15840 15841 15842 15843 15844 15845 15846 15847 15848 15849 15850 15851 15852 15853 15854 15855
15856 15857 15858 15859 15860 15861 15862 15863 15864 15865 15866 15867 15868 15869 15870 15871 15872 15873 15874 15875
15876 15877 15878 15879 15880 15881 15882 15883 15884 15885 15886 15887 15888 15889 15890 15891 15892 15893 15894 15895
15896 15897 15898 15899 15900 15901 15902 15903 15904 15905 15906 15907 15908 15909 15910 15911 15912 15913 15914 15915
15916 15917 15918 15919 15920 15921 15922 15923 15924 15925 15926 15927 15928 15929 15930 15931 15932 15933 15934 15935
15936 15937 15938 15939 15940 15941 15942 15943 15944 15945 15946 15947 15948 15949 15950 15951 15952 15953 15954 15955
15956 15957 15958 15959 15960 15961 15962 15963 15964 15965 15966 15967 15968 15969 15970 15971 15972 15973 15974 15976
15977 15978 15979 15980 15982 15984 15988 15989 15990 15991 15992 15993 15994 15995 15996 15997 15999 16000 16001 16002
16004 16012 16015 16016 16017 16021 16022 16027 16028 16029 16030 16032 16034 16035 16037 16042 16043 16044 16045 16046
16049 16050 16051 16052 16752 16753 17139 17142 17150 17151 17152 17153 17154 17155 17156 17157 17158 17159 17160 17161
17162 17163 17164 17165 17166 17167 17169 17170 17173 17174 17232 17773 18091 18092 18601 18602 18603 18604 18605 18606
18607 18608 18609 18610 18612 18613 18614 18615 18616 18617 18618 18619 18620 18621 18622 18623 18624 18625 18626 18627
18628 18629 18630 18631 18632 18633 18634 18635 18636 20776 20777 20778 20779 20780 20781 20782 20783 20784 20785 20786
20787 20788 20789 20790 20791 20792 20793 20794 20795 20796 20797 20798 20799 20800 20801 20802 20803 20804 20805 20806
20807 20808 20809 20810 20811 20812 20813 20814 20815 20816 20817 20818 20819 20820 20821 20822 20823 20824 20825 20826
20827 20828 20829 20830 20831 20832 20833 20834 20835 20836 20837 20838 20839 20840 20841 20842 20843 20844 20845 20846
20847 20848 20849 20850 20851 20852 20853 20854 20855 20856 20857 20858 20859 20860 20861 20862 20863 20864 20865 20866
20867 20868 20869 20870 20871 20872 20873 20874 20875 20876 20877 20878 20879 20880 20881 20882 20883 20884 20885 20886
20887 20888 20889 20890 20891 20892 20893 20894 20895 20896 20897 20898 20899 20900 20901 20902 20903 20904 20905 20906
20907 20908 20909 20910 20911 20912 20913 20914 20915 20916 20917 20918 20919 20920 20921 20922 20923 20924 20925 20926
20927 20928 20929 20930 20931 20932 20933 20934 20935 20936 20937 20938 20939 20940 20941 20942 20943 20944 20945 20946
20947 20948 20949 20950 20951 20952 20953 20954 20955 20956 21893 21920 21921 21922 21923 21924 21925 21926 21927 21928
21929 21930 21931 21932 21933 21934 21935 21936 21937 21938 21939 21940 21941 21942 21943 21944 21945 21946 21947 21948
21949 21950 21951 21952 21953 21954 21955 21956 21957 21958 21959 21960 21961 21962 21963 21964 21965 21966 21967 21968
21969 21970 21971 21972 21973 21974 21975 21976 21977 21978 21979 21980 21981 21982 21983 21984 21985 21986 21987 21988
21989 21990 21991 21992 21993 21994 21995 21996 21997 21998 21999 22000 22001 22002 22003 22004 22005 22006 22007 22008
22009 22010 22011 22012 22013 22014 22015 22016 22017 22018 22019 22020 22021 22022 22023 22024 22025 22026 22027 22028
22029 22030 22031 22032 22033 22034 22035 22036 22037 22038 22039 22040 22041 22042 22043 22044 22045 22046 22047 22048
22049 22050 22051 22052 22053 22054 22055 22056 22057 22058 22059 22060 22061 22062 22063 22064 22065 22066 22067 22068
22069 22070 22071 22072 22073 22074 22075 22076 22077 22078 22079 22080 22081 22082 22083 22084 22085 22086 22087 22088
22089 22090 22091 22092 22093 22094 22095 22096 22097 22098 22099 22100 22101 22104 22105 22107 22108 22115 22989 22990
22991 22992 22993 22994 22995 22996 22997 22998 22999 23000 23001 23002 23003 23004 23005 23006 23007 23008 23009 23010
23011 23012 23013 23014 23015 23016 23017 23018 23019 23020 23021 23022 23023 23024 23025 23026 23027 23028 23029 23030
23031 23032 23033 23034 23035 23036 23037 23038 23039 23040 23041 23042 23043 23044 23045 23046 23047 23048 23049 23050
23051 23052 23053 23054 23055 23056 23057 23058 23059 23060 23061 23062 23063 23064 23065 23066 23067 23068 23069 23070
23071 23072 23073 23074 23075 23076 23077 23078 23079 23080 23081 23082 23083 23084 23085 23086 23087 23088 23089 23090
23091 23092 23093 23094 23095 23096 23097 23098 23099 23100 23101 23102 23103 23104 23105 23106 23107 23108 23109 23110
23111 23112 23113 23114 23115 23116 23117 23118 23119 23120 23121 23122 23123 23124 23125 23126 23127 23128 23129 23130
23131 23132 23133 23134 23135 23136 23137 23138 23139 23140 23141 23142 23143 23144 23145 23146 23147 23148 23149 23150
23151 23152 23153 23154 23155 23156 23157 23158 23159 23160 23163 23164 23165 23166 23167 23168 23169 23170 23171 23172
23173 23174 23175 23176 23177 23178 23179 23180 23181 23182 23183 23184 23185 23186 23187 23202 23208 23209 23215 23223
23296 23833 23835 24174 24194 24195 24298 24306 24338 24342 24343 24344 24354 24355 24356 24375 24376 24377 24378 24379
24380 24381 24382 24383 24384 24385 24386 24387 24388 24389 24390 24391 24392 24393 24394 24395 24396 24397 24398 24399
24400 24401 24402 24403 24404 24405 24406 24407 24408 24409 24410 24411 24412 24413 24414 24415 24416 24417 24418 24419
24420 24421 24422 24423 24424 24425 24426 24427 24428 24429 24430 24431 24432 24433 24434 24435 24436 24437 24438 24439
24440 24441 24442 24443 24444 24445 24446 24447 24448 24449 24450 24451 24452 24453 24454 24455 24456 24457 24458 24459
24460 24461 24462 24463 24464 24465 24466 24467 24468 24469 24470 24471 24472 24473 24474 24914 24915 25106 25134 25150
25152 25153 25161 25163 25294 25295 25296 25378 25381 25747 25906 25907 25908 25909 25910 25911 25912 25913 25914 25915
25916 25917 25918 25919 25920 25921 25922 25923 25924 25925 25926 25927 25928 25929 25930 25931 25932 25933 25934 25935
25936 25937 25938 25939 25940 25941 25942 25943 25944 25945 25946 25947 25948 25949 25950 25951 25952 25954 26551 26552
26553 26797 26798 26799 26800 26801 26802 26803 26804 26805 26806 26807 26808 26809 26810 26811 26812 26813 26814 26815
26816 26817 26818 26819 26820 26821 26822 26823 27138 27403 27404 27405 27406 27407 27408 27409 27410 27411 27412 27413
27414 27415 27416 27417 27418 27419 27420 27421 27422 27423 27438 27439 27449 27450 27783 27784 27785 27786 27787 27788
27789 27790 27791 27792 27793 27794 27795 27796 27797 27798 27799 27800 27801 27802 27803 27804 27805 27806 27807 27808
27809 27810 27811 27812 27813 27814 27815 27816 27817 27818 27819 27820 27821 27822 27823 27824 27825 27826 27827 27828
27829 27830 27831 27832 27833 27835 27838 27839 27840 27841 27842 27843 27844 27845 27846 27848 27849 27850 27851 27853
27854 27855 27858 27859 27860 27861 27862 27865 27866 27867 27868 27869 27878 27886 27893 27894 27895 27896 27897 27898
27899 27900 27901 27902 27903 27904 27905 27906 27907 27908 27909 27910 27911 27912 28218 28251 28252 28253 28254 28255
28256 28257 28258 28259 28260 28261 28262 28263 28264 28265 28266 28267 28268 28269 28270 28271 28272 28273 28274 28275
28276 28277 28278 28279 28280 28281 28282 28283 28284 28285 28286 28287 28288 28289 28290 28291 28292 28293 28294 28295
28296 28297 28298 28299 28300 28301 28302 28303 28304 28305 28306 28307 28308 28309 28310 28311 28312 28313 28314 28315
28316 28317 28318 28319 28320 28321 28322 28323 28324 28325 28326 28327 28328 28329 28330 28331 28332 28333 28334 28335
28336 28337 28338 28339 28340 28341 28342 28343 28344 28345 28346 28347 28348 28349 28350 28351 28352 28353 28354 28355
28356 28357 28360 28381 28382 28383 28384 28385 28386 28387 28388 28389 28390 28391 28392 28393 28394 28395 28396 28397
28398 28399 28400 28401 28402 28403 28434 28436 28437 28438 28760 28761 28762 28763 28764 28765 28766 28767 28768 28769
28770 28771 28772 28773 28774 28775 28776 28777 28778 28779 28780 28781 28782 28783 28784 28785 28786 28787 28788 29132
29133 29255 29256 29347 29348 29349 29350 29351 29352 29353 29354 29355 29356 29357 29358 29359 29360 29361 29362 29363
29364 29365 29366 29367 29368 29369 29370 29371 29372 29373 29374 29375 29376 29377 29378 29379 29380 29381 29382 29383
29384 29385 29478 29489 29783 30303 30319 30320 30321 30322 30323 30324 30325 30326 30327 30328 30329 30330 30331 30332
30333 30334 30335 30336 30337 30338 30339 30340 30341 30342 30343 30344 30345 30346 30347 30348 30349 30350 30351 30352
30353 30354 30355 30356 30357 30358 30359 30360 30361 30362 30363 30364 30365 30366 30367 30368 30369 30370 30371 30372
30373 30374 30375 30376 30377 30378 30379 30380 30381 30382 30383 30384 30385 30386 30387 30388 30389 30390 30391 30392
30393 30394 30395 30396 30397 30398 30399 30400 30401 30402 30403 30404 30405 30406 30407 30408 30409 30410 30411 30412
30413 30414 30415 30416 30417 30418 30419 30420 30421 30425 30426 30427 30428 30429 30430 30431 30432 30437 30438 30439
30440 30441 30442 30443 30444 30445 30446 30447 30448 30449 30450 30451 30452 30453 30454 30455 30456 30457 30458 30459
30460 30461 30462 30463 30464 30465 30466 30467 30468 30469 30470 30471 30472 30473 30474 30475 30476 30477 30478 30479
30480 30481 30482 30483 30484 30485 30486 30487 30488 30489 30490 30491 30492 30493 30494 30495 30496 30497 30498 30499
30500 30501 30502 30503 30504 30505 30506 30507 30508 30509 30510 30511 30512 30513 30514 30515 30516 30517 30518 30519
30520 30521 30522 30523 30524 30525 30526 30527 30528 30529 30530 30531 30532 30533 30534 30535 30536 30537 30538 30539
30540 30541 30542 30543 30544 30545 30546 30547 30548 30549 30550 30551 30552 30553 30554 30555 30556 30557 30558 30559
30560 30561 30562 30563 30564 30565 30566 30567 30568 30569 30570 30571 30572 30573 30574 30575 30576 30577 30578 30579
30580 30581 30582 30583 30584 30585 30586 30587 30588 30589 30590 30591 30592 30593 30594 30595 30596 30597 30598 30599
30600 30601 30602 30603 30604 30605 30606 30607 30608 30609 30610 30611 30612 30613 30614 30615 30616 30617 30618 30619
30620 30621 30622 30623 30624 30625 30626 30627 30628 30629 30630 30631 30632 30633 30634 30635 30636 30637 30638 30639
30640 30641 30642 30643 30644 30645 30646 30647 30648 30649 30650 30651 30652 30653 30654 30655 30656 30657 30658 30659
30660 30661 30662 30663 30664 30665 30666 30667 30668 30669 30670 30671 30672 30673 30674 30675 30676 30677 30678 30679
30680 30681 30682 30683 30684 30685 30686 30687 30688 30689 30690 30691 30692 30693 30694 30695 30696 30697 30698 30699
30700 30701 30702 30703 30704 30705 30706 30707 30708 30709 30710 30711 30712 30713 30714 30715 30716 30717 30718 30719
30720 30721 30722 30723 30724 30725 30726 30727 30728 30729 30730 30731 30732 30733 30734 30735 30736 30737 30738 30739
30740 30741 30742 30743 30744 30745 30746 30747 30748 30749 30750 30751 30752 30753 30754 30755 30756 30757 30758 30759
30760 30761 30762 30763 30764 30765 30766 30767 30768 30769 30770 30771 30772 30773 30774 30775 30776 30777 30778 30779
30780 30781 30782 30783 30784 30785 30786 30787 30788 30789 30790 30791 30792 30793 30794 30795 30796 30797 30798 30799
30800 30801 30802 30803 30804 30805 30806 30807 30808 30809 30810 30811 30812 30813 30814 30815 30816 30817 30818 30819
30820 30821 30822 30823 30824 30825 30826 30827 30828 30829 30830 30831 30832 30833 30834 30835 30836 30837 30838 30839
30840 30841 30842 30843 30844 30845 30846 30847 30848 30849 30850 30851 30852 30853 30854 30855 30856 30857 30858 30859
30860 30861 30862 30863 30864 30865 30866 30867 30868 30869 30870 30871 30872 30873 30874 30875 30876 30877 30878 30879
30880 30881 30882 30883 30884 30885 30886 30887 30888 30889 30891 30892 30893 30894
WARNING: removing columns having less than 79.667 counts:
1 2 3 4 5 6 7 8 10 13 14 28 130 1203 1207 1208 1209 1213 1214 1217
1218 1219 1220 1221 1222 1223 1224 1225 1226 1227 1228 1229 1230 1231 1232 1233 1234 1235 1236 1237
1238 1239 1240 1241 1242 1243 1244 1245 1246 1247 1248 1249 1253 1254 1255 1256 1257 1258 1259 1260
1261 1262 1263 1264 1265 1266 1267 1268 1269 1270 1271 1272 1273 1274 1275 1276 1277 1278 1279 1280
1281 1282 1283 1284 1285 1286 1287 1288 1289 1290 1291 1292 1293 1294 1295 1296 1297 1298 1299 1300
1301 1302 1303 1304 1305 1306 1307 1308 1309 1310 1311 1312 1313 1314 1315 1316 1317 1318 1319 1320
1321 1322 1323 1324 1325 1326 1327 1328 1329 1330 1331 1332 1333 1334 1335 1336 1337 1338 1339 1340
1341 1342 1343 1344 1345 1346 1347 1348 1349 1350 1351 1352 1353 1354 1355 1356 1357 1358 1359 1360
1361 1362 1363 1364 1365 1366 1367 1368 1369 1370 1371 1372 1373 1374 1375 1376 1377 1378 1379 1380
1381 1382 1383 1384 1385 1386 1387 1388 1389 1390 1391 1392 1393 1394 1395 1396 1397 1398 1399 1400
1401 1402 1403 1404 1405 1406 1407 1408 1409 1410 1411 1412 1413 1414 1415 1416 1417 1418 1419 1420
1421 1422 1423 1424 1425 1426 1427 1428 1429 1430 1431 1432 1433 1434 1435 1438 1440 1441 1442 1443
1444 1446 1447 1448 1449 1450 1451 1452 1454 1455 1462 1463 1464 1465 1486 1487 1488 1489 1492 1493
1494 1495 1496 1498 3360 3364 3365 3368 3369 3384 3385 3386 3388 3394 3395 3396 3397 3398 3399 3400
3401 3402 3403 3404 3412 3427 3431 3434 3435 3465 3589 3590 3591 3594 3595 3596 3795 3796 4912 5820
5825 5828 5829 5830 5831 5832 5833 5834 5835 5836 5837 5838 5839 5840 5841 5842 5843 5844 5845 5846
5848 5849 6895 6988 6989 7387 7389 7392 7393 7394 7395 7396 7397 7398 7399 7400 7401 7402 7403 7404
7405 7406 7407 7408 7409 7410 7411 7412 7413 8797 8798 9264 9265 9266 9267 9268 9271 9272 9273 9274
9275 9276 9277 9278 9279 9280 9281 9282 9283 9284 9285 9286 9287 9288 9289 9290 9291 9292 9293 9294
9295 9296 9297 9298 9299 9495 9496 9497 9498 9499 9500 9501 9502 9503 9504 9505 9506 9507 9508 9509
9510 9511 10612 10613 10614 10615 10882 11200 11201 11202 11203 11204 11205 11206 11207 11208 11209 11210 11211 11212
11213 11214 11215 11216 12322 12323 12383 12892 12893 12905 12906 12907 12908 12909 12910 12911 12912 12913 12914 12915
12916 12917 12918 12919 12920 12921 12922 12923 12924 12925 12926 12927 12928 12929 12930 12931 12932 12936 12937 12938
12939 12944 12957 12958 12959 12960 12975 12980 13054 13055 13056 13072 13073 13075 13076 13077 13349 13350 13918 13919
13942 13990 13993 13994 13995 13997 13998 13999 14040 14312 14358 14359 14360 14361 14362 14363 14364 14365 14366 14367
14368 14369 14370 14371 14372 14373 14374 14375 14376 14775 15357 15361 15369 15370 15578 15579 15580 15581 15582 15583
15584 15585 15586 15587 15588 15589 15590 15591 15592 15593 15594 15595 15596 15597 15598 15599 15600 15601 15602 15603
15604 15605 15606 15607 15608 15609 15610 15611 15612 15613 15614 15615 15616 15617 15618 15619 15620 15621 15622 15623
15624 15625 15626 15627 15628 15629 15630 15631 15632 15633 15634 15635 15656 15657 15658 15659 15660 15661 15662 15663
15664 15665 15666 15667 15668 15669 15670 15671 15672 15673 15674 15675 15676 15677 15678 15679 15680 15681 15682 15683
15684 15685 15758 15759 15762 15763 15764 15766 15767 15769 15770 15771 15772 15773 15774 15775 15776 15777 15781 15784
15785 15786 15788 15789 15790 15791 15795 15796 15798 15800 15802 15803 15804 15805 15806 15807 15808 15809 15810 15811
15812 15813 15814 15815 15816 15817 15818 15819 15820 15821 15822 15823 15824 15825 15826 15827 15828 15829 15830 15831
15832 15833 15834 15835 15836 15837 15838 15839 15840 15841 15842 15843 15844 15845 15846 15847 15848 15849 15850 15851
15852 15853 15854 15855 15856 15857 15858 15859 15860 15861 15862 15863 15864 15865 15866 15867 15868 15869 15870 15871
15872 15873 15874 15875 15876 15877 15878 15879 15880 15881 15882 15883 15884 15885 15886 15887 15888 15889 15890 15891
15892 15893 15894 15895 15896 15897 15898 15899 15900 15901 15902 15903 15904 15905 15906 15907 15908 15909 15910 15911
15912 15913 15914 15915 15916 15917 15918 15919 15920 15921 15922 15923 15924 15925 15926 15927 15928 15929 15930 15931
15932 15933 15934 15935 15936 15937 15938 15939 15940 15941 15942 15943 15944 15945 15946 15947 15948 15949 15950 15951
15952 15953 15954 15955 15956 15957 15958 15959 15960 15961 15962 15963 15964 15965 15966 15967 15968 15969 15970 15971
15972 15973 15974 15975 15976 15977 15978 15979 15980 15981 15982 15983 15984 15985 15987 15988 15989 15990 15991 15992
15993 15994 15995 15996 15997 15999 16000 16001 16002 16003 16004 16005 16006 16009 16010 16012 16013 16014 16015 16016
16017 16018 16019 16021 16022 16023 16025 16026 16027 16028 16029 16030 16031 16032 16033 16034 16035 16036 16037 16038
16039 16040 16041 16042 16043 16044 16045 16046 16047 16048 16049 16050 16051 16052 16752 16753 17139 17142 17150 17151
17152 17153 17154 17155 17156 17157 17158 17159 17160 17161 17162 17163 17164 17165 17166 17167 17168 17169 17170 17173
17174 17231 17232 17234 17773 18091 18092 18093 18097 18600 18601 18602 18603 18604 18605 18606 18607 18608 18609 18610
18611 18612 18613 18614 18615 18616 18617 18618 18619 18620 18621 18622 18623 18624 18625 18626 18627 18628 18629 18630
18631 18632 18633 18634 18635 18636 18991 19538 19792 19793 19794 19795 19796 19798 19799 19800 19801 19802 19803 19804
19805 19807 19808 19809 19810 19811 19813 19814 20768 20776 20777 20778 20779 20780 20781 20782 20783 20784 20785 20786
20787 20788 20789 20790 20791 20792 20793 20794 20795 20796 20797 20798 20799 20800 20801 20802 20803 20804 20805 20806
20807 20808 20809 20810 20811 20812 20813 20814 20815 20816 20817 20818 20819 20820 20821 20822 20823 20824 20825 20826
20827 20828 20829 20830 20831 20832 20833 20834 20835 20836 20837 20838 20839 20840 20841 20842 20843 20844 20845 20846
20847 20848 20849 20850 20851 20852 20853 20854 20855 20856 20857 20858 20859 20860 20861 20862 20863 20864 20865 20866
20867 20868 20869 20870 20871 20872 20873 20874 20875 20876 20877 20878 20879 20880 20881 20882 20883 20884 20885 20886
20887 20888 20889 20890 20891 20892 20893 20894 20895 20896 20897 20898 20899 20900 20901 20902 20903 20904 20905 20906
20907 20908 20909 20910 20911 20912 20913 20914 20915 20916 20917 20918 20919 20920 20921 20922 20923 20924 20925 20926
20927 20928 20929 20930 20931 20932 20933 20934 20935 20936 20937 20938 20939 20940 20941 20942 20943 20944 20945 20946
20947 20948 20949 20950 20951 20952 20953 20954 20955 20956 20960 20961 21893 21920 21921 21922 21923 21924 21925 21926
21927 21928 21929 21930 21931 21932 21933 21934 21935 21936 21937 21938 21939 21940 21941 21942 21943 21944 21945 21946
21947 21948 21949 21950 21951 21952 21953 21954 21955 21956 21957 21958 21959 21960 21961 21962 21963 21964 21965 21966
21967 21968 21969 21970 21971 21972 21973 21974 21975 21976 21977 21978 21979 21980 21981 21982 21983 21984 21985 21986
21987 21988 21989 21990 21991 21992 21993 21994 21995 21996 21997 21998 21999 22000 22001 22002 22003 22004 22005 22006
22007 22008 22009 22010 22011 22012 22013 22014 22015 22016 22017 22018 22019 22020 22021 22022 22023 22024 22025 22026
22027 22028 22029 22030 22031 22032 22033 22034 22035 22036 22037 22038 22039 22040 22041 22042 22043 22044 22045 22046
22047 22048 22049 22050 22051 22052 22053 22054 22055 22056 22057 22058 22059 22060 22061 22062 22063 22064 22065 22066
22067 22068 22069 22070 22071 22072 22073 22074 22075 22076 22077 22078 22079 22080 22081 22082 22083 22084 22085 22086
22087 22088 22089 22090 22091 22092 22093 22094 22095 22096 22097 22098 22099 22100 22101 22102 22103 22104 22105 22107
22108 22109 22110 22114 22115 22989 22990 22991 22992 22993 22994 22995 22996 22997 22998 22999 23000 23001 23002 23003
23004 23005 23006 23007 23008 23009 23010 23011 23012 23013 23014 23015 23016 23017 23018 23019 23020 23021 23022 23023
23024 23025 23026 23027 23028 23029 23030 23031 23032 23033 23034 23035 23036 23037 23038 23039 23040 23041 23042 23043
23044 23045 23046 23047 23048 23049 23050 23051 23052 23053 23054 23055 23056 23057 23058 23059 23060 23061 23062 23063
23064 23065 23066 23067 23068 23069 23070 23071 23072 23073 23074 23075 23076 23077 23078 23079 23080 23081 23082 23083
23084 23085 23086 23087 23088 23089 23090 23091 23092 23093 23094 23095 23096 23097 23098 23099 23100 23101 23102 23103
23104 23105 23106 23107 23108 23109 23110 23111 23112 23113 23114 23115 23116 23117 23118 23119 23120 23121 23122 23123
23124 23125 23126 23127 23128 23129 23130 23131 23132 23133 23134 23135 23136 23137 23138 23139 23140 23141 23142 23143
23144 23145 23146 23147 23148 23149 23150 23151 23152 23153 23154 23155 23156 23157 23158 23159 23160 23162 23163 23164
23165 23166 23167 23168 23169 23170 23171 23172 23173 23174 23175 23176 23177 23178 23179 23180 23181 23182 23183 23184
23185 23186 23187 23191 23196 23198 23199 23201 23202 23204 23205 23206 23207 23208 23209 23210 23211 23212 23214 23215
23216 23222 23223 23274 23275 23276 23277 23292 23294 23295 23296 23313 23314 23315 23316 23833 23835 24021 24032 24158
24174 24175 24177 24193 24194 24195 24225 24295 24298 24305 24306 24307 24313 24331 24332 24338 24340 24341 24342 24343
24344 24345 24354 24355 24356 24359 24360 24375 24376 24377 24378 24379 24380 24381 24382 24383 24384 24385 24386 24387
24388 24389 24390 24391 24392 24393 24394 24395 24396 24397 24398 24399 24400 24401 24402 24403 24404 24405 24406 24407
24408 24409 24410 24411 24412 24413 24414 24415 24416 24417 24418 24419 24420 24421 24422 24423 24424 24425 24426 24427
24428 24429 24430 24431 24432 24433 24434 24435 24436 24437 24438 24439 24440 24441 24442 24443 24444 24445 24446 24447
24448 24449 24450 24451 24452 24453 24454 24455 24456 24457 24458 24459 24460 24461 24462 24463 24464 24465 24466 24467
24468 24469 24470 24471 24472 24473 24474 24907 24912 24913 24914 24915 25099 25100 25105 25106 25119 25134 25143 25144
25145 25146 25147 25148 25149 25150 25151 25152 25153 25155 25157 25158 25159 25160 25161 25162 25163 25164 25165 25167
25168 25169 25170 25171 25172 25173 25175 25177 25178 25179 25278 25294 25295 25296 25370 25378 25379 25380 25381 25732
25734 25747 25906 25907 25908 25909 25910 25911 25912 25913 25914 25915 25916 25917 25918 25919 25920 25921 25922 25923
25924 25925 25926 25927 25928 25929 25930 25931 25932 25933 25934 25935 25936 25937 25938 25939 25940 25941 25942 25943
25944 25945 25946 25947 25948 25949 25950 25951 25952 25954 25955 25956 26551 26552 26553 26556 26558 26559 26560 26561
26562 26563 26564 26572 26576 26586 26597 26609 26619 26656 26663 26681 26683 26692 26724 26732 26735 26797 26798 26799
26800 26801 26802 26803 26804 26805 26806 26807 26808 26809 26810 26811 26812 26813 26814 26815 26816 26817 26818 26819
26820 26821 26822 26823 26940 27029 27030 27037 27089 27138 27399 27403 27404 27405 27406 27407 27408 27409 27410 27411
27412 27413 27414 27415 27416 27417 27418 27419 27420 27421 27422 27423 27438 27439 27449 27450 27782 27783 27784 27785
27786 27787 27788 27789 27790 27791 27792 27793 27794 27795 27796 27797 27798 27799 27800 27801 27802 27803 27804 27805
27806 27807 27808 27809 27810 27811 27812 27813 27814 27815 27816 27817 27818 27819 27820 27821 27822 27823 27824 27825
27826 27827 27828 27829 27830 27831 27832 27833 27834 27835 27838 27839 27840 27841 27842 27843 27844 27845 27846 27848
27849 27850 27851 27852 27853 27854 27855 27858 27859 27860 27861 27862 27864 27865 27866 27867 27868 27869 27870 27878
27886 27893 27894 27895 27896 27897 27898 27899 27900 27901 27902 27903 27904 27905 27906 27907 27908 27909 27910 27911
27912 27921 28218 28225 28251 28252 28253 28254 28255 28256 28257 28258 28259 28260 28261 28262 28263 28264 28265 28266
28267 28268 28269 28270 28271 28272 28273 28274 28275 28276 28277 28278 28279 28280 28281 28282 28283 28284 28285 28286
28287 28288 28289 28290 28291 28292 28293 28294 28295 28296 28297 28298 28299 28300 28301 28302 28303 28304 28305 28306
28307 28308 28309 28310 28311 28312 28313 28314 28315 28316 28317 28318 28319 28320 28321 28322 28323 28324 28325 28326
28327 28328 28329 28330 28331 28332 28333 28334 28335 28336 28337 28338 28339 28340 28341 28342 28343 28344 28345 28346
28347 28348 28349 28350 28351 28352 28353 28354 28355 28356 28357 28358 28360 28363 28369 28374 28381 28382 28383 28384
28385 28386 28387 28388 28389 28390 28391 28392 28393 28394 28395 28396 28397 28398 28399 28400 28401 28402 28403 28404
28405 28406 28409 28410 28434 28435 28436 28437 28438 28439 28664 28711 28754 28760 28761 28762 28763 28764 28765 28766
28767 28768 28769 28770 28771 28772 28773 28774 28775 28776 28777 28778 28779 28780 28781 28782 28783 28784 28785 28786
28787 28788 28799 28998 29132 29133 29255 29256 29288 29346 29347 29348 29349 29350 29351 29352 29353 29354 29355 29356
29357 29358 29359 29360 29361 29362 29363 29364 29365 29366 29367 29368 29369 29370 29371 29372 29373 29374 29375 29376
29377 29378 29379 29380 29381 29382 29383 29384 29385 29478 29489 29783 29784 29918 29919 29970 30119 30171 30176 30202
30303 30306 30318 30319 30320 30321 30322 30323 30324 30325 30326 30327 30328 30329 30330 30331 30332 30333 30334 30335
30336 30337 30338 30339 30340 30341 30342 30343 30344 30345 30346 30347 30348 30349 30350 30351 30352 30353 30354 30355
30356 30357 30358 30359 30360 30361 30362 30363 30364 30365 30366 30367 30368 30369 30370 30371 30372 30373 30374 30375
30376 30377 30378 30379 30380 30381 30382 30383 30384 30385 30386 30387 30388 30389 30390 30391 30392 30393 30394 30395
30396 30397 30398 30399 30400 30401 30402 30403 30404 30405 30406 30407 30408 30409 30410 30411 30412 30413 30414 30415
30416 30417 30418 30419 30420 30421 30422 30423 30424 30425 30426 30427 30428 30429 30430 30431 30432 30436 30437 30438
30439 30440 30441 30442 30443 30444 30445 30446 30447 30448 30449 30450 30451 30452 30453 30454 30455 30456 30457 30458
30459 30460 30461 30462 30463 30464 30465 30466 30467 30468 30469 30470 30471 30472 30473 30474 30475 30476 30477 30478
30479 30480 30481 30482 30483 30484 30485 30486 30487 30488 30489 30490 30491 30492 30493 30494 30495 30496 30497 30498
30499 30500 30501 30502 30503 30504 30505 30506 30507 30508 30509 30510 30511 30512 30513 30514 30515 30516 30517 30518
30519 30520 30521 30522 30523 30524 30525 30526 30527 30528 30529 30530 30531 30532 30533 30534 30535 30536 30537 30538
30539 30540 30541 30542 30543 30544 30545 30546 30547 30548 30549 30550 30551 30552 30553 30554 30555 30556 30557 30558
30559 30560 30561 30562 30563 30564 30565 30566 30567 30568 30569 30570 30571 30572 30573 30574 30575 30576 30577 30578
30579 30580 30581 30582 30583 30584 30585 30586 30587 30588 30589 30590 30591 30592 30593 30594 30595 30596 30597 30598
30599 30600 30601 30602 30603 30604 30605 30606 30607 30608 30609 30610 30611 30612 30613 30614 30615 30616 30617 30618
30619 30620 30621 30622 30623 30624 30625 30626 30627 30628 30629 30630 30631 30632 30633 30634 30635 30636 30637 30638
30639 30640 30641 30642 30643 30644 30645 30646 30647 30648 30649 30650 30651 30652 30653 30654 30655 30656 30657 30658
30659 30660 30661 30662 30663 30664 30665 30666 30667 30668 30669 30670 30671 30672 30673 30674 30675 30676 30677 30678
30679 30680 30681 30682 30683 30684 30685 30686 30687 30688 30689 30690 30691 30692 30693 30694 30695 30696 30697 30698
30699 30700 30701 30702 30703 30704 30705 30706 30707 30708 30709 30710 30711 30712 30713 30714 30715 30716 30717 30718
30719 30720 30721 30722 30723 30724 30725 30726 30727 30728 30729 30730 30731 30732 30733 30734 30735 30736 30737 30738
30739 30740 30741 30742 30743 30744 30745 30746 30747 30748 30749 30750 30751 30752 30753 30754 30755 30756 30757 30758
30759 30760 30761 30762 30763 30764 30765 30766 30767 30768 30769 30770 30771 30772 30773 30774 30775 30776 30777 30778
30779 30780 30781 30782 30783 30784 30785 30786 30787 30788 30789 30790 30791 30792 30793 30794 30795 30796 30797 30798
30799 30800 30801 30802 30803 30804 30805 30806 30807 30808 30809 30810 30811 30812 30813 30814 30815 30816 30817 30818
30819 30820 30821 30822 30823 30824 30825 30826 30827 30828 30829 30830 30831 30832 30833 30834 30835 30836 30837 30838
30839 30840 30841 30842 30843 30844 30845 30846 30847 30848 30849 30850 30851 30852 30853 30854 30855 30856 30857 30858
30859 30860 30861 30862 30863 30864 30865 30866 30867 30868 30869 30870 30871 30872 30873 30874 30875 30876 30877 30878
30879 30880 30881 30882 30883 30884 30885 30886 30887 30888 30889 30891 30892 30893 30894
Found 2815 of 30894 columns with poor signal
iterative correction
- copying matrix
- computing baises
rescaling to factor 1
- getting the sum of the matrix
=> 31879.088
- rescaling biases
In [ ]:
Content source: 4DGenome/Chromosomal-Conformation-Course
Similar notebooks: