parte 1

library(htmltab)
democracia= htmltab(doc = "https://es.wikipedia.org/wiki/%C3%8Dndice_de_democracia", 
               which ="//*/div/table[2]",
               encoding = "UTF-8")
names(democracia)
[1] "Puesto"                        "País"                          "Puntuación"                   
[4] "Proceso electoraly pluralismo" "Funcionamientodel gobierno"    "Participaciónpolítica"        
[7] "Culturapolítica"               "Derechosciviles"               "Categoría"                    

me quedo con pais y score

democracia=democracia[,c(2,3)]
names(democracia)
[1] "País"       "Puntuación"
str(democracia)
'data.frame':   167 obs. of  2 variables:
 $ País      : chr  "Noruega" " Islandia" "Suecia" " Nueva Zelanda" ...
 $ Puntuación: chr  "9.87" "9.58" "9.39" "9.26" ...
democracia$Puntuación=as.numeric(democracia$Puntuación)
str(democracia)
'data.frame':   167 obs. of  2 variables:
 $ País      : chr  "Noruega" " Islandia" "Suecia" " Nueva Zelanda" ...
 $ Puntuación: num  9.87 9.58 9.39 9.26 9.22 9.15 9.15 9.14 9.09 9.03 ...
row.names(democracia)=NULL

Historiograma

library(ggplot2)
histNum=ggplot(democracia,aes(x=Puntuación))+ geom_histogram(bins=5) 
histNum

Del gráfico anterior, podemos decir que la distribución no es simétrica por lo tanto la MEDIA no será una medida representativa (la MEDIANA será mejor opción Boxplot Esto es mejor para revisar los valores perdidos

box=ggplot(democracia,aes(y=Puntuación)) + geom_boxplot() + coord_flip()
box 

Estadigrafos

summary(democracia$Puntuación)
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  1.080   3.545   5.690   5.479   7.175   9.870 
library(ggplot2)

estadigrafos=round(as.vector(summary(democracia$Puntuación)),10)

box + scale_y_continuous(breaks = estadigrafos)

En este caso, observamos que la media es mayor a la mediana, eso significa que se tiene una asimetría positiva.

Ahora podemos confirmarla con el siguiente codigo

library(DescTools)
Skew(democracia$Puntuación,conf.level = 0.05)
       skew      lwr.ci      upr.ci 
-0.07107195 -0.08745086 -0.07241390 

distancia intercuartilica

IQR(democracia$Puntuación)
[1] 3.63

Es decir, entre el primer y tercer cuartil hay 3.63 valores. Así, podemos proponer que un atípico es aquel que está a una distancia lejana de estos valores centrales.

Ahora podemos representarlos por GINI

Gini(democracia$Puntuación,conf.level=0.95)
     gini    lwr.ci    upr.ci 
0.2316285 0.2089967 0.2574779 

Grafiquemos en la curva de Lorentz

library(ggplot2)
library(gglorenz)
ggplot(democracia,aes(x=Puntuación))+ gglorenz::stat_lorenz(color='purple') +
    geom_abline(linetype = "dashed") + coord_fixed() +
    labs(x = "% Paises ordenados por puntuación",
         y = " Puntuación",
         title = "Relación pais/puntuación",
         caption = "Fuente: Wikipedia")

Add a new chunk by clicking the Insert Chunk button on the toolbar or by pressing Ctrl+Alt+I.

When you save the notebook, an HTML file containing the code and output will be saved alongside it (click the Preview button or press Ctrl+Shift+K to preview the HTML file).

The preview shows you a rendered HTML copy of the contents of the editor. Consequently, unlike Knit, Preview does not run any R code chunks. Instead, the output of the chunk when it was last run in the editor is displayed.

LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQpwYXJ0ZSAxCmBgYHtyfQpsaWJyYXJ5KGh0bWx0YWIpCmRlbW9jcmFjaWE9IGh0bWx0YWIoZG9jID0gImh0dHBzOi8vZXMud2lraXBlZGlhLm9yZy93aWtpLyVDMyU4RG5kaWNlX2RlX2RlbW9jcmFjaWEiLCAKICAgICAgICAgICAgICAgd2hpY2ggPSIvLyovZGl2L3RhYmxlWzJdIiwKICAgICAgICAgICAgICAgZW5jb2RpbmcgPSAiVVRGLTgiKQpgYGAKCgpgYGB7cn0KbmFtZXMoZGVtb2NyYWNpYSkKYGBgCgptZSBxdWVkbyBjb24gcGFpcyB5IHNjb3JlCmBgYHtyfQpkZW1vY3JhY2lhPWRlbW9jcmFjaWFbLGMoMiwzKV0KYGBgCgpgYGB7cn0KbmFtZXMoZGVtb2NyYWNpYSkKYGBgCgpgYGB7cn0Kc3RyKGRlbW9jcmFjaWEpCmBgYApgYGB7cn0KZGVtb2NyYWNpYSRQdW50dWFjacOzbj1hcy5udW1lcmljKGRlbW9jcmFjaWEkUHVudHVhY2nDs24pCmBgYAoKYGBge3J9CnN0cihkZW1vY3JhY2lhKQpgYGAKCmBgYHtyfQpyb3cubmFtZXMoZGVtb2NyYWNpYSk9TlVMTApgYGAKCgoKSGlzdG9yaW9ncmFtYQoKYGBge3J9CmxpYnJhcnkoZ2dwbG90MikKaGlzdE51bT1nZ3Bsb3QoZGVtb2NyYWNpYSxhZXMoeD1QdW50dWFjacOzbikpKyBnZW9tX2hpc3RvZ3JhbShiaW5zPTUpIApoaXN0TnVtCgpgYGAKCkRlbCBncsOhZmljbyBhbnRlcmlvciwgcG9kZW1vcyBkZWNpciBxdWUgbGEgZGlzdHJpYnVjacOzbiBubyBlcyBzaW3DqXRyaWNhIHBvciBsbyB0YW50byBsYSBNRURJQSBubyBzZXLDoSB1bmEgbWVkaWRhIHJlcHJlc2VudGF0aXZhIChsYSBNRURJQU5BIHNlcsOhIG1lam9yIG9wY2nDs24KQm94cGxvdCBFc3RvIGVzIG1lam9yIHBhcmEgcmV2aXNhciBsb3MgdmFsb3JlcyBwZXJkaWRvcwoKCmBgYHtyfQpib3g9Z2dwbG90KGRlbW9jcmFjaWEsYWVzKHk9UHVudHVhY2nDs24pKSArIGdlb21fYm94cGxvdCgpICsgY29vcmRfZmxpcCgpCmJveCAKYGBgCgoKRXN0YWRpZ3JhZm9zCgpgYGB7cn0Kc3VtbWFyeShkZW1vY3JhY2lhJFB1bnR1YWNpw7NuKQpgYGAKCmBgYHtyfQpsaWJyYXJ5KGdncGxvdDIpCgplc3RhZGlncmFmb3M9cm91bmQoYXMudmVjdG9yKHN1bW1hcnkoZGVtb2NyYWNpYSRQdW50dWFjacOzbikpLDEwKQoKYm94ICsgc2NhbGVfeV9jb250aW51b3VzKGJyZWFrcyA9IGVzdGFkaWdyYWZvcykKYGBgCgpFbiBlc3RlIGNhc28sIG9ic2VydmFtb3MgcXVlIGxhIG1lZGlhIGVzIG1heW9yIGEgbGEgbWVkaWFuYSwgZXNvIHNpZ25pZmljYSBxdWUgc2UgdGllbmUgdW5hIGFzaW1ldHLDrWEgcG9zaXRpdmEuCgpBaG9yYSBwb2RlbW9zIGNvbmZpcm1hcmxhIGNvbiBlbCBzaWd1aWVudGUgY29kaWdvCgoKYGBge3J9CmxpYnJhcnkoRGVzY1Rvb2xzKQpTa2V3KGRlbW9jcmFjaWEkUHVudHVhY2nDs24sY29uZi5sZXZlbCA9IDAuMDUpCmBgYAoKZGlzdGFuY2lhIGludGVyY3VhcnRpbGljYQoKYGBge3J9CklRUihkZW1vY3JhY2lhJFB1bnR1YWNpw7NuKQpgYGAKRXMgZGVjaXIsIGVudHJlIGVsIHByaW1lciB5IHRlcmNlciBjdWFydGlsIGhheSAzLjYzIHZhbG9yZXMuIEFzw60sIHBvZGVtb3MgcHJvcG9uZXIgcXVlIHVuIGF0w61waWNvIGVzIGFxdWVsIHF1ZSBlc3TDoSBhIHVuYSBkaXN0YW5jaWEgbGVqYW5hIGRlIGVzdG9zIHZhbG9yZXMgY2VudHJhbGVzLgoKQWhvcmEgcG9kZW1vcyByZXByZXNlbnRhcmxvcyBwb3IgR0lOSQoKYGBge3J9CkdpbmkoZGVtb2NyYWNpYSRQdW50dWFjacOzbixjb25mLmxldmVsPTAuOTUpCmBgYApHcmFmaXF1ZW1vcyBlbiBsYSBjdXJ2YSBkZSBMb3JlbnR6CgoKCgpgYGB7cn0KbGlicmFyeShnZ3Bsb3QyKQpsaWJyYXJ5KGdnbG9yZW56KQpgYGAKCmBgYHtyfQpnZ3Bsb3QoZGVtb2NyYWNpYSxhZXMoeD1QdW50dWFjacOzbikpKyBnZ2xvcmVuejo6c3RhdF9sb3JlbnooY29sb3I9J3B1cnBsZScpICsKICAgIGdlb21fYWJsaW5lKGxpbmV0eXBlID0gImRhc2hlZCIpICsgY29vcmRfZml4ZWQoKSArCiAgICBsYWJzKHggPSAiJSBQYWlzZXMgb3JkZW5hZG9zIHBvciBwdW50dWFjacOzbiIsCiAgICAgICAgIHkgPSAiIFB1bnR1YWNpw7NuIiwKICAgICAgICAgdGl0bGUgPSAiUmVsYWNpw7NuIHBhaXMvcHVudHVhY2nDs24iLAogICAgICAgICBjYXB0aW9uID0gIkZ1ZW50ZTogV2lraXBlZGlhIikKYGBgCgoKCgoKCgoKCgoKCgoKCgoKCkFkZCBhIG5ldyBjaHVuayBieSBjbGlja2luZyB0aGUgKkluc2VydCBDaHVuayogYnV0dG9uIG9uIHRoZSB0b29sYmFyIG9yIGJ5IHByZXNzaW5nICpDdHJsK0FsdCtJKi4KCldoZW4geW91IHNhdmUgdGhlIG5vdGVib29rLCBhbiBIVE1MIGZpbGUgY29udGFpbmluZyB0aGUgY29kZSBhbmQgb3V0cHV0IHdpbGwgYmUgc2F2ZWQgYWxvbmdzaWRlIGl0IChjbGljayB0aGUgKlByZXZpZXcqIGJ1dHRvbiBvciBwcmVzcyAqQ3RybCtTaGlmdCtLKiB0byBwcmV2aWV3IHRoZSBIVE1MIGZpbGUpLgoKVGhlIHByZXZpZXcgc2hvd3MgeW91IGEgcmVuZGVyZWQgSFRNTCBjb3B5IG9mIHRoZSBjb250ZW50cyBvZiB0aGUgZWRpdG9yLiBDb25zZXF1ZW50bHksIHVubGlrZSAqS25pdCosICpQcmV2aWV3KiBkb2VzIG5vdCBydW4gYW55IFIgY29kZSBjaHVua3MuIEluc3RlYWQsIHRoZSBvdXRwdXQgb2YgdGhlIGNodW5rIHdoZW4gaXQgd2FzIGxhc3QgcnVuIGluIHRoZSBlZGl0b3IgaXMgZGlzcGxheWVkLgo=