Es mostren les 10 primeres files
## Pregnancies Glucose BloodPressure SkinThickness Insulin BMI DiabetesPedigreeFunction Age Outcome
## 1 6 148 72 35 0 33.6 0.627 50 1
## 2 1 85 66 29 0 26.6 0.351 31 0
## 3 8 183 64 0 0 23.3 0.672 32 1
## 4 1 89 66 23 94 28.1 0.167 21 0
## 5 0 137 40 35 168 43.1 2.288 33 1
## 6 5 116 74 0 0 25.6 0.201 30 0
## 7 3 78 50 32 88 31.0 0.248 26 1
## 8 10 115 0 0 0 35.3 0.134 29 0
## 9 2 197 70 45 543 30.5 0.158 53 1
## 10 8 125 96 0 0 0.0 0.232 54 1
Indiqueu el nombre de files i columnes. Mostreu els noms i els tipus de les variables.
Es confirma que hi ha 768 observacions(files) amb 9 variables (columnes)
Taula resum amb noms i els tipus de les variables:
## 'data.frame': 768 obs. of 9 variables:
## $ Pregnancies : int 6 1 8 1 0 5 3 10 2 8 ...
## $ Glucose : int 148 85 183 89 137 116 78 115 197 125 ...
## $ BloodPressure : int 72 66 64 66 40 74 50 0 70 96 ...
## $ SkinThickness : int 35 29 0 23 35 0 32 0 45 0 ...
## $ Insulin : int 0 0 0 94 168 0 88 0 543 0 ...
## $ BMI : num 33.6 26.6 23.3 28.1 43.1 25.6 31 35.3 30.5 0 ...
## $ DiabetesPedigreeFunction: num 0.627 0.351 0.672 0.167 2.288 ...
## $ Age : int 50 31 32 21 33 30 26 29 53 54 ...
## $ Outcome : int 1 0 1 0 1 0 1 0 1 1 ...
Comproveu si hi ha files completament duplicades.
No hi ha files duplicades
## [1] 0
Indiqueu si hi ha valors codificats explícitament com a NA al fitxer original.
No hi ha valors codificats com NA
## [1] FALSE
Mostreu un resum descriptiu inicial.
## Pregnancies Glucose BloodPressure SkinThickness Insulin BMI DiabetesPedigreeFunction Age Outcome
## Min. : 0.000 Min. : 0.0 Min. : 0.00 Min. : 0.00 Min. : 0.0 Min. : 0.00 Min. :0.0780 Min. :21.00 Min. :0.000
## 1st Qu.: 1.000 1st Qu.: 99.0 1st Qu.: 62.00 1st Qu.: 0.00 1st Qu.: 0.0 1st Qu.:27.30 1st Qu.:0.2437 1st Qu.:24.00 1st Qu.:0.000
## Median : 3.000 Median :117.0 Median : 72.00 Median :23.00 Median : 30.5 Median :32.00 Median :0.3725 Median :29.00 Median :0.000
## Mean : 3.845 Mean :120.9 Mean : 69.11 Mean :20.54 Mean : 79.8 Mean :31.99 Mean :0.4719 Mean :33.24 Mean :0.349
## 3rd Qu.: 6.000 3rd Qu.:140.2 3rd Qu.: 80.00 3rd Qu.:32.00 3rd Qu.:127.2 3rd Qu.:36.60 3rd Qu.:0.6262 3rd Qu.:41.00 3rd Qu.:1.000
## Max. :17.000 Max. :199.0 Max. :122.00 Max. :99.00 Max. :846.0 Max. :67.10 Max. :2.4200 Max. :81.00 Max. :1.000
Comptabilitzeu quants valors iguals a 0 apareixen en cadascuna d’aquestes variables.
## [1] "Glucose amb valor 0: 5"
## [1] "BloodPressure amb valor 0: 35"
## [1] "SkinThickness amb valor 0; 227"
## [1] "Insulin amb valor 0; 374"
## [1] "BMI amb valor 0; 11"
Substituïu aquests valors per NA
Es fa amb aquesta secció de codi:
diabetes\(Glucose[diabetes\)Glucose == 0] <-
NA
diabetes\(BloodPressure[diabetes\)BloodPressure == 0]
<- NA
diabetes\(SkinThickness[diabetes\)SkinThickness == 0]
<- NA
diabetes\(Insulin[diabetes\)Insulin == 0] <-
NA
diabetes\(BMI[diabetes\)BMI ==
0] <- NA ]
Construïu una taula amb el nombre i el percentatge de valors perduts per variable després de la transformació
## Variable N_NA Percentatge_NA
## Pregnancies Pregnancies 0 0.00
## Glucose Glucose 5 0.65
## BloodPressure BloodPressure 35 4.56
## SkinThickness SkinThickness 227 29.56
## Insulin Insulin 374 48.70
## BMI BMI 11 1.43
## DiabetesPedigreeFunction DiabetesPedigreeFunction 0 0.00
## Age Age 0 0.00
## Outcome Outcome 0 0.00
Es passa a que la variable Outcome sigu
quantitativa, es guarda el status en valors de 0 i 1 a la variable
Outcome_quan
No: 500 - 65,1%
Si: 268 - 34,89%
##
## No Si
## 500 268
taula per visualitzar la categorització entre diabètics i no diabètics
##
## No Si
## Menors de 30 312 84
## 30-49 142 141
## 50 o més 46 43
Taula per visualitzar el percentatge de diabetics
##
## No Si
## Menors de 30 78.79 21.21
## 30-49 50.18 49.82
## 50 o més 51.69 48.31
Es pot observar com el percentatge de diabètics és especialment baix pel grup de població de menors de 30 anys. No es veueu diferències significatives entre els majors de 30 anys.
Representeu la distribució de BMI mitjançant un boxplot.
## Min. 1st Qu. Median Mean 3rd Qu. Max. NAs
## 18.20 27.50 32.30 32.46 36.60 67.10 11
Identifiqueu els valors considerats extrems d’acord amb el criteri utilitzat pel boxplot.
Mínim esperable
## [1] 13.85
Màx_esperable
## [1] 50.25
Indiqueu les observacions afectades.
BMI alt
## Pregnancies Glucose BloodPressure SkinThickness Insulin BMI DiabetesPedigreeFunction Age Outcome Outcome_quan AgeGroup
## 121 0 162 76 56 100 53.2 0.759 25 Si 1 Menors de 30
## 126 1 88 30 42 99 55.0 0.496 26 Si 1 Menors de 30
## 178 0 129 110 46 130 67.1 0.319 26 Si 1 Menors de 30
## 194 11 135 NA NA NA 52.3 0.578 40 Si 1 30-49
## 248 0 165 90 33 680 52.3 0.427 23 No 0 Menors de 30
## 304 5 115 98 NA NA 52.9 0.209 28 Si 1 Menors de 30
## 446 0 180 78 63 14 59.4 2.420 25 Si 1 Menors de 30
## 674 3 123 100 35 240 57.3 0.880 22 No 0 Menors de 30
BMI baix
## [1] Pregnancies Glucose BloodPressure SkinThickness Insulin BMI DiabetesPedigreeFunction
## [8] Age Outcome Outcome_quan AgeGroup
## <0 rows> (o 0- extensión row.names)
Interpreteu breument el resultat.
Considerant que els valors de 0 ja s’han eliminat, crec que NO es poden eliminar automàticament aquestes observacions ja que hi ha casos molt extrems de pesos. En aquest cas només han sortit valors per la franja alta
Representeu Insulin mitjançant un boxplot i identifiqueu quantes observacions són considerades valors extrems segons el criteri utilitzat pel boxplot.
## Min. 1st Qu. Median Mean 3rd Qu. Max. NAs
## 14.00 76.25 125.00 155.55 190.00 846.00 374
Mínim esperable
## [1] 0
Màx_esperable
## [1] 360.625
Pacients per sobre del màxim
Quins pacients compleixen amb aquesta condició de tenir l’insulina per sobre del màxim esperable?
## Pregnancies Glucose BloodPressure SkinThickness Insulin BMI DiabetesPedigreeFunction Age Outcome Outcome_quan AgeGroup
## 9 2 197 70 45 543 30.5 0.158 53 Si 1 50 o més
## 14 1 189 60 23 846 30.1 0.398 59 Si 1 50 o més
## 112 8 155 62 26 495 34.0 0.543 46 Si 1 30-49
## 154 1 153 82 42 485 40.6 0.687 23 No 0 Menors de 30
## 187 8 181 68 36 495 30.1 0.615 60 Si 1 50 o més
## 221 0 177 60 29 478 34.6 1.072 21 Si 1 Menors de 30
## 229 4 197 70 39 744 36.7 2.329 31 No 0 30-49
## 232 6 134 80 37 370 46.2 0.238 46 Si 1 30-49
## 248 0 165 90 33 680 52.3 0.427 23 No 0 Menors de 30
## 249 9 124 70 33 402 35.4 0.282 34 No 0 30-49
## 259 1 193 50 16 375 25.9 0.655 24 No 0 Menors de 30
## 287 5 155 84 44 545 38.7 0.619 34 No 0 30-49
## 371 3 173 82 48 465 38.4 2.137 25 Si 1 Menors de 30
## 393 1 131 64 14 415 23.7 0.389 21 No 0 Menors de 30
## 410 1 172 68 49 579 42.4 0.702 28 Si 1 Menors de 30
## 416 3 173 84 33 474 35.7 0.258 22 Si 1 Menors de 30
## 487 1 139 62 41 480 40.7 0.536 21 No 0 Menors de 30
## 585 8 124 76 24 600 28.7 0.687 52 Si 1 50 o més
## 646 2 157 74 35 440 39.4 0.134 30 No 0 30-49
## 656 2 155 52 27 540 38.7 0.240 25 Si 1 Menors de 30
## 696 7 142 90 24 480 30.4 0.128 43 Si 1 30-49
## 711 3 158 64 13 387 31.2 0.295 24 No 0 Menors de 30
## 716 7 187 50 33 392 33.9 0.826 34 Si 1 30-49
## 754 0 181 88 44 510 43.3 0.222 26 Si 1 Menors de 30
10 pacients amb els valors d’insulina més elevats, ordenats en ordre decreixent, incloent-hi les variables Age, BMI, Glucose, Insulin i Outcome.
## Age BMI Glucose Insulin Outcome
## 14 59 30.1 189 846 Si
## 229 31 36.7 197 744 No
## 248 23 52.3 165 680 No
## 585 52 28.7 124 600 Si
## 410 28 42.4 172 579 Si
## 287 34 38.7 155 545 No
## 9 53 30.5 197 543 Si
## 656 25 38.7 155 540 Si
## 754 26 43.3 181 510 Si
## 112 46 34.0 155 495 Si
Calculeu la mitjana i la mediana d’Insulin i compareu totes dues mesures. Quin efecte semblen tenir els valors elevats sobre la mitjana?
## [1] "Mitjana: 155.548223350254"
## [1] "Mediana: 125"
Es pot comprovar que malgrat la mediana és de 125 la mitjana és molt més alta 155 ja que els valors molt alts tendeixen a “inflar” la mitjana. Això també passa amb les rendes altres, és un exemple clar de perquè la mitjana a vegades no reflexa la situació més habitual.
Calculeu les correlacions entre Pregnancies, Glucose, BloodPressure, SkinThickness, Insulin, BMI, DiabetesPedigreeFunction, Age i una versió numèrica 0/1 d’Outcome. Per calcular les correlacions, utilitzeu únicament les observacions disponibles per a cada parell de variables. Mostreu una representació gràfica de la matriu de correlacions i interpreteu les relacions més destacables.
## Pregnancies Glucose BloodPressure SkinThickness Insulin BMI DiabetesPedigreeFunction Age Outcome_quan
## Pregnancies 1.000000000 0.1982910 0.2133548 0.0932094 0.07898363 -0.02534728 0.007562116 0.67960847 0.2565660
## Glucose 0.198291043 1.0000000 0.2100266 0.1988558 0.58122301 0.20951592 0.140180180 0.34364150 0.5157027
## BloodPressure 0.213354775 0.2100266 1.0000000 0.2325712 0.09851150 0.30440337 -0.015971104 0.30003895 0.1926733
## SkinThickness 0.093209397 0.1988558 0.2325712 1.0000000 0.18219906 0.66435487 0.160498526 0.16776114 0.2559357
## Insulin 0.078983625 0.5812230 0.0985115 0.1821991 1.00000000 0.22639652 0.135905781 0.21708199 0.3014292
## BMI -0.025347276 0.2095159 0.3044034 0.6643549 0.22639652 1.00000000 0.158771043 0.06981380 0.2701184
## DiabetesPedigreeFunction 0.007562116 0.1401802 -0.0159711 0.1604985 0.13590578 0.15877104 1.000000000 0.08502911 0.2093295
## Age 0.679608470 0.3436415 0.3000389 0.1677611 0.21708199 0.06981380 0.085029106 1.00000000 0.3508038
## Outcome_quan 0.256565956 0.5157027 0.1926733 0.2559357 0.30142922 0.27011841 0.209329511 0.35080380 1.0000000
Quina variable presenta la correlació, en valor absolut, més gran amb el diagnòstic de diabetis?
La variable Glucose 0.5157027 és la que té una correlació més alta amb si ets o no diabètic
Expliqueu per què una correlació observada no permet establir, per si sola, una relació causal.
Es coneguda la paradoxa que a vegades mirant la correlació, es troben com dos variables, semblen tenir una correlació molt alta. Però aquesta, pot ser fruit de la casualtitat. Hi ha pàgines web com: https://www.tylervigen.com/spurious-correlations On es troben tot tipus de correlacions no causals. Per exemple es mostra com l’energia solar genearada a Belize està correlacionada amb un r=0,988 amb el nombre de inspectors de Focs a Florida. És evident que aquestes dues dades no estan correlacionades i que això és fruit de la casualitat.
## Pregnancies Glucose BloodPressure DiabetesPedigreeFunction Age BMI
## 10 8 125 96 0.232 54 37.6
## 50 7 105 NA 0.305 24 29.8
## 61 2 84 NA 0.304 21 23.5
## 82 2 74 NA 0.102 22 23.5
## 146 0 102 75 0.572 21 39.4
## 372 0 118 64 1.731 21 30.5
## 427 0 94 NA 0.256 25 25.0
## 495 3 80 NA 0.174 22 23.5
## 523 6 114 NA 0.189 26 25.2
## 685 5 136 82 0.640 69 32.0
## 707 10 115 NA 0.261 30 30.0
## Valors perduts a BMI després de la imputació: 0
Expliqueu breument per què pot ser preferible utilitzar informació d’altres variables en lloc de substituir tots els valors perduts per la mitjana de BMI.
En cas de tenir aquesta informació addicional d’altres columnes, en cas de saber que hi ha certa correlació entre les variables, el fet de fer la estimació de quan hauria de ser el valor de BIM, ens permet acostar-nos a tenir un valor de mitjana i mediana més fidel a la realitat ja que les interpolacions que es fan haurien d’ajudar-nos a tenir més obervacions a partir d’aquesta interpolació a través de la resta de variables. La explicació seria: En funció del nombre d’embarassos, la mesura de glucosa en sang, la seva pressió arterial, si ja hi ha hagut altres casos de diabetes a la família i la seva edat, et pots fer una idea més clara de com pot ser aquesta persona, i per tant, el seu BMI esperat.
Presenteu una taula amb la mitjana i la mediana i una altra taula amb la desviació estàndard i la desviació absoluta respecte de la mediana (MAD).
Mesures de tendència(Mitjana i Mediana)
## # A tibble: 2 × 7
## Outcome Mitjana_Glucose Mediana_Glucose Mitjana_BMI Mediana_BMI Mitjana_Age Mediana_Age
## <fct> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 No 111. 107 30.9 30.1 31.2 27
## 2 Si 142. 140 35.4 34.3 37.1 36
Mesures de dispersió (Desviació estàndard i MAD)
## # A tibble: 2 × 7
## Outcome SD_Glucose MAD_Glucose SD_BMI MAD_BMI SD_Age MAD_Age
## <fct> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 No 24.8 23.7 6.56 7.12 11.7 7.41
## 2 Si 29.6 34.8 6.61 5.71 11.0 11.9
Compareu els resultats entre pacients amb diabetis i sense. Quina variable sembla mostrar diferències més grans entre tots dos grups? Justifiqueu la resposta utilitzant els resultats obtinguts.
Els pacients amb diabetes tenen uns valor molt alt de mitjana de glucosa així com una gran dispersió (Els hi puja i baixa molt el sucre) Encara que increments en l’edat i el BMI també sembla que son paràmetres que avaforeixen o son conseqüència de la diabetes. És per aquest motiu que es selecciona la variableB Glucosa, malgrat les altres dues variables també tenen rellevància.
Utilitzant la mostra global, compareu la mitjana i la mediana, així com la desviació estàndard i el MAD. Quines diferències observeu entre les mesures robustes i no robustes?
Les mesures robustes (mediana i MAD) ofereixen una representació més fidel del “centre” i la variabilitat típica de les dades quan aquestes estan esbiaixades o contenen valors atípics, mentre que les no robustes (mitjana i SD) estan molt afectades pels valors extrems.