CARGA DE LIBRERÍAS
##
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
##
## filter, lag
## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, union
CARGA DE DATOS
# Cargar datos
# Cargar datos
ruta_completa <- "C:/Users/klaus/Downloads/5000_Datos.xlsx"
datos <- read_excel(ruta_completa, sheet = "Hoja1")
# Convertir DAYS_LOST a numérico
datos$DAYS_LOST <- as.numeric(datos$DAYS_LOST)## Warning: NAs introduced by coercion
# Generar DAYS_RESTRICT correlacionado con DAYS_LOST
set.seed(123)
datos$DAYS_RESTRICT <- ifelse(
is.na(datos$DAYS_LOST),
NA,
ifelse(
datos$DAYS_LOST == 0,
sample(0:3, nrow(datos), replace = TRUE, prob = c(0.7, 0.2, 0.07, 0.03)),
round(datos$DAYS_LOST * 0.85 + rnorm(nrow(datos), 0, 0.25 * datos$DAYS_LOST))
)
)## Warning in rnorm(nrow(datos), 0, 0.25 * datos$DAYS_LOST): NAs produced
# Definir variables
x <- as.numeric(datos$DAYS_RESTRICT)
y <- as.numeric(datos$DAYS_LOST)
# Tabla de pares de valores
TPV <- data.frame(
DAYS_RESTRICT = x,
DAYS_LOST = y
)
# Eliminar NA, ceros y negativos
TPV <- na.omit(TPV)
TPV <- TPV[TPV$DAYS_RESTRICT >= 0 & TPV$DAYS_LOST >= 0, ]# Eliminación de Outliers
# Outliers Días restringidos
Q1_x <- quantile(TPV$DAYS_RESTRICT, 0.25, na.rm = TRUE)
Q3_x <- quantile(TPV$DAYS_RESTRICT, 0.75, na.rm = TRUE)
IQR_x <- Q3_x - Q1_x
lim_inf_x <- max(0, Q1_x - 1.5 * IQR_x)
lim_sup_x <- Q3_x + 1.5 * IQR_x
# Outliers Días perdidos
Q1_y <- quantile(TPV$DAYS_LOST, 0.25, na.rm = TRUE)
Q3_y <- quantile(TPV$DAYS_LOST, 0.75, na.rm = TRUE)
IQR_y <- Q3_y - Q1_y
lim_inf_y <- max(0, Q1_y - 1.5 * IQR_y)
lim_sup_y <- Q3_y + 1.5 * IQR_y
# Tabla sin outliers
TPV_limpio <- TPV[
TPV$DAYS_RESTRICT >= lim_inf_x & TPV$DAYS_RESTRICT <= lim_sup_x &
TPV$DAYS_LOST >= lim_inf_y & TPV$DAYS_LOST <= lim_sup_y,
]
TPV_limpio## DAYS_RESTRICT DAYS_LOST
## 1 0 0
## 2 1 0
## 3 0 0
## 4 1 0
## 5 2 0
## 7 0 0
## 8 9 10
## 10 0 0
## 11 13 13
## 12 0 0
## 13 0 0
## 14 0 0
## 15 0 0
## 16 1 0
## 17 0 0
## 19 5 5
## 20 2 0
## 21 1 0
## 22 12 12
## 23 0 0
## 25 0 0
## 26 1 0
## 27 0 0
## 29 0 0
## 30 0 0
## 31 2 0
## 32 3 3
## 33 0 0
## 34 1 0
## 35 5 6
## 36 0 0
## 37 10 11
## 38 0 0
## 40 0 0
## 41 0 0
## 42 13 24
## 44 0 0
## 46 0 0
## 47 0 0
## 50 1 0
## 52 0 0
## 53 16 25
## 54 0 0
## 55 0 0
## 56 0 0
## 57 0 0
## 59 11 11
## 63 0 0
## 64 0 0
## 65 4 5
## 66 1 2
## 67 1 0
## 68 1 0
## 71 1 2
## 72 7 11
## 73 0 1
## 75 0 0
## 78 0 0
## 82 12 9
## 88 1 0
## 89 1 0
## 92 5 5
## 93 1 1
## 94 1 1
## 96 0 0
## 97 1 0
## 98 7 6
## 99 0 0
## 100 0 0
## 102 5 14
## 103 0 0
## 104 2 0
## 106 1 0
## 108 1 2
## 109 0 0
## 110 0 0
## 111 2 0
## 112 0 0
## 114 5 4
## 115 1 0
## 116 0 0
## 117 0 0
## 118 2 0
## 119 0 0
## 120 0 0
## 121 0 0
## 122 0 0
## 123 0 0
## 124 0 0
## 126 2 2
## 127 0 0
## 128 0 0
## 129 0 0
## 131 19 13
## 132 1 0
## 133 0 0
## 135 0 0
## 136 0 0
## 137 1 0
## 139 2 2
## 140 0 0
## 141 0 0
## 142 0 0
## 143 0 0
## 144 0 0
## 146 0 0
## 147 0 0
## 148 0 0
## 149 0 0
## 150 4 6
## 151 10 9
## 152 0 0
## 153 0 0
## 154 0 0
## 155 0 0
## 157 0 0
## 159 0 0
## 160 13 16
## 162 0 0
## 163 0 0
## 164 0 0
## 165 0 0
## 166 2 2
## 168 2 4
## 169 0 0
## 170 0 0
## 171 0 0
## 172 0 0
## 173 1 0
## 174 1 0
## 175 0 0
## 176 0 0
## 177 0 0
## 178 1 1
## 179 1 0
## 180 0 0
## 181 1 0
## 182 0 0
## 184 0 0
## 185 0 0
## 186 0 0
## 187 0 0
## 189 2 0
## 190 2 0
## 193 2 2
## 194 0 0
## 195 2 0
## 196 0 0
## 197 0 0
## 198 0 0
## 199 0 0
## 200 0 0
## 201 0 0
## 202 2 0
## 203 0 0
## 204 0 0
## 205 0 0
## 206 3 4
## 207 0 0
## 208 17 21
## 209 0 0
## 211 0 0
## 213 0 0
## 214 2 2
## 215 0 0
## 217 0 0
## 219 1 0
## 220 2 0
## 221 0 0
## 222 2 0
## 223 1 0
## 224 2 2
## 225 0 0
## 226 0 0
## 227 0 0
## 228 0 0
## 229 0 0
## 231 0 0
## 232 18 17
## 234 1 1
## 235 0 0
## 236 0 0
## 237 0 0
## 238 1 1
## 239 2 2
## 240 1 0
## 241 0 0
## 242 0 0
## 243 0 0
## 245 0 0
## 246 1 0
## 247 0 0
## 249 2 0
## 250 1 0
## 251 0 0
## 252 0 0
## 253 0 0
## 254 0 0
## 255 1 1
## 256 1 0
## 257 10 12
## 258 0 0
## 259 0 0
## 260 1 0
## 261 2 0
## 262 1 0
## 264 5 7
## 265 0 0
## 266 0 0
## 267 1 1
## 268 19 26
## 269 0 0
## 270 0 0
## 271 1 0
## 272 0 0
## 273 0 0
## 274 9 11
## 275 6 5
## 277 3 0
## 278 0 0
## 279 0 0
## 281 1 0
## 283 0 0
## 284 0 0
## 285 0 0
## 287 0 0
## 288 5 3
## 289 0 0
## 290 0 0
## 292 0 0
## 294 1 0
## 295 12 16
## 296 1 0
## 297 3 0
## 298 0 0
## 299 2 4
## 300 1 0
## 302 0 0
## 303 1 0
## 305 0 0
## 306 3 6
## 307 0 0
## 308 0 0
## 309 0 0
## 310 0 0
## 311 0 0
## 313 1 0
## 314 11 15
## 315 0 0
## 316 1 0
## 317 1 0
## 318 0 0
## 319 0 0
## 320 1 0
## 321 3 5
## 322 0 0
## 324 2 3
## 325 0 0
## 326 0 0
## 328 0 0
## 332 0 0
## 334 1 0
## 336 0 0
## 337 0 0
## 338 1 1
## 339 0 0
## 340 3 0
## 341 0 0
## 342 0 0
## 343 0 0
## 344 0 0
## 345 0 0
## 346 0 0
## 349 0 0
## 350 0 0
## 351 0 0
## 352 2 0
## 353 0 0
## 354 1 1
## 355 16 23
## 356 3 0
## 357 0 0
## 358 0 0
## 359 0 0
## 361 10 17
## 362 0 0
## 363 1 0
## 364 0 0
## 365 0 0
## 367 0 0
## 368 3 2
## 369 0 0
## 370 0 0
## 372 0 0
## 373 1 0
## 374 0 0
## 375 0 0
## 377 1 0
## 378 0 0
## 379 12 12
## 380 7 10
## 381 0 0
## 382 1 1
## 383 0 0
## 385 1 1
## 387 0 0
## 388 6 8
## 389 1 1
## 392 0 0
## 393 1 0
## 394 1 0
## 396 0 0
## 397 5 7
## 398 0 0
## 401 18 17
## 403 14 17
## 404 0 0
## 405 0 0
## 406 0 0
## 408 0 0
## 410 0 0
## 412 1 0
## 413 0 0
## 414 0 0
## 415 13 19
## 416 0 0
## 417 1 1
## 418 14 23
## 419 0 0
## 420 0 0
## 421 12 15
## 422 0 0
## 423 0 0
## 424 0 0
## 425 1 1
## 426 9 18
## 427 0 0
## 428 0 0
## 429 0 0
## 430 1 0
## 431 2 0
## 432 0 0
## 433 0 0
## 434 7 12
## 435 0 0
## 436 0 0
## 437 0 0
## 438 0 0
## 439 0 0
## 440 5 6
## 441 10 11
## 442 0 0
## 443 1 0
## 444 0 0
## 445 1 0
## 446 1 0
## 448 0 0
## 449 0 0
## 450 0 0
## 451 0 0
## 452 0 0
## 453 0 0
## 454 8 9
## 455 3 3
## 456 1 0
## 458 1 0
## 459 0 0
## 460 0 0
## 461 3 0
## 462 1 2
## 463 0 0
## 464 0 0
## 465 0 0
## 466 0 0
## 467 0 0
## 468 0 0
## 469 21 20
## 470 1 0
## 471 8 9
## 472 1 0
## 473 0 0
## 474 1 0
## 476 0 0
## 478 0 0
## 479 0 0
## 480 1 0
## 481 0 0
## 482 1 0
## 483 0 0
## 484 0 0
## 486 0 0
## 487 0 0
## 488 1 1
## 489 0 0
## 490 1 0
## 491 3 0
## 492 0 0
## 493 0 0
## 494 1 0
## 495 8 11
## 496 3 0
## 497 0 0
## 498 0 0
## 499 1 1
## 500 3 3
## 501 0 0
## 503 0 0
## 505 4 6
## 506 1 1
## 507 0 0
## 508 7 10
## 509 2 0
## 510 0 0
## 511 4 4
## 512 0 0
## 514 9 14
## 517 5 5
## 518 1 0
## 519 2 4
## 520 1 0
## 521 0 0
## 522 0 0
## 523 0 0
## 524 0 0
## 525 0 0
## 526 0 0
## 527 2 0
## 529 22 17
## 530 0 0
## 531 2 0
## 532 0 0
## 533 0 0
## 535 0 0
## 536 1 0
## 537 0 0
## 538 1 1
## 539 0 0
## 541 1 0
## 542 0 0
## 543 1 1
## 544 0 0
## 546 2 0
## 547 2 2
## 548 1 0
## 549 1 0
## 550 0 0
## 551 14 19
## 553 0 0
## 554 8 10
## 555 0 0
## 556 0 0
## 557 0 0
## 558 0 0
## 559 0 0
## 560 0 0
## 561 0 0
## 563 1 0
## 565 0 0
## 566 0 0
## 567 0 0
## 570 1 0
## 571 0 0
## 572 2 0
## 573 0 0
## 575 2 0
## 576 1 0
## 578 0 0
## 579 0 0
## 580 0 0
## 582 1 0
## 583 1 0
## 584 0 0
## 585 1 0
## 586 0 0
## 588 0 0
## 589 3 0
## 592 5 6
## 593 2 0
## 594 0 0
## 595 10 11
## 596 1 0
## 597 0 0
## 598 8 8
## 599 6 6
## 601 3 5
## 603 0 0
## 604 0 0
## 605 5 8
## 606 1 0
## 609 7 14
## 610 4 6
## 611 18 24
## 612 18 17
## 614 8 7
## 616 1 0
## 617 4 4
## 620 0 0
## 621 2 0
## 622 2 2
## 623 1 0
## 624 0 0
## 625 0 0
## 626 0 0
## 627 0 0
## 629 0 0
## 631 1 0
## 632 1 0
## 633 20 20
## 636 0 0
## 637 2 0
## 638 1 0
## 640 0 0
## 641 0 0
## 642 1 0
## 644 9 7
## 645 2 2
## 648 0 0
## 649 7 9
## 650 0 0
## 654 0 0
## 657 2 3
## 658 0 0
## 659 0 0
## 660 10 11
## 662 13 14
## 663 0 0
## 667 22 19
## 668 1 1
## 669 2 3
## 670 0 0
## 671 2 2
## 673 14 14
## 674 0 0
## 676 0 0
## 679 0 0
## 682 8 11
## 685 6 7
## 686 0 0
## 687 0 0
## 688 0 0
## 689 0 0
## 691 0 0
## 693 0 0
## 694 0 0
## 695 0 0
## 697 0 0
## 700 0 0
## 702 0 0
## 703 2 2
## 704 2 6
## 706 0 0
## 707 0 0
## 708 0 0
## 709 0 0
## 710 0 0
## 711 11 15
## 712 0 0
## 713 0 0
## 714 1 0
## 715 0 0
## 716 0 0
## 717 0 0
## 718 0 0
## 719 9 8
## 720 0 0
## 721 0 0
## 722 1 0
## 723 1 0
## 724 1 0
## 725 1 1
## 726 1 0
## 727 1 0
## 728 0 0
## 730 0 0
## 731 0 0
## 732 1 0
## 733 1 2
## 734 3 4
## 735 2 2
## 736 0 0
## 737 1 0
## 739 0 0
## 742 0 0
## 743 0 0
## 746 0 0
## 747 1 0
## 749 4 6
## 750 2 0
## 751 0 0
## 752 0 0
## 753 0 0
## 754 0 0
## 755 0 0
## 756 0 0
## 757 0 0
## 758 1 0
## 759 0 0
## 761 0 0
## 762 0 0
## 763 0 0
## 764 0 0
## 766 0 0
## 767 6 7
## 768 6 6
## 769 1 1
## 770 0 0
## 771 3 4
## 772 7 7
## 773 0 0
## 774 1 0
## 775 15 16
## 780 1 0
## 781 0 0
## 782 1 0
## 783 0 0
## 784 7 8
## 785 1 0
## 786 0 0
## 787 0 0
## 788 0 0
## 789 0 0
## 790 2 0
## 791 2 0
## 792 0 0
## 794 0 0
## 797 0 0
## 798 17 20
## 799 9 13
## 800 0 0
## 801 0 0
## 802 3 3
## 803 0 0
## 804 0 0
## 805 0 0
## 808 0 0
## 812 5 6
## 813 0 0
## 815 1 0
## 816 0 0
## 819 0 0
## 820 0 0
## 821 0 0
## 823 0 0
## 824 15 18
## 825 1 0
## 826 2 2
## 827 0 0
## 828 0 0
## 829 1 1
## 831 0 0
## 833 0 0
## 834 1 0
## 835 0 0
## 836 18 19
## 837 1 0
## 838 0 0
## 839 0 0
## 840 0 0
## 841 2 0
## 843 21 18
## 844 0 0
## 845 0 0
## 846 1 1
## 847 1 0
## 848 1 1
## 849 0 0
## 850 0 0
## 852 5 4
## 853 3 0
## 854 2 3
## 855 7 7
## 857 2 0
## 858 3 4
## 859 0 0
## 860 1 0
## 862 0 0
## 863 0 0
## 864 0 0
## 867 0 0
## 868 0 0
## 871 3 4
## 873 1 0
## 874 0 0
## 875 1 0
## 876 0 0
## 877 0 0
## 878 0 0
## 880 0 0
## 881 0 0
## 882 7 6
## 884 0 0
## 886 2 3
## 889 0 0
## 891 0 0
## 893 0 0
## 894 18 14
## 895 1 0
## 896 1 1
## 897 0 0
## 898 0 0
## 899 0 0
## 900 2 2
## 901 2 0
## 902 0 0
## 903 1 0
## 904 1 1
## 905 0 0
## 906 0 0
## 910 0 0
## 911 0 0
## 913 1 0
## 914 0 0
## 915 0 0
## 916 2 3
## 917 0 0
## 919 0 0
## 920 12 16
## 921 1 0
## 923 0 0
## 924 1 1
## 926 0 0
## 927 0 0
## 929 0 0
## 931 1 0
## 933 2 2
## 935 7 6
## 937 0 0
## 939 20 26
## 941 3 0
## 942 0 0
## 943 0 0
## 945 2 3
## 946 12 15
## 949 0 0
## 951 0 0
## 952 0 0
## 954 0 0
## 955 1 1
## 956 2 0
## 957 0 0
## 959 0 0
## 960 0 0
## 961 0 0
## 962 0 0
## 963 0 0
## 964 0 0
## 965 2 0
## 966 0 0
## 967 0 0
## 968 0 0
## 969 1 0
## 970 0 0
## 971 1 0
## 972 4 6
## 975 12 14
## 977 2 0
## 978 0 0
## 980 1 1
## 982 1 0
## 983 20 20
## 984 0 0
## 985 2 0
## 986 0 1
## 987 0 0
## 989 0 0
## 990 0 0
## 991 0 0
## 992 0 0
## 993 0 0
## 994 0 0
## 995 0 0
## 996 1 0
## 997 5 6
## 999 1 0
## 1000 2 4
## 1001 1 1
## 1002 0 0
## 1003 0 0
## 1004 3 5
## 1006 0 0
## 1009 9 7
## 1010 0 0
## 1011 0 0
## 1012 0 0
## 1013 0 0
## 1014 6 7
## 1015 4 8
## 1016 11 20
## 1017 0 0
## 1019 0 0
## 1020 0 0
## 1021 1 2
## 1022 1 0
## 1024 0 0
## 1025 1 0
## 1026 0 0
## 1027 0 0
## 1028 0 0
## 1029 0 0
## 1030 11 21
## 1031 2 2
## 1032 20 22
## 1033 4 4
## 1035 0 0
## 1037 1 2
## 1038 15 17
## 1040 1 1
## 1041 1 0
## 1042 1 0
## 1043 0 0
## 1044 0 0
## 1045 0 0
## 1046 0 0
## 1047 0 0
## 1048 0 0
## 1049 1 0
## 1051 1 0
## 1052 1 0
## 1053 1 0
## 1054 0 0
## 1055 0 0
## 1056 0 0
## 1057 0 0
## 1059 16 19
## 1061 1 0
## 1062 0 0
## 1063 1 1
## 1064 1 1
## 1065 0 0
## 1066 1 0
## 1067 2 2
## 1068 3 4
## 1069 3 0
## 1070 1 1
## 1071 1 1
## 1073 1 0
## 1074 1 0
## 1075 6 5
## 1076 2 2
## 1077 0 0
## 1078 0 0
## 1079 22 21
## 1080 4 3
## 1085 1 0
## 1087 1 1
## 1088 0 0
## 1089 1 0
## 1090 0 0
## 1092 0 0
## 1094 0 0
## 1095 1 0
## 1096 1 0
## 1100 8 8
## 1101 3 4
## 1103 3 3
## 1104 1 0
## 1106 2 2
## 1107 0 0
## 1108 1 1
## 1109 1 1
## 1110 11 11
## 1111 0 0
## 1112 1 0
## 1113 0 0
## 1115 2 2
## 1116 0 0
## 1117 0 0
## 1118 15 12
## 1120 0 0
## 1122 4 3
## 1125 0 0
## 1128 0 0
## 1129 1 1
## 1131 0 0
## 1132 1 0
## 1135 1 0
## 1136 1 0
## 1137 4 3
## 1139 3 4
## 1140 0 0
## 1141 2 2
## 1143 12 16
## 1144 2 0
## 1147 0 0
## 1151 5 7
## 1152 2 3
## 1153 1 0
## 1155 0 0
## 1157 9 12
## 1159 1 1
## 1160 16 23
## 1161 1 1
## 1163 0 0
## 1164 0 0
## 1165 0 0
## 1166 0 0
## 1167 0 0
## 1168 3 4
## 1171 15 15
## 1172 2 2
## 1173 0 0
## 1174 2 0
## 1175 0 0
## 1176 0 0
## 1177 1 0
## 1178 1 1
## 1179 0 0
## 1180 0 0
## 1181 7 6
## 1182 1 10
## 1183 0 0
## 1184 1 0
## 1186 0 0
## 1187 0 0
## 1188 1 0
## 1189 0 0
## 1190 1 1
## 1191 0 0
## 1192 0 1
## 1194 5 6
## 1195 17 24
## 1196 0 0
## 1197 1 0
## 1198 6 6
## 1199 2 0
## 1200 0 0
## 1201 21 18
## 1202 10 13
## 1203 0 0
## 1207 0 0
## 1208 1 1
## 1209 22 27
## 1210 0 0
## 1211 0 0
## 1212 5 7
## 1215 0 0
## 1216 1 0
## 1217 0 0
## 1218 1 0
## 1219 0 0
## 1221 0 0
## 1223 0 0
## 1224 1 0
## 1225 2 0
## 1226 2 2
## 1227 0 0
## 1228 4 3
## 1229 1 0
## 1230 3 5
## 1231 1 1
## 1232 0 0
## 1233 0 0
## 1234 16 13
## 1235 0 0
## 1236 2 2
## 1237 0 0
## 1238 1 0
## 1239 2 0
## 1240 9 9
## 1241 0 0
## 1242 0 0
## 1243 0 0
## 1244 0 0
## 1245 1 0
## 1246 0 0
## 1247 2 2
## 1248 0 0
## 1250 2 2
## 1251 5 6
## 1252 0 0
## 1253 10 12
## 1254 0 0
## 1255 0 0
## 1256 1 0
## 1257 0 0
## 1258 1 1
## 1259 3 0
## 1260 9 10
## 1261 0 0
## 1262 0 0
## 1263 3 10
## 1264 0 0
## 1265 12 17
## 1268 2 0
## 1270 3 0
## 1273 0 0
## 1275 4 5
## 1277 2 4
## 1278 0 0
## 1279 0 0
## 1281 10 13
## 1282 1 0
## 1283 0 0
## 1284 19 24
## 1285 1 1
## 1286 0 0
## 1287 1 0
## 1288 0 0
## 1289 0 0
## 1290 0 0
## 1292 0 0
## 1293 1 1
## 1294 1 0
## 1295 0 0
## 1296 0 0
## 1297 1 0
## 1298 1 0
## 1299 1 0
## 1300 0 0
## 1301 12 27
## 1302 0 0
## 1303 0 0
## 1304 1 0
## 1305 1 0
## 1306 0 0
## 1307 3 3
## 1308 1 1
## 1309 0 0
## 1310 15 26
## 1311 0 0
## 1312 0 0
## 1313 3 3
## 1314 10 9
## 1315 8 11
## 1316 0 0
## 1318 0 0
## 1319 0 0
## 1320 0 0
## 1321 1 0
## 1323 2 0
## 1324 0 0
## 1325 3 0
## 1326 0 0
## 1327 7 7
## 1328 0 0
## 1329 0 0
## 1330 0 0
## 1331 0 0
## 1332 0 0
## 1333 0 0
## 1334 0 0
## 1336 18 15
## 1337 1 0
## 1338 1 0
## 1341 0 0
## 1342 0 0
## 1343 0 0
## 1345 0 0
## 1347 0 0
## 1348 0 0
## 1349 0 0
## 1350 1 2
## 1351 0 0
## 1352 2 2
## 1354 2 3
## 1355 0 0
## 1356 0 0
## 1358 1 0
## 1359 0 0
## 1360 0 0
## 1361 0 0
## 1362 0 0
## 1363 1 0
## 1364 13 11
## 1365 0 0
## 1366 0 0
## 1367 0 0
## 1368 0 0
## 1369 0 0
## 1372 0 0
## 1373 0 0
## 1374 0 0
## 1375 5 9
## 1376 1 0
## 1377 0 0
## 1378 2 0
## 1379 0 0
## 1381 0 0
## 1382 0 0
## 1383 2 2
## 1384 3 3
## 1386 0 0
## 1387 4 3
## 1388 1 0
## 1390 12 14
## 1391 3 0
## 1392 0 0
## 1394 0 0
## 1395 1 1
## 1396 1 0
## 1397 0 0
## 1398 11 21
## 1399 8 8
## 1400 2 2
## 1402 16 22
## 1403 0 0
## 1404 0 0
## 1405 0 0
## 1407 0 0
## 1408 2 0
## 1411 1 1
## 1412 2 0
## 1413 5 6
## 1414 3 2
## 1415 0 0
## 1416 1 0
## 1417 0 0
## 1418 0 0
## 1419 1 0
## 1420 0 0
## 1421 1 0
## 1423 3 4
## 1424 2 0
## 1425 0 0
## 1426 0 0
## 1427 0 0
## 1429 2 0
## 1431 1 0
## 1432 2 2
## 1433 1 0
## 1434 1 0
## 1435 0 0
## 1436 17 20
## 1437 0 0
## 1438 1 0
## 1440 5 8
## 1441 0 0
## 1442 0 0
## 1444 0 0
## 1445 1 0
## 1447 1 2
## 1448 0 0
## 1449 0 0
## 1450 0 0
## 1451 0 0
## 1452 1 1
## 1453 1 0
## 1454 0 0
## 1455 1 2
## 1456 1 0
## 1457 0 0
## 1458 1 2
## 1459 2 2
## 1460 1 1
## 1461 2 0
## 1464 4 5
## 1465 0 0
## 1466 1 0
## 1467 2 0
## 1468 0 0
## 1470 3 2
## 1471 1 0
## 1472 0 0
## 1473 1 0
## 1474 15 15
## 1475 0 0
## 1476 5 7
## 1477 0 0
## 1478 3 0
## 1479 1 0
## 1480 1 0
## 1481 1 0
## 1482 0 0
## 1483 1 0
## 1484 1 0
## 1485 3 0
## 1486 0 0
## 1487 1 0
## 1488 11 16
## 1492 0 0
## 1493 3 0
## 1494 1 0
## 1495 0 0
## 1496 1 0
## 1497 3 0
## 1498 0 0
## 1499 0 0
## 1502 3 0
## 1503 0 0
## 1504 1 1
## 1505 0 0
## 1507 0 0
## 1508 0 0
## 1509 0 0
## 1510 2 3
## 1511 0 0
## 1512 2 0
## 1513 1 0
## 1516 0 0
## 1517 0 0
## 1518 0 0
## 1520 0 0
## 1521 1 0
## 1522 0 0
## 1523 3 3
## 1524 3 3
## 1525 1 0
## 1526 0 0
## 1527 0 0
## 1528 6 7
## 1529 0 0
## 1530 0 0
## 1532 0 0
## 1534 0 0
## 1535 0 0
## 1536 0 0
## 1537 0 0
## 1538 0 0
## 1539 2 0
## 1541 0 0
## 1542 1 0
## 1543 4 8
## 1544 1 1
## 1545 0 0
## 1546 1 2
## 1547 0 0
## 1548 0 0
## 1549 0 0
## 1550 3 0
## 1551 2 0
## 1552 0 0
## 1553 0 0
## 1554 0 0
## 1555 0 0
## 1557 0 0
## 1558 0 0
## 1559 0 0
## 1560 3 0
## 1562 0 0
## 1563 1 0
## 1565 0 0
## 1566 0 0
## 1568 7 7
## 1570 1 0
## 1572 0 0
## 1573 1 0
## 1574 0 0
## 1575 0 0
## 1576 6 8
## 1577 0 0
## 1578 0 0
## 1579 0 0
## 1580 0 0
## 1582 0 0
## 1583 0 0
## 1584 1 1
## 1585 1 1
## 1587 2 3
## 1588 0 0
## 1589 0 0
## 1590 7 9
## 1591 3 7
## 1592 2 3
## 1593 0 0
## 1596 0 0
## 1597 0 0
## 1598 0 0
## 1599 4 4
## 1600 5 5
## 1601 5 6
## 1602 1 1
## 1603 0 0
## 1604 1 0
## 1606 0 0
## 1607 4 4
## 1608 0 0
## 1609 0 0
## 1610 0 0
## 1611 0 0
## 1612 1 0
## 1613 1 1
## 1614 3 0
## 1615 3 0
## 1616 0 0
## 1617 0 0
## 1618 10 20
## 1619 0 0
## 1620 2 2
## 1621 0 0
## 1622 0 0
## 1623 0 0
## 1624 6 10
## 1627 0 0
## 1628 1 1
## 1629 0 0
## 1631 1 0
## 1632 0 0
## 1633 2 2
## 1634 1 1
## 1635 1 0
## 1636 1 0
## 1637 0 0
## 1640 1 0
## 1642 1 0
## 1643 1 0
## 1645 0 0
## 1646 7 7
## 1647 1 0
## 1649 0 0
## 1650 0 0
## 1651 0 0
## 1652 1 0
## 1653 0 0
## 1654 0 0
## 1655 0 0
## 1659 3 4
## 1661 1 0
## 1662 3 4
## 1664 0 0
## 1665 1 0
## 1668 0 0
## 1669 14 11
## 1670 0 0
## 1671 0 0
## 1674 10 14
## 1676 0 0
## 1677 0 0
## 1678 0 0
## 1679 0 0
## 1680 2 2
## 1681 0 0
## 1682 0 0
## 1683 3 0
## 1684 2 2
## 1687 0 0
## 1690 2 0
## 1691 3 3
## 1692 1 1
## 1694 13 14
## 1695 4 3
## 1696 1 0
## 1698 0 0
## 1699 0 0
## 1700 0 0
## 1701 1 0
## 1702 0 0
## 1703 2 1
## 1705 1 0
## 1706 1 0
## 1707 0 0
## 1708 5 4
## 1709 1 1
## 1710 0 0
## 1711 0 0
## 1712 0 0
## 1714 1 0
## 1717 1 1
## 1718 0 0
## 1720 11 17
## 1721 1 0
## 1722 0 0
## 1723 0 0
## 1724 0 0
## 1725 14 15
## 1726 1 0
## 1728 10 8
## 1730 0 0
## 1731 3 0
## 1732 0 0
## 1733 0 0
## 1735 1 0
## 1739 2 0
## 1740 1 0
## 1741 0 0
## 1742 2 0
## 1743 0 0
## 1745 0 0
## 1746 0 0
## 1747 1 0
## 1749 4 5
## 1750 0 0
## 1752 1 1
## 1753 5 12
## 1754 0 0
## 1755 2 0
## 1756 0 0
## 1758 0 0
## 1760 2 0
## 1761 0 0
## 1762 0 0
## 1763 0 1
## 1764 0 0
## 1765 0 0
## 1766 1 1
## 1768 2 3
## 1769 8 9
## 1770 2 0
## 1772 0 0
## 1773 0 0
## 1774 1 1
## 1775 0 0
## 1776 0 0
## 1777 0 0
## 1779 4 5
## 1780 0 0
## 1781 0 0
## 1782 0 0
## 1783 1 0
## 1784 0 0
## 1785 0 0
## 1787 2 0
## 1788 0 0
## 1791 0 0
## 1792 0 0
## 1793 0 0
## 1794 2 2
## 1796 0 0
## 1797 0 0
## 1798 2 0
## 1799 0 0
## 1800 0 0
## 1801 0 0
## 1802 1 0
## 1803 0 0
## 1804 1 1
## 1806 11 17
## 1807 2 0
## 1808 0 0
## 1810 0 0
## 1813 1 0
## 1815 0 0
## 1816 12 9
## 1818 15 25
## 1819 0 0
## 1821 11 13
## 1822 0 0
## 1825 0 0
## 1827 0 0
## 1830 1 0
## 1837 0 0
## 1842 0 0
## 1843 22 27
## 1844 1 0
## 1846 0 0
## 1849 1 0
## 1851 3 8
## 1852 1 0
## 1853 13 18
## 1854 0 0
## 1858 0 0
## 1861 0 0
## 1863 15 16
## 1866 0 0
## 1868 1 0
## 1870 0 0
## 1871 0 0
## 1872 1 0
## 1873 0 0
## 1874 1 1
## 1877 1 1
## 1879 2 0
## 1880 1 2
## 1883 0 0
## 1887 11 13
## 1889 0 0
## 1890 0 0
## 1891 1 2
## 1892 5 7
## 1893 0 0
## 1894 1 0
## 1896 0 0
## 1897 0 0
## 1898 1 2
## 1906 0 0
## 1908 3 3
## 1909 0 0
## 1910 0 0
## 1911 0 0
## 1912 1 2
## 1913 1 1
## 1914 1 2
## 1915 0 0
## 1916 0 0
## 1918 1 2
## 1919 1 0
## 1922 1 1
## 1923 3 0
## 1924 0 0
## 1925 0 0
## 1926 0 0
## 1927 2 2
## 1928 0 0
## 1930 2 0
## 1931 0 0
## 1932 2 5
## 1933 13 10
## 1935 2 3
## 1936 10 9
## 1940 1 0
## 1941 9 12
## 1942 0 0
## 1943 0 0
## 1944 6 7
## 1948 3 3
## 1950 0 0
## 1952 1 0
## 1954 0 0
## 1955 21 15
## 1957 2 2
## 1958 1 0
## 1959 1 1
## 1960 0 0
## 1961 3 0
## 1962 0 0
## 1963 22 20
## 1965 0 0
## 1966 1 3
## 1967 0 0
## 1969 2 0
## 1970 2 2
## 1971 9 11
## 1972 1 0
## 1974 0 0
## 1975 1 0
## 1977 2 2
## 1978 0 0
## 1980 0 0
## 1981 8 18
## 1983 0 0
## 1984 4 4
## 1986 0 0
## 1988 0 0
## 1989 20 20
## 1990 0 0
## 1991 0 0
## 1992 0 0
## 1993 2 3
## 1994 0 0
## 1995 0 0
## 1996 1 0
## 1997 2 2
## 1998 1 0
## 2000 0 0
## 2001 0 0
## 2005 3 4
## 2008 0 0
## 2009 12 15
## 2010 0 0
## 2011 1 0
## 2012 0 0
## 2013 0 0
## 2014 0 0
## 2015 3 0
## 2016 9 10
## 2018 0 0
## 2019 0 0
## 2020 0 0
## 2022 10 12
## 2023 3 0
## 2024 5 5
## 2026 0 0
## 2027 0 0
## 2028 0 0
## 2029 0 0
## 2030 1 0
## 2031 0 0
## 2032 2 5
## 2033 1 1
## 2034 0 0
## 2035 0 0
## 2036 1 2
## 2037 0 0
## 2039 2 0
## 2040 1 0
## 2041 0 0
## 2042 2 2
## 2043 0 0
## 2044 1 0
## 2045 2 0
## 2047 1 0
## 2048 0 0
## 2050 0 0
## 2051 0 0
## 2052 0 0
## 2053 10 10
## 2054 0 0
## 2055 2 0
## 2056 2 2
## 2057 5 6
## 2061 0 0
## 2062 1 0
## 2063 0 0
## 2064 9 8
## 2065 0 0
## 2066 0 0
## 2067 0 0
## 2069 1 0
## 2070 2 0
## 2071 0 0
## 2072 2 0
## 2073 0 0
## 2075 0 0
## 2076 1 0
## 2077 0 0
## 2078 0 0
## 2079 2 0
## 2080 0 0
## 2081 0 0
## 2082 1 1
## 2083 1 0
## 2084 0 0
## 2085 0 0
## 2086 0 0
## 2087 0 0
## 2088 3 0
## 2089 0 0
## 2090 0 0
## 2091 0 0
## 2092 0 0
## 2093 0 0
## 2094 1 1
## 2095 0 0
## 2098 0 0
## 2100 0 0
## 2101 0 0
## 2104 2 0
## 2106 1 0
## 2107 0 0
## 2108 0 0
## 2109 1 0
## 2110 1 0
## 2111 1 0
## 2112 0 0
## 2113 0 0
## 2114 0 0
## 2116 0 0
## 2117 1 0
## 2119 1 1
## 2121 8 17
## 2122 0 0
## 2124 0 0
## 2126 1 1
## 2127 1 0
## 2128 0 0
## 2129 0 0
## 2130 1 1
## 2131 1 1
## 2132 0 0
## 2133 0 0
## 2134 5 4
## 2135 0 0
## 2136 0 0
## 2137 0 0
## 2138 1 0
## 2139 2 0
## 2140 0 0
## 2141 1 0
## 2142 0 0
## 2143 0 0
## 2144 0 0
## 2145 0 0
## 2146 0 0
## 2147 0 0
## 2148 2 3
## 2149 3 3
## 2150 1 0
## 2151 1 1
## 2152 2 0
## 2153 0 0
## 2155 0 0
## 2156 0 0
## 2157 0 1
## 2158 1 0
## 2159 1 1
## 2161 1 1
## 2162 0 0
## 2163 0 0
## 2164 1 0
## 2165 7 8
## 2166 7 8
## 2167 0 0
## 2169 1 0
## 2171 0 0
## 2173 0 0
## 2174 0 0
## 2175 4 4
## 2176 0 0
## 2177 0 0
## 2178 8 14
## 2179 1 2
## 2180 1 0
## 2181 2 0
## 2184 0 0
## 2186 0 0
## 2187 1 0
## 2188 1 0
## 2190 1 0
## 2191 0 0
## 2192 3 2
## 2193 0 0
## 2195 5 4
## 2196 1 0
## 2197 0 0
## 2199 0 0
## 2200 0 0
## 2201 1 0
## 2202 1 0
## 2204 5 6
## 2205 12 16
## 2206 0 0
## 2207 0 0
## 2210 0 0
## 2212 0 0
## 2214 0 0
## 2215 1 1
## 2216 12 18
## 2217 10 13
## 2218 0 0
## 2219 0 0
## 2220 0 0
## 2221 0 0
## 2222 0 0
## 2223 1 0
## 2224 0 0
## 2225 9 10
## 2226 0 0
## 2228 1 1
## 2229 15 13
## 2230 0 0
## 2231 1 0
## 2232 1 0
## 2233 0 0
## 2234 0 0
## 2235 0 0
## 2236 0 0
## 2237 2 0
## 2238 0 0
## 2239 0 0
## 2240 2 0
## 2241 4 3
## 2242 0 0
## 2243 1 1
## 2244 2 0
## 2245 2 0
## 2246 1 0
## 2248 1 0
## 2249 0 0
## 2250 1 2
## 2251 3 0
## 2252 3 0
## 2254 18 17
## 2255 0 0
## 2256 0 0
## 2257 0 0
## 2258 0 0
## 2259 0 0
## 2260 8 6
## 2261 13 12
## 2262 1 0
## 2263 6 7
## 2264 0 0
## 2265 0 0
## 2266 0 0
## 2269 0 0
## 2270 0 0
## 2271 20 18
## 2274 0 0
## 2275 2 0
## 2277 2 0
## 2278 0 0
## 2279 0 0
## 2280 0 0
## 2282 3 3
## 2283 0 0
## 2284 0 0
## 2285 0 0
## 2286 0 0
## 2287 0 0
## 2289 0 0
## 2290 0 0
## 2291 14 22
## 2293 14 17
## 2294 2 2
## 2295 1 0
## 2296 7 7
## 2298 1 0
## 2299 0 0
## 2300 1 0
## 2301 0 0
## 2304 0 0
## 2305 11 12
## 2306 0 0
## 2307 0 0
## 2308 0 0
## 2309 2 0
## 2311 0 0
## 2313 4 7
## 2314 0 0
## 2315 2 3
## 2316 14 14
## 2318 13 16
## 2319 0 0
## 2320 0 0
## 2322 0 0
## 2323 1 1
## 2324 0 0
## 2325 0 0
## 2326 13 20
## 2327 0 0
## 2328 0 0
## 2329 0 0
## 2330 0 0
## 2331 0 0
## 2332 0 0
## 2333 0 0
## 2334 1 0
## 2335 1 0
## 2337 0 0
## 2338 2 0
## 2339 2 0
## 2340 2 0
## 2342 19 22
## 2343 2 0
## 2344 1 1
## 2346 5 5
## 2347 0 0
## 2348 0 0
## 2349 0 0
## 2350 1 1
## 2351 1 0
## 2352 5 5
## 2353 0 0
## 2354 1 0
## 2355 0 0
## 2356 13 17
## 2358 0 0
## 2359 5 5
## 2361 2 3
## 2363 0 0
## 2364 1 0
## 2365 0 0
## 2366 0 0
## 2367 0 0
## 2369 1 0
## 2370 0 0
## 2371 1 0
## 2375 0 0
## 2376 0 0
## 2377 2 0
## 2378 0 0
## 2380 2 0
## 2381 1 1
## 2382 0 0
## 2383 3 0
## 2384 0 0
## 2386 0 0
## 2388 1 0
## 2389 1 1
## 2390 1 0
## 2392 0 0
## 2393 0 0
## 2394 0 0
## 2395 2 0
## 2396 3 0
## 2397 6 9
## 2398 4 7
## 2399 1 2
## 2400 2 0
## 2401 1 0
## 2402 9 10
## 2403 0 0
## 2404 2 0
## 2406 0 0
## 2409 1 0
## 2410 8 15
## 2412 0 0
## 2413 11 12
## 2414 0 0
## 2416 5 7
## 2422 1 0
## 2423 11 10
## 2424 0 0
## 2425 7 6
## 2427 15 23
## 2428 1 0
## 2429 0 0
## 2431 1 0
## 2436 0 1
## 2437 5 5
## 2439 0 0
## 2441 3 6
## 2444 12 11
## 2445 17 18
## 2449 16 19
## 2450 4 9
## 2452 9 10
## 2453 0 0
## 2454 14 19
## 2455 1 0
## 2457 0 0
## 2460 20 25
## 2464 12 14
## 2466 14 11
## 2467 20 27
## 2470 2 11
## 2479 1 1
## 2481 13 15
## 2485 3 3
## 2486 8 10
## 2487 2 7
## 2488 0 0
## 2490 1 0
## 2491 1 0
## 2494 2 0
## 2495 1 0
## 2496 1 0
## 2497 2 3
## 2499 1 1
## 2500 0 0
## 2501 8 7
## 2502 0 0
## 2503 2 2
## 2505 0 0
## 2507 0 0
## 2508 1 0
## 2510 0 0
## 2511 0 0
## 2512 0 0
## 2514 0 0
## 2517 0 0
## 2518 10 17
## 2519 5 4
## 2523 1 0
## 2524 0 0
## 2525 0 0
## 2526 6 6
## 2528 0 0
## 2530 2 0
## 2531 1 0
## 2533 10 11
## 2534 0 0
## 2535 2 0
## 2536 0 0
## 2541 1 0
## 2542 0 0
## 2543 1 0
## 2547 0 0
## 2550 1 1
## 2551 0 0
## 2552 3 7
## 2553 5 5
## 2555 1 0
## 2556 0 0
## 2558 0 0
## 2560 0 0
## 2562 1 0
## 2563 0 0
## 2565 1 2
## 2566 1 18
## 2567 0 0
## 2568 0 0
## 2569 11 15
## 2570 1 0
## 2571 0 0
## 2577 0 0
## 2578 0 0
## 2580 1 0
## 2582 0 0
## 2583 3 8
## 2584 0 0
## 2585 2 5
## 2587 2 0
## 2588 1 0
## 2590 5 6
## 2593 0 0
## 2595 0 0
## 2596 1 1
## 2599 1 0
## 2600 0 0
## 2603 0 0
## 2604 16 13
## 2605 0 0
## 2606 2 2
## 2613 3 3
## 2614 0 0
## 2616 3 4
## 2617 0 0
## 2618 0 0
## 2619 0 0
## 2620 2 0
## 2621 14 22
## 2622 2 2
## 2623 0 0
## 2624 7 8
## 2625 1 0
## 2627 0 0
## 2628 2 2
## 2629 0 0
## 2630 4 6
## 2631 3 3
## 2632 1 0
## 2633 0 0
## 2634 3 3
## 2636 0 0
## 2637 20 27
## 2638 0 0
## 2640 0 0
## 2641 0 0
## 2642 4 3
## 2643 10 10
## 2644 0 0
## 2645 7 19
## 2646 1 0
## 2647 3 0
## 2649 0 0
## 2650 0 0
## 2651 1 1
## 2653 4 5
## 2654 0 0
## 2656 1 1
## 2657 0 0
## 2658 1 0
## 2659 2 2
## 2661 0 0
## 2663 1 0
## 2665 0 0
## 2666 0 0
## 2667 3 3
## 2668 0 0
## 2669 7 8
## 2670 0 0
## 2671 3 3
## 2672 6 6
## 2673 1 0
## 2675 0 0
## 2676 0 0
## 2677 0 0
## 2678 2 4
## 2679 0 0
## 2680 0 0
## 2681 3 4
## 2682 0 0
## 2683 1 2
## 2686 1 0
## 2689 13 11
## 2690 3 0
## 2691 0 0
## 2692 0 0
## 2693 0 0
## 2694 0 0
## 2695 1 1
## 2696 3 0
## 2697 0 0
## 2699 0 0
## 2700 3 0
## 2701 0 0
## 2705 0 0
## 2706 2 0
## 2707 1 0
## 2708 1 0
## 2709 1 0
## 2710 10 17
## 2711 3 0
## 2712 0 0
## 2715 0 0
## 2717 0 0
## 2718 0 0
## 2719 3 4
## 2720 0 0
## 2721 3 0
## 2722 1 1
## 2723 0 0
## 2727 1 0
## 2729 0 0
## 2730 0 0
## 2731 1 1
## 2732 1 2
## 2734 0 0
## 2735 1 0
## 2736 0 0
## 2737 1 1
## 2738 1 1
## 2739 0 0
## 2742 0 0
## 2743 4 4
## 2744 0 0
## 2745 2 0
## 2746 1 0
## 2747 1 2
## 2750 0 0
## 2753 0 0
## 2754 12 21
## 2755 0 0
## 2760 0 0
## 2762 0 0
## 2764 6 7
## 2765 17 16
## 2766 4 3
## 2767 4 5
## 2768 4 6
## 2770 2 2
## 2771 0 0
## 2772 3 3
## 2773 0 0
## 2774 0 0
## 2775 0 0
## 2781 3 6
## 2783 0 0
## 2784 1 1
## 2786 1 2
## 2787 10 9
## 2788 5 8
## 2789 1 0
## 2790 2 0
## 2791 5 8
## 2792 0 0
## 2793 0 0
## 2794 0 0
## 2795 1 2
## 2796 0 0
## 2799 0 0
## 2800 1 0
## 2801 0 0
## 2802 0 0
## 2803 0 0
## 2804 11 14
## 2805 2 2
## 2808 18 26
## 2809 16 14
## 2810 1 1
## 2813 1 0
## 2814 1 0
## 2816 0 0
## 2817 0 0
## 2818 11 11
## 2819 0 0
## 2820 1 0
## 2821 0 0
## 2822 0 0
## 2826 0 0
## 2827 1 0
## 2828 0 0
## 2830 0 0
## 2831 0 0
## 2832 2 2
## 2836 1 1
## 2837 1 0
## 2838 0 0
## 2839 1 0
## 2840 5 5
## 2842 0 0
## 2844 1 0
## 2845 0 0
## 2847 1 0
## 2848 0 0
## 2849 0 0
## 2850 0 0
## 2851 2 0
## 2852 0 0
## 2853 0 0
## 2854 2 3
## 2855 1 0
## 2856 0 0
## 2857 6 5
## 2858 2 0
## 2859 6 7
## 2860 1 1
## 2861 1 0
## 2863 15 17
## 2865 0 0
## 2866 9 8
## 2871 0 0
## 2873 0 0
## 2876 1 0
## 2877 0 0
## 2878 0 0
## 2879 0 0
## 2880 2 0
## 2881 0 0
## 2882 0 0
## 2883 4 4
## 2884 0 0
## 2885 0 0
## 2886 1 1
## 2887 1 0
## 2888 8 14
## 2889 0 0
## 2890 22 25
## 2891 0 0
## 2892 2 3
## 2894 1 1
## 2895 1 3
## 2898 0 0
## 2900 0 0
## 2901 2 0
## 2902 5 5
## 2903 1 0
## 2904 0 0
## 2905 1 0
## 2906 0 0
## 2907 2 0
## 2908 0 0
## 2909 0 0
## 2910 1 0
## 2911 0 0
## 2913 1 0
## 2914 0 0
## 2915 2 0
## 2916 1 1
## 2917 1 0
## 2918 1 0
## 2919 0 0
## 2920 1 0
## 2921 0 0
## 2922 6 5
## 2923 1 0
## 2924 15 20
## 2925 12 13
## 2926 0 0
## 2927 1 0
## 2929 11 22
## 2932 0 0
## 2933 0 0
## 2934 0 0
## 2936 0 0
## 2938 18 21
## 2939 5 5
## 2940 4 6
## 2942 1 0
## 2943 0 0
## 2944 0 0
## 2945 5 9
## 2946 4 4
## 2947 0 0
## 2950 0 0
## 2951 0 0
## 2953 6 6
## 2954 1 1
## 2955 0 0
## 2956 1 0
## 2957 0 0
## 2958 1 0
## 2959 0 0
## 2960 0 0
## 2961 2 0
## 2962 0 0
## 2963 1 0
## 2965 0 0
## 2966 0 0
## 2967 1 0
## 2968 0 0
## 2972 0 0
## 2975 0 0
## 2976 6 5
## 2980 0 0
## 2982 0 0
## 2983 0 0
## 2988 0 0
## 2989 0 0
## 2991 0 0
## 2992 0 0
## 2994 3 3
## 2995 1 0
## 2996 1 2
## 2997 1 0
## 2998 2 2
## 3001 0 0
## 3002 2 0
## 3003 0 0
## 3005 0 0
## 3010 1 0
## 3011 3 0
## 3012 0 0
## 3013 0 0
## 3014 0 0
## 3015 0 0
## 3018 7 5
## 3020 2 3
## 3021 0 0
## 3022 5 5
## 3026 1 0
## 3027 9 8
## 3028 2 0
## 3029 0 0
## 3030 2 0
## 3031 1 0
## 3032 18 17
## 3033 0 0
## 3034 1 1
## 3035 1 0
## 3036 0 0
## 3037 0 0
## 3038 0 0
## 3039 0 0
## 3040 1 0
## 3041 0 0
## 3042 0 0
## 3043 0 0
## 3044 0 0
## 3045 0 0
## 3046 1 0
## 3047 0 0
## 3049 0 0
## 3051 0 0
## 3055 0 0
## 3056 1 0
## 3058 12 16
## 3061 2 2
## 3062 4 5
## 3063 12 13
## 3067 0 0
## 3068 2 4
## 3072 1 0
## 3073 8 8
## 3074 1 1
## 3075 9 9
## 3076 1 0
## 3077 15 15
## 3078 4 4
## 3079 15 13
## 3080 4 4
## 3081 1 0
## 3083 0 0
## 3084 6 7
## 3087 10 21
## 3088 18 23
## 3089 0 0
## 3090 14 14
## 3091 15 24
## 3092 3 4
## 3093 1 1
## 3094 0 0
## 3095 1 0
## 3096 2 4
## 3098 0 0
## 3103 4 4
## 3105 1 1
## 3106 0 0
## 3107 1 0
## 3108 1 1
## 3109 1 0
## 3110 0 0
## 3113 0 0
## 3115 0 0
## 3116 3 0
## 3117 0 0
## 3118 2 0
## 3119 21 23
## 3120 0 0
## 3121 0 0
## 3122 0 0
## 3123 0 2
## 3125 0 0
## 3126 1 0
## 3127 1 0
## 3130 2 0
## 3131 0 0
## 3132 0 0
## 3133 1 0
## 3134 1 1
## 3136 3 2
## 3137 0 0
## 3138 0 0
## 3139 0 0
## 3140 1 1
## 3141 2 0
## 3143 0 0
## 3144 0 0
## 3145 7 5
## 3146 0 0
## 3147 0 0
## 3148 2 0
## 3149 0 0
## 3150 0 0
## 3151 1 0
## 3153 0 0
## 3154 4 6
## 3155 0 0
## 3156 0 0
## 3157 0 0
## 3158 21 20
## 3160 0 0
## 3161 0 0
## 3162 0 0
## 3163 1 0
## 3164 0 0
## 3166 1 0
## 3167 0 0
## 3168 0 0
## 3169 0 0
## 3170 0 0
## 3171 2 2
## 3173 0 0
## 3174 16 14
## 3176 0 0
## 3177 1 1
## 3178 1 1
## 3179 1 0
## 3180 1 0
## 3181 0 0
## 3182 1 0
## 3183 0 0
## 3184 0 1
## 3186 10 8
## 3188 0 0
## 3191 0 0
## 3192 0 0
## 3193 0 0
## 3194 2 0
## 3195 14 17
## 3196 3 4
## 3197 0 0
## 3199 1 0
## 3201 0 0
## 3202 0 0
## 3203 1 0
## 3205 0 0
## 3206 0 0
## 3207 0 0
## 3209 1 1
## 3210 0 0
## 3211 3 4
## 3212 0 0
## 3213 3 0
## 3214 0 0
## 3215 1 0
## 3216 0 0
## 3217 2 4
## 3219 2 4
## 3220 7 8
## 3221 1 0
## 3222 0 0
## 3223 1 0
## 3224 0 0
## 3225 0 0
## 3226 1 0
## 3227 0 0
## 3228 0 0
## 3229 0 0
## 3230 0 0
## 3231 0 0
## 3232 1 1
## 3233 0 0
## 3234 0 0
## 3235 1 0
## 3236 8 10
## 3237 1 0
## 3239 0 0
## 3241 0 0
## 3242 0 0
## 3243 3 0
## 3244 2 2
## 3245 1 0
## 3246 1 0
## 3247 2 0
## 3248 0 0
## 3249 0 0
## 3250 1 0
## 3251 1 1
## 3252 0 0
## 3254 0 0
## 3255 2 0
## 3256 0 0
## 3257 0 0
## 3258 1 0
## 3259 1 0
## 3260 2 2
## 3261 0 0
## 3262 0 0
## 3263 17 16
## 3264 0 0
## 3265 0 0
## 3266 2 0
## 3268 0 0
## 3269 3 3
## 3270 0 0
## 3271 1 2
## 3272 0 0
## 3273 2 0
## 3274 7 6
## 3276 1 0
## 3277 0 0
## 3280 0 0
## 3282 0 0
## 3283 0 0
## 3284 0 0
## 3285 2 0
## 3286 0 0
## 3287 2 0
## 3288 5 6
## 3289 1 0
## 3291 0 0
## 3292 0 0
## 3293 0 0
## 3295 0 0
## 3296 0 0
## 3299 0 0
## 3300 2 0
## 3301 1 1
## 3302 0 0
## 3303 0 0
## 3305 2 0
## 3306 0 0
## 3307 1 0
## 3308 0 0
## 3309 0 0
## 3310 0 0
## 3312 1 0
## 3314 1 0
## 3316 0 0
## 3317 1 0
## 3318 0 0
## 3319 0 0
## 3320 0 0
## 3322 5 8
## 3324 6 7
## 3327 3 4
## 3328 1 0
## 3329 0 0
## 3330 0 0
## 3331 1 2
## 3332 0 0
## 3334 0 0
## 3335 1 0
## 3337 0 0
## 3338 0 0
## 3340 2 2
## 3342 2 2
## 3343 1 0
## 3344 0 0
## 3345 0 0
## 3346 1 0
## 3347 0 0
## 3351 2 0
## 3352 0 0
## 3353 1 0
## 3355 0 0
## 3356 0 0
## 3357 0 0
## 3358 0 0
## 3360 0 0
## 3361 0 0
## 3362 0 0
## 3364 1 0
## 3365 0 0
## 3366 9 10
## 3367 10 13
## 3368 7 7
## 3371 1 0
## 3373 0 0
## 3374 2 0
## 3375 0 0
## 3377 0 0
## 3379 0 0
## 3380 0 0
## 3382 0 0
## 3383 0 0
## 3384 5 4
## 3385 0 0
## 3386 0 0
## 3387 2 0
## 3388 0 0
## 3389 7 9
## 3390 0 0
## 3391 0 0
## 3392 0 0
## 3393 0 0
## 3394 0 0
## 3395 0 0
## 3396 2 0
## 3397 0 0
## 3398 0 0
## 3400 0 0
## 3401 1 1
## 3402 2 2
## 3403 0 0
## 3406 0 0
## 3407 2 0
## 3408 0 0
## 3409 2 0
## 3410 0 0
## 3411 0 0
## 3412 0 0
## 3413 0 0
## 3414 1 0
## 3417 1 0
## 3418 0 0
## 3420 0 0
## 3421 0 0
## 3422 6 10
## 3423 3 0
## 3424 2 3
## 3425 5 6
## 3426 0 0
## 3427 1 0
## 3428 0 0
## 3429 0 0
## 3430 6 6
## 3432 0 0
## 3433 0 0
## 3434 8 8
## 3435 0 0
## 3436 0 0
## 3437 2 0
## 3438 3 0
## 3439 0 0
## 3440 0 0
## 3441 0 0
## 3442 1 0
## 3443 0 0
## 3444 0 0
## 3445 0 0
## 3446 0 0
## 3447 0 0
## 3448 0 0
## 3449 0 0
## 3450 0 0
## 3451 0 0
## 3452 0 0
## 3453 0 0
## 3454 0 0
## 3455 0 0
## 3456 0 0
## 3457 0 0
## 3458 0 0
## 3459 0 0
## 3460 13 14
## 3461 1 1
## 3462 0 0
## 3463 0 0
## 3464 0 0
## 3465 0 0
## 3466 12 20
## 3467 1 0
## 3468 0 0
## 3469 0 0
## 3470 0 0
## 3471 0 0
## 3472 0 0
## 3473 1 0
## 3474 1 0
## 3475 0 0
## 3476 22 21
## 3477 0 0
## 3478 0 0
## 3479 0 0
## 3480 2 0
## 3481 0 0
## 3482 3 4
## 3483 0 0
## 3484 15 18
## 3485 19 27
## 3486 0 0
## 3487 1 1
## 3488 1 0
## 3489 1 0
## 3490 1 0
## 3491 0 0
## 3492 0 0
## 3493 0 0
## 3498 0 0
## 3499 0 0
## 3500 0 0
## 3501 6 7
## 3502 0 0
## 3503 1 0
## 3506 0 0
## 3508 0 0
## 3509 0 0
## 3513 0 0
## 3514 0 0
## 3515 2 3
## 3516 0 0
## 3517 0 0
## 3518 0 1
## 3519 3 4
## 3520 2 3
## 3521 0 0
## 3522 0 0
## 3524 0 0
## 3525 1 1
## 3526 1 1
## 3527 0 0
## 3529 2 3
## 3530 5 6
## 3531 0 0
## 3532 1 0
## 3533 0 0
## 3534 0 0
## 3535 1 0
## 3537 8 8
## 3538 9 8
## 3539 0 0
## 3541 16 15
## 3542 2 0
## 3543 0 0
## 3544 0 0
## 3546 3 0
## 3547 0 0
## 3548 2 3
## 3549 2 2
## 3552 0 0
## 3553 12 17
## 3554 10 26
## 3555 5 6
## 3556 1 0
## 3557 0 0
## 3558 0 0
## 3560 0 0
## 3561 0 0
## 3562 0 0
## 3564 10 14
## 3565 0 0
## 3566 4 8
## 3567 0 0
## 3568 3 0
## 3569 0 0
## 3572 0 0
## 3576 1 0
## 3577 0 0
## 3579 0 0
## 3581 1 0
## 3583 7 9
## 3584 13 18
## 3585 2 0
## 3587 0 0
## 3588 1 0
## 3589 0 0
## 3591 0 0
## 3592 0 0
## 3593 0 0
## 3596 2 0
## 3597 4 8
## 3598 2 0
## 3599 0 0
## 3604 2 2
## 3605 0 0
## 3606 11 15
## 3607 1 0
## 3608 0 0
## 3609 2 0
## 3611 0 0
## 3612 0 0
## 3613 0 0
## 3615 19 27
## 3616 0 0
## 3617 0 0
## 3618 0 0
## 3619 0 0
## 3620 1 0
## 3623 0 0
## 3624 0 0
## 3628 0 0
## 3629 0 0
## 3631 0 0
## 3632 0 0
## 3635 1 0
## 3636 0 0
## 3638 0 0
## 3640 0 0
## 3641 0 0
## 3644 1 2
## 3647 0 0
## 3648 1 1
## 3649 0 0
## 3650 1 0
## 3654 1 0
## 3655 2 0
## 3656 0 0
## 3657 0 0
## 3658 0 0
## 3660 7 7
## 3662 1 0
## 3664 2 0
## 3665 10 17
## 3666 0 0
## 3668 1 1
## 3669 0 0
## 3671 1 0
## 3672 1 1
## 3673 0 0
## 3674 2 0
## 3675 2 0
## 3676 9 8
## 3677 2 0
## 3678 0 0
## 3679 0 0
## 3680 3 4
## 3681 0 0
## 3682 2 0
## 3684 1 0
## 3685 0 0
## 3686 0 0
## 3687 1 0
## 3689 1 0
## 3690 0 0
## 3691 0 0
## 3692 0 0
## 3694 0 0
## 3695 0 0
## 3696 0 0
## 3697 1 1
## 3699 0 0
## 3700 0 0
## 3701 1 2
## 3702 1 1
## 3703 2 0
## 3704 0 0
## 3705 3 0
## 3706 1 0
## 3707 0 0
## 3708 0 0
## 3709 1 0
## 3710 0 0
## 3711 5 7
## 3714 0 0
## 3715 1 0
## 3716 0 0
## 3717 0 0
## 3718 2 0
## 3719 0 0
## 3720 0 0
## 3721 1 0
## 3722 0 0
## 3723 0 0
## 3724 2 0
## 3725 5 6
## 3726 0 0
## 3727 0 0
## 3728 0 0
## 3730 8 6
## 3731 5 5
## 3732 0 0
## 3733 0 0
## 3734 2 3
## 3735 0 0
## 3736 2 0
## 3737 2 0
## 3738 0 0
## 3739 4 4
## 3740 0 0
## 3742 19 17
## 3743 2 0
## 3745 3 5
## 3746 2 2
## 3748 15 17
## 3749 1 0
## 3750 0 0
## 3752 0 0
## 3754 0 0
## 3755 1 1
## 3756 2 3
## 3757 4 10
## 3758 0 0
## 3759 1 0
## 3760 0 0
## 3762 2 0
## 3763 0 0
## 3764 2 0
## 3765 1 0
## 3766 1 0
## 3768 3 0
## 3769 0 0
## 3770 4 4
## 3771 1 0
## 3772 2 0
## 3774 1 1
## 3775 17 18
## 3776 1 0
## 3777 1 0
## 3778 6 7
## 3779 0 0
## 3780 1 0
## 3781 1 1
## 3782 6 9
## 3783 2 0
## 3784 5 7
## 3785 4 3
## 3786 0 0
## 3787 0 0
## 3788 0 0
## 3789 1 0
## 3791 0 0
## 3793 9 10
## 3794 2 0
## 3795 4 5
## 3797 0 0
## 3798 1 0
## 3799 1 0
## 3800 0 0
## 3803 0 0
## 3804 0 0
## 3805 0 0
## 3806 1 0
## 3807 2 0
## 3808 0 0
## 3809 0 0
## 3810 0 0
## 3812 1 0
## 3813 1 2
## 3817 20 24
## 3819 9 16
## 3820 0 0
## 3821 2 3
## 3822 2 3
## 3823 0 0
## 3825 1 0
## 3827 0 0
## 3828 0 0
## 3829 0 0
## 3830 10 9
## 3832 2 3
## 3833 0 0
## 3836 3 0
## 3839 1 0
## 3840 0 0
## 3842 0 0
## 3845 0 0
## 3846 21 22
## 3848 5 5
## 3850 2 0
## 3851 1 1
## 3853 0 0
## 3856 19 15
## 3857 8 19
## 3859 0 0
## 3862 0 0
## 3865 1 0
## 3866 0 0
## 3867 0 0
## 3869 1 0
## 3872 4 3
## 3873 0 0
## 3874 0 0
## 3875 1 0
## 3876 0 0
## 3877 0 0
## 3878 7 11
## 3879 0 0
## 3880 12 15
## 3881 2 0
## 3882 6 8
## 3883 0 0
## 3886 14 20
## 3887 2 4
## 3889 0 0
## 3890 4 3
## 3891 1 0
## 3892 0 0
## 3893 7 18
## 3894 0 0
## 3895 19 25
## 3896 0 0
## 3897 0 0
## 3898 0 0
## 3899 3 0
## 3900 2 2
## 3901 0 0
## 3902 0 0
## 3903 0 0
## 3905 0 0
## 3906 0 0
## 3907 0 0
## 3908 1 0
## 3909 2 0
## 3910 1 0
## 3911 0 0
## 3912 1 1
## 3914 9 8
## 3915 5 7
## 3916 0 0
## 3917 0 0
## 3918 0 0
## 3919 2 2
## 3922 0 0
## 3926 0 0
## 3927 2 0
## 3929 0 0
## 3930 1 0
## 3931 1 0
## 3932 1 0
## 3933 1 0
## 3934 0 0
## 3938 1 1
## 3939 0 0
## 3940 1 0
## 3941 1 0
## 3943 0 0
## 3944 4 5
## 3945 0 0
## 3947 13 13
## 3948 3 0
## 3949 0 0
## 3950 0 0
## 3951 1 0
## 3954 2 2
## 3955 18 20
## 3957 0 0
## 3959 1 0
## 3963 1 2
## 3966 5 4
## 3968 0 0
## 3969 1 0
## 3970 0 0
## 3971 0 0
## 3972 1 0
## 3973 0 0
## 3974 1 0
## 3975 2 2
## 3976 0 0
## 3977 0 0
## 3978 20 21
## 3979 2 0
## 3980 0 0
## 3981 0 0
## 3982 0 0
## 3983 1 0
## 3985 0 0
## 3988 0 0
## 3990 0 0
## 3991 0 0
## 3992 0 0
## 3993 0 0
## 3994 3 0
## 3995 0 0
## 3996 0 0
## 3997 9 15
## 3998 0 0
## 3999 0 0
## 4000 0 0
## 4001 0 0
## 4002 1 0
## 4003 0 0
## 4005 0 0
## 4007 1 2
## 4008 0 0
## 4009 7 11
## 4010 2 0
## 4011 0 0
## 4012 3 0
## 4013 0 0
## 4014 7 14
## 4015 0 0
## 4016 2 0
## 4017 3 0
## 4019 0 0
## 4020 0 0
## 4023 1 0
## 4024 0 0
## 4025 0 0
## 4026 11 12
## 4027 1 0
## 4028 0 0
## 4029 0 0
## 4030 0 0
## 4031 1 1
## 4032 0 0
## 4033 0 0
## 4034 0 0
## 4036 2 0
## 4038 1 0
## 4039 1 0
## 4040 0 0
## 4041 0 0
## 4042 0 0
## 4043 1 0
## 4044 16 19
## 4045 3 0
## 4046 0 0
## 4047 19 26
## 4048 0 0
## 4049 0 0
## 4050 1 0
## 4051 4 3
## 4052 0 0
## 4053 5 4
## 4054 2 3
## 4056 1 0
## 4057 3 0
## 4058 2 2
## 4059 0 0
## 4060 4 3
## 4061 2 2
## 4062 2 6
## 4063 0 0
## 4064 2 7
## 4065 19 21
## 4066 0 0
## 4067 1 0
## 4068 7 6
## 4070 0 0
## 4072 4 3
## 4073 0 0
## 4075 13 15
## 4078 2 3
## 4079 8 12
## 4080 2 5
## 4081 0 0
## 4082 0 0
## 4084 1 0
## 4087 0 0
## 4088 0 0
## 4089 1 1
## 4090 0 0
## 4091 0 0
## 4092 0 0
## 4093 4 4
## 4095 0 0
## 4097 0 0
## 4098 2 0
## 4099 6 7
## 4100 0 0
## 4101 0 0
## 4102 13 15
## 4103 1 1
## 4104 14 17
## 4105 1 0
## 4106 0 0
## 4107 2 0
## 4108 2 3
## 4109 0 0
## 4112 0 0
## 4114 0 0
## 4115 1 0
## 4116 0 0
## 4117 0 0
## 4118 0 0
## 4119 0 0
## 4120 1 0
## 4122 1 0
## 4124 14 18
## 4125 4 3
## 4126 2 0
## 4127 0 0
## 4128 21 21
## 4129 0 0
## 4132 0 0
## 4135 14 13
## 4136 8 9
## 4137 0 0
## 4138 0 0
## 4139 0 0
## 4141 1 0
## 4142 4 4
## 4143 0 0
## 4145 0 0
## 4146 1 0
## 4147 0 0
## 4148 0 0
## 4150 1 1
## 4152 0 0
## 4153 1 0
## 4154 1 0
## 4155 0 0
## 4156 0 0
## 4157 0 0
## 4158 0 0
## 4159 0 0
## 4160 2 0
## 4162 0 0
## 4163 0 0
## 4166 5 12
## 4167 0 0
## 4168 0 0
## 4169 0 1
## 4170 3 0
## 4171 0 0
## 4172 4 5
## 4173 0 0
## 4174 10 13
## 4175 0 0
## 4176 1 1
## 4178 0 0
## 4179 0 0
## 4180 0 0
## 4181 1 0
## 4182 1 0
## 4183 0 0
## 4184 1 0
## 4185 0 0
## 4186 0 0
## 4187 0 0
## 4188 0 0
## 4189 0 0
## 4191 0 0
## 4192 0 0
## 4193 0 1
## 4194 18 20
## 4195 0 0
## 4196 0 0
## 4197 1 0
## 4199 6 6
## 4200 0 0
## 4201 0 0
## 4202 0 0
## 4206 0 0
## 4207 0 0
## 4208 1 0
## 4209 0 0
## 4210 14 17
## 4211 0 0
## 4213 0 0
## 4214 0 0
## 4215 7 10
## 4216 0 0
## 4217 0 0
## 4218 1 1
## 4219 0 0
## 4220 0 0
## 4222 15 25
## 4223 0 0
## 4224 0 0
## 4225 0 0
## 4226 1 1
## 4227 0 0
## 4228 1 1
## 4230 0 0
## 4231 0 0
## 4232 0 0
## 4233 0 0
## 4234 2 0
## 4235 4 5
## 4236 1 1
## 4237 7 11
## 4238 0 0
## 4241 4 7
## 4242 1 0
## 4243 4 5
## 4244 0 0
## 4245 1 0
## 4246 0 0
## 4247 1 0
## 4248 0 0
## 4249 1 0
## 4250 1 0
## 4251 1 0
## 4252 0 0
## 4253 0 0
## 4255 0 0
## 4256 2 0
## 4257 1 1
## 4258 0 0
## 4261 1 1
## 4262 0 0
## 4263 0 0
## 4265 4 4
## 4266 1 0
## 4267 12 13
## 4268 0 0
## 4269 0 0
## 4270 2 0
## 4271 0 0
## 4272 1 0
## 4273 0 0
## 4275 0 0
## 4276 0 0
## 4277 0 0
## 4278 0 0
## 4279 0 0
## 4280 0 0
## 4281 1 0
## 4283 7 9
## 4285 0 0
## 4287 0 0
## 4288 0 0
## 4289 1 1
## 4294 0 0
## 4295 0 0
## 4296 0 0
## 4297 0 0
## 4298 0 0
## 4299 0 0
## 4301 0 0
## 4302 3 0
## 4303 2 3
## 4304 2 0
## 4305 0 0
## 4306 1 0
## 4307 2 0
## 4308 0 0
## 4309 1 1
## 4310 8 7
## 4312 1 0
## 4314 0 0
## 4315 0 0
## 4316 1 0
## 4317 0 0
## 4318 0 0
## 4319 1 0
## 4320 4 3
## 4321 0 0
## 4322 1 1
## 4324 0 0
## 4325 0 0
## 4326 0 0
## 4327 0 0
## 4328 1 3
## 4329 0 0
## 4331 10 15
## 4332 2 0
## 4333 0 0
## 4334 5 5
## 4335 0 0
## 4336 0 0
## 4337 0 0
## 4338 0 0
## 4339 0 0
## 4340 0 0
## 4341 1 2
## 4342 2 0
## 4343 0 0
## 4344 1 1
## 4345 0 0
## 4346 0 0
## 4347 1 2
## 4348 0 0
## 4350 0 0
## 4351 0 0
## 4352 0 0
## 4353 0 0
## 4354 3 0
## 4355 0 0
## 4356 3 3
## 4357 1 0
## 4358 0 0
## 4359 0 0
## 4360 3 3
## 4361 1 0
## 4362 0 0
## 4363 0 0
## 4364 2 3
## 4365 0 0
## 4366 4 5
## 4368 0 0
## 4369 0 0
## 4370 1 0
## 4371 0 0
## 4372 0 0
## 4373 3 9
## 4374 1 0
## 4375 1 0
## 4376 14 27
## 4377 0 0
## 4378 3 3
## 4379 0 0
## 4380 6 6
## 4381 0 0
## 4382 0 0
## 4383 0 0
## 4384 0 0
## 4385 0 0
## 4386 0 0
## 4387 20 13
## 4388 0 0
## 4389 0 0
## 4391 1 0
## 4392 2 0
## 4393 0 0
## 4394 0 0
## 4395 0 0
## 4396 0 0
## 4397 0 0
## 4398 1 0
## 4399 1 0
## 4400 0 0
## 4401 4 9
## 4402 3 0
## 4403 0 0
## 4404 2 2
## 4407 1 0
## 4408 3 4
## 4409 1 0
## 4410 1 0
## 4411 0 0
## 4412 0 0
## 4413 1 0
## 4414 0 0
## 4416 1 0
## 4417 0 0
## 4418 0 0
## 4419 1 0
## 4420 0 0
## 4421 0 0
## 4422 1 0
## 4425 0 0
## 4426 3 0
## 4427 2 0
## 4428 0 0
## 4429 1 0
## 4431 11 18
## 4432 0 0
## 4434 10 15
## 4435 0 0
## 4437 0 1
## 4438 0 0
## 4439 11 18
## 4440 1 0
## 4443 2 2
## 4446 0 0
## 4447 13 21
## 4448 1 1
## 4449 0 0
## 4451 2 2
## 4452 0 0
## 4453 2 0
## 4454 1 0
## 4456 15 20
## 4457 1 2
## 4458 0 0
## 4459 2 3
## 4460 0 0
## 4461 0 0
## 4464 1 1
## 4466 12 20
## 4467 1 0
## 4470 1 0
## 4471 1 0
## 4472 1 0
## 4473 0 0
## 4474 0 0
## 4475 1 0
## 4476 0 0
## 4477 0 0
## 4479 0 0
## 4480 1 2
## 4481 0 0
## 4482 2 3
## 4483 0 0
## 4484 2 0
## 4485 0 0
## 4486 1 1
## 4487 1 3
## 4490 1 1
## 4492 0 0
## 4493 0 0
## 4494 0 0
## 4495 0 0
## 4496 0 0
## 4497 22 24
## 4498 0 0
## 4499 0 1
## 4500 3 3
## 4503 0 0
## 4504 0 0
## 4505 0 0
## 4508 0 0
## 4509 1 0
## 4510 1 0
## 4511 1 0
## 4513 1 1
## 4514 0 0
## 4515 0 0
## 4516 2 5
## 4517 11 15
## 4518 0 0
## 4522 3 3
## 4523 5 6
## 4524 1 0
## 4526 0 0
## 4527 0 0
## 4529 2 3
## 4530 2 3
## 4531 1 0
## 4533 0 0
## 4536 1 2
## 4537 0 0
## 4539 1 0
## 4540 1 0
## 4541 0 0
## 4542 14 14
## 4543 1 2
## 4544 0 0
## 4546 0 0
## 4547 2 4
## 4548 1 1
## 4549 1 1
## 4550 5 7
## 4552 1 0
## 4553 0 0
## 4554 0 0
## 4555 0 0
## 4556 0 0
## 4557 0 0
## 4558 0 0
## 4559 0 0
## 4561 0 0
## 4562 0 0
## 4563 0 0
## 4564 6 13
## 4565 4 4
## 4566 1 0
## 4567 0 0
## 4569 0 0
## 4570 1 0
## 4571 3 3
## 4572 2 0
## 4573 7 10
## 4574 0 0
## 4575 4 5
## 4576 3 3
## 4578 0 0
## 4579 1 0
## 4580 0 0
## 4581 1 1
## 4582 0 0
## 4584 2 0
## 4586 0 0
## 4587 0 0
## 4589 7 10
## 4590 1 0
## 4592 0 0
## 4593 1 0
## 4594 1 1
## 4596 2 2
## 4598 0 0
## 4599 0 0
## 4600 0 0
## 4604 3 4
## 4607 0 0
## 4608 1 0
## 4609 0 0
## 4610 19 25
## 4616 1 0
## 4617 0 0
## 4618 1 2
## 4619 0 0
## 4620 1 0
## 4621 0 0
## 4622 0 0
## 4623 1 0
## 4624 1 0
## 4625 1 0
## 4626 0 0
## 4627 0 0
## 4629 0 0
## 4630 1 0
## 4631 0 0
## 4632 0 0
## 4636 1 0
## 4637 1 0
## 4638 0 0
## 4639 0 0
## 4642 0 0
## 4643 0 0
## 4644 0 0
## 4645 0 0
## 4646 0 0
## 4647 1 0
## 4648 1 0
## 4649 0 0
## 4650 2 0
## 4652 0 0
## 4653 2 0
## 4654 21 25
## 4655 0 0
## 4657 2 0
## 4658 14 12
## 4659 4 5
## 4660 0 0
## 4661 1 0
## 4662 1 0
## 4663 0 0
## 4664 0 0
## 4665 0 0
## 4666 0 0
## 4667 0 0
## 4668 1 0
## 4669 0 0
## 4670 0 0
## 4671 0 0
## 4672 0 0
## 4673 0 0
## 4675 0 0
## 4677 0 0
## 4678 3 4
## 4681 1 1
## 4682 0 0
## 4683 8 9
## 4684 1 0
## 4685 0 0
## 4686 3 0
## 4687 1 0
## 4688 1 1
## 4689 0 0
## 4690 0 0
## 4691 1 1
## 4692 0 0
## 4693 2 2
## 4694 0 0
## 4695 0 0
## 4696 0 0
## 4697 0 0
## 4700 0 0
## 4701 0 0
## 4702 1 0
## 4703 0 1
## 4704 3 4
## 4705 1 1
## 4706 0 0
## 4707 1 1
## 4710 2 0
## 4711 1 0
## 4712 0 0
## 4713 1 0
## 4715 0 0
## 4716 4 16
## 4717 0 0
## 4718 0 0
## 4719 0 0
## 4720 1 0
## 4723 1 0
## 4724 0 0
## 4725 1 1
## 4726 3 5
## 4731 0 0
## 4732 2 0
## 4733 3 0
## 4734 1 2
## 4735 0 0
## 4736 0 0
## 4737 13 15
## 4738 1 1
## 4739 0 0
## 4740 19 26
## 4741 0 0
## 4743 0 0
## 4744 0 0
## 4745 0 0
## 4747 0 3
## 4749 9 9
## 4751 0 0
## 4752 3 3
## 4754 1 0
## 4755 2 2
## 4756 0 0
## 4757 1 0
## 4758 0 0
## 4759 0 0
## 4761 0 0
## 4762 0 0
## 4763 0 0
## 4764 15 20
## 4765 3 4
## 4768 1 1
## 4770 5 6
## 4772 0 0
## 4773 1 0
## 4776 0 0
## 4782 0 0
## 4785 1 0
## 4786 0 0
## 4788 1 0
## 4789 0 0
## 4790 11 19
## 4793 3 8
## 4794 0 0
## 4795 2 0
## 4797 2 0
## 4798 20 19
## 4799 1 0
## 4800 6 5
## 4801 0 0
## 4803 6 5
## 4804 7 8
## 4805 0 0
## 4806 0 0
## 4807 0 0
## 4808 6 5
## 4810 0 0
## 4811 0 0
## 4812 0 0
## 4813 0 0
## 4814 9 13
## 4817 0 0
## 4818 1 0
## 4819 0 0
## 4820 0 0
## 4821 8 6
## 4822 0 0
## 4823 0 0
## 4824 1 0
## 4825 5 5
## 4829 0 0
## 4830 0 0
## 4831 19 26
## 4832 0 0
## 4833 2 2
## 4835 0 0
## 4837 0 0
## 4842 0 0
## 4843 12 15
## 4845 1 1
## 4846 1 1
## 4847 0 0
## 4848 3 4
## 4851 0 0
## 4852 1 0
## 4853 0 0
## 4855 0 0
## 4859 1 1
## 4861 0 0
## 4862 1 0
## 4863 17 18
## 4865 0 0
## 4866 2 2
## 4868 0 0
## 4870 0 0
## 4871 2 0
## 4872 20 27
## 4873 1 0
## 4874 0 0
## 4875 0 0
## 4876 0 0
## 4877 0 0
## 4879 0 0
## 4882 3 3
## 4883 1 0
## 4884 0 0
## 4885 0 0
## 4886 0 0
## 4887 0 0
## 4888 0 0
## 4889 1 0
## 4894 0 0
## 4895 0 0
## 4896 0 0
## 4897 3 3
## 4899 0 0
## 4900 0 0
## 4901 0 0
## 4902 19 20
## 4903 2 0
## 4905 0 0
## 4906 3 4
## 4907 9 14
## 4909 0 0
## 4910 3 2
## 4911 15 18
## 4912 1 0
## 4914 0 0
## 4916 0 0
## 4919 0 0
## 4920 0 0
## 4921 0 0
## 4922 0 0
## 4925 0 0
## 4926 1 0
## 4927 2 2
## 4928 0 0
## 4929 0 0
## 4930 1 0
## 4931 0 0
## 4932 17 26
## 4933 1 0
## 4934 1 0
## 4935 1 0
## 4936 0 0
## 4937 3 3
## 4938 0 0
## 4939 5 7
## 4940 0 0
## 4941 0 0
## 4942 4 4
## 4943 0 0
## 4944 0 0
## 4945 1 0
# Correspondencia una X a una sola y
tabla_inf <- aggregate(
DAYS_LOST ~ DAYS_RESTRICT,
data = TPV_limpio,
FUN = min
)
tabla_inf## DAYS_RESTRICT DAYS_LOST
## 1 0 0
## 2 1 0
## 3 2 0
## 4 3 0
## 5 4 3
## 6 5 3
## 7 6 5
## 8 7 5
## 9 8 6
## 10 9 7
## 11 10 8
## 12 11 10
## 13 12 9
## 14 13 10
## 15 14 11
## 16 15 12
## 17 16 13
## 18 17 16
## 19 18 14
## 20 19 13
## 21 20 13
## 22 21 15
## 23 22 17
# Diagrama de dispersión para una x a una sola y (mínimos)
plot(tabla_inf$DAYS_RESTRICT,
tabla_inf$DAYS_LOST,
pch = 16,
col = "orange",
main = "Gráfica N°1: Diagrama de dispersión entre días restringidos y días
perdidos por accidentes, lesiones y enfermedades en operaciones mineras",
xlab = "Días con trabajo restringido",
ylab = "Días perdidos",
xlim = c(0, max(tabla_inf$DAYS_RESTRICT, na.rm = TRUE)),
ylim = c(0, max(tabla_inf$DAYS_LOST, na.rm = TRUE)))Debido a la forma creciente de la nube de puntos, se conjetura un modelo exponencial.
# Ajuste del modelo de regresión exponencial
x <- tabla_inf$DAYS_RESTRICT
y <- tabla_inf$DAYS_LOST
# Transformación logarítmica
y_log <- log(y + 1)
# Ajuste del modelo
modelo_exp <- lm(y_log ~ x)
summary(modelo_exp)##
## Call:
## lm(formula = y_log ~ x)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.8629 -0.4398 0.1441 0.4000 0.5375
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.47160 0.18145 2.599 0.0168 *
## x 0.13044 0.01413 9.234 7.67e-09 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4494 on 21 degrees of freedom
## Multiple R-squared: 0.8024, Adjusted R-squared: 0.793
## F-statistic: 85.27 on 1 and 21 DF, p-value: 7.674e-09
# Parámetros del modelo
beta0 <- coef(modelo_exp)[1]
beta1 <- coef(modelo_exp)[2]
a <- exp(beta0)
b <- beta1
cat("Parámetro a:", round(a, 2), "\n")## Parámetro a: 1.6
## Parámetro b: 0.1304
DIAGRAMA DE DISPERSIÓN
# Diagrama de dispersión
plot(x, y,
pch = 16,
col = "blue",
main = "Gráfica N°4: Comparación realidad vs modelo exponencial entre días restringidos
y días perdidos por accidentes, lesiones y enfermedades en operaciones mineras",
xlab = "Días con trabajo restringido",
ylab = "Días perdidos",
xlim = c(0, max(x, na.rm = TRUE)),
ylim = c(0, max(y, na.rm = TRUE)))
curve(a * exp(b * x),
from = 0,
to = max(x, na.rm = TRUE),
add = TRUE,
col = "red",
lwd = 2)
# Leyenda
legend("topleft",
legend = c("Datos reales", "Modelo exponencial"),
col = c("blue", "red"),
pch = c(16, NA),
lty = c(NA, 1),
lwd = c(NA, 2))# ECUACIÓN DEL MODELO
plot(1, type = "n", axes = FALSE, xlab = "", ylab = "")
text(1, 1,
labels = paste0(
"Ecuación Exponencial\n",
"Y = ae^bx\n",
"Y = ", round(a, 2), " e^(", round(b, 4), "x)"
),
cex = 2,
font = 6)TEST DE PEARSON
## Porcentaje de aprobación: 89.58 %
if(r > 0.8) {
cat("✓ APRUEBA el test de Pearson ( >80% )\n")
} else {
cat("✗ NO aprueba el test de Pearson\n")
}## ✓ APRUEBA el test de Pearson ( >80% )
RESTRICCIONES
## Rango de DAYS_RESTRICT observado:
## Mínimo: 0
## Máximo: 22
## ✓ El modelo es válido para valores positivos de días restringidos
## Pronósticos:
cat("Si hay 5 días con trabajo restringido, se esperan aproximadamente",
round(pronostico_5, 2), "días perdidos.\n")## Si hay 5 días con trabajo restringido, se esperan aproximadamente 3.08 días perdidos.
# Visualización de pronósticos
plot(1, type = "n", axes = FALSE, xlab = "", ylab = "")
text(1, 1,
labels = paste0(
"Pronóstico del Modelo:\n\n",
"5 días restringidos → ", round(pronostico_5, 1), " días perdidos\n"
),
cex = 1.8,
font = 6)Entre los días restringidos y los días perdidos existe una relación exponencial donde el modelo f(x)=1.6e^(0.1304x) siendo “x” los días restringidos y “y” los días perdidos. Sí existen restricciones, ya que el modelo solo es válido para valores positivos (días).
Ejemplo: Si hay 5 días restringidos se esperan aproximádamente 3.1 días perdidos.