Parte 1 : Indice de felicidad
link1="https://es.wikipedia.org/wiki/%C3%8Dndice_global_de_felicidad"
path1='//*/div/table/tbody'
library(htmltab)
felicidad=htmltab(doc = link1,which = path1)
str(felicidad)
## 'data.frame': 156 obs. of 9 variables:
## $ № : chr "1" "2" "3" "4" ...
## $ País : chr "Finlandia" "Colombia" "Noruega" "Dinamarca" ...
## $ Puntuación : chr "7.633" "7.594" "7.560" "7.555" ...
## $ PIB per cápita : chr "1.305" "1.456" "1.372" "1.351" ...
## $ Apoyo social : chr "1.592" "1.582" "1.595" "1.590" ...
## $ Esperanza de años de vida saludable : chr "0.874" "0.873" "0.870" "0.868" ...
## $ Libertad para tomar decisiones vitales: chr "0.681" "0.686" "0.685" "0.683" ...
## $ Generosidad : chr "0.192" "0.286" "0.285" "0.284" ...
## $ Percepción de la corrupción : chr "0.393" "0.130" "0.410" "0.408" ...
names(felicidad)
## [1] "№"
## [2] "País"
## [3] "Puntuación"
## [4] "PIB per cápita"
## [5] "Apoyo social"
## [6] "Esperanza de años de vida saludable"
## [7] "Libertad para tomar decisiones vitales"
## [8] "Generosidad"
## [9] "Percepción de la corrupción"
newN=c("n","pais","puntuacion","pbi","apoyosoc","esperanza","libertad","generosidad","percepcioncorrupcion")
names(felicidad)=newN
felicidad$n =NULL
felicidad[,c(2:8)]=lapply(felicidad[,c(2:8)],as.numeric)
felicidad[!complete.cases(felicidad),]
## [1] pais puntuacion pbi
## [4] apoyosoc esperanza libertad
## [7] generosidad percepcioncorrupcion
## <0 rows> (or 0-length row.names)
felicidad$pais=as.character(felicidad$pais)
str(felicidad)
## 'data.frame': 156 obs. of 8 variables:
## $ pais : chr "Finlandia" "Colombia" "Noruega" "Dinamarca" ...
## $ puntuacion : num 7.63 7.59 7.56 7.55 7.5 ...
## $ pbi : num 1.3 1.46 1.37 1.35 1.34 ...
## $ apoyosoc : num 1.59 1.58 1.59 1.59 1.64 ...
## $ esperanza : num 0.874 0.873 0.87 0.868 0.914 0.927 0.878 0.896 0.876 0.913 ...
## $ libertad : num 0.681 0.686 0.685 0.683 0.677 0.66 0.638 0.653 0.669 0.659 ...
## $ generosidad : num 0.192 0.286 0.285 0.284 0.353 0.256 0.333 0.321 0.365 0.285 ...
## $ percepcioncorrupcion: num 0.393 0.13 0.41 0.408 0.138 0.357 0.295 0.291 0.389 0.383 ...
summary(felicidad)
## pais puntuacion pbi apoyosoc
## Length:156 Min. :2.905 Min. :0.0000 Min. :0.000
## Class :character 1st Qu.:4.454 1st Qu.:0.6162 1st Qu.:1.067
## Mode :character Median :5.378 Median :0.9495 Median :1.255
## Mean :5.384 Mean :0.8941 Mean :1.214
## 3rd Qu.:6.168 3rd Qu.:1.2025 3rd Qu.:1.466
## Max. :7.633 Max. :2.0960 Max. :1.644
## esperanza libertad generosidad percepcioncorrupcion
## Min. :0.0000 Min. :0.0000 Min. :0.0000 Min. :0.00000
## 1st Qu.:0.4223 1st Qu.:0.3560 1st Qu.:0.1108 1st Qu.:0.05175
## Median :0.6510 Median :0.4870 Median :0.1750 Median :0.08200
## Mean :0.5989 Mean :0.4555 Mean :0.1821 Mean :0.11306
## 3rd Qu.:0.7820 3rd Qu.:0.5800 3rd Qu.:0.2422 3rd Qu.:0.13650
## Max. :1.0300 Max. :0.7240 Max. :0.5980 Max. :0.45700
Analisis:
shapiro.test(felicidad$puntuacion)
##
## Shapiro-Wilk normality test
##
## data: felicidad$puntuacion
## W = 0.9838, p-value = 0.06475
el p-value es mayor a 0.05 Es normmal
library(DescTools)
Skew(felicidad$puntuacion)
## [1] 0.03563908
es asimetria positiva
No hay atipicos
Correlacion entre PIB y generosidad
shapiro.test(felicidad$pbi)
##
## Shapiro-Wilk normality test
##
## data: felicidad$pbi
## W = 0.97748, p-value = 0.01176
cor.test(felicidad$pbi,felicidad$generosidad,method = "spearman")
## Warning in cor.test.default(felicidad$pbi, felicidad$generosidad, method =
## "spearman"): Cannot compute exact p-value with ties
##
## Spearman's rank correlation rho
##
## data: felicidad$pbi and felicidad$generosidad
## S = 628978, p-value = 0.9417
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
## rho
## 0.005898693
p-value es de 0.9417 es mayor a 0.05 No correlacion
Tomando a esperanza como variable dep
shapiro.test(felicidad$esperanza)
##
## Shapiro-Wilk normality test
##
## data: felicidad$esperanza
## W = 0.95408, p-value = 5.151e-05
el p-value es menor a 0.05 no es normal
todocontrol2=lm(esperanza~pbi+apoyosoc+percepcioncorrupcion+libertad+generosidad,data = felicidad)
summary(todocontrol2)
##
## Call:
## lm(formula = esperanza ~ pbi + apoyosoc + percepcioncorrupcion +
## libertad + generosidad, data = felicidad)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.44391 -0.06484 0.01431 0.07939 0.24825
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.003605 0.047854 -0.075 0.94004
## pbi 0.445871 0.036908 12.081 < 2e-16 ***
## apoyosoc 0.139466 0.048666 2.866 0.00476 **
## percepcioncorrupcion 0.122647 0.128416 0.955 0.34108
## libertad 0.046434 0.078767 0.590 0.55641
## generosidad -0.002610 0.116482 -0.022 0.98216
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.1295 on 150 degrees of freedom
## Multiple R-squared: 0.7373, Adjusted R-squared: 0.7285
## F-statistic: 84.18 on 5 and 150 DF, p-value: < 2.2e-16
Solo generosidad tiene un efecto inverso
Tomando a percepcion como variable dep
shapiro.test(felicidad$percepcioncorrupcion)
##
## Shapiro-Wilk normality test
##
## data: felicidad$percepcioncorrupcion
## W = 0.81431, p-value = 8.486e-13
El p-value es menor a 0.05 NO es normal
todocontrol3=lm(percepcioncorrupcion~pbi+apoyosoc+esperanza+libertad+generosidad,data = felicidad)
summary(todocontrol3)
##
## Call:
## lm(formula = percepcioncorrupcion ~ pbi + apoyosoc + esperanza +
## libertad + generosidad, data = felicidad)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.18413 -0.05413 -0.01276 0.04097 0.31134
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.04356 0.03013 -1.446 0.150266
## pbi 0.04141 0.03269 1.267 0.207190
## apoyosoc -0.03841 0.03153 -1.218 0.225052
## esperanza 0.04928 0.05160 0.955 0.341075
## libertad 0.19913 0.04727 4.213 4.34e-05 ***
## generosidad 0.25261 0.07090 3.563 0.000492 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.08212 on 150 degrees of freedom
## Multiple R-squared: 0.3098, Adjusted R-squared: 0.2868
## F-statistic: 13.47 on 5 and 150 DF, p-value: 7.529e-11
link2="https://docs.google.com/spreadsheets/d/e/2PACX-1vTFLhCO2eqAth81eQsaT03RQx32n9GRxx6ixAR2trHCt4rWR6QIHk-Ig-b9VmXpBg/pub?output=csv"
ide=read.csv(link2,stringsAsFactors = F)
str(ide)
## 'data.frame': 195 obs. of 13 variables:
## $ regionUbigeo : int 10000 10000 10000 10000 10000 10000 10000 20000 20000 20000 ...
## $ provinciaUbigeo: int 10100 10200 10300 10400 10500 10600 10700 20100 20200 20300 ...
## $ PROVINCIA : chr "CHACHAPOYAS" "BAGUA" "BONGARA" "CONDORCANQUI" ...
## $ IDE : num 0.774 0.662 0.632 0.46 0.605 ...
## $ identidad : num 98.6 94.6 97.5 86.2 96.2 ...
## $ salud : num 25.45 14.61 9.01 8.56 12.42 ...
## $ educacion : num 91.5 79.8 76.4 52.2 74.7 ...
## $ saneamiento : num 70.3 64.5 54.8 37.7 43.3 ...
## $ electrificacion: num 84 67.9 72.2 39.5 67.4 ...
## $ poblacion : int 54783 77438 32317 51802 52185 30236 118747 161003 7974 16879 ...
## $ costa : chr "NO" "NO" "NO" "NO" ...
## $ capital : chr "SI" "NO" "NO" "NO" ...
## $ tamano : chr "Pequena" "Pequena" "Muy pequena" "Pequena" ...
names(ide)
## [1] "regionUbigeo" "provinciaUbigeo" "PROVINCIA"
## [4] "IDE" "identidad" "salud"
## [7] "educacion" "saneamiento" "electrificacion"
## [10] "poblacion" "costa" "capital"
## [13] "tamano"
ide[,c(3,11:13)]=lapply(ide[,c(3,11:13)],as.factor)
ide=ide[complete.cases(ide),]
str(ide)
## 'data.frame': 195 obs. of 13 variables:
## $ regionUbigeo : int 10000 10000 10000 10000 10000 10000 10000 20000 20000 20000 ...
## $ provinciaUbigeo: int 10100 10200 10300 10400 10500 10600 10700 20100 20200 20300 ...
## $ PROVINCIA : Factor w/ 195 levels "ABANCAY","ACOBAMBA",..: 46 19 24 59 116 156 186 90 4 11 ...
## $ IDE : num 0.774 0.662 0.632 0.46 0.605 ...
## $ identidad : num 98.6 94.6 97.5 86.2 96.2 ...
## $ salud : num 25.45 14.61 9.01 8.56 12.42 ...
## $ educacion : num 91.5 79.8 76.4 52.2 74.7 ...
## $ saneamiento : num 70.3 64.5 54.8 37.7 43.3 ...
## $ electrificacion: num 84 67.9 72.2 39.5 67.4 ...
## $ poblacion : int 54783 77438 32317 51802 52185 30236 118747 161003 7974 16879 ...
## $ costa : Factor w/ 2 levels "NO","SI": 1 1 1 1 1 1 1 1 1 1 ...
## $ capital : Factor w/ 2 levels "NO","SI": 2 1 1 1 1 1 1 2 1 1 ...
## $ tamano : Factor w/ 5 levels "Grande","Mediana",..: 5 5 4 5 5 4 2 2 4 4 ...
summary(ide)
## regionUbigeo provinciaUbigeo PROVINCIA IDE
## Min. : 10000 Min. : 10100 ABANCAY : 1 Min. :0.4245
## 1st Qu.: 50000 1st Qu.: 50750 ACOBAMBA : 1 1st Qu.:0.5806
## Median :100000 Median :101100 ACOMAYO : 1 Median :0.6367
## Mean :113795 Mean :114358 AIJA : 1 Mean :0.6519
## 3rd Qu.:170000 3rd Qu.:170250 ALTO AMAZONAS: 1 3rd Qu.:0.7164
## Max. :250000 Max. :250400 AMBO : 1 Max. :0.9104
## (Other) :189
## identidad salud educacion saneamiento
## Min. :81.97 Min. : 2.598 Min. :44.03 Min. : 0.2351
## 1st Qu.:96.54 1st Qu.: 6.485 1st Qu.:74.06 1st Qu.:40.9610
## Median :97.90 Median :10.380 Median :82.49 Median :54.2229
## Mean :97.30 Mean :11.919 Mean :80.42 Mean :56.2641
## 3rd Qu.:98.87 3rd Qu.:14.744 3rd Qu.:89.02 3rd Qu.:72.4629
## Max. :99.50 Max. :44.741 Max. :99.50 Max. :99.5000
##
## electrificacion poblacion costa capital tamano
## Min. :33.83 Min. : 4251 NO:163 NO:170 Grande :15
## 1st Qu.:61.31 1st Qu.: 32188 SI: 32 SI: 25 Mediana :46
## Median :74.29 Median : 63039 Muy grande : 1
## Mean :72.11 Mean : 154543 Muy pequena:69
## 3rd Qu.:85.02 3rd Qu.: 121804 Pequena :64
## Max. :99.50 Max. :8481415
##
Analisis de Variables:
Tomando a identidad como variable dep
shapiro.test(ide$identidad)
##
## Shapiro-Wilk normality test
##
## data: ide$identidad
## W = 0.72664, p-value < 2.2e-16
Es menor a 0.05 , es decir es No Normal
names(ide)
## [1] "regionUbigeo" "provinciaUbigeo" "PROVINCIA"
## [4] "IDE" "identidad" "salud"
## [7] "educacion" "saneamiento" "electrificacion"
## [10] "poblacion" "costa" "capital"
## [13] "tamano"
todocontrol=lm(identidad~regionUbigeo+provinciaUbigeo+PROVINCIA+IDE+salud+educacion+saneamiento+electrificacion+poblacion+costa+capital+tamano,data = ide)
summary(todocontrol)
##
## Call:
## lm(formula = identidad ~ regionUbigeo + provinciaUbigeo + PROVINCIA +
## IDE + salud + educacion + saneamiento + electrificacion +
## poblacion + costa + capital + tamano, data = ide)
##
## Residuals:
## ALL 195 residuals are 0: no residual degrees of freedom!
##
## Coefficients: (14 not defined because of singularities)
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 9.791e+01 NA NA NA
## regionUbigeo -8.213e-04 NA NA NA
## provinciaUbigeo 8.185e-04 NA NA NA
## PROVINCIAACOBAMBA 9.348e-01 NA NA NA
## PROVINCIAACOMAYO -6.460e-01 NA NA NA
## PROVINCIAAIJA 2.201e-01 NA NA NA
## PROVINCIAALTO AMAZONAS -3.579e+00 NA NA NA
## PROVINCIAAMBO 1.319e-01 NA NA NA
## PROVINCIAANDAHUAYLAS 3.800e-01 NA NA NA
## PROVINCIAANGARAES 5.404e-01 NA NA NA
## PROVINCIAANTA 4.868e-02 NA NA NA
## PROVINCIAANTABAMBA -1.907e+00 NA NA NA
## PROVINCIAANTONIO RAYMONDI -1.280e+00 NA NA NA
## PROVINCIAAREQUIPA 1.542e+00 NA NA NA
## PROVINCIAASCOPE 1.796e+00 NA NA NA
## PROVINCIAASUNCION -1.030e+00 NA NA NA
## PROVINCIAATALAYA -1.308e+01 NA NA NA
## PROVINCIAAYABACA -2.685e+00 NA NA NA
## PROVINCIAAYMARAES -6.395e-01 NA NA NA
## PROVINCIAAZANGARO 1.278e+00 NA NA NA
## PROVINCIABAGUA -3.435e+00 NA NA NA
## PROVINCIABARRANCA 1.648e+00 NA NA NA
## PROVINCIABELLAVISTA -3.251e+00 NA NA NA
## PROVINCIABOLIVAR -1.001e+00 NA NA NA
## PROVINCIABOLOGNESI -1.886e+00 NA NA NA
## PROVINCIABONGARA -6.562e-01 NA NA NA
## PROVINCIACAJABAMBA -2.631e+00 NA NA NA
## PROVINCIACAJAMARCA 2.482e-01 NA NA NA
## PROVINCIACAJATAMBO 4.367e-02 NA NA NA
## PROVINCIACALCA -1.055e+00 NA NA NA
## PROVINCIACALLAO 1.703e+00 NA NA NA
## PROVINCIACAMANA 7.702e-01 NA NA NA
## PROVINCIACANAS 2.498e-01 NA NA NA
## PROVINCIACANCHIS -3.344e-01 NA NA NA
## PROVINCIACANDARAVE 2.078e+00 NA NA NA
## PROVINCIACANETE 1.610e+00 NA NA NA
## PROVINCIACANGALLO 1.570e+00 NA NA NA
## PROVINCIACANTA 8.824e-01 NA NA NA
## PROVINCIACARABAYA 6.689e-01 NA NA NA
## PROVINCIACARAVELI 4.532e-01 NA NA NA
## PROVINCIACARHUAZ 1.308e-01 NA NA NA
## PROVINCIACARLOS F. FITZCARRALD -1.358e+00 NA NA NA
## PROVINCIACASMA -1.553e+00 NA NA NA
## PROVINCIACASTILLA 4.374e-01 NA NA NA
## PROVINCIACASTROVIRREYNA 1.047e-01 NA NA NA
## PROVINCIACAYLLOMA 3.746e-01 NA NA NA
## PROVINCIACELENDIN -1.593e+00 NA NA NA
## PROVINCIACHACHAPOYAS 6.573e-01 NA NA NA
## PROVINCIACHANCHAMAYO -2.622e-01 NA NA NA
## PROVINCIACHEPEN 1.228e+00 NA NA NA
## PROVINCIACHICLAYO 8.626e-01 NA NA NA
## PROVINCIACHINCHA 1.739e+00 NA NA NA
## PROVINCIACHINCHEROS -3.474e-01 NA NA NA
## PROVINCIACHOTA -2.199e+00 NA NA NA
## PROVINCIACHUCUITO 1.040e+00 NA NA NA
## PROVINCIACHUMBIVILCAS -3.934e-01 NA NA NA
## PROVINCIACHUPACA -4.826e-01 NA NA NA
## PROVINCIACHURCAMPA -1.259e-01 NA NA NA
## PROVINCIACONCEPCION 8.848e-02 NA NA NA
## PROVINCIACONDESUYOS -3.408e-01 NA NA NA
## PROVINCIACONDORCANQUI -1.197e+01 NA NA NA
## PROVINCIACONTRALMIRANTE VILLAR 1.794e+00 NA NA NA
## PROVINCIACONTUMAZA -1.197e+00 NA NA NA
## PROVINCIACORONEL PORTILLO 3.087e-01 NA NA NA
## PROVINCIACORONGO -2.323e+00 NA NA NA
## PROVINCIACOTABAMBAS -2.587e+00 NA NA NA
## PROVINCIACUSCO 1.223e-01 NA NA NA
## PROVINCIACUTERVO -1.104e+00 NA NA NA
## PROVINCIADANIEL ALCIDES CARRION 1.965e+00 NA NA NA
## PROVINCIADATEM DEL MARANON -1.606e+01 NA NA NA
## PROVINCIADOS DE MAYO -1.068e+00 NA NA NA
## PROVINCIAEL COLLAO 1.683e+00 NA NA NA
## PROVINCIAEL DORADO -6.969e-01 NA NA NA
## PROVINCIAESPINAR 3.956e-01 NA NA NA
## PROVINCIAFERRENAFE -7.932e-01 NA NA NA
## PROVINCIAGENERAL SANCHEZ CERRO 1.308e+00 NA NA NA
## PROVINCIAGRAN CHIMU -1.431e+00 NA NA NA
## PROVINCIAGRAU -2.757e+00 NA NA NA
## PROVINCIAHUACAYBAMBA -1.431e+00 NA NA NA
## PROVINCIAHUALGAYOC -2.087e+00 NA NA NA
## PROVINCIAHUALLAGA -8.829e-01 NA NA NA
## PROVINCIAHUAMALIES -3.041e+00 NA NA NA
## PROVINCIAHUAMANGA 1.594e+00 NA NA NA
## PROVINCIAHUANCA SANCOS 1.489e+00 NA NA NA
## PROVINCIAHUANCABAMBA -3.402e+00 NA NA NA
## PROVINCIAHUANCANE 7.844e-01 NA NA NA
## PROVINCIAHUANCAVELICA 1.188e+00 NA NA NA
## PROVINCIAHUANCAYO 1.220e+00 NA NA NA
## PROVINCIAHUANTA 2.338e-01 NA NA NA
## PROVINCIAHUANUCO -2.398e-01 NA NA NA
## PROVINCIAHUARAL 1.078e+00 NA NA NA
## PROVINCIAHUARAZ 5.749e-01 NA NA NA
## PROVINCIAHUARI -1.357e+00 NA NA NA
## PROVINCIAHUARMEY -7.458e-01 NA NA NA
## PROVINCIAHUAROCHIRI 9.755e-01 NA NA NA
## PROVINCIAHUAURA 1.325e+00 NA NA NA
## PROVINCIAHUAYLAS -1.440e+00 NA NA NA
## PROVINCIAHUAYTARA -7.017e-01 NA NA NA
## PROVINCIAICA 1.825e+00 NA NA NA
## PROVINCIAILO 1.855e+00 NA NA NA
## PROVINCIAISLAY 1.369e-01 NA NA NA
## PROVINCIAJAEN -2.579e+00 NA NA NA
## PROVINCIAJAUJA 1.179e-01 NA NA NA
## PROVINCIAJORGE BASADRE 1.996e+00 NA NA NA
## PROVINCIAJULCAN -2.218e-01 NA NA NA
## PROVINCIAJUNIN 9.399e-01 NA NA NA
## PROVINCIALA CONVENCION -2.095e+00 NA NA NA
## PROVINCIALA MAR 6.124e-02 NA NA NA
## PROVINCIALA UNION -1.684e+00 NA NA NA
## PROVINCIALAMAS -1.139e+00 NA NA NA
## PROVINCIALAMBAYEQUE 1.415e-01 NA NA NA
## PROVINCIALAMPA 1.422e+00 NA NA NA
## PROVINCIALAURICOCHA -2.136e+00 NA NA NA
## PROVINCIALEONCIO PRADO -2.637e+00 NA NA NA
## PROVINCIALIMA 1.793e+00 NA NA NA
## PROVINCIALORETO -7.817e+00 NA NA NA
## PROVINCIALUCANAS 9.325e-01 NA NA NA
## PROVINCIALUYA -2.095e+00 NA NA NA
## PROVINCIAMANU -4.056e+00 NA NA NA
## PROVINCIAMARANON -1.362e+00 NA NA NA
## PROVINCIAMARISCAL CACERES -1.751e+00 NA NA NA
## PROVINCIAMARISCAL LUZURIAGA -2.260e+00 NA NA NA
## PROVINCIAMARISCAL NIETO 1.956e+00 NA NA NA
## PROVINCIAMARISCAL RAMON CASTILLA -6.480e+00 NA NA NA
## PROVINCIAMAYNAS -1.557e+00 NA NA NA
## PROVINCIAMELGAR 1.129e+00 NA NA NA
## PROVINCIAMOHO 6.078e-02 NA NA NA
## PROVINCIAMORROPON -5.520e-01 NA NA NA
## PROVINCIAMOYOBAMBA -1.611e+00 NA NA NA
## PROVINCIANAZCA 1.542e+00 NA NA NA
## PROVINCIAOCROS -1.039e+00 NA NA NA
## PROVINCIAOTUZCO -1.004e+00 NA NA NA
## PROVINCIAOXAPAMPA -2.855e+00 NA NA NA
## PROVINCIAOYON 7.878e-01 NA NA NA
## PROVINCIAPACASMAYO 1.188e+00 NA NA NA
## PROVINCIAPACHITEA -3.571e+00 NA NA NA
## PROVINCIAPADRE ABAD -1.385e-01 NA NA NA
## PROVINCIAPAITA 1.512e+00 NA NA NA
## PROVINCIAPALLASCA -4.542e+00 NA NA NA
## PROVINCIAPALPA 1.382e+00 NA NA NA
## PROVINCIAPARINACOCHAS 1.068e+00 NA NA NA
## PROVINCIAPARURO -8.900e-01 NA NA NA
## PROVINCIAPASCO 2.047e+00 NA NA NA
## PROVINCIAPATAZ -2.441e+00 NA NA NA
## PROVINCIAPAUCAR DEL SARA SARA 2.443e-01 NA NA NA
## PROVINCIAPAUCARTAMBO -1.797e+00 NA NA NA
## PROVINCIAPICOTA -1.693e+00 NA NA NA
## PROVINCIAPISCO 1.251e+00 NA NA NA
## PROVINCIAPIURA 1.138e+00 NA NA NA
## PROVINCIAPOMABAMBA -2.521e+00 NA NA NA
## PROVINCIAPUERTO INCA -8.285e+00 NA NA NA
## PROVINCIAPUNO 1.587e+00 NA NA NA
## PROVINCIAPURUS -5.562e+00 NA NA NA
## PROVINCIAQUISPICANCHI -1.055e+00 NA NA NA
## PROVINCIARECUAY -1.522e+00 NA NA NA
## PROVINCIAREQUENA -5.367e+00 NA NA NA
## PROVINCIARIOJA -2.533e+00 NA NA NA
## PROVINCIARODRIGUEZ DE MENDOZA -1.027e+00 NA NA NA
## PROVINCIASAN ANTONIO DE PUTINA 6.834e-01 NA NA NA
## PROVINCIASAN IGNACIO -2.720e+00 NA NA NA
## PROVINCIASAN MARCOS -2.961e+00 NA NA NA
## PROVINCIASAN MARTIN 3.029e-01 NA NA NA
## PROVINCIASAN MIGUEL -2.112e+00 NA NA NA
## PROVINCIASAN PABLO -1.729e+00 NA NA NA
## PROVINCIASAN ROMAN 8.119e-01 NA NA NA
## PROVINCIASANCHEZ CARRION -2.065e+00 NA NA NA
## PROVINCIASANDIA 2.969e-01 NA NA NA
## PROVINCIASANTA -1.237e+00 NA NA NA
## PROVINCIASANTA CRUZ -2.723e+00 NA NA NA
## PROVINCIASANTIAGO DE CHUCO -2.818e-01 NA NA NA
## PROVINCIASATIPO -1.948e+00 NA NA NA
## PROVINCIASECHURA -2.628e-02 NA NA NA
## PROVINCIASIHUAS -2.701e+00 NA NA NA
## PROVINCIASUCRE 7.612e-01 NA NA NA
## PROVINCIASULLANA 9.331e-02 NA NA NA
## PROVINCIATACNA 2.153e+00 NA NA NA
## PROVINCIATAHUAMANU 1.827e+00 NA NA NA
## PROVINCIATALARA 1.088e+00 NA NA NA
## PROVINCIATAMBOPATA 1.695e+00 NA NA NA
## PROVINCIATARATA 1.730e+00 NA NA NA
## PROVINCIATARMA 7.093e-01 NA NA NA
## PROVINCIATAYACAJA -9.555e-01 NA NA NA
## PROVINCIATOCACHE -3.159e+00 NA NA NA
## PROVINCIATRUJILLO 1.724e+00 NA NA NA
## PROVINCIATUMBES 1.807e+00 NA NA NA
## PROVINCIAUCAYALI -7.310e+00 NA NA NA
## PROVINCIAURUBAMBA -2.678e+00 NA NA NA
## PROVINCIAUTCUBAMBA -3.277e+00 NA NA NA
## PROVINCIAVICTOR FAJARDO 5.525e-01 NA NA NA
## PROVINCIAVILCAS HUAMAN 4.725e-01 NA NA NA
## PROVINCIAVIRU -1.520e-01 NA NA NA
## PROVINCIAYAROWILCA -2.646e-01 NA NA NA
## PROVINCIAYAULI 1.277e+00 NA NA NA
## PROVINCIAYAUYOS 1.177e-01 NA NA NA
## PROVINCIAYUNGAY -2.212e+00 NA NA NA
## PROVINCIAYUNGUYO NA NA NA NA
## PROVINCIAZARUMILLA NA NA NA NA
## IDE NA NA NA NA
## salud NA NA NA NA
## educacion NA NA NA NA
## saneamiento NA NA NA NA
## electrificacion NA NA NA NA
## poblacion NA NA NA NA
## costaSI NA NA NA NA
## capitalSI NA NA NA NA
## tamanoMediana NA NA NA NA
## tamanoMuy grande NA NA NA NA
## tamanoMuy pequena NA NA NA NA
## tamanoPequena NA NA NA NA
##
## Residual standard error: NaN on 0 degrees of freedom
## Multiple R-squared: 1, Adjusted R-squared: NaN
## F-statistic: NaN on 194 and 0 DF, p-value: NA