1 1. Cargar Paquetes

library(readr)
library(dplyr)
library(ggplot2)
library(stringr)
library(tidyr)
library(BSDA)

2 2. Cargar y Preparar la Base de Datos

setwd("C:/Users/Terumo/Downloads")
insurance <- read_csv("smoking_health_data_final.csv")
insurance <- as.data.frame(unclass(insurance), stringsAsFactors = TRUE)

str(insurance)
## 'data.frame':    3900 obs. of  7 variables:
##  $ age           : num  54 45 58 42 42 57 43 42 37 49 ...
##  $ sex           : Factor w/ 2 levels "female","male": 2 2 2 2 2 2 2 2 2 2 ...
##  $ current_smoker: Factor w/ 2 levels "no","yes": 2 2 2 2 2 2 2 2 2 2 ...
##  $ heart_rate    : num  95 64 81 90 62 62 75 66 65 93 ...
##  $ blood_pressure: Factor w/ 2317 levels "100.5/62","100.5/66",..: 230 642 869 671 572 143 198 723 709 874 ...
##  $ cigs_per_day  : num  NA NA NA NA NA NA NA NA NA NA ...
##  $ chol          : num  219 248 235 225 226 223 222 196 188 256 ...
dim(insurance)
## [1] 3900    7

3 3. Limpieza de Datos

insurance <- insurance %>% distinct()
colSums(is.na(insurance))
##            age            sex current_smoker     heart_rate blood_pressure 
##              0              0              0              0              0 
##   cigs_per_day           chol 
##             14              7
summary(insurance)
##       age            sex       current_smoker   heart_rate     blood_pressure
##  Min.   :32.00   female:2081   no :1968       Min.   : 44.00   130/80 :  18  
##  1st Qu.:42.00   male  :1819   yes:1932       1st Qu.: 68.00   120/80 :  17  
##  Median :49.00                                Median : 75.00   110/70 :  15  
##  Mean   :49.54                                Mean   : 75.69   125/80 :  15  
##  3rd Qu.:56.00                                3rd Qu.: 82.00   105/70 :   9  
##  Max.   :70.00                                Max.   :143.00   107/73 :   9  
##                                                                (Other):3817  
##   cigs_per_day         chol      
##  Min.   : 0.000   Min.   :113.0  
##  1st Qu.: 0.000   1st Qu.:206.0  
##  Median : 0.000   Median :234.0  
##  Mean   : 9.169   Mean   :236.6  
##  3rd Qu.:20.000   3rd Qu.:263.0  
##  Max.   :70.000   Max.   :696.0  
##  NA's   :14       NA's   :7
unique(insurance$sex)
## [1] male   female
## Levels: female male
unique(insurance$smoker)
## NULL
unique(insurance$region)
## NULL
age_sex <- insurance %>% arrange(age, sex)

4 4. Enriquecer la Base con Nuevas Variables

insurance <- insurance %>%
  mutate(
    grupo_edad = case_when(
      age < 25 ~ "Joven",
      age < 45 ~ "Adulto joven",
      age < 65 ~ "Adulto maduro",
      TRUE ~ "Mayor"
    )) %>%
  separate(blood_pressure, into = c("systolic", "diastolic"), sep = "/") %>%
  mutate(across(c(systolic, diastolic), as.numeric))

5 5. Visualización Exploratoria

resumen_heatmap <- insurance %>%
  group_by(sex, current_smoker) %>%
  summarise(promedio_charges = mean(heart_rate), .groups = "drop")

ggplot(resumen_heatmap, aes(x = sex, y = current_smoker, fill = promedio_charges)) +
  geom_tile(color = "white") +
  geom_text(aes(label = round(promedio_charges, 0)), size = 4) +
  scale_fill_gradient(low = "lightblue", high = "darkred") +
  labs(title = "Promedio de ritmo cardiaco por Sexo y fumador", fill = "Promedio") +
  theme_minimal()

resumen_heatmap <- insurance %>%
  group_by(sex, current_smoker) %>%
  summarise(promedio_charges = mean(diastolic), .groups = "drop")

ggplot(resumen_heatmap, aes(x = sex, y = current_smoker, fill = promedio_charges)) +
  geom_tile(color = "white") +
  geom_text(aes(label = round(promedio_charges, 0)), size = 4) +
  scale_fill_gradient(low = "lightblue", high = "darkred") +
  labs(title = "Promedio de diastole por Sexo y fumador", fill = "Promedio") +
  theme_minimal()

resumen_heatmap <- insurance %>%
  group_by(sex, current_smoker) %>%
  summarise(promedio_charges = mean(systolic), .groups = "drop")

ggplot(resumen_heatmap, aes(x = sex, y = current_smoker, fill = promedio_charges)) +
  geom_tile(color = "white") +
  geom_text(aes(label = round(promedio_charges, 0)), size = 4) +
  scale_fill_gradient(low = "lightblue", high = "darkred") +
  labs(title = "Promedio de sistole por Sexo y fumador", fill = "Promedio") +
  theme_minimal()

6 6. Estadísticas Descriptivas Adicionales

insurance %>% summarise(media_hr = mean(heart_rate), sd_hr = sd(heart_rate))
##   media_hr    sd_hr
## 1 75.68897 12.01538
insurance %>% group_by(grupo_edad) %>% 
  summarise(media_hr = mean(heart_rate), sd_hr = sd(heart_rate))
## # A tibble: 3 × 3
##   grupo_edad    media_hr sd_hr
##   <chr>            <dbl> <dbl>
## 1 Adulto joven      75.7  11.9
## 2 Adulto maduro     75.7  12.0
## 3 Mayor             76.5  12.9
insurance %>% group_by(grupo_edad, current_smoker) %>%
  summarise(media_hr = mean(heart_rate), sd_hr = sd(heart_rate))
## # A tibble: 6 × 4
## # Groups:   grupo_edad [3]
##   grupo_edad    current_smoker media_hr sd_hr
##   <chr>         <fct>             <dbl> <dbl>
## 1 Adulto joven  no                 74.8  12.5
## 2 Adulto joven  yes                76.2  11.4
## 3 Adulto maduro no                 75.0  12.1
## 4 Adulto maduro yes                76.5  11.9
## 5 Mayor         no                 76.5  12.8
## 6 Mayor         yes                76.5  13.2
insurance %>% group_by(grupo_edad, current_smoker) %>%
  summarise(media_systolic = mean(systolic), sd_systolic = sd(systolic))
## # A tibble: 6 × 4
## # Groups:   grupo_edad [3]
##   grupo_edad    current_smoker media_systolic sd_systolic
##   <chr>         <fct>                   <dbl>       <dbl>
## 1 Adulto joven  no                       124.        15.3
## 2 Adulto joven  yes                      123.        16.1
## 3 Adulto maduro no                       139.        23.5
## 4 Adulto maduro yes                      134.        21.5
## 5 Mayor         no                       149.        28.6
## 6 Mayor         yes                      146.        22.5

7 7. Exploración Rápida

glimpse(insurance)
## Rows: 3,900
## Columns: 9
## $ age            <dbl> 54, 45, 58, 42, 42, 57, 43, 42, 37, 49, 55, 39, 53, 45,…
## $ sex            <fct> male, male, male, male, male, male, male, male, male, m…
## $ current_smoker <fct> yes, yes, yes, yes, yes, yes, yes, yes, yes, yes, yes, …
## $ heart_rate     <dbl> 95, 64, 81, 90, 62, 62, 75, 66, 65, 93, 70, 85, 58, 83,…
## $ systolic       <dbl> 110.0, 121.0, 127.5, 122.5, 119.0, 107.5, 109.5, 123.0,…
## $ diastolic      <dbl> 72.0, 72.0, 76.0, 80.0, 80.0, 72.5, 69.0, 73.0, 77.0, 8…
## $ cigs_per_day   <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,…
## $ chol           <dbl> 219, 248, 235, 225, 226, 223, 222, 196, 188, 256, 214, …
## $ grupo_edad     <chr> "Adulto maduro", "Adulto maduro", "Adulto maduro", "Adu…
insurance %>% select(sex, age)
##         sex age
## 1      male  54
## 2      male  45
## 3      male  58
## 4      male  42
## 5      male  42
## 6      male  57
## 7      male  43
## 8      male  42
## 9      male  37
## 10     male  49
## 11     male  55
## 12     male  39
## 13     male  53
## 14     male  45
## 15     male  48
## 16   female  58
## 17   female  46
## 18   female  45
## 19   female  37
## 20   female  46
## 21   female  40
## 22   female  56
## 23   female  38
## 24   female  58
## 25   female  56
## 26     male  65
## 27   female  41
## 28     male  58
## 29   female  39
## 30   female  38
## 31     male  44
## 32   female  43
## 33   female  44
## 34   female  36
## 35   female  41
## 36   female  65
## 37     male  64
## 38   female  55
## 39     male  68
## 40   female  39
## 41   female  38
## 42   female  54
## 43   female  57
## 44   female  41
## 45   female  52
## 46   female  41
## 47   female  60
## 48   female  38
## 49   female  67
## 50     male  41
## 51   female  53
## 52   female  43
## 53   female  35
## 54     male  64
## 55   female  39
## 56   female  46
## 57   female  39
## 58   female  42
## 59   female  34
## 60     male  53
## 61     male  51
## 62   female  34
## 63     male  49
## 64   female  57
## 65     male  48
## 66   female  40
## 67     male  63
## 68     male  42
## 69     male  35
## 70     male  40
## 71     male  62
## 72     male  39
## 73   female  40
## 74   female  59
## 75   female  49
## 76   female  53
## 77   female  39
## 78     male  57
## 79   female  40
## 80     male  45
## 81   female  43
## 82   female  42
## 83   female  45
## 84     male  54
## 85   female  49
## 86   female  40
## 87   female  59
## 88   female  46
## 89   female  52
## 90   female  59
## 91   female  38
## 92   female  57
## 93   female  38
## 94   female  42
## 95   female  48
## 96   female  38
## 97   female  61
## 98   female  38
## 99   female  50
## 100  female  61
## 101  female  44
## 102  female  54
## 103  female  55
## 104  female  50
## 105  female  36
## 106  female  58
## 107  female  56
## 108    male  53
## 109  female  46
## 110  female  53
## 111    male  67
## 112  female  53
## 113    male  64
## 114  female  65
## 115    male  52
## 116  female  41
## 117    male  66
## 118  female  38
## 119  female  48
## 120  female  41
## 121    male  56
## 122    male  57
## 123    male  60
## 124  female  46
## 125  female  58
## 126  female  50
## 127  female  42
## 128    male  51
## 129  female  44
## 130    male  52
## 131  female  49
## 132    male  55
## 133  female  48
## 134  female  38
## 135  female  55
## 136  female  36
## 137  female  66
## 138  female  67
## 139    male  41
## 140    male  50
## 141    male  61
## 142    male  52
## 143    male  41
## 144  female  62
## 145  female  42
## 146    male  39
## 147    male  63
## 148  female  46
## 149    male  37
## 150  female  39
## 151  female  50
## 152  female  57
## 153    male  45
## 154  female  55
## 155    male  67
## 156  female  64
## 157    male  51
## 158  female  48
## 159  female  64
## 160  female  45
## 161    male  52
## 162  female  55
## 163  female  48
## 164    male  37
## 165    male  64
## 166  female  56
## 167  female  59
## 168    male  53
## 169    male  39
## 170  female  52
## 171  female  66
## 172    male  68
## 173    male  63
## 174  female  38
## 175  female  67
## 176  female  39
## 177    male  59
## 178  female  40
## 179  female  51
## 180    male  57
## 181  female  46
## 182    male  56
## 183  female  57
## 184    male  53
## 185  female  55
## 186  female  36
## 187  female  58
## 188    male  57
## 189  female  40
## 190  female  57
## 191  female  38
## 192    male  44
## 193  female  45
## 194  female  50
## 195  female  55
## 196    male  55
## 197    male  38
## 198    male  62
## 199  female  53
## 200    male  43
## 201    male  63
## 202    male  43
## 203  female  60
## 204    male  60
## 205    male  53
## 206    male  54
## 207    male  48
## 208  female  48
## 209  female  48
## 210    male  46
## 211  female  56
## 212    male  41
## 213    male  57
## 214    male  48
## 215  female  51
## 216  female  54
## 217  female  63
## 218  female  50
## 219  female  46
## 220    male  44
## 221  female  37
## 222    male  42
## 223    male  61
## 224  female  63
## 225  female  51
## 226    male  39
## 227  female  53
## 228    male  48
## 229  female  59
## 230    male  57
## 231  female  51
## 232    male  62
## 233  female  39
## 234    male  56
## 235  female  39
## 236  female  54
## 237  female  45
## 238    male  42
## 239    male  48
## 240  female  49
## 241    male  41
## 242    male  53
## 243    male  42
## 244  female  48
## 245  female  37
## 246    male  45
## 247  female  45
## 248    male  55
## 249  female  55
## 250    male  62
## 251  female  42
## 252  female  55
## 253    male  44
## 254  female  46
## 255  female  56
## 256  female  51
## 257  female  61
## 258  female  48
## 259    male  53
## 260  female  42
## 261  female  46
## 262    male  62
## 263    male  65
## 264    male  57
## 265  female  60
## 266    male  60
## 267  female  58
## 268  female  53
## 269  female  65
## 270  female  60
## 271  female  51
## 272  female  64
## 273  female  54
## 274  female  42
## 275    male  43
## 276  female  63
## 277  female  49
## 278  female  57
## 279    male  34
## 280  female  60
## 281  female  64
## 282  female  64
## 283  female  55
## 284    male  45
## 285  female  59
## 286    male  55
## 287  female  39
## 288  female  50
## 289  female  51
## 290  female  47
## 291  female  47
## 292    male  64
## 293  female  45
## 294    male  59
## 295    male  38
## 296  female  48
## 297    male  41
## 298  female  54
## 299  female  46
## 300  female  39
## 301  female  57
## 302  female  53
## 303  female  41
## 304    male  53
## 305  female  59
## 306  female  52
## 307    male  48
## 308  female  70
## 309  female  52
## 310  female  55
## 311  female  55
## 312    male  50
## 313  female  54
## 314    male  41
## 315  female  55
## 316    male  37
## 317    male  55
## 318  female  55
## 319  female  48
## 320    male  49
## 321    male  51
## 322  female  64
## 323    male  38
## 324    male  59
## 325    male  58
## 326  female  43
## 327  female  41
## 328    male  41
## 329  female  60
## 330  female  66
## 331    male  65
## 332  female  50
## 333  female  49
## 334  female  52
## 335  female  54
## 336    male  42
## 337    male  48
## 338    male  46
## 339    male  50
## 340  female  43
## 341    male  55
## 342  female  43
## 343    male  36
## 344  female  57
## 345  female  56
## 346    male  66
## 347  female  49
## 348    male  38
## 349  female  41
## 350  female  50
## 351  female  41
## 352  female  42
## 353  female  41
## 354  female  58
## 355  female  63
## 356  female  63
## 357  female  42
## 358  female  59
## 359  female  67
## 360  female  59
## 361  female  55
## 362  female  45
## 363    male  44
## 364  female  49
## 365  female  45
## 366    male  52
## 367    male  48
## 368  female  44
## 369  female  61
## 370    male  58
## 371  female  54
## 372  female  55
## 373  female  55
## 374  female  63
## 375    male  58
## 376    male  50
## 377  female  53
## 378  female  67
## 379  female  51
## 380  female  57
## 381    male  53
## 382    male  42
## 383  female  55
## 384  female  64
## 385  female  44
## 386  female  49
## 387    male  38
## 388  female  55
## 389  female  49
## 390    male  54
## 391  female  55
## 392    male  46
## 393    male  41
## 394    male  67
## 395  female  52
## 396  female  55
## 397    male  39
## 398  female  40
## 399    male  37
## 400    male  50
## 401    male  59
## 402  female  49
## 403    male  46
## 404  female  51
## 405    male  51
## 406    male  43
## 407  female  62
## 408  female  43
## 409  female  50
## 410    male  60
## 411  female  63
## 412  female  61
## 413    male  44
## 414    male  38
## 415  female  62
## 416    male  56
## 417    male  52
## 418    male  60
## 419  female  65
## 420    male  48
## 421    male  56
## 422  female  63
## 423    male  61
## 424  female  47
## 425    male  47
## 426    male  51
## 427    male  43
## 428  female  66
## 429  female  47
## 430  female  50
## 431    male  58
## 432  female  67
## 433  female  66
## 434    male  47
## 435  female  68
## 436    male  61
## 437    male  60
## 438    male  57
## 439    male  48
## 440  female  60
## 441    male  50
## 442  female  67
## 443    male  40
## 444    male  55
## 445  female  48
## 446  female  61
## 447    male  65
## 448  female  64
## 449  female  46
## 450    male  65
## 451  female  53
## 452    male  54
## 453    male  50
## 454    male  54
## 455  female  63
## 456  female  65
## 457  female  57
## 458    male  46
## 459  female  42
## 460  female  66
## 461  female  64
## 462    male  57
## 463    male  58
## 464    male  63
## 465  female  58
## 466    male  57
## 467    male  49
## 468    male  48
## 469  female  52
## 470  female  59
## 471    male  47
## 472  female  62
## 473  female  57
## 474    male  51
## 475  female  58
## 476  female  46
## 477    male  50
## 478  female  61
## 479  female  59
## 480    male  61
## 481  female  51
## 482    male  45
## 483    male  56
## 484  female  54
## 485    male  40
## 486  female  42
## 487    male  45
## 488  female  46
## 489  female  47
## 490  female  55
## 491    male  64
## 492    male  44
## 493  female  50
## 494    male  64
## 495  female  41
## 496  female  35
## 497    male  37
## 498    male  58
## 499  female  44
## 500    male  53
## 501  female  57
## 502  female  45
## 503  female  61
## 504    male  43
## 505  female  65
## 506    male  39
## 507    male  52
## 508    male  46
## 509  female  55
## 510  female  40
## 511  female  61
## 512  female  58
## 513    male  38
## 514  female  42
## 515  female  59
## 516  female  45
## 517  female  40
## 518  female  38
## 519    male  43
## 520  female  52
## 521    male  48
## 522    male  48
## 523  female  63
## 524  female  61
## 525  female  46
## 526    male  47
## 527  female  54
## 528  female  37
## 529  female  59
## 530    male  37
## 531    male  51
## 532  female  47
## 533  female  39
## 534    male  54
## 535  female  51
## 536    male  41
## 537  female  45
## 538  female  47
## 539  female  63
## 540    male  54
## 541  female  50
## 542    male  54
## 543    male  60
## 544  female  57
## 545    male  48
## 546  female  63
## 547    male  43
## 548  female  51
## 549    male  40
## 550  female  58
## 551  female  59
## 552  female  63
## 553  female  55
## 554  female  45
## 555  female  62
## 556    male  55
## 557    male  40
## 558  female  52
## 559    male  53
## 560  female  61
## 561    male  46
## 562  female  49
## 563    male  62
## 564  female  59
## 565    male  61
## 566  female  37
## 567  female  42
## 568    male  37
## 569  female  47
## 570  female  55
## 571    male  39
## 572  female  53
## 573    male  36
## 574    male  52
## 575  female  42
## 576    male  44
## 577    male  39
## 578    male  51
## 579  female  59
## 580  female  52
## 581  female  36
## 582  female  36
## 583  female  41
## 584  female  48
## 585  female  53
## 586  female  35
## 587  female  52
## 588  female  42
## 589    male  42
## 590  female  37
## 591  female  51
## 592  female  35
## 593    male  64
## 594  female  63
## 595    male  50
## 596  female  62
## 597  female  44
## 598    male  36
## 599  female  41
## 600  female  65
## 601  female  55
## 602  female  60
## 603  female  48
## 604    male  42
## 605  female  59
## 606  female  35
## 607    male  52
## 608  female  47
## 609  female  41
## 610    male  42
## 611    male  38
## 612  female  53
## 613    male  54
## 614    male  41
## 615    male  41
## 616  female  45
## 617    male  64
## 618  female  35
## 619  female  49
## 620  female  43
## 621  female  52
## 622  female  44
## 623  female  41
## 624  female  45
## 625    male  53
## 626  female  67
## 627    male  36
## 628  female  45
## 629  female  69
## 630    male  38
## 631  female  61
## 632  female  46
## 633    male  65
## 634  female  65
## 635  female  48
## 636    male  44
## 637    male  38
## 638  female  64
## 639  female  47
## 640    male  60
## 641  female  51
## 642  female  43
## 643  female  67
## 644    male  52
## 645  female  51
## 646  female  56
## 647  female  64
## 648  female  42
## 649    male  65
## 650  female  48
## 651    male  57
## 652    male  64
## 653  female  51
## 654  female  54
## 655    male  54
## 656    male  57
## 657    male  48
## 658  female  57
## 659  female  45
## 660  female  47
## 661    male  39
## 662  female  46
## 663    male  44
## 664    male  61
## 665  female  51
## 666    male  59
## 667  female  46
## 668  female  54
## 669    male  37
## 670  female  50
## 671  female  35
## 672  female  46
## 673    male  56
## 674  female  46
## 675  female  46
## 676  female  43
## 677    male  64
## 678  female  46
## 679    male  61
## 680  female  45
## 681  female  46
## 682  female  55
## 683    male  44
## 684    male  55
## 685  female  62
## 686  female  40
## 687    male  38
## 688  female  49
## 689    male  48
## 690    male  49
## 691  female  54
## 692  female  61
## 693  female  40
## 694  female  53
## 695    male  47
## 696  female  57
## 697  female  48
## 698  female  44
## 699  female  45
## 700  female  63
## 701    male  37
## 702    male  63
## 703  female  50
## 704  female  55
## 705  female  40
## 706  female  53
## 707    male  38
## 708  female  51
## 709    male  63
## 710    male  44
## 711    male  47
## 712    male  60
## 713    male  65
## 714  female  38
## 715    male  65
## 716    male  53
## 717    male  37
## 718  female  55
## 719    male  50
## 720  female  55
## 721  female  55
## 722  female  60
## 723    male  48
## 724    male  40
## 725  female  41
## 726  female  48
## 727    male  42
## 728  female  63
## 729  female  58
## 730    male  48
## 731  female  69
## 732    male  58
## 733  female  59
## 734  female  59
## 735  female  44
## 736  female  43
## 737  female  55
## 738  female  50
## 739    male  68
## 740  female  40
## 741    male  56
## 742    male  50
## 743    male  59
## 744    male  50
## 745  female  48
## 746  female  55
## 747    male  46
## 748  female  48
## 749    male  42
## 750  female  49
## 751    male  52
## 752  female  51
## 753  female  60
## 754  female  57
## 755  female  57
## 756    male  61
## 757  female  52
## 758  female  38
## 759  female  54
## 760  female  43
## 761  female  46
## 762  female  64
## 763  female  55
## 764  female  64
## 765  female  62
## 766  female  50
## 767  female  63
## 768    male  62
## 769    male  36
## 770    male  49
## 771  female  60
## 772  female  53
## 773    male  36
## 774  female  45
## 775  female  52
## 776  female  60
## 777    male  38
## 778  female  67
## 779  female  38
## 780  female  52
## 781  female  57
## 782  female  40
## 783    male  51
## 784  female  44
## 785    male  46
## 786  female  63
## 787  female  58
## 788  female  54
## 789  female  56
## 790    male  40
## 791  female  55
## 792  female  47
## 793  female  50
## 794  female  60
## 795  female  59
## 796  female  59
## 797  female  40
## 798  female  54
## 799  female  57
## 800    male  52
## 801  female  67
## 802  female  62
## 803    male  53
## 804    male  52
## 805  female  54
## 806    male  37
## 807  female  52
## 808  female  55
## 809  female  49
## 810  female  50
## 811    male  45
## 812  female  63
## 813    male  60
## 814    male  36
## 815  female  53
## 816    male  52
## 817  female  61
## 818  female  60
## 819  female  62
## 820  female  54
## 821  female  65
## 822  female  49
## 823  female  40
## 824  female  61
## 825  female  39
## 826  female  53
## 827  female  56
## 828  female  51
## 829  female  63
## 830    male  61
## 831  female  63
## 832  female  60
## 833  female  67
## 834  female  60
## 835  female  53
## 836  female  61
## 837    male  45
## 838  female  54
## 839  female  52
## 840    male  40
## 841  female  53
## 842    male  43
## 843  female  54
## 844  female  41
## 845  female  61
## 846  female  52
## 847  female  50
## 848  female  55
## 849  female  54
## 850  female  52
## 851  female  39
## 852    male  62
## 853    male  57
## 854    male  48
## 855  female  67
## 856  female  51
## 857  female  54
## 858  female  63
## 859  female  62
## 860    male  65
## 861  female  41
## 862  female  43
## 863  female  54
## 864  female  39
## 865  female  40
## 866  female  41
## 867  female  64
## 868  female  64
## 869  female  49
## 870    male  54
## 871  female  49
## 872    male  42
## 873  female  65
## 874  female  47
## 875  female  45
## 876  female  44
## 877    male  62
## 878  female  52
## 879  female  55
## 880  female  46
## 881    male  56
## 882    male  57
## 883  female  64
## 884  female  53
## 885  female  63
## 886    male  47
## 887  female  56
## 888  female  57
## 889  female  62
## 890    male  68
## 891  female  50
## 892    male  46
## 893  female  58
## 894  female  48
## 895  female  57
## 896  female  61
## 897  female  51
## 898    male  58
## 899  female  52
## 900    male  68
## 901  female  45
## 902  female  61
## 903    male  41
## 904  female  58
## 905  female  62
## 906  female  58
## 907  female  53
## 908    male  41
## 909  female  58
## 910    male  55
## 911  female  48
## 912  female  59
## 913  female  63
## 914  female  53
## 915    male  44
## 916  female  57
## 917  female  64
## 918  female  60
## 919  female  67
## 920    male  60
## 921    male  46
## 922  female  36
## 923  female  50
## 924  female  59
## 925  female  39
## 926    male  55
## 927    male  53
## 928  female  54
## 929  female  45
## 930    male  51
## 931  female  62
## 932  female  56
## 933  female  63
## 934  female  63
## 935  female  53
## 936  female  61
## 937    male  49
## 938  female  47
## 939    male  63
## 940    male  55
## 941    male  56
## 942  female  53
## 943    male  57
## 944  female  50
## 945  female  61
## 946  female  61
## 947  female  61
## 948  female  55
## 949  female  62
## 950  female  50
## 951  female  41
## 952  female  35
## 953    male  64
## 954  female  38
## 955  female  36
## 956    male  46
## 957  female  48
## 958  female  36
## 959    male  43
## 960    male  63
## 961  female  36
## 962  female  49
## 963  female  44
## 964    male  39
## 965    male  46
## 966    male  61
## 967  female  51
## 968  female  58
## 969  female  48
## 970  female  64
## 971    male  51
## 972    male  54
## 973  female  41
## 974    male  59
## 975  female  34
## 976  female  42
## 977    male  60
## 978  female  37
## 979    male  38
## 980  female  46
## 981    male  40
## 982  female  52
## 983  female  45
## 984    male  58
## 985  female  48
## 986  female  44
## 987    male  36
## 988    male  37
## 989  female  55
## 990  female  63
## 991  female  59
## 992  female  58
## 993    male  41
## 994  female  59
## 995    male  44
## 996  female  49
## 997  female  42
## 998  female  42
## 999    male  41
## 1000 female  36
## 1001 female  40
## 1002 female  62
## 1003 female  50
## 1004   male  44
## 1005 female  61
## 1006 female  58
## 1007   male  42
## 1008   male  60
## 1009 female  57
## 1010   male  61
## 1011 female  51
## 1012 female  43
## 1013 female  39
## 1014   male  65
## 1015 female  65
## 1016 female  42
## 1017 female  58
## 1018 female  64
## 1019   male  62
## 1020   male  43
## 1021   male  60
## 1022 female  62
## 1023   male  47
## 1024   male  62
## 1025 female  39
## 1026 female  48
## 1027 female  61
## 1028 female  54
## 1029 female  55
## 1030 female  58
## 1031 female  38
## 1032 female  60
## 1033   male  40
## 1034   male  50
## 1035 female  61
## 1036 female  67
## 1037 female  37
## 1038 female  40
## 1039   male  64
## 1040 female  53
## 1041 female  56
## 1042   male  64
## 1043 female  49
## 1044 female  40
## 1045   male  62
## 1046 female  39
## 1047 female  63
## 1048   male  51
## 1049   male  56
## 1050 female  34
## 1051   male  66
## 1052   male  39
## 1053   male  59
## 1054 female  45
## 1055 female  63
## 1056   male  36
## 1057 female  37
## 1058 female  43
## 1059 female  54
## 1060 female  47
## 1061   male  51
## 1062 female  63
## 1063   male  39
## 1064   male  61
## 1065 female  48
## 1066 female  60
## 1067 female  59
## 1068 female  39
## 1069 female  52
## 1070 female  60
## 1071 female  39
## 1072 female  58
## 1073   male  64
## 1074 female  65
## 1075 female  54
## 1076   male  68
## 1077 female  45
## 1078   male  60
## 1079 female  55
## 1080 female  60
## 1081 female  42
## 1082 female  59
## 1083 female  51
## 1084 female  39
## 1085 female  55
## 1086 female  57
## 1087 female  46
## 1088 female  59
## 1089 female  56
## 1090 female  41
## 1091 female  49
## 1092 female  38
## 1093 female  50
## 1094   male  64
## 1095   male  39
## 1096 female  41
## 1097 female  40
## 1098   male  45
## 1099 female  43
## 1100 female  40
## 1101 female  50
## 1102 female  53
## 1103   male  37
## 1104 female  61
## 1105 female  50
## 1106 female  63
## 1107 female  61
## 1108 female  44
## 1109   male  45
## 1110 female  43
## 1111 female  55
## 1112   male  42
## 1113 female  41
## 1114 female  37
## 1115 female  44
## 1116   male  61
## 1117   male  67
## 1118 female  48
## 1119 female  64
## 1120   male  49
## 1121 female  47
## 1122   male  38
## 1123 female  57
## 1124   male  62
## 1125   male  36
## 1126 female  40
## 1127 female  46
## 1128 female  34
## 1129   male  54
## 1130 female  48
## 1131   male  61
## 1132 female  54
## 1133   male  60
## 1134 female  50
## 1135 female  56
## 1136 female  59
## 1137 female  45
## 1138   male  57
## 1139 female  37
## 1140 female  45
## 1141 female  35
## 1142   male  58
## 1143   male  43
## 1144   male  58
## 1145 female  66
## 1146 female  40
## 1147 female  44
## 1148   male  57
## 1149 female  60
## 1150   male  52
## 1151   male  38
## 1152   male  45
## 1153   male  50
## 1154 female  52
## 1155   male  39
## 1156   male  42
## 1157 female  51
## 1158   male  43
## 1159 female  48
## 1160   male  64
## 1161   male  59
## 1162   male  60
## 1163 female  64
## 1164 female  40
## 1165 female  43
## 1166   male  50
## 1167   male  52
## 1168 female  49
## 1169 female  65
## 1170   male  55
## 1171 female  63
## 1172 female  40
## 1173 female  40
## 1174 female  38
## 1175 female  58
## 1176 female  61
## 1177 female  39
## 1178 female  59
## 1179   male  55
## 1180   male  52
## 1181 female  63
## 1182   male  59
## 1183   male  63
## 1184 female  60
## 1185 female  39
## 1186 female  54
## 1187 female  47
## 1188   male  57
## 1189 female  44
## 1190   male  55
## 1191 female  44
## 1192 female  39
## 1193   male  45
## 1194 female  49
## 1195 female  36
## 1196   male  62
## 1197   male  44
## 1198   male  36
## 1199   male  61
## 1200   male  64
## 1201 female  40
## 1202 female  42
## 1203   male  57
## 1204 female  62
## 1205 female  63
## 1206   male  50
## 1207 female  55
## 1208 female  40
## 1209 female  64
## 1210   male  66
## 1211 female  64
## 1212 female  45
## 1213   male  44
## 1214 female  45
## 1215   male  60
## 1216   male  45
## 1217   male  36
## 1218 female  55
## 1219   male  55
## 1220   male  42
## 1221 female  47
## 1222   male  47
## 1223   male  50
## 1224   male  48
## 1225   male  40
## 1226   male  62
## 1227 female  46
## 1228 female  55
## 1229 female  55
## 1230   male  47
## 1231 female  60
## 1232 female  44
## 1233 female  60
## 1234   male  60
## 1235 female  42
## 1236 female  56
## 1237 female  52
## 1238   male  40
## 1239 female  51
## 1240 female  46
## 1241 female  64
## 1242 female  42
## 1243 female  64
## 1244 female  52
## 1245   male  38
## 1246   male  52
## 1247 female  63
## 1248   male  48
## 1249   male  60
## 1250   male  52
## 1251 female  64
## 1252 female  44
## 1253 female  39
## 1254 female  43
## 1255   male  49
## 1256 female  64
## 1257 female  44
## 1258   male  54
## 1259 female  59
## 1260 female  49
## 1261 female  42
## 1262   male  49
## 1263   male  53
## 1264   male  59
## 1265   male  38
## 1266 female  37
## 1267   male  55
## 1268 female  63
## 1269   male  60
## 1270 female  55
## 1271 female  64
## 1272   male  38
## 1273   male  41
## 1274   male  36
## 1275   male  56
## 1276   male  37
## 1277 female  44
## 1278 female  47
## 1279 female  38
## 1280   male  58
## 1281   male  66
## 1282 female  66
## 1283   male  56
## 1284 female  61
## 1285 female  59
## 1286 female  41
## 1287 female  41
## 1288 female  42
## 1289 female  52
## 1290   male  58
## 1291   male  46
## 1292 female  55
## 1293 female  43
## 1294 female  39
## 1295 female  59
## 1296 female  39
## 1297   male  57
## 1298   male  44
## 1299   male  45
## 1300 female  64
## 1301 female  61
## 1302   male  52
## 1303 female  48
## 1304 female  41
## 1305 female  39
## 1306 female  47
## 1307 female  46
## 1308 female  38
## 1309   male  65
## 1310 female  56
## 1311 female  48
## 1312   male  53
## 1313   male  44
## 1314   male  36
## 1315   male  42
## 1316 female  44
## 1317 female  50
## 1318 female  66
## 1319   male  49
## 1320   male  43
## 1321   male  39
## 1322 female  50
## 1323   male  44
## 1324 female  59
## 1325 female  40
## 1326 female  59
## 1327 female  63
## 1328   male  56
## 1329 female  62
## 1330   male  58
## 1331   male  45
## 1332   male  38
## 1333 female  63
## 1334   male  41
## 1335 female  41
## 1336 female  47
## 1337 female  62
## 1338 female  43
## 1339   male  39
## 1340 female  58
## 1341 female  40
## 1342 female  59
## 1343 female  46
## 1344   male  57
## 1345 female  42
## 1346   male  48
## 1347 female  48
## 1348 female  64
## 1349 female  51
## 1350 female  50
## 1351   male  45
## 1352 female  42
## 1353 female  57
## 1354 female  51
## 1355 female  42
## 1356 female  58
## 1357 female  62
## 1358   male  63
## 1359 female  51
## 1360 female  50
## 1361   male  57
## 1362 female  49
## 1363   male  63
## 1364 female  60
## 1365 female  48
## 1366   male  48
## 1367 female  69
## 1368 female  59
## 1369 female  41
## 1370 female  38
## 1371 female  48
## 1372 female  65
## 1373 female  64
## 1374 female  66
## 1375 female  64
## 1376 female  55
## 1377 female  66
## 1378 female  58
## 1379 female  66
## 1380   male  50
## 1381 female  56
## 1382 female  52
## 1383 female  39
## 1384 female  40
## 1385 female  60
## 1386 female  54
## 1387 female  62
## 1388   male  57
## 1389 female  58
## 1390   male  63
## 1391 female  54
## 1392   male  51
## 1393   male  52
## 1394 female  52
## 1395   male  44
## 1396 female  62
## 1397 female  56
## 1398 female  47
## 1399 female  57
## 1400 female  53
## 1401   male  48
## 1402   male  67
## 1403   male  50
## 1404   male  44
## 1405 female  63
## 1406   male  52
## 1407 female  47
## 1408 female  52
## 1409 female  66
## 1410   male  56
## 1411   male  36
## 1412   male  41
## 1413 female  54
## 1414 female  53
## 1415 female  61
## 1416   male  59
## 1417 female  62
## 1418 female  58
## 1419   male  42
## 1420 female  63
## 1421 female  65
## 1422 female  50
## 1423 female  51
## 1424   male  48
## 1425 female  55
## 1426 female  62
## 1427 female  50
## 1428   male  48
## 1429 female  46
## 1430   male  54
## 1431   male  49
## 1432   male  58
## 1433   male  63
## 1434   male  45
## 1435   male  40
## 1436   male  56
## 1437   male  40
## 1438   male  56
## 1439   male  55
## 1440   male  54
## 1441   male  36
## 1442 female  67
## 1443 female  46
## 1444 female  56
## 1445   male  58
## 1446 female  65
## 1447 female  54
## 1448   male  47
## 1449   male  49
## 1450 female  65
## 1451 female  56
## 1452 female  59
## 1453   male  52
## 1454   male  41
## 1455   male  63
## 1456 female  64
## 1457 female  60
## 1458 female  45
## 1459 female  57
## 1460   male  67
## 1461 female  64
## 1462 female  55
## 1463   male  64
## 1464 female  58
## 1465 female  43
## 1466   male  39
## 1467 female  60
## 1468   male  44
## 1469   male  61
## 1470 female  62
## 1471 female  54
## 1472 female  55
## 1473   male  60
## 1474   male  53
## 1475 female  65
## 1476   male  63
## 1477 female  62
## 1478   male  59
## 1479 female  52
## 1480 female  47
## 1481 female  53
## 1482 female  39
## 1483 female  49
## 1484   male  61
## 1485 female  48
## 1486 female  49
## 1487   male  62
## 1488 female  49
## 1489   male  68
## 1490 female  63
## 1491 female  61
## 1492 female  49
## 1493 female  38
## 1494 female  42
## 1495   male  57
## 1496 female  39
## 1497   male  42
## 1498   male  40
## 1499 female  59
## 1500 female  42
## 1501 female  61
## 1502 female  35
## 1503 female  61
## 1504   male  54
## 1505 female  46
## 1506 female  39
## 1507 female  40
## 1508   male  58
## 1509 female  62
## 1510   male  53
## 1511   male  38
## 1512   male  45
## 1513   male  36
## 1514 female  44
## 1515   male  51
## 1516 female  65
## 1517 female  49
## 1518 female  60
## 1519 female  42
## 1520   male  43
## 1521 female  56
## 1522 female  55
## 1523   male  56
## 1524   male  43
## 1525   male  60
## 1526 female  62
## 1527 female  53
## 1528 female  47
## 1529 female  39
## 1530 female  40
## 1531 female  51
## 1532 female  43
## 1533   male  55
## 1534   male  55
## 1535   male  46
## 1536 female  48
## 1537   male  36
## 1538   male  42
## 1539 female  60
## 1540   male  52
## 1541 female  65
## 1542 female  36
## 1543   male  57
## 1544   male  54
## 1545   male  46
## 1546 female  65
## 1547 female  56
## 1548 female  53
## 1549   male  52
## 1550 female  51
## 1551 female  56
## 1552 female  54
## 1553   male  50
## 1554 female  53
## 1555 female  50
## 1556 female  53
## 1557   male  49
## 1558 female  61
## 1559   male  61
## 1560   male  66
## 1561   male  43
## 1562   male  41
## 1563 female  36
## 1564   male  64
## 1565   male  59
## 1566 female  41
## 1567 female  43
## 1568 female  46
## 1569 female  58
## 1570 female  33
## 1571   male  61
## 1572 female  49
## 1573 female  64
## 1574   male  52
## 1575 female  37
## 1576   male  41
## 1577 female  63
## 1578   male  66
## 1579 female  34
## 1580 female  38
## 1581 female  40
## 1582 female  53
## 1583   male  53
## 1584 female  42
## 1585 female  56
## 1586 female  47
## 1587 female  48
## 1588 female  36
## 1589   male  52
## 1590 female  42
## 1591 female  51
## 1592 female  54
## 1593   male  61
## 1594 female  68
## 1595 female  41
## 1596 female  56
## 1597 female  61
## 1598 female  57
## 1599 female  52
## 1600 female  55
## 1601 female  45
## 1602 female  51
## 1603 female  62
## 1604   male  66
## 1605 female  59
## 1606 female  39
## 1607   male  59
## 1608 female  46
## 1609   male  58
## 1610 female  39
## 1611 female  41
## 1612 female  40
## 1613 female  42
## 1614 female  63
## 1615 female  63
## 1616 female  63
## 1617   male  65
## 1618 female  64
## 1619 female  64
## 1620 female  37
## 1621 female  40
## 1622 female  43
## 1623 female  40
## 1624 female  55
## 1625   male  56
## 1626   male  57
## 1627   male  63
## 1628 female  50
## 1629 female  52
## 1630   male  38
## 1631 female  63
## 1632 female  46
## 1633 female  61
## 1634   male  38
## 1635 female  65
## 1636 female  58
## 1637 female  63
## 1638 female  38
## 1639 female  49
## 1640 female  47
## 1641   male  51
## 1642 female  62
## 1643   male  56
## 1644   male  33
## 1645 female  37
## 1646 female  42
## 1647   male  54
## 1648 female  42
## 1649 female  62
## 1650   male  50
## 1651 female  62
## 1652   male  61
## 1653 female  65
## 1654 female  62
## 1655   male  52
## 1656 female  37
## 1657 female  42
## 1658   male  48
## 1659 female  60
## 1660   male  39
## 1661   male  64
## 1662 female  65
## 1663 female  59
## 1664   male  64
## 1665   male  41
## 1666 female  46
## 1667   male  48
## 1668 female  46
## 1669 female  45
## 1670 female  61
## 1671   male  43
## 1672 female  42
## 1673   male  40
## 1674 female  42
## 1675 female  49
## 1676   male  60
## 1677 female  59
## 1678 female  43
## 1679   male  45
## 1680   male  52
## 1681 female  54
## 1682   male  35
## 1683   male  41
## 1684   male  44
## 1685 female  51
## 1686 female  60
## 1687 female  46
## 1688   male  56
## 1689 female  52
## 1690   male  44
## 1691 female  39
## 1692   male  62
## 1693   male  47
## 1694   male  57
## 1695 female  62
## 1696 female  43
## 1697 female  40
## 1698 female  56
## 1699   male  63
## 1700 female  48
## 1701 female  63
## 1702 female  52
## 1703   male  68
## 1704 female  47
## 1705 female  46
## 1706 female  41
## 1707   male  53
## 1708 female  52
## 1709   male  62
## 1710 female  41
## 1711 female  53
## 1712 female  43
## 1713   male  37
## 1714 female  45
## 1715   male  58
## 1716   male  44
## 1717 female  60
## 1718   male  59
## 1719 female  44
## 1720 female  56
## 1721   male  41
## 1722 female  39
## 1723   male  40
## 1724 female  55
## 1725   male  54
## 1726   male  53
## 1727   male  56
## 1728 female  63
## 1729   male  40
## 1730   male  67
## 1731   male  56
## 1732 female  62
## 1733   male  67
## 1734 female  57
## 1735   male  48
## 1736   male  58
## 1737 female  57
## 1738 female  67
## 1739 female  41
## 1740 female  54
## 1741 female  46
## 1742   male  52
## 1743   male  56
## 1744   male  61
## 1745 female  54
## 1746 female  41
## 1747 female  45
## 1748 female  57
## 1749 female  49
## 1750 female  58
## 1751   male  37
## 1752 female  53
## 1753 female  51
## 1754   male  60
## 1755   male  40
## 1756   male  52
## 1757   male  38
## 1758 female  61
## 1759 female  52
## 1760 female  63
## 1761 female  42
## 1762   male  64
## 1763 female  46
## 1764 female  55
## 1765 female  56
## 1766 female  60
## 1767 female  40
## 1768 female  50
## 1769   male  50
## 1770 female  56
## 1771 female  51
## 1772   male  47
## 1773 female  63
## 1774 female  46
## 1775 female  51
## 1776   male  50
## 1777 female  55
## 1778 female  61
## 1779   male  38
## 1780 female  49
## 1781 female  48
## 1782   male  41
## 1783 female  60
## 1784 female  54
## 1785   male  45
## 1786 female  62
## 1787   male  38
## 1788 female  39
## 1789 female  65
## 1790   male  52
## 1791 female  44
## 1792   male  41
## 1793   male  54
## 1794 female  52
## 1795 female  53
## 1796 female  52
## 1797 female  49
## 1798 female  42
## 1799 female  55
## 1800 female  54
## 1801 female  54
## 1802   male  57
## 1803 female  39
## 1804 female  64
## 1805   male  62
## 1806 female  57
## 1807 female  60
## 1808 female  51
## 1809 female  56
## 1810 female  39
## 1811   male  55
## 1812 female  64
## 1813 female  52
## 1814 female  39
## 1815 female  56
## 1816 female  56
## 1817 female  52
## 1818   male  36
## 1819   male  44
## 1820 female  53
## 1821 female  42
## 1822 female  43
## 1823   male  44
## 1824   male  50
## 1825 female  54
## 1826 female  49
## 1827   male  63
## 1828   male  47
## 1829 female  60
## 1830 female  62
## 1831 female  42
## 1832   male  42
## 1833   male  46
## 1834   male  39
## 1835 female  63
## 1836   male  48
## 1837 female  53
## 1838   male  66
## 1839   male  59
## 1840   male  62
## 1841 female  65
## 1842 female  66
## 1843 female  49
## 1844 female  56
## 1845   male  47
## 1846 female  57
## 1847   male  61
## 1848 female  58
## 1849   male  35
## 1850 female  55
## 1851 female  54
## 1852   male  33
## 1853 female  52
## 1854 female  66
## 1855   male  38
## 1856 female  49
## 1857   male  50
## 1858 female  40
## 1859   male  62
## 1860 female  51
## 1861   male  42
## 1862 female  66
## 1863   male  43
## 1864 female  53
## 1865 female  62
## 1866   male  60
## 1867   male  46
## 1868 female  53
## 1869 female  58
## 1870 female  59
## 1871 female  46
## 1872 female  54
## 1873   male  43
## 1874   male  48
## 1875 female  62
## 1876   male  50
## 1877 female  63
## 1878 female  51
## 1879   male  52
## 1880 female  47
## 1881 female  46
## 1882   male  39
## 1883 female  52
## 1884 female  49
## 1885 female  59
## 1886   male  51
## 1887 female  48
## 1888   male  53
## 1889 female  54
## 1890   male  59
## 1891   male  48
## 1892 female  39
## 1893 female  40
## 1894 female  45
## 1895 female  44
## 1896 female  46
## 1897 female  59
## 1898 female  46
## 1899 female  57
## 1900   male  56
## 1901 female  59
## 1902   male  44
## 1903 female  53
## 1904   male  52
## 1905 female  61
## 1906 female  48
## 1907 female  64
## 1908   male  44
## 1909 female  48
## 1910 female  52
## 1911   male  67
## 1912   male  40
## 1913 female  49
## 1914   male  62
## 1915 female  54
## 1916 female  54
## 1917 female  61
## 1918 female  51
## 1919   male  48
## 1920 female  53
## 1921 female  45
## 1922 female  53
## 1923 female  54
## 1924 female  60
## 1925 female  52
## 1926 female  57
## 1927 female  44
## 1928 female  47
## 1929 female  41
## 1930 female  38
## 1931 female  56
## 1932 female  46
## 1933   male  62
## 1934   male  36
## 1935 female  62
## 1936 female  43
## 1937 female  39
## 1938 female  48
## 1939   male  54
## 1940 female  55
## 1941 female  47
## 1942 female  55
## 1943   male  56
## 1944 female  51
## 1945 female  45
## 1946 female  54
## 1947 female  54
## 1948   male  52
## 1949 female  46
## 1950 female  64
## 1951 female  47
## 1952   male  43
## 1953 female  49
## 1954 female  48
## 1955   male  57
## 1956 female  55
## 1957   male  59
## 1958   male  51
## 1959   male  44
## 1960 female  45
## 1961   male  48
## 1962   male  51
## 1963   male  57
## 1964 female  61
## 1965 female  56
## 1966   male  59
## 1967 female  46
## 1968   male  40
## 1969   male  35
## 1970   male  48
## 1971   male  46
## 1972 female  43
## 1973 female  52
## 1974 female  47
## 1975 female  51
## 1976 female  60
## 1977   male  49
## 1978 female  56
## 1979 female  58
## 1980 female  54
## 1981 female  64
## 1982 female  65
## 1983 female  44
## 1984 female  43
## 1985 female  55
## 1986 female  38
## 1987 female  38
## 1988 female  67
## 1989   male  45
## 1990 female  43
## 1991 female  40
## 1992 female  56
## 1993   male  39
## 1994 female  56
## 1995   male  50
## 1996 female  61
## 1997 female  39
## 1998 female  57
## 1999 female  46
## 2000   male  57
## 2001 female  35
## 2002   male  47
## 2003   male  59
## 2004   male  69
## 2005   male  58
## 2006 female  52
## 2007 female  48
## 2008 female  57
## 2009 female  58
## 2010 female  57
## 2011 female  66
## 2012 female  42
## 2013 female  47
## 2014 female  37
## 2015 female  40
## 2016 female  61
## 2017 female  45
## 2018   male  57
## 2019 female  41
## 2020 female  36
## 2021 female  55
## 2022 female  49
## 2023 female  49
## 2024 female  39
## 2025 female  41
## 2026 female  46
## 2027 female  44
## 2028 female  43
## 2029 female  59
## 2030 female  54
## 2031 female  43
## 2032 female  57
## 2033   male  56
## 2034   male  39
## 2035 female  53
## 2036   male  60
## 2037   male  52
## 2038 female  51
## 2039 female  48
## 2040 female  49
## 2041 female  49
## 2042 female  40
## 2043 female  57
## 2044   male  37
## 2045   male  48
## 2046 female  65
## 2047   male  40
## 2048   male  67
## 2049 female  41
## 2050 female  45
## 2051 female  42
## 2052 female  58
## 2053 female  43
## 2054   male  49
## 2055   male  52
## 2056 female  51
## 2057   male  50
## 2058 female  38
## 2059 female  47
## 2060 female  54
## 2061 female  45
## 2062 female  47
## 2063 female  64
## 2064   male  49
## 2065 female  41
## 2066 female  50
## 2067   male  50
## 2068   male  65
## 2069 female  56
## 2070 female  47
## 2071 female  45
## 2072 female  57
## 2073   male  62
## 2074 female  64
## 2075 female  41
## 2076   male  61
## 2077 female  45
## 2078 female  41
## 2079 female  38
## 2080   male  53
## 2081   male  47
## 2082   male  35
## 2083   male  58
## 2084 female  56
## 2085   male  66
## 2086 female  46
## 2087   male  49
## 2088 female  51
## 2089 female  56
## 2090   male  60
## 2091   male  57
## 2092 female  67
## 2093 female  60
## 2094   male  55
## 2095 female  63
## 2096 female  59
## 2097 female  40
## 2098 female  58
## 2099 female  41
## 2100 female  39
## 2101   male  40
## 2102 female  62
## 2103 female  42
## 2104 female  62
## 2105   male  45
## 2106 female  44
## 2107 female  60
## 2108   male  39
## 2109   male  54
## 2110 female  47
## 2111 female  45
## 2112 female  55
## 2113 female  37
## 2114 female  46
## 2115 female  36
## 2116 female  51
## 2117   male  41
## 2118 female  50
## 2119   male  60
## 2120 female  48
## 2121 female  47
## 2122 female  67
## 2123 female  49
## 2124   male  49
## 2125 female  63
## 2126   male  45
## 2127 female  48
## 2128   male  45
## 2129 female  56
## 2130   male  52
## 2131   male  59
## 2132   male  45
## 2133 female  52
## 2134 female  55
## 2135 female  35
## 2136 female  58
## 2137   male  55
## 2138 female  56
## 2139 female  48
## 2140 female  55
## 2141 female  53
## 2142 female  44
## 2143 female  46
## 2144 female  40
## 2145 female  40
## 2146 female  40
## 2147   male  54
## 2148 female  54
## 2149   male  58
## 2150 female  66
## 2151 female  65
## 2152   male  44
## 2153 female  59
## 2154 female  55
## 2155   male  44
## 2156 female  48
## 2157 female  47
## 2158 female  51
## 2159   male  69
## 2160 female  41
## 2161   male  44
## 2162 female  49
## 2163 female  35
## 2164 female  37
## 2165 female  36
## 2166 female  44
## 2167 female  41
## 2168 female  42
## 2169 female  38
## 2170 female  46
## 2171 female  40
## 2172 female  56
## 2173 female  44
## 2174 female  36
## 2175   male  60
## 2176 female  58
## 2177 female  43
## 2178 female  46
## 2179 female  53
## 2180   male  42
## 2181   male  46
## 2182   male  39
## 2183 female  60
## 2184   male  62
## 2185 female  41
## 2186 female  37
## 2187 female  36
## 2188 female  38
## 2189   male  48
## 2190 female  43
## 2191 female  38
## 2192 female  40
## 2193 female  54
## 2194 female  46
## 2195 female  46
## 2196 female  41
## 2197 female  45
## 2198 female  36
## 2199   male  55
## 2200 female  48
## 2201 female  38
## 2202 female  47
## 2203 female  51
## 2204   male  46
## 2205   male  42
## 2206 female  39
## 2207 female  48
## 2208 female  47
## 2209 female  44
## 2210   male  41
## 2211 female  44
## 2212 female  40
## 2213 female  39
## 2214 female  55
## 2215   male  38
## 2216   male  34
## 2217   male  49
## 2218 female  57
## 2219 female  43
## 2220 female  45
## 2221 female  46
## 2222 female  33
## 2223 female  35
## 2224 female  38
## 2225 female  44
## 2226 female  38
## 2227 female  56
## 2228   male  45
## 2229 female  41
## 2230 female  67
## 2231 female  58
## 2232 female  40
## 2233   male  45
## 2234 female  54
## 2235 female  36
## 2236 female  45
## 2237   male  45
## 2238 female  52
## 2239   male  58
## 2240   male  42
## 2241 female  56
## 2242 female  56
## 2243   male  61
## 2244 female  48
## 2245 female  49
## 2246   male  40
## 2247 female  60
## 2248 female  62
## 2249   male  49
## 2250   male  61
## 2251 female  39
## 2252 female  35
## 2253 female  54
## 2254 female  51
## 2255 female  38
## 2256 female  55
## 2257 female  59
## 2258   male  64
## 2259   male  45
## 2260   male  40
## 2261 female  51
## 2262 female  54
## 2263   male  53
## 2264   male  55
## 2265   male  42
## 2266   male  51
## 2267 female  44
## 2268 female  64
## 2269 female  46
## 2270 female  51
## 2271   male  63
## 2272   male  65
## 2273 female  42
## 2274   male  54
## 2275   male  65
## 2276   male  51
## 2277   male  55
## 2278   male  62
## 2279 female  62
## 2280   male  45
## 2281   male  63
## 2282 female  56
## 2283 female  61
## 2284 female  47
## 2285 female  41
## 2286 female  46
## 2287 female  38
## 2288   male  62
## 2289 female  42
## 2290   male  38
## 2291 female  45
## 2292   male  51
## 2293   male  61
## 2294 female  46
## 2295   male  56
## 2296 female  64
## 2297 female  44
## 2298 female  43
## 2299 female  39
## 2300 female  49
## 2301 female  50
## 2302   male  64
## 2303 female  45
## 2304   male  40
## 2305 female  40
## 2306 female  45
## 2307   male  67
## 2308 female  50
## 2309   male  54
## 2310   male  37
## 2311 female  40
## 2312 female  54
## 2313 female  66
## 2314 female  41
## 2315 female  37
## 2316 female  59
## 2317   male  38
## 2318 female  42
## 2319   male  38
## 2320 female  55
## 2321 female  45
## 2322 female  49
## 2323 female  65
## 2324   male  48
## 2325   male  50
## 2326 female  47
## 2327 female  48
## 2328   male  46
## 2329 female  53
## 2330 female  37
## 2331 female  44
## 2332 female  39
## 2333 female  38
## 2334 female  39
## 2335   male  42
## 2336 female  47
## 2337 female  44
## 2338   male  59
## 2339 female  47
## 2340 female  61
## 2341   male  60
## 2342 female  47
## 2343 female  40
## 2344 female  50
## 2345 female  42
## 2346 female  51
## 2347   male  46
## 2348 female  48
## 2349 female  41
## 2350 female  50
## 2351   male  49
## 2352 female  48
## 2353   male  46
## 2354   male  60
## 2355 female  44
## 2356   male  58
## 2357   male  58
## 2358 female  52
## 2359 female  38
## 2360 female  49
## 2361 female  38
## 2362 female  48
## 2363 female  43
## 2364 female  59
## 2365 female  42
## 2366 female  44
## 2367   male  39
## 2368 female  40
## 2369 female  39
## 2370 female  39
## 2371 female  40
## 2372   male  63
## 2373   male  45
## 2374 female  53
## 2375   male  50
## 2376   male  46
## 2377 female  57
## 2378 female  55
## 2379   male  51
## 2380   male  57
## 2381   male  47
## 2382 female  55
## 2383 female  40
## 2384 female  42
## 2385   male  47
## 2386 female  41
## 2387   male  51
## 2388 female  49
## 2389 female  48
## 2390 female  64
## 2391 female  40
## 2392 female  52
## 2393 female  42
## 2394   male  38
## 2395 female  64
## 2396 female  42
## 2397   male  49
## 2398   male  63
## 2399   male  44
## 2400 female  36
## 2401 female  53
## 2402   male  43
## 2403   male  51
## 2404 female  61
## 2405 female  55
## 2406 female  42
## 2407   male  45
## 2408   male  53
## 2409   male  39
## 2410 female  56
## 2411 female  38
## 2412 female  55
## 2413 female  45
## 2414 female  49
## 2415 female  44
## 2416 female  46
## 2417 female  38
## 2418 female  47
## 2419 female  41
## 2420   male  51
## 2421 female  38
## 2422 female  57
## 2423 female  44
## 2424   male  48
## 2425 female  36
## 2426 female  43
## 2427 female  44
## 2428 female  43
## 2429 female  48
## 2430 female  49
## 2431 female  50
## 2432   male  60
## 2433   male  45
## 2434 female  54
## 2435 female  57
## 2436   male  35
## 2437   male  62
## 2438   male  67
## 2439 female  58
## 2440   male  63
## 2441   male  63
## 2442   male  51
## 2443 female  63
## 2444 female  38
## 2445 female  43
## 2446   male  44
## 2447 female  51
## 2448   male  39
## 2449   male  42
## 2450 female  47
## 2451   male  62
## 2452 female  40
## 2453 female  46
## 2454 female  40
## 2455   male  68
## 2456 female  54
## 2457 female  55
## 2458 female  46
## 2459 female  47
## 2460   male  34
## 2461   male  51
## 2462   male  45
## 2463 female  53
## 2464   male  38
## 2465 female  56
## 2466 female  45
## 2467   male  60
## 2468   male  50
## 2469 female  44
## 2470   male  38
## 2471   male  63
## 2472 female  38
## 2473 female  39
## 2474 female  39
## 2475 female  44
## 2476 female  38
## 2477 female  34
## 2478   male  43
## 2479 female  38
## 2480   male  61
## 2481   male  68
## 2482   male  44
## 2483   male  40
## 2484 female  34
## 2485 female  42
## 2486 female  42
## 2487 female  46
## 2488   male  41
## 2489 female  44
## 2490 female  48
## 2491   male  40
## 2492   male  42
## 2493   male  53
## 2494   male  46
## 2495   male  50
## 2496 female  59
## 2497 female  45
## 2498   male  49
## 2499 female  53
## 2500   male  40
## 2501 female  44
## 2502 female  43
## 2503   male  48
## 2504 female  50
## 2505 female  40
## 2506 female  46
## 2507   male  60
## 2508 female  60
## 2509 female  48
## 2510 female  52
## 2511 female  37
## 2512   male  41
## 2513   male  45
## 2514 female  50
## 2515 female  40
## 2516 female  46
## 2517 female  40
## 2518   male  38
## 2519 female  42
## 2520 female  55
## 2521 female  42
## 2522 female  52
## 2523 female  48
## 2524 female  45
## 2525 female  58
## 2526 female  46
## 2527 female  42
## 2528 female  41
## 2529   male  57
## 2530 female  51
## 2531 female  37
## 2532   male  57
## 2533 female  51
## 2534   male  47
## 2535 female  53
## 2536 female  63
## 2537 female  50
## 2538 female  42
## 2539 female  37
## 2540 female  38
## 2541   male  48
## 2542 female  44
## 2543 female  61
## 2544 female  47
## 2545   male  53
## 2546 female  49
## 2547 female  41
## 2548   male  59
## 2549 female  39
## 2550 female  38
## 2551   male  60
## 2552   male  51
## 2553   male  61
## 2554   male  58
## 2555 female  41
## 2556   male  68
## 2557 female  45
## 2558   male  58
## 2559 female  41
## 2560 female  46
## 2561 female  36
## 2562   male  43
## 2563 female  41
## 2564 female  39
## 2565 female  42
## 2566   male  47
## 2567   male  57
## 2568   male  49
## 2569 female  41
## 2570   male  38
## 2571   male  43
## 2572 female  50
## 2573   male  64
## 2574   male  39
## 2575   male  51
## 2576 female  61
## 2577 female  38
## 2578 female  46
## 2579   male  54
## 2580 female  41
## 2581   male  41
## 2582   male  61
## 2583   male  52
## 2584 female  50
## 2585 female  39
## 2586   male  41
## 2587   male  48
## 2588   male  39
## 2589 female  41
## 2590   male  61
## 2591 female  49
## 2592   male  51
## 2593 female  50
## 2594   male  45
## 2595   male  48
## 2596 female  51
## 2597 female  40
## 2598   male  53
## 2599 female  46
## 2600 female  40
## 2601   male  42
## 2602   male  46
## 2603   male  49
## 2604   male  67
## 2605   male  63
## 2606 female  50
## 2607 female  59
## 2608   male  44
## 2609 female  39
## 2610   male  61
## 2611   male  42
## 2612 female  38
## 2613 female  43
## 2614 female  36
## 2615   male  39
## 2616 female  39
## 2617 female  36
## 2618 female  45
## 2619 female  50
## 2620   male  55
## 2621 female  41
## 2622 female  43
## 2623 female  36
## 2624   male  41
## 2625 female  37
## 2626 female  40
## 2627   male  49
## 2628 female  53
## 2629   male  52
## 2630   male  57
## 2631 female  33
## 2632 female  43
## 2633 female  42
## 2634   male  64
## 2635   male  51
## 2636   male  40
## 2637   male  41
## 2638 female  64
## 2639 female  48
## 2640   male  43
## 2641 female  67
## 2642 female  44
## 2643   male  51
## 2644   male  52
## 2645   male  47
## 2646   male  47
## 2647   male  47
## 2648 female  38
## 2649 female  46
## 2650 female  60
## 2651 female  39
## 2652   male  43
## 2653 female  59
## 2654   male  55
## 2655   male  38
## 2656 female  39
## 2657 female  46
## 2658 female  41
## 2659   male  38
## 2660   male  38
## 2661   male  59
## 2662   male  56
## 2663 female  44
## 2664 female  42
## 2665 female  35
## 2666 female  42
## 2667 female  47
## 2668   male  56
## 2669 female  43
## 2670 female  45
## 2671 female  32
## 2672 female  46
## 2673   male  45
## 2674 female  57
## 2675   male  52
## 2676 female  49
## 2677 female  56
## 2678 female  52
## 2679   male  44
## 2680   male  36
## 2681 female  36
## 2682 female  58
## 2683 female  43
## 2684 female  60
## 2685 female  40
## 2686   male  51
## 2687 female  43
## 2688   male  40
## 2689   male  35
## 2690 female  40
## 2691   male  44
## 2692   male  43
## 2693   male  58
## 2694   male  65
## 2695   male  67
## 2696   male  46
## 2697   male  53
## 2698 female  61
## 2699   male  40
## 2700   male  50
## 2701 female  64
## 2702 female  54
## 2703   male  52
## 2704   male  35
## 2705   male  36
## 2706 female  43
## 2707 female  45
## 2708   male  53
## 2709   male  36
## 2710 female  46
## 2711 female  40
## 2712   male  52
## 2713   male  48
## 2714   male  46
## 2715 female  42
## 2716   male  35
## 2717 female  51
## 2718 female  38
## 2719   male  68
## 2720 female  36
## 2721 female  57
## 2722 female  37
## 2723   male  38
## 2724 female  37
## 2725 female  56
## 2726 female  47
## 2727   male  38
## 2728 female  41
## 2729 female  42
## 2730   male  40
## 2731   male  39
## 2732   male  37
## 2733 female  43
## 2734 female  51
## 2735 female  40
## 2736   male  44
## 2737   male  61
## 2738   male  59
## 2739   male  47
## 2740   male  50
## 2741   male  37
## 2742 female  51
## 2743   male  41
## 2744   male  48
## 2745   male  44
## 2746 female  51
## 2747 female  43
## 2748 female  43
## 2749   male  39
## 2750   male  58
## 2751 female  59
## 2752 female  45
## 2753 female  51
## 2754   male  37
## 2755   male  50
## 2756 female  41
## 2757   male  48
## 2758 female  48
## 2759   male  49
## 2760 female  41
## 2761 female  54
## 2762   male  56
## 2763   male  48
## 2764   male  66
## 2765   male  45
## 2766 female  46
## 2767 female  42
## 2768 female  56
## 2769   male  47
## 2770 female  51
## 2771   male  60
## 2772   male  40
## 2773 female  34
## 2774 female  47
## 2775 female  46
## 2776 female  46
## 2777   male  39
## 2778 female  57
## 2779 female  43
## 2780 female  56
## 2781 female  39
## 2782   male  50
## 2783 female  46
## 2784 female  51
## 2785   male  57
## 2786 female  40
## 2787 female  48
## 2788   male  52
## 2789   male  39
## 2790   male  47
## 2791 female  43
## 2792 female  46
## 2793   male  52
## 2794   male  45
## 2795 female  37
## 2796 female  64
## 2797 female  37
## 2798   male  45
## 2799 female  39
## 2800 female  45
## 2801   male  51
## 2802   male  48
## 2803 female  39
## 2804   male  51
## 2805 female  40
## 2806 female  45
## 2807 female  50
## 2808   male  60
## 2809 female  38
## 2810   male  45
## 2811 female  47
## 2812   male  45
## 2813   male  55
## 2814   male  46
## 2815   male  51
## 2816 female  40
## 2817 female  54
## 2818   male  36
## 2819 female  42
## 2820 female  46
## 2821   male  49
## 2822   male  45
## 2823   male  40
## 2824   male  49
## 2825   male  43
## 2826   male  42
## 2827 female  45
## 2828   male  50
## 2829   male  58
## 2830   male  40
## 2831   male  48
## 2832 female  63
## 2833 female  39
## 2834   male  45
## 2835   male  57
## 2836   male  52
## 2837   male  38
## 2838   male  42
## 2839   male  52
## 2840 female  51
## 2841   male  37
## 2842   male  63
## 2843   male  44
## 2844   male  42
## 2845 female  38
## 2846 female  52
## 2847   male  53
## 2848   male  62
## 2849   male  39
## 2850   male  41
## 2851 female  43
## 2852   male  37
## 2853   male  58
## 2854   male  38
## 2855   male  37
## 2856   male  40
## 2857   male  45
## 2858 female  55
## 2859 female  52
## 2860 female  49
## 2861 female  43
## 2862 female  45
## 2863   male  63
## 2864   male  59
## 2865   male  46
## 2866 female  49
## 2867 female  47
## 2868 female  40
## 2869   male  44
## 2870   male  45
## 2871 female  59
## 2872   male  37
## 2873   male  56
## 2874   male  50
## 2875   male  48
## 2876 female  48
## 2877   male  55
## 2878 female  46
## 2879   male  41
## 2880   male  62
## 2881 female  35
## 2882   male  54
## 2883   male  50
## 2884   male  62
## 2885 female  54
## 2886 female  41
## 2887 female  45
## 2888 female  44
## 2889   male  58
## 2890   male  49
## 2891   male  44
## 2892   male  66
## 2893   male  46
## 2894   male  44
## 2895   male  42
## 2896 female  43
## 2897   male  38
## 2898   male  50
## 2899   male  43
## 2900   male  43
## 2901   male  52
## 2902   male  47
## 2903   male  43
## 2904   male  38
## 2905 female  42
## 2906 female  43
## 2907   male  57
## 2908   male  61
## 2909   male  44
## 2910   male  43
## 2911 female  51
## 2912   male  44
## 2913   male  39
## 2914   male  60
## 2915 female  44
## 2916   male  54
## 2917   male  36
## 2918   male  53
## 2919   male  50
## 2920 female  51
## 2921   male  52
## 2922   male  45
## 2923   male  46
## 2924   male  65
## 2925   male  43
## 2926   male  60
## 2927   male  57
## 2928   male  50
## 2929   male  46
## 2930 female  50
## 2931 female  41
## 2932   male  44
## 2933   male  55
## 2934   male  49
## 2935   male  50
## 2936   male  55
## 2937   male  55
## 2938   male  56
## 2939   male  56
## 2940   male  66
## 2941   male  43
## 2942 female  56
## 2943   male  35
## 2944 female  52
## 2945   male  41
## 2946   male  37
## 2947   male  44
## 2948   male  60
## 2949   male  61
## 2950   male  50
## 2951 female  42
## 2952   male  34
## 2953   male  40
## 2954   male  35
## 2955   male  41
## 2956   male  42
## 2957   male  52
## 2958   male  36
## 2959   male  36
## 2960   male  41
## 2961   male  46
## 2962 female  57
## 2963   male  52
## 2964 female  44
## 2965   male  46
## 2966   male  53
## 2967   male  61
## 2968 female  37
## 2969 female  40
## 2970   male  50
## 2971 female  40
## 2972 female  46
## 2973   male  39
## 2974   male  52
## 2975 female  42
## 2976 female  39
## 2977 female  46
## 2978   male  50
## 2979 female  43
## 2980 female  43
## 2981   male  49
## 2982   male  38
## 2983 female  50
## 2984   male  59
## 2985 female  50
## 2986   male  51
## 2987 female  44
## 2988   male  43
## 2989   male  49
## 2990   male  48
## 2991 female  45
## 2992   male  51
## 2993 female  43
## 2994   male  48
## 2995   male  46
## 2996 female  40
## 2997   male  38
## 2998 female  36
## 2999 female  61
## 3000   male  40
## 3001   male  48
## 3002 female  60
## 3003 female  43
## 3004   male  46
## 3005   male  39
## 3006   male  40
## 3007   male  42
## 3008 female  58
## 3009   male  62
## 3010   male  40
## 3011   male  47
## 3012 female  38
## 3013   male  57
## 3014   male  35
## 3015   male  42
## 3016   male  48
## 3017   male  46
## 3018   male  48
## 3019 female  57
## 3020   male  42
## 3021   male  64
## 3022   male  53
## 3023 female  47
## 3024 female  58
## 3025   male  36
## 3026 female  44
## 3027   male  44
## 3028   male  44
## 3029 female  58
## 3030   male  59
## 3031 female  39
## 3032   male  57
## 3033 female  48
## 3034   male  43
## 3035 female  42
## 3036   male  64
## 3037   male  48
## 3038   male  45
## 3039   male  43
## 3040   male  47
## 3041 female  48
## 3042 female  47
## 3043   male  50
## 3044 female  60
## 3045   male  40
## 3046 female  61
## 3047   male  46
## 3048   male  43
## 3049   male  37
## 3050   male  59
## 3051   male  37
## 3052 female  43
## 3053 female  39
## 3054 female  47
## 3055   male  40
## 3056 female  39
## 3057   male  54
## 3058   male  43
## 3059   male  57
## 3060   male  43
## 3061 female  39
## 3062   male  40
## 3063 female  44
## 3064   male  42
## 3065 female  45
## 3066 female  42
## 3067 female  64
## 3068   male  38
## 3069   male  61
## 3070   male  41
## 3071 female  47
## 3072   male  62
## 3073   male  53
## 3074 female  51
## 3075   male  58
## 3076 female  53
## 3077   male  49
## 3078   male  48
## 3079 female  49
## 3080 female  43
## 3081   male  40
## 3082   male  41
## 3083   male  36
## 3084   male  54
## 3085   male  55
## 3086 female  59
## 3087 female  37
## 3088   male  47
## 3089   male  59
## 3090   male  63
## 3091   male  47
## 3092   male  53
## 3093   male  51
## 3094 female  57
## 3095   male  49
## 3096   male  58
## 3097   male  48
## 3098   male  38
## 3099   male  50
## 3100 female  57
## 3101   male  42
## 3102   male  60
## 3103   male  56
## 3104   male  38
## 3105 female  34
## 3106   male  42
## 3107 female  37
## 3108   male  46
## 3109 female  36
## 3110   male  62
## 3111   male  37
## 3112 female  40
## 3113   male  44
## 3114   male  36
## 3115   male  44
## 3116   male  39
## 3117   male  57
## 3118   male  52
## 3119 female  44
## 3120   male  58
## 3121   male  58
## 3122   male  55
## 3123   male  56
## 3124 female  59
## 3125   male  43
## 3126   male  43
## 3127 female  58
## 3128 female  56
## 3129 female  49
## 3130 female  42
## 3131 female  35
## 3132 female  43
## 3133 female  48
## 3134 female  40
## 3135 female  40
## 3136   male  49
## 3137   male  62
## 3138 female  42
## 3139 female  41
## 3140   male  40
## 3141 female  63
## 3142   male  42
## 3143 female  36
## 3144 female  66
## 3145   male  40
## 3146   male  44
## 3147   male  42
## 3148 female  37
## 3149   male  59
## 3150 female  43
## 3151   male  53
## 3152   male  40
## 3153 female  42
## 3154 female  39
## 3155   male  56
## 3156   male  56
## 3157   male  64
## 3158   male  58
## 3159   male  58
## 3160   male  42
## 3161   male  39
## 3162 female  50
## 3163 female  47
## 3164   male  63
## 3165   male  48
## 3166 female  39
## 3167 female  38
## 3168 female  38
## 3169   male  54
## 3170 female  50
## 3171   male  54
## 3172 female  37
## 3173   male  39
## 3174 female  46
## 3175   male  42
## 3176   male  51
## 3177   male  46
## 3178   male  37
## 3179 female  52
## 3180   male  51
## 3181   male  47
## 3182 female  40
## 3183   male  39
## 3184   male  53
## 3185 female  43
## 3186   male  40
## 3187   male  43
## 3188 female  48
## 3189   male  54
## 3190 female  46
## 3191 female  47
## 3192 female  46
## 3193   male  46
## 3194   male  41
## 3195   male  56
## 3196 female  46
## 3197   male  52
## 3198   male  38
## 3199 female  50
## 3200 female  51
## 3201   male  51
## 3202   male  40
## 3203 female  41
## 3204 female  49
## 3205 female  39
## 3206   male  46
## 3207 female  43
## 3208   male  45
## 3209   male  40
## 3210 female  52
## 3211 female  53
## 3212 female  50
## 3213 female  46
## 3214   male  44
## 3215   male  46
## 3216   male  58
## 3217   male  39
## 3218   male  41
## 3219   male  48
## 3220   male  60
## 3221 female  57
## 3222   male  52
## 3223   male  41
## 3224   male  53
## 3225   male  55
## 3226   male  63
## 3227 female  53
## 3228   male  49
## 3229   male  35
## 3230   male  41
## 3231 female  53
## 3232   male  41
## 3233   male  45
## 3234 female  52
## 3235   male  47
## 3236   male  55
## 3237 female  42
## 3238   male  41
## 3239   male  41
## 3240   male  49
## 3241 female  49
## 3242   male  56
## 3243   male  42
## 3244   male  38
## 3245   male  40
## 3246   male  42
## 3247   male  41
## 3248 female  42
## 3249   male  39
## 3250 female  38
## 3251   male  36
## 3252 female  37
## 3253   male  44
## 3254 female  35
## 3255 female  62
## 3256   male  39
## 3257 female  67
## 3258   male  63
## 3259   male  44
## 3260 female  45
## 3261 female  55
## 3262   male  48
## 3263   male  57
## 3264   male  48
## 3265   male  43
## 3266   male  42
## 3267 female  53
## 3268   male  60
## 3269 female  55
## 3270   male  44
## 3271 female  46
## 3272   male  63
## 3273 female  68
## 3274   male  52
## 3275   male  44
## 3276 female  47
## 3277 female  41
## 3278 female  63
## 3279   male  52
## 3280   male  58
## 3281 female  53
## 3282   male  44
## 3283   male  47
## 3284   male  45
## 3285 female  40
## 3286   male  51
## 3287   male  46
## 3288 female  41
## 3289   male  59
## 3290 female  62
## 3291   male  47
## 3292 female  52
## 3293   male  45
## 3294 female  60
## 3295 female  46
## 3296 female  36
## 3297 female  55
## 3298 female  39
## 3299   male  43
## 3300 female  38
## 3301   male  39
## 3302   male  47
## 3303   male  45
## 3304   male  65
## 3305   male  36
## 3306   male  63
## 3307   male  65
## 3308   male  45
## 3309 female  47
## 3310 female  44
## 3311   male  62
## 3312   male  42
## 3313   male  51
## 3314   male  51
## 3315 female  43
## 3316   male  36
## 3317 female  58
## 3318   male  49
## 3319 female  60
## 3320   male  39
## 3321 female  37
## 3322   male  41
## 3323 female  36
## 3324   male  41
## 3325 female  39
## 3326 female  37
## 3327   male  59
## 3328   male  44
## 3329 female  40
## 3330   male  36
## 3331   male  56
## 3332   male  39
## 3333 female  45
## 3334 female  45
## 3335   male  51
## 3336 female  44
## 3337   male  46
## 3338   male  55
## 3339   male  41
## 3340   male  47
## 3341   male  42
## 3342 female  42
## 3343 female  46
## 3344 female  46
## 3345   male  45
## 3346   male  53
## 3347 female  42
## 3348   male  45
## 3349 female  52
## 3350   male  57
## 3351 female  45
## 3352 female  61
## 3353   male  51
## 3354   male  48
## 3355   male  60
## 3356   male  56
## 3357   male  36
## 3358   male  64
## 3359   male  51
## 3360 female  40
## 3361   male  38
## 3362   male  54
## 3363   male  52
## 3364   male  56
## 3365 female  40
## 3366 female  37
## 3367 female  53
## 3368   male  52
## 3369 female  54
## 3370 female  45
## 3371   male  40
## 3372   male  38
## 3373   male  62
## 3374 female  54
## 3375   male  42
## 3376   male  53
## 3377   male  46
## 3378   male  48
## 3379 female  42
## 3380   male  41
## 3381 female  56
## 3382   male  61
## 3383   male  52
## 3384 female  61
## 3385   male  38
## 3386   male  51
## 3387 female  40
## 3388 female  45
## 3389   male  45
## 3390   male  42
## 3391   male  54
## 3392   male  42
## 3393 female  50
## 3394 female  51
## 3395 female  50
## 3396   male  37
## 3397 female  46
## 3398 female  47
## 3399 female  45
## 3400 female  47
## 3401   male  62
## 3402   male  42
## 3403 female  45
## 3404 female  43
## 3405   male  38
## 3406 female  42
## 3407   male  55
## 3408 female  53
## 3409   male  35
## 3410   male  42
## 3411   male  45
## 3412   male  39
## 3413 female  61
## 3414 female  40
## 3415 female  48
## 3416   male  35
## 3417 female  37
## 3418   male  47
## 3419   male  39
## 3420 female  45
## 3421   male  40
## 3422   male  46
## 3423 female  59
## 3424 female  47
## 3425 female  68
## 3426   male  44
## 3427   male  39
## 3428   male  49
## 3429   male  61
## 3430   male  58
## 3431   male  59
## 3432   male  64
## 3433   male  49
## 3434   male  41
## 3435   male  40
## 3436   male  60
## 3437   male  56
## 3438   male  46
## 3439   male  56
## 3440   male  57
## 3441 female  54
## 3442 female  38
## 3443   male  39
## 3444 female  34
## 3445   male  61
## 3446 female  52
## 3447   male  53
## 3448 female  46
## 3449   male  59
## 3450   male  45
## 3451   male  37
## 3452   male  43
## 3453   male  42
## 3454   male  38
## 3455   male  38
## 3456 female  47
## 3457   male  62
## 3458   male  46
## 3459   male  60
## 3460 female  41
## 3461   male  44
## 3462 female  42
## 3463   male  57
## 3464   male  44
## 3465   male  45
## 3466   male  38
## 3467 female  46
## 3468   male  62
## 3469   male  69
## 3470 female  39
## 3471   male  56
## 3472   male  50
## 3473 female  42
## 3474   male  42
## 3475   male  36
## 3476   male  40
## 3477   male  67
## 3478 female  43
## 3479 female  48
## 3480   male  43
## 3481   male  52
## 3482 female  41
## 3483   male  34
## 3484   male  45
## 3485 female  36
## 3486   male  39
## 3487   male  47
## 3488   male  52
## 3489   male  63
## 3490   male  47
## 3491 female  46
## 3492 female  51
## 3493   male  47
## 3494   male  42
## 3495   male  55
## 3496   male  44
## 3497   male  44
## 3498   male  40
## 3499   male  39
## 3500 female  52
## 3501   male  55
## 3502   male  50
## 3503   male  43
## 3504 female  56
## 3505   male  44
## 3506   male  45
## 3507   male  43
## 3508 female  48
## 3509   male  54
## 3510   male  36
## 3511   male  37
## 3512 female  44
## 3513   male  48
## 3514   male  49
## 3515   male  41
## 3516   male  38
## 3517   male  65
## 3518   male  48
## 3519   male  48
## 3520   male  54
## 3521   male  41
## 3522   male  42
## 3523   male  44
## 3524   male  41
## 3525 female  38
## 3526   male  36
## 3527   male  54
## 3528   male  37
## 3529   male  46
## 3530 female  40
## 3531   male  39
## 3532 female  50
## 3533   male  46
## 3534   male  40
## 3535   male  40
## 3536   male  37
## 3537   male  53
## 3538   male  51
## 3539 female  43
## 3540   male  50
## 3541   male  54
## 3542   male  37
## 3543 female  35
## 3544   male  41
## 3545   male  50
## 3546   male  57
## 3547   male  57
## 3548   male  44
## 3549   male  55
## 3550   male  61
## 3551   male  43
## 3552   male  42
## 3553   male  52
## 3554   male  59
## 3555 female  41
## 3556 female  39
## 3557   male  54
## 3558   male  39
## 3559   male  59
## 3560   male  42
## 3561   male  38
## 3562   male  51
## 3563 female  39
## 3564   male  40
## 3565   male  60
## 3566   male  43
## 3567   male  38
## 3568   male  40
## 3569 female  44
## 3570 female  38
## 3571 female  42
## 3572   male  47
## 3573   male  49
## 3574   male  55
## 3575   male  66
## 3576 female  50
## 3577   male  43
## 3578 female  44
## 3579   male  57
## 3580   male  43
## 3581 female  47
## 3582   male  46
## 3583   male  40
## 3584   male  64
## 3585   male  41
## 3586   male  48
## 3587 female  48
## 3588 female  58
## 3589   male  46
## 3590   male  37
## 3591 female  40
## 3592   male  40
## 3593 female  47
## 3594   male  46
## 3595 female  45
## 3596   male  46
## 3597   male  39
## 3598   male  44
## 3599   male  39
## 3600 female  64
## 3601   male  55
## 3602   male  38
## 3603   male  60
## 3604   male  37
## 3605   male  62
## 3606   male  41
## 3607   male  59
## 3608   male  45
## 3609   male  42
## 3610   male  56
## 3611   male  56
## 3612 female  53
## 3613   male  45
## 3614   male  50
## 3615 female  52
## 3616   male  39
## 3617 female  50
## 3618   male  41
## 3619   male  50
## 3620   male  49
## 3621   male  40
## 3622   male  53
## 3623   male  48
## 3624   male  45
## 3625   male  40
## 3626   male  42
## 3627   male  58
## 3628   male  49
## 3629   male  47
## 3630   male  56
## 3631   male  42
## 3632   male  54
## 3633   male  47
## 3634 female  56
## 3635   male  51
## 3636   male  37
## 3637   male  52
## 3638   male  40
## 3639   male  44
## 3640   male  56
## 3641   male  54
## 3642   male  41
## 3643   male  37
## 3644   male  47
## 3645   male  50
## 3646   male  45
## 3647   male  50
## 3648   male  43
## 3649 female  46
## 3650   male  63
## 3651 female  50
## 3652   male  63
## 3653   male  37
## 3654   male  52
## 3655   male  42
## 3656   male  52
## 3657   male  39
## 3658 female  49
## 3659   male  48
## 3660   male  54
## 3661   male  40
## 3662   male  50
## 3663   male  47
## 3664   male  42
## 3665 female  39
## 3666   male  44
## 3667   male  39
## 3668   male  58
## 3669 female  42
## 3670 female  49
## 3671   male  46
## 3672   male  48
## 3673   male  40
## 3674   male  45
## 3675   male  61
## 3676   male  46
## 3677   male  38
## 3678   male  52
## 3679   male  64
## 3680   male  50
## 3681   male  57
## 3682 female  42
## 3683   male  54
## 3684   male  58
## 3685   male  51
## 3686   male  41
## 3687   male  39
## 3688 female  42
## 3689   male  41
## 3690   male  49
## 3691   male  43
## 3692   male  48
## 3693   male  39
## 3694 female  46
## 3695   male  48
## 3696   male  60
## 3697   male  47
## 3698   male  51
## 3699 female  43
## 3700   male  39
## 3701   male  40
## 3702   male  37
## 3703   male  57
## 3704 female  42
## 3705   male  44
## 3706   male  43
## 3707   male  36
## 3708   male  44
## 3709 female  47
## 3710   male  47
## 3711   male  42
## 3712   male  54
## 3713   male  38
## 3714 female  61
## 3715   male  45
## 3716   male  58
## 3717 female  63
## 3718   male  42
## 3719   male  52
## 3720   male  43
## 3721   male  39
## 3722   male  53
## 3723   male  53
## 3724   male  53
## 3725   male  62
## 3726 female  44
## 3727 female  50
## 3728   male  54
## 3729   male  36
## 3730   male  36
## 3731 female  38
## 3732   male  40
## 3733   male  55
## 3734 female  48
## 3735 female  55
## 3736   male  43
## 3737   male  50
## 3738   male  50
## 3739   male  52
## 3740   male  45
## 3741   male  42
## 3742   male  42
## 3743   male  44
## 3744   male  40
## 3745   male  48
## 3746   male  54
## 3747   male  48
## 3748   male  44
## 3749 female  59
## 3750   male  43
## 3751   male  49
## 3752   male  43
## 3753   male  39
## 3754   male  59
## 3755 female  37
## 3756   male  36
## 3757   male  55
## 3758   male  41
## 3759   male  44
## 3760   male  57
## 3761   male  48
## 3762 female  53
## 3763   male  40
## 3764   male  41
## 3765   male  45
## 3766   male  50
## 3767   male  62
## 3768   male  41
## 3769   male  51
## 3770   male  43
## 3771   male  39
## 3772   male  35
## 3773   male  44
## 3774   male  35
## 3775 female  56
## 3776 female  63
## 3777   male  51
## 3778   male  53
## 3779   male  39
## 3780   male  40
## 3781   male  43
## 3782   male  40
## 3783   male  49
## 3784 female  47
## 3785   male  50
## 3786   male  55
## 3787   male  53
## 3788   male  44
## 3789   male  44
## 3790   male  43
## 3791   male  48
## 3792   male  50
## 3793   male  41
## 3794   male  67
## 3795   male  48
## 3796   male  42
## 3797 female  46
## 3798   male  56
## 3799 female  37
## 3800   male  39
## 3801   male  49
## 3802   male  39
## 3803   male  64
## 3804   male  51
## 3805   male  46
## 3806   male  54
## 3807   male  54
## 3808   male  44
## 3809   male  43
## 3810   male  54
## 3811   male  42
## 3812   male  44
## 3813   male  41
## 3814   male  49
## 3815   male  62
## 3816   male  41
## 3817   male  41
## 3818   male  44
## 3819   male  43
## 3820 female  47
## 3821   male  59
## 3822   male  46
## 3823   male  60
## 3824   male  58
## 3825   male  44
## 3826   male  36
## 3827   male  38
## 3828   male  63
## 3829   male  53
## 3830 female  40
## 3831   male  43
## 3832 female  37
## 3833   male  51
## 3834   male  47
## 3835   male  40
## 3836   male  51
## 3837   male  46
## 3838   male  48
## 3839   male  54
## 3840   male  53
## 3841   male  39
## 3842   male  46
## 3843   male  50
## 3844   male  54
## 3845 female  36
## 3846   male  59
## 3847 female  47
## 3848   male  38
## 3849   male  45
## 3850   male  46
## 3851   male  40
## 3852   male  53
## 3853   male  43
## 3854   male  59
## 3855 female  46
## 3856   male  42
## 3857   male  63
## 3858   male  63
## 3859   male  41
## 3860   male  48
## 3861   male  41
## 3862   male  40
## 3863   male  38
## 3864 female  41
## 3865   male  47
## 3866 female  57
## 3867   male  44
## 3868   male  57
## 3869   male  40
## 3870   male  45
## 3871   male  53
## 3872   male  62
## 3873   male  43
## 3874   male  42
## 3875   male  46
## 3876   male  41
## 3877   male  40
## 3878 female  40
## 3879   male  48
## 3880   male  56
## 3881   male  46
## 3882   male  55
## 3883   male  43
## 3884   male  47
## 3885   male  38
## 3886   male  51
## 3887   male  53
## 3888   male  49
## 3889   male  67
## 3890   male  39
## 3891   male  48
## 3892   male  58
## 3893   male  59
## 3894   male  48
## 3895   male  46
## 3896   male  37
## 3897   male  49
## 3898   male  56
## 3899   male  50
## 3900   male  40
insurance %>% select(diastolic, systolic, heart_rate) %>% glimpse()
## Rows: 3,900
## Columns: 3
## $ diastolic  <dbl> 72.0, 72.0, 76.0, 80.0, 80.0, 72.5, 69.0, 73.0, 77.0, 81.5,…
## $ systolic   <dbl> 110.0, 121.0, 127.5, 122.5, 119.0, 107.5, 109.5, 123.0, 123…
## $ heart_rate <dbl> 95, 64, 81, 90, 62, 62, 75, 66, 65, 93, 70, 85, 58, 83, 75,…
insurance %>% top_n(10, heart_rate) %>% arrange(desc(heart_rate))
##    age    sex current_smoker heart_rate systolic diastolic cigs_per_day chol
## 1   57 female             no        143    133.0        94            0  257
## 2   64 female             no        140    160.0        82            0  312
## 3   40 female            yes        130    118.0        79           20  210
## 4   44 female             no        125    107.0        69            0  160
## 5   36   male             no        125    155.0        74            0  167
## 6   40 female             no        125    158.0       103            0  202
## 7   48 female             no        122    165.0        84            0  246
## 8   67 female            yes        122    169.0        82            5  245
## 9   55   male             no        120    134.0        80            0  199
## 10  50 female             no        120    152.0        92            0  236
## 11  59 female             no        120    149.0        85            0  233
## 12  42 female             no        120    127.5        90            0  158
## 13  59 female            yes        120    175.0        82           10  312
## 14  48 female            yes        120    130.0        84           20  271
## 15  46   male            yes        120    132.5        69           20  273
##       grupo_edad
## 1  Adulto maduro
## 2  Adulto maduro
## 3   Adulto joven
## 4   Adulto joven
## 5   Adulto joven
## 6   Adulto joven
## 7  Adulto maduro
## 8          Mayor
## 9  Adulto maduro
## 10 Adulto maduro
## 11 Adulto maduro
## 12  Adulto joven
## 13 Adulto maduro
## 14 Adulto maduro
## 15 Adulto maduro

8 8. Estadística Inferencial

8.1 Intervalo de Confianza para la Media

intervalo_confianza_media <- function(data, variable, alpha = 0.05) {
  x <- data[[variable]]
  media_x <- mean(x, na.rm = TRUE)
  sd_x <- sd(x, na.rm = TRUE)
  n <- sum(!is.na(x))
  z <- qnorm(1 - alpha / 2)
  inf <- media_x - z * (sd_x / sqrt(n))
  sup <- media_x + z * (sd_x / sqrt(n))
  data.frame(Variable = variable, Media = media_x, Desviacion = sd_x,
             n = n, Z = z, Limite_Inferior = inf, Limite_Superior = sup)
}

intervalo_confianza_media(insurance, "heart_rate")
##     Variable    Media Desviacion    n        Z Limite_Inferior Limite_Superior
## 1 heart_rate 75.68897   12.01538 3900 1.959964        75.31188        76.06607

8.2 Distribución Normal con ggplot

grafico_normal <- function(media = 0, sd = 1, z1 = -1.95, z2 = 1.95, region = "outside") {
  x_vals <- seq(-4, 4, 0.01)
  y_vals <- dnorm(x_vals, mean = media, sd = sd)
  df <- data.frame(x = x_vals, y = y_vals)

  gg <- ggplot(df, aes(x = x, y = y)) +
    geom_line(color = "blue") +
    labs(title = "Distribución Normal Estándar", x = "Z", y = "Densidad") +
    theme_minimal()

  if (region == "outside") {
    gg <- gg + geom_area(data = subset(df, x < z1 | x > z2), aes(y = y), fill = "red", alpha = 0.4)
  } else if (region == "between") {
    gg <- gg + geom_area(data = subset(df, x > z1 & x < z2), aes(y = y), fill = "green", alpha = 0.4)
  } else if (region == "below") {
    gg <- gg + geom_area(data = subset(df, x < z2), aes(y = y), fill = "blue", alpha = 0.4)
  } else if (region == "above") {
    gg <- gg + geom_area(data = subset(df, x > z2), aes(y = y), fill = "blue", alpha = 0.4)
  }

  gg + geom_vline(xintercept = c(z1, z2), linetype = "dashed", color = "black")
}

grafico_normal(region = "outside")

grafico_normal(z1 = -1, z2 = 1, region = "between")

grafico_normal(z1 = -4, z2 = 1, region = "below")

grafico_normal(z1 = 1, z2 = 4, region = "above")

8.3 Pruebas de Hipótesis (una muestra)

z.test(x = insurance$age, sigma.x = sd(insurance$age), mu = 30, alternative = "two.sided", conf.level = 0.95)
## 
##  One-sample z-Test
## 
## data:  insurance$age
## z = 142.55, p-value < 2.2e-16
## alternative hypothesis: true mean is not equal to 30
## 95 percent confidence interval:
##  49.26932 49.80658
## sample estimates:
## mean of x 
##  49.53795
z.test(x = insurance$age, sigma.x = sd(insurance$age), mu = 30, alternative = "less", conf.level = 0.95)
## 
##  One-sample z-Test
## 
## data:  insurance$age
## z = 142.55, p-value = 1
## alternative hypothesis: true mean is less than 30
## 95 percent confidence interval:
##        NA 49.76339
## sample estimates:
## mean of x 
##  49.53795
z.test(x = insurance$age, sigma.x = sd(insurance$age), mu = 30, alternative = "greater", conf.level = 0.95)
## 
##  One-sample z-Test
## 
## data:  insurance$age
## z = 142.55, p-value < 2.2e-16
## alternative hypothesis: true mean is greater than 30
## 95 percent confidence interval:
##  49.31251       NA
## sample estimates:
## mean of x 
##  49.53795

8.4 Comparación de Medias por Sexo

sex.heart_rate <- split(insurance, insurance$sex)
z.test(x = sex.heart_rate$female$heart_rate,
       y = sex.heart_rate$male$heart_rate,
       sigma.x = sd(sex.heart_rate$female$heart_rate),
       sigma.y = sd(sex.heart_rate$male$heart_rate),
       mu = 0,
       alternative = "greater",
       conf.level = 0.95)
## 
##  Two-sample z-Test
## 
## data:  sex.heart_rate$female$heart_rate and sex.heart_rate$male$heart_rate
## z = 7.0142, p-value = 1.156e-12
## alternative hypothesis: true difference in means is greater than 0
## 95 percent confidence interval:
##  2.053926       NA
## sample estimates:
## mean of x mean of y 
##  76.94041  74.25728

8.5 Pruebas de Proporciones e Independencia

fa <- table(insurance$sex)
prop.test(x = fa[1], n = sum(fa), p = 0.5, alternative = "two.sided", conf.level = 0.95)
## 
##  1-sample proportions test with continuity correction
## 
## data:  fa[1] out of sum(fa), null probability 0.5
## X-squared = 17.467, df = 1, p-value = 2.924e-05
## alternative hypothesis: true p is not equal to 0.5
## 95 percent confidence interval:
##  0.5177792 0.5493336
## sample estimates:
##         p 
## 0.5335897
prop.test(c(fa[1], fa[2]), c(sum(fa), sum(fa)), alternative = "two.sided", conf.level = 0.95)
## 
##  2-sample test for equality of proportions with continuity correction
## 
## data:  c(fa[1], fa[2]) out of c(sum(fa), sum(fa))
## X-squared = 34.934, df = 1, p-value = 3.411e-09
## alternative hypothesis: two.sided
## 95 percent confidence interval:
##  0.04478099 0.08957799
## sample estimates:
##    prop 1    prop 2 
## 0.5335897 0.4664103
# Prueba de independencia
chisq.test(fa)
## 
##  Chi-squared test for given probabilities
## 
## data:  fa
## X-squared = 17.601, df = 1, p-value = 2.724e-05