diff --git a/R/ANEVA_DOT.R b/R/ANEVA_DOT.R index 99a8010..08a75c1 100644 --- a/R/ANEVA_DOT.R +++ b/R/ANEVA_DOT.R @@ -70,7 +70,6 @@ ANEVA_DOT<-function(ASEdat, output_columns = c("refCount","altCount"), eh1 = "re return(output) } - #' @param Eg_std Standard deviation in the log2 transformed total gene expression in a healthy #' population. #' @param eh1 Integer value for of expression of haplotype 1. @@ -154,6 +153,27 @@ Binom_test_ctm_dbl<-function(X,N,p1,p2,log_BinCoeff,r0){ return(p.val) } +#' This function calculates Binomial coefficients for choosing 0:n out of n, +#' in log scale. It is useful for precalculating the coefficients and +#' providing them to "Binom_test_fast.R" or "pdf_binom_fast.R" when they +#' are redone on the same N, over and over again. This function was +#' designed by Pejman Mohammadi, 2018, Scripps Research, San Diego, CA. + +log_BinCoeffs<-function(n){ + log_z<-numeric(n+1) #note first and last element will remain 0 + + lgamm<-lgamma(0:n+1) #precalculate gamma function + + m1<-floor(n/2) #this is the middle. binom coeffs symmetric so we calculate only half + m2<-ceiling(n/2) #this is the "other" middle + x<- 1:m1 + xc<- n-x + log_z[2:(m1+1)]<- (lgamm[n+1]-lgamm[x+1]-lgamm[xc+1]) + + log_z[n:(m2+1)]<-log_z[2:(m1+1)] #fill in the other side + return(log_z) +} + #' ################################################################################################# #' Generate plot of ASE data showing significant outliers. (This Function not currently implemented) #' diff --git a/R/ANEVAdot.R b/R/ANEVAdot.R index 6eae1dc..002b938 100644 --- a/R/ANEVAdot.R +++ b/R/ANEVAdot.R @@ -58,10 +58,10 @@ ANEVAdot<-function(filepath, output_columns = c("refCount","altCount"), eh1 = "r output$p.val[i]<-NA next } - if (dat[i,eh1]==dat[i,eh2]){ + if (dat[i,eh1]==dat[i,eh2] & (dat[i,eh1]+dat[i,eh2])>10){ output$p.val[i]<-1 } - else if ((dat[i,eh1]+dat[i,eh2])>8){ + else if ((dat[i,eh1]+dat[i,eh2])>10){ Hh<-max(dat[i,eh1],dat[i,eh2]) #higher expressed haplotype Lh<-min(dat[i,eh1],dat[i,eh2]) #lower expressed haplotype rad<-Eg_std[i]*4 #integration radius diff --git a/R/Evaluate_ANEVAdot.R b/R/Evaluate_ANEVAdot.R index bff144f..018533a 100644 --- a/R/Evaluate_ANEVAdot.R +++ b/R/Evaluate_ANEVAdot.R @@ -23,14 +23,9 @@ test.filepath<-"testdata.txt" ####add column titles manually#### -<<<<<<< HEAD start_time <- Sys.time() test.result<-ANEVAdot(filepath = test.filepath, output_columns = c("eh1","eh2"),eh1 = "eh1", eh2 = "eh2",Eg_std=Sg) end_time <- Sys.time() runtime<-end_time - start_time qqplot(runif(10000),test.result$p.val) -======= -test.result<-ANEVAdot(filepath = test.filepath, output_columns = c("eh1","eh2"), eh1 = "eh1", eh2 = "eh2",Eg_std=Sg) -qqplot(runif(1000),test.result$p.val) ->>>>>>> e615cc169111775374b560ddc35e4c3f964d07b2 diff --git a/testdata.txt b/testdata.txt index 1dd46ea..db28fb4 100644 --- a/testdata.txt +++ b/testdata.txt @@ -10001,834 +10001,1008 @@ 343 256 ======= "eh1" "eh2" -1686 1174 -669 495 -447 352 -1023 986 -544 435 -736 739 -981 911 -757 596 -724 819 -500 459 -666 596 -1219 1266 -270 292 -332 398 -375 322 -348 718 -269 292 -380 437 -322 379 -109 79 -271 356 -447 436 -1394 1087 -903 951 -1663 1760 -343 384 -361 275 -189 213 -846 865 -324 308 -285 292 -293 293 -639 639 -509 447 -1004 1208 -395 358 -601 581 -291 258 -375 473 -240 221 -407 367 -260 303 -249 229 -656 696 -461 391 -437 373 -644 736 -145 167 -104 108 -391 311 -304 334 -288 286 -692 727 -418 354 -684 789 -441 523 -1041 1028 -510 594 -683 695 -300 344 -855 771 -307 652 -918 856 -385 358 -461 437 -530 789 -341 513 -939 861 -512 372 -592 598 -930 1181 -328 260 -501 543 -717 561 -308 280 -387 368 -341 369 -1374 1130 -498 481 -280 320 -773 500 -520 575 -223 271 -1885 2241 -500 436 -695 594 -155 119 -326 174 -509 515 -313 383 -445 694 -298 400 -813 732 -270 252 -481 486 -279 238 -368 163 -278 254 -314 346 -850 909 -551 557 -582 449 -408 378 -533 585 -538 683 -444 460 -321 500 -2114 2257 -379 473 -684 648 -675 590 -289 332 -898 714 -290 185 -673 551 -690 553 -483 492 -695 636 -842 979 -467 341 -280 313 -1292 1270 -471 482 -303 175 -287 249 -219 247 -185 217 -344 363 -663 634 -615 666 -749 712 -797 976 -226 191 -466 504 -850 687 -343 516 -2160 1961 -270 245 -547 444 -859 906 -580 453 -483 668 -602 643 -324 290 -1042 881 -512 434 -273 288 -454 444 -1203 1472 -457 458 -796 734 -495 660 -426 413 -843 819 -245 175 -395 387 -1157 1043 -707 731 -670 530 -515 540 -738 646 -200 183 -738 580 -1696 1311 -507 407 -341 434 -537 511 -384 329 -412 510 -485 485 -546 457 -606 752 -929 1354 -401 510 -668 660 -522 605 -231 348 -414 253 -394 371 -1068 935 -414 566 -755 706 -87 126 -638 770 -370 586 -318 303 -341 665 -145 150 -562 690 -156 178 -700 755 -378 507 -822 885 -752 667 -859 676 -740 871 -593 592 -656 535 -585 583 -498 792 -1109 920 -225 324 -572 737 -545 388 -436 550 -468 454 -327 370 -439 428 -416 253 -626 582 -418 489 -485 602 -435 459 -268 212 -723 712 -186 243 -1039 1129 -2126 1970 -680 740 -408 511 -287 283 -433 575 -711 654 -734 473 -305 223 -234 233 -167 151 -415 452 -563 525 -408 366 -453 553 -99 172 -631 633 -711 564 -807 774 -341 282 -492 585 -813 825 -491 597 -320 364 -278 299 -529 431 -493 553 -530 566 -152 260 -1006 872 -441 460 -716 616 -647 679 -136 159 -478 518 -722 593 -1020 1080 -453 439 -245 259 -407 749 -652 621 -717 584 -207 312 -415 450 -484 482 -151 175 -984 670 -500 423 -458 559 -1458 1063 -417 465 -212 199 -568 609 -573 588 -416 289 -167 250 -828 590 -745 802 -575 627 -442 559 -594 609 -503 509 -333 417 -663 645 -578 682 -248 220 -294 286 -312 393 -257 224 -962 1037 -274 278 -557 1090 -790 1047 -1557 1473 -524 511 -607 908 -457 494 -325 248 -276 277 -1018 1012 -107 135 -389 419 -844 1153 -564 590 -794 793 -299 331 -362 343 -578 548 -788 669 -575 675 -1140 1389 -999 1151 -289 397 -414 487 -473 447 -2375 2491 -1228 1313 -1160 1247 -584 616 -898 791 -998 680 -172 160 -303 388 -619 770 -583 507 -838 680 -287 287 -224 205 -593 585 -267 245 -473 427 -253 607 -1003 923 -483 428 -894 953 -891 584 -488 341 -101 103 -1341 1100 -600 472 -187 198 -658 743 -627 636 -359 316 -219 217 -289 271 -400 456 -487 476 -373 517 -450 599 -163 255 -430 443 -1020 884 -383 299 -669 646 -389 412 -153 118 -116 161 -409 394 -679 563 -499 530 -919 837 -1117 988 -498 530 -986 926 -243 229 -1868 1542 -2233 1755 -880 860 -464 511 -1088 1183 -760 1025 -550 606 -918 887 -490 451 -733 659 -994 1092 -480 555 -649 662 -658 415 -538 556 -921 1028 -1023 1041 -1384 1396 -157 228 -237 215 -387 396 -413 398 -243 179 -359 327 -543 361 -325 299 -752 828 -713 747 -290 329 -1248 1091 -221 273 -568 384 -289 326 -957 933 -519 544 -709 670 -298 285 -786 611 -862 765 -373 610 -905 699 -348 334 -502 484 -1535 1578 -574 611 -1295 1864 +493 793 +491 485 +233 217 +194 231 +889 882 +203 226 +780 790 +231 228 +289 334 +310 343 +852 732 +178 284 +208 237 +408 482 +1211 808 +704 718 +1656 1496 +471 490 +1141 1136 +670 733 +549 505 +453 528 +763 944 +193 155 +543 586 +469 412 +166 173 +297 287 +365 349 +413 356 +1656 1283 +137 102 +3102 2255 +780 747 +1423 1315 +708 577 +1675 1545 +638 771 +316 301 +342 534 +401 439 +659 792 +573 563 +432 372 +444 504 +501 448 +218 375 +709 602 +663 489 +762 970 +317 334 +467 486 +315 423 +936 766 +447 450 +625 698 +275 235 +825 554 +695 750 +461 554 +266 364 +194 169 +522 578 +686 826 +339 381 +629 528 +391 451 +1347 1554 +492 345 +109 151 +328 315 +540 517 +420 465 +396 487 +580 869 +1338 1331 +631 1245 +433 440 +345 559 +1441 1847 +554 585 +953 1068 +599 367 +704 776 +537 513 +322 204 +387 292 +460 442 +303 273 +562 905 +704 718 +1752 1786 +434 388 +266 327 +1274 1288 +1248 1325 +758 775 +445 403 +270 276 +703 686 +709 639 +303 272 +310 249 +477 484 +419 514 +275 258 +509 517 +384 319 +286 271 +478 484 +149 144 +1413 1416 +911 812 +649 767 +383 436 +619 498 522 606 -464 573 -448 442 -753 508 -177 170 -468 415 -340 348 -554 550 -877 835 -3451 3442 -311 329 -172 247 -286 457 -516 444 -501 475 -639 612 -1032 1062 -847 1281 -621 679 -378 384 -1650 1273 -535 533 -570 595 -420 342 -719 885 -213 339 -634 686 -691 969 -891 756 -579 642 -589 419 -175 154 -291 328 -330 298 -792 760 -586 537 -565 574 -1996 2014 -537 485 -516 502 -892 680 -607 1019 -681 674 -223 311 -413 388 -969 1015 -570 613 -246 1048 -356 325 -398 359 -327 267 -1010 1103 -401 429 -413 495 -411 437 -333 388 -315 283 -270 261 -706 688 -492 352 -642 659 -597 590 -772 693 -167 145 -1090 898 -377 350 -721 551 -434 432 -727 643 -237 238 -248 239 -1748 1687 -773 697 -415 675 -277 268 -447 417 -348 248 -493 585 -499 875 -545 374 -1261 1543 -586 508 -274 288 -211 174 -1198 786 -649 653 -551 544 -564 442 -577 587 -606 656 -628 526 -921 880 -703 671 -380 368 -530 637 -834 677 -173 197 -1581 1575 -659 756 -555 563 -1413 1055 -532 501 -203 199 -652 831 -247 243 -1269 1314 -1215 1461 -230 227 -789 825 -1515 1443 -478 588 -769 565 -471 754 -181 186 -749 552 -350 452 -337 340 -545 556 -521 541 -736 588 -725 757 -517 857 -529 595 -188 169 -361 403 -383 381 -480 370 -1045 926 -1491 1229 -366 346 -3561 3198 -913 860 -295 202 -487 476 -440 412 -557 430 -530 445 -211 225 -759 738 -1234 1442 -558 542 -203 216 -607 604 -171 151 -201 205 -241 293 -812 1014 -574 419 -557 399 -383 343 -622 658 -806 965 -1013 797 -592 971 -360 447 -574 492 -369 520 -280 292 -1262 1482 -1010 1328 -1154 1128 -213 198 -851 644 -490 641 -756 800 -823 545 -498 451 -541 465 -391 616 -206 243 -466 397 -385 310 -416 427 -892 863 -301 331 -383 389 -382 299 -618 654 -437 369 -504 657 -504 406 -398 463 -425 421 -365 462 -218 221 -207 205 -459 524 -402 478 -511 423 -474 450 -312 318 -271 262 -806 660 -1296 1212 -1109 1128 -644 655 -516 381 -192 185 -618 616 -330 362 -556 518 -895 807 -695 743 -224 194 -372 517 -658 686 -910 1028 -400 761 -561 563 -52 55 -1241 1058 -429 445 -922 1018 -267 293 -634 549 -289 319 -174 207 -577 429 -446 391 -330 415 -349 222 -657 928 -299 215 -407 418 -356 388 -448 224 -399 375 -194 201 -203 194 -775 692 -612 855 -418 499 -501 805 -942 1008 -832 838 -754 842 -344 338 -988 970 -666 565 -294 338 -412 394 -406 520 -1441 1403 -290 256 -155 297 -246 206 -308 259 -248 361 -387 395 -685 442 -229 286 -788 722 -438 610 -1291 1151 -841 935 -1631 1196 -156 381 -318 313 -1075 1202 -691 660 -445 400 -164 178 -695 799 -399 208 -404 298 -1408 1572 -310 242 -975 965 -626 501 -459 415 -549 380 -393 352 -2444 2062 -1121 1288 -520 573 -1282 1118 -554 829 -415 408 -1091 1012 -702 826 -667 548 -918 1060 -419 481 -627 604 -988 1029 -285 326 -145 116 -802 801 -369 278 -2374 2386 -231 384 -521 347 -216 231 -338 279 -623 358 -370 455 -732 842 -386 354 -704 701 -858 953 -196 208 -1303 1252 -347 264 -265 359 -296 229 -389 455 -857 803 -751 768 -483 558 -164 176 -345 345 -798 737 -740 718 -1782 2031 -800 738 -526 453 -442 418 -743 802 -430 500 -520 525 -836 585 -994 949 -518 428 -622 664 -980 1255 -719 726 -598 602 -674 708 -226 195 -302 298 -809 783 -260 307 -768 794 -311 301 -1449 1046 -294 316 -1030 1011 -552 544 -1472 1902 -252 243 -589 341 -590 704 -691 483 -967 1156 -367 404 -1002 862 -502 472 -640 761 -349 405 -567 588 -184 201 -461 441 -347 310 -873 897 -116 100 -385 382 -280 264 -107 94 -450 374 -795 550 -629 645 -441 335 -150 293 -871 874 -240 229 -524 550 -286 300 -311 355 -659 576 -434 430 -285 322 -538 836 -660 628 -212 229 -643 668 -366 415 -267 237 -254 265 -143 147 -822 991 -388 458 -294 314 -346 306 -1089 891 -642 619 -1102 1127 -434 582 -1018 1151 -248 245 -619 929 -905 903 -398 487 -339 367 -371 333 -432 424 -206 191 -724 694 -439 433 -306 313 -309 301 -663 683 -666 706 -286 413 -885 813 -174 165 -361 339 -542 607 -452 428 -650 709 -226 264 -954 1052 -628 708 -367 434 -465 425 -905 787 -535 681 -2069 1826 -1035 999 -1313 1145 +238 394 +529 453 +325 365 +449 490 +470 509 +1119 1072 +355 290 +522 483 +114 104 +803 472 +730 692 +228 298 +538 654 +1074 1071 +506 476 +472 497 +123 175 +711 622 +862 694 +1047 1107 +588 611 +1143 1022 +554 601 +591 609 +1298 1227 +1061 914 +564 538 +683 866 +233 269 +321 358 +639 928 +545 508 +388 350 +861 1129 +1500 1631 +420 433 +501 537 +516 570 +1232 1316 +595 494 +276 401 +1230 1275 +896 988 +634 856 +619 476 +523 565 +119 105 +1181 930 +307 382 +311 545 +142 239 +266 286 +403 583 +141 148 +271 280 +489 387 +480 520 +935 921 +640 535 +400 490 +528 615 +863 549 +448 468 +786 609 +1230 1291 +366 344 +638 531 +1053 1055 +592 616 +620 523 +959 936 +362 442 +455 627 +787 823 +432 400 +474 430 +724 679 +705 549 +382 278 +819 450 +1194 1308 +755 690 +740 1023 +749 739 +922 1062 +604 625 +293 302 +537 540 +1171 1394 +351 533 +179 140 +352 213 +408 328 +588 450 +276 721 +674 721 +383 330 +470 440 +585 371 +947 802 +325 472 +422 421 +681 693 +678 668 +967 1173 +255 330 +232 279 +708 775 +233 238 +583 624 +208 219 +521 516 +429 599 +435 404 +839 823 +294 278 +917 971 +875 1045 +461 514 +400 399 +817 868 +90 96 +526 806 +1124 1178 +742 637 +1450 1641 +577 687 +349 344 +307 353 +252 227 +369 320 +671 693 +503 441 +576 579 +796 1387 +356 401 +223 201 +237 350 +621 587 +760 653 +378 288 +240 234 +981 866 +207 214 +1047 1055 +335 323 +767 731 +313 310 +353 360 +321 327 +379 301 +1209 1147 +1264 1094 +348 354 +678 649 +313 316 +663 810 +305 334 +470 418 +234 274 +298 264 +541 395 +567 477 +782 705 +886 876 +893 949 +388 417 +487 564 +607 594 +903 876 +745 786 +365 422 +389 359 +1044 791 +808 681 +746 845 +863 779 +886 868 +288 179 +1237 1095 +315 466 +222 202 +1129 1249 +500 452 +296 386 +718 826 +383 407 +370 319 +280 245 +800 982 +416 354 +1084 928 +553 577 +948 741 +943 984 +380 372 +727 703 +127 113 +1686 1293 +815 468 +371 500 +798 738 +561 509 +723 959 +822 830 +375 355 +492 493 +401 452 +1039 878 +788 893 +841 537 +698 670 +694 936 +351 399 +1046 1056 +533 547 +357 386 +1236 654 +1003 1009 +481 441 +1922 1938 +193 222 +963 919 +958 693 +1725 1813 +374 327 +761 497 +472 503 +655 899 +723 620 +403 407 +862 861 +913 1041 +432 399 +688 584 +342 310 +423 223 +537 686 +1869 1732 +1521 1512 +959 1403 +507 476 +836 510 +305 301 +1091 832 +497 320 +1258 945 +926 632 +67 90 +549 485 +551 930 +1275 1418 +309 274 +401 461 +557 530 +349 341 +374 356 +1572 1413 +492 508 +1280 1602 +304 360 +764 791 +398 660 +220 216 +518 361 +749 883 +790 811 +116 293 +185 191 +464 557 +1611 1165 +843 891 +465 567 +513 628 +396 575 +448 462 +618 539 +1024 907 +613 677 +341 375 +221 300 +799 799 +244 267 +501 656 +517 559 +294 260 +363 478 +580 653 +624 575 +1605 1432 +927 911 +154 265 +380 396 +421 406 +1022 923 +323 308 +476 299 +298 305 +797 773 +382 361 +234 316 +1414 1460 +226 231 +302 249 +314 370 +504 504 +482 452 +773 749 +762 594 +537 525 +483 508 +209 200 +302 308 +303 290 +706 661 +339 470 +881 640 +460 517 +567 435 +330 356 +186 187 +550 748 +244 256 +582 801 +588 641 +538 414 +198 208 +402 424 +308 276 +596 663 +285 725 +469 691 +328 393 +742 625 +572 368 +563 931 +271 205 +223 280 +298 286 +288 284 +492 435 +443 406 +526 654 +765 825 +420 433 +734 998 +472 481 +859 840 +601 550 +130 126 +104 143 +235 534 +670 888 +851 555 +258 265 +453 397 +1214 1274 +351 360 +274 281 +147 220 +1255 1346 +1130 1335 +567 436 +170 134 +495 523 +399 437 +766 722 +119 139 +796 850 +438 605 +775 794 +386 352 +581 499 +843 917 +268 289 +707 650 +515 705 +822 885 +234 326 +491 842 +231 221 +201 224 +332 325 +1019 902 +440 462 +264 250 +1442 1408 +651 661 +428 944 +847 827 +292 291 +287 309 +429 479 +217 158 +634 573 +409 429 +718 890 +547 444 +522 528 +548 460 +545 537 +260 193 +576 642 +452 431 +922 828 +302 330 +376 814 +849 879 +705 1576 +315 338 +389 393 +415 479 +767 684 +560 669 +179 168 +1528 1458 +404 398 +235 288 +364 441 +420 280 +246 227 +259 338 +1108 976 +1319 1464 +558 522 +290 445 +228 305 +489 388 +389 515 +444 439 +351 414 +364 340 +486 758 +940 954 +339 378 +571 572 +1273 1334 +1273 1052 +673 634 +708 506 +639 671 +464 392 +259 285 +854 869 +362 417 +753 741 +410 363 +759 843 +459 316 +393 401 +446 468 +210 288 316 306 +<<<<<<< HEAD +637 550 +264 301 +601 593 +466 570 +830 782 +672 483 +774 666 +393 367 +296 370 +505 613 +924 795 +711 714 +349 461 +270 271 +1193 1005 +380 380 +142 140 +1354 1129 +540 419 +145 99 +1416 1274 +903 1034 +404 434 +980 1064 +490 488 +1125 925 +375 425 +151 133 +278 241 +635 576 +1511 1827 +399 376 +204 216 +668 686 +755 681 +668 689 +232 226 +553 649 +718 1349 +684 614 +421 436 +714 589 +673 686 +1286 1112 +1174 814 +428 498 +1020 1037 +267 286 +671 665 +614 626 +541 477 +441 461 +254 274 +354 233 +975 848 +520 466 +393 318 +427 407 +358 357 +633 678 +1076 802 +304 209 +418 433 +359 371 +1148 1173 +1765 1780 +368 349 +741 585 +532 386 +748 752 +1318 1295 +421 428 +231 219 +370 367 +522 602 +296 230 +330 342 +239 232 +844 701 +720 831 +408 404 +1119 1241 +404 388 +476 745 +524 577 +362 666 +295 315 +292 289 +1272 1809 +358 366 +901 598 +845 943 +612 466 +459 531 +389 430 +367 368 +300 298 +607 539 +1635 1433 +288 180 +444 740 +319 350 +447 562 +1083 989 +672 835 +171 116 +261 265 +527 539 +748 689 +737 831 +706 851 +1162 983 +366 560 +675 502 +850 877 +352 386 +320 335 +174 157 +3011 2598 +1001 1258 +487 406 +917 924 +299 298 +665 446 +388 391 +522 563 +579 545 +571 562 +207 240 +920 767 +526 574 +268 296 +638 549 +1025 1215 +262 270 +838 815 +261 309 +431 401 +423 404 +158 142 +484 742 +505 561 +260 362 +744 612 +1063 1096 +477 426 +613 553 +701 707 +953 821 +514 518 +301 287 +535 420 +743 594 +917 830 +434 369 +827 807 +577 704 +374 331 +263 294 +403 307 +145 183 +463 477 +433 430 +397 427 +301 269 +828 919 +1112 1311 +722 765 +386 313 +482 507 +255 185 +322 351 +1035 1139 +1430 1402 +697 513 +226 234 +337 387 +406 280 +605 515 +541 731 +725 653 +276 252 +321 280 +1175 1341 +389 413 +549 397 +796 846 +512 544 +418 394 +279 249 +686 654 +907 1044 +735 672 +255 305 +871 724 +338 248 +329 427 +416 362 +880 882 +262 246 +531 496 +248 280 +712 868 +4515 4240 +1561 1297 +367 391 +181 270 +1448 1580 +175 194 +311 491 +1075 1363 +629 1126 +886 1055 +385 481 +1022 1130 +264 196 +793 715 +576 431 +654 563 +398 287 +2408 2098 +1193 1196 +402 375 +708 599 +724 776 +299 355 +164 186 +345 382 +625 526 +251 208 +508 462 +847 731 +172 203 +182 228 +231 267 +674 484 +694 631 +676 820 +312 299 +879 871 +722 751 +619 680 +394 458 +362 538 +1292 2185 +674 698 +480 420 +899 425 +737 652 +1320 1534 +252 219 +1110 1335 +209 265 +535 373 +666 647 +403 322 +246 270 +134 122 +451 402 +737 689 +1026 971 +512 661 +273 265 +987 859 +334 396 +349 345 +965 598 +536 417 +808 722 +839 696 +508 481 +485 558 +238 213 +485 567 +202 204 +1454 1481 +1028 1433 +149 144 +929 906 +1065 1242 +331 293 +413 454 +308 290 +622 754 +614 724 +2281 2320 +177 822 +305 213 +360 306 +168 228 +385 350 +453 498 +396 361 +306 322 +998 882 +462 493 +345 344 +984 1077 +306 375 +291 237 +417 358 +378 345 +377 493 +704 527 +896 898 +721 761 +1123 1246 +982 894 +753 636 +329 311 +442 509 +984 1330 +577 636 +453 439 +652 617 +232 343 +578 434 +663 631 +481 498 +418 391 +512 990 +1023 831 +1490 1322 +903 755 +187 147 +884 1100 +492 487 +310 337 +460 426 +565 1326 +210 216 +393 612 +179 181 +529 423 +487 515 +932 933 +279 311 +312 407 +600 431 +781 662 +811 847 +269 308 +988 939 +663 796 +560 535 +278 242 +277 338 +444 453 +189 177 +236 154 +443 431 +406 374 +2094 2043 +715 917 +147 153 +967 782 +426 411 +1226 1223 +498 538 +369 325 +524 687 +810 672 +348 356 +468 442 +616 835 +864 752 +529 549 +1193 1190 +201 216 +380 434 +342 476 +418 391 +171 174 +777 841 +264 236 +307 250 +1343 1296 +556 525 +347 301 +653 509 +408 290 +488 380 +326 289 +217 190 +917 618 +260 259 +211 214 +1101 1194 +834 460 +737 693 +1041 863 +919 881 +229 184 +194 239 +1676 1258 +285 268 +154 131 +306 323 +203 138 +317 344 +898 1034 +278 337 +477 428 +463 545 +141 153 +1211 869 +1217 1242 +487 427 +373 303 +284 251 +417 428 +424 581 +367 313 +389 383 +425 431 +192 189 +415 364 +312 265 +319 485 +789 780 +300 403 +2006 1968 +430 513 +724 565 +985 747 +1117 1083 +1011 750 +364 379 +340 350 +348 354 +======= 280 265 191 201 222 445 @@ -11002,3 +11176,4 @@ 574 532 571 576 >>>>>>> e615cc169111775374b560ddc35e4c3f964d07b2 +>>>>>>> upstream/master