Parte 2.

Para esta parte descargue la información de cada país sobre el índice de felicidad

library(htmltab)

link="https://es.wikipedia.org/wiki/%C3%8Dndice_global_de_felicidad"

path='//*/div/table/tbody'

indice=htmltab(doc = link,which = path)

summary(indice$Puntuación)

##    Length     Class      Mode 
##       156 character character

names(indice)

## [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"

names(indice)=c("N","pais","puntuacion","PBI","apoyosoc","esperanzavida","libertadtomadecisiones","generosidad","percepcioncorrupcion")
names(indice)

## [1] "N"                      "pais"                  
## [3] "puntuacion"             "PBI"                   
## [5] "apoyosoc"               "esperanzavida"         
## [7] "libertadtomadecisiones" "generosidad"           
## [9] "percepcioncorrupcion"

str(indice)

## 'data.frame':    156 obs. of  9 variables:
##  $ N                     : chr  "1" "2" "3" "4" ...
##  $ pais                  : chr  "Finlandia" "Noruega" "Dinamarca" "Islandia" ...
##  $ puntuacion            : chr  "7.633" "7.594" "7.555" "7.495" ...
##  $ PBI                   : chr  "1.305" "1.456" "1.351" "1.343" ...
##  $ apoyosoc              : chr  "1.592" "1.582" "1.590" "1.644" ...
##  $ esperanzavida         : chr  "0.874" "0.861" "0.868" "0.914" ...
##  $ libertadtomadecisiones: chr  "0.681" "0.686" "0.683" "0.677" ...
##  $ generosidad           : chr  "0.192" "0.286" "0.284" "0.353" ...
##  $ percepcioncorrupcion  : chr  "0.393" "0.340" "0.408" "0.138" ...

indice$puntuacion=as.numeric(indice$puntuacion)

str(indice)

## 'data.frame':    156 obs. of  9 variables:
##  $ N                     : chr  "1" "2" "3" "4" ...
##  $ pais                  : chr  "Finlandia" "Noruega" "Dinamarca" "Islandia" ...
##  $ puntuacion            : num  7.63 7.59 7.55 7.5 7.49 ...
##  $ PBI                   : chr  "1.305" "1.456" "1.351" "1.343" ...
##  $ apoyosoc              : chr  "1.592" "1.582" "1.590" "1.644" ...
##  $ esperanzavida         : chr  "0.874" "0.861" "0.868" "0.914" ...
##  $ libertadtomadecisiones: chr  "0.681" "0.686" "0.683" "0.677" ...
##  $ generosidad           : chr  "0.192" "0.286" "0.284" "0.353" ...
##  $ percepcioncorrupcion  : chr  "0.393" "0.340" "0.408" "0.138" ...

La puntuación de la felicidad se distribuye asimetricamente?

library(ggplot2)
base=ggplot(indice,aes(x=puntuacion))
hist=base+geom_histogram(bins = 7)
hist

library(ggplot2)
basee=ggplot(indice,aes(y=puntuacion))
box=basee+geom_boxplot()+coord_flip()
box

Segun el grafico de boxplot, no figuran atipicos, se nota una ligera simetria

¿El valor representativo de indice de felicidad es robusto?

summary(indice$puntuacion)

##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   2.905   4.454   5.378   5.376   6.168   7.633

basen=ggplot(data=indice,aes(x=puntuacion))
basen+geom_histogram(bins = 10)

¿Puedes concluir que hay mucha desigualdad en este indice?

No hay desigualdad notoria

library(DescTools)
Gini(indice$puntuacion)

## [1] 0.1196207

La desigualdad no es notoria ¿Se puede afirmar que aproximadamente el 80% de los paises del mundo suman el 75% del puntaje acumulado de felicidad?

library(ggplot2)
library(gglorenz)

## Registered S3 methods overwritten by 'ineq':
##   method   from     
##   plot.Lc  DescTools
##   lines.Lc DescTools

ggplot(indice,aes(x=puntuacion))+gglorenz::stat_lorenz(color='red')+geom_abline(linetype="dashed")+coord_fixed()+labs(x="% Paises Indice felicidad", y="% Acumulado puntuacion",title = "Relacion pais/indice felicidad",caption = "Fuente:Wikipedia")+scale_y_continuous(breaks = seq(0,1,0.15))+scale_x_continuous(breaks = seq(0,1,0.2))

¿Hay atipicos en el puntaje de felicidad?

Por el boxplot notamos que no hay atipicos para el puntaje de felicidad

# Cuartiles

q3=quantile(indice$puntuacion,0.75)
q1=quantile(indice$puntuacion,0.25)

#Umbrales calculando

umbralA=q3+1.5*IQR(indice$puntuacion)
umbralB=q1-1.5*IQR(indice$puntuacion)