Estadisticas - Recibo de Luz

Exportar base de datos

#file.choose()
poblacion <- read.csv ("/Users/ricardogc/Desktop/R - Analisis de datos para la toma de decisiones. /poblacion.csv")

#file.choose()
muestra <- read.csv ("/Users/ricardogc/Desktop/R - Analisis de datos para la toma de decisiones. /muestra.csv")

Tamaño de la población (N)

N <- length(poblacion$Pago)
N
## [1] 12

Tamaño de la muestra (n)

n <- length(muestra$Pago)
n
## [1] 5

Medidas de tendencia Central: Permite conocer el valor al que tiende el conjunto de datos

Media o Promedio: Valor que se obtiene al sumar todos los datos y dividirlos entre la cantidad total de datos.

Media poblacional (x barra)
media_poblacional <- mean(poblacion$Pago)
media_poblacional
## [1] 245.0167
Media muestral
media_muestral <- mean(muestra$Pago)
media_muestral
## [1] 249.432

####Mediana: Valor que ocupa el valor central de todos los datos cuando estos estan ordenados de menor a mayor ##### Mediana poblacional

mediana_poblacional <- median(poblacion$Pago)
mediana_poblacional
## [1] 228.63
Mediana muestral
   mediana_muestral <- median(muestra$Pago)
   mediana_muestral
## [1] 230.46
  #Moda: Valor que aparece con mayor frecuencia en un conjunto de datos
   #Funcion para calcular la Moda
   mode <-  function(x) {
     ux <- unique(x) 
     ux[which.max(tabulate(match(x,ux)))] 
     }

Nota: Si ningun dato se repite la función coloca el primer valor en lugar de marcar error

Moda poblacional
moda_poblacional <- mode(poblacion$Pago)
   moda_poblacional
## [1] 266.63
Moda muestral
moda_muestral <- mode(muestra$Pago)
   moda_muestral
## [1] 266.63

Relación entre la media, mediana y moda

Si la media = mediana = moda, los datos tienen una distribución SIMETRICA
Si la media < mediana < moda , los datos tienen SESGO NEGATIVO

Histograma de Población (X- población) (Y - Frecuencia)

La población tiene un sesgo positivo
hist(poblacion$Pago)

Histograma de Muestra (X- Muestra) (Y - Frecuencia)

hist(muestra$Pago) 

Medidas de dispersión: Miden que tan esparcidos se encuentran los datos

Rango: intervalo o diferencia entre el valor maximo y el minimo de un conjunto de datos

Rango poblacional
rango_poblacional <- max(poblacion$Pago) - min(poblacion$Pago)
   rango_poblacional
## [1] 180.86
    r <-    range(poblacion$Pago)
    r
## [1] 162.64 343.50
La función range() devuelve el valor minimo y maximo pero no su diferencia, que es el valor que buscamos
Rango muestral
rango_muestral <-  max(muestra$Pago) - min(muestra$Pago)
rango_muestral
## [1] 156.34

Varianza: Es el promedio elevado al cuadrado de las desviaciones individuales de cada observacion con respecto a la media de una distribución.

Nota: Si es población se divide entre N: Si es muestra se divide entre n-1.

Varianza Poblacional (Sigma cuadrada)
 varianza_poblacional <- var(poblacion$Pago)*(N-1)/N
    varianza_poblacional
## [1] 3614.659
Varianza Muestral (Sigma cuadrada)
varianza_muestral <- var(muestra$Pago)
    varianza_muestral
## [1] 3570.905

Desviación estandar: Raiz cuadrada de la varianza

Desviación estandar poblacional (sigma)
 desviacion_estandar_poblacional <- sqrt(varianza_poblacional)
    desviacion_estandar_poblacional
## [1] 60.12203
Desviación estandar muestral (sigma)
desviacion_estandar_muestral <- sqrt(varianza_muestral)
    desviacion_estandar_muestral
## [1] 59.75705

Conclusiones

Para obtener información más especifica y Acertiva, bien se sabe que Rstudio es una excelente herramienta que nos brinda un estudio poblaciona con estadisticas automaticas que genera una mejor respuesta a la obtención de información. Es por esto que Rstudio es muy recomendado para las estadisticas para evitar errores que podrían ocasionar mala interpretación de los datos y además de que R studio genera de manera automatica las formulas o numero con simples comandos sencillos de utilizar y poner en practica.

LS0tCnRpdGxlOiA8c3BhbiBzdHlsZT0iQ29sb3I6R3JleSI+ICJFc3RhZMOtc3RpY2FzIgphdXRob3I6ICJSaWNhcmRvIEdhbGljaWEgLSBBMDE2NTMyNzkiCmRhdGU6ICcyMDIyLTA5LTE5JwpvdXRwdXQ6IAogIGh0bWxfZG9jdW1lbnQ6CiAgICB0b2M6IHRydWUKICAgIHRvY19mbG9hdDogdHJ1ZQogICAgY29kZV9kb3dubG9hZDogdHJ1ZSAKLS0tCiMjICoqRXN0YWRpc3RpY2FzIC0gUmVjaWJvIGRlIEx1eioqCgohW10oL1VzZXJzL3JpY2FyZG9nYy9EZXNrdG9wL0FjdGl2aWRhZCAzLjIvRXN0YWRpc3RpY2FzLWRlLXZlbnRhcy5QTkcpCgojIyAqRXhwb3J0YXIgYmFzZSBkZSBkYXRvcyoKYGBge3J9CiNmaWxlLmNob29zZSgpCnBvYmxhY2lvbiA8LSByZWFkLmNzdiAoIi9Vc2Vycy9yaWNhcmRvZ2MvRGVza3RvcC9SIC0gQW5hbGlzaXMgZGUgZGF0b3MgcGFyYSBsYSB0b21hIGRlIGRlY2lzaW9uZXMuIC9wb2JsYWNpb24uY3N2IikKCiNmaWxlLmNob29zZSgpCm11ZXN0cmEgPC0gcmVhZC5jc3YgKCIvVXNlcnMvcmljYXJkb2djL0Rlc2t0b3AvUiAtIEFuYWxpc2lzIGRlIGRhdG9zIHBhcmEgbGEgdG9tYSBkZSBkZWNpc2lvbmVzLiAvbXVlc3RyYS5jc3YiKQpgYGAKCiMjIFRhbWHDsW8gZGUgbGEgcG9ibGFjacOzbiAoTikKYGBge3J9Ck4gPC0gbGVuZ3RoKHBvYmxhY2lvbiRQYWdvKQpOCmBgYAoKIyMgVGFtYcOxbyBkZSBsYSBtdWVzdHJhIChuKQpgYGB7cn0KbiA8LSBsZW5ndGgobXVlc3RyYSRQYWdvKQpuCmBgYAoKIyMgTWVkaWRhcyBkZSB0ZW5kZW5jaWEgQ2VudHJhbDogUGVybWl0ZSBjb25vY2VyIGVsIHZhbG9yIGFsIHF1ZSB0aWVuZGUgZWwgY29uanVudG8gZGUgZGF0b3MgCiMjIyMgTWVkaWEgbyBQcm9tZWRpbzogVmFsb3IgcXVlIHNlIG9idGllbmUgYWwgc3VtYXIgdG9kb3MgbG9zIGRhdG9zIHkgZGl2aWRpcmxvcyBlbnRyZSBsYSBjYW50aWRhZCB0b3RhbCBkZSBkYXRvcy4gCiMjIyMjIE1lZGlhIHBvYmxhY2lvbmFsICh4IGJhcnJhKQpgYGB7cn0KbWVkaWFfcG9ibGFjaW9uYWwgPC0gbWVhbihwb2JsYWNpb24kUGFnbykKbWVkaWFfcG9ibGFjaW9uYWwKYGBgCgojIyMjIyBNZWRpYSBtdWVzdHJhbApgYGB7cn0KbWVkaWFfbXVlc3RyYWwgPC0gbWVhbihtdWVzdHJhJFBhZ28pCm1lZGlhX211ZXN0cmFsCmBgYAoKIyMjI01lZGlhbmE6IFZhbG9yIHF1ZSBvY3VwYSBlbCB2YWxvciBjZW50cmFsIGRlIHRvZG9zIGxvcyBkYXRvcyBjdWFuZG8gZXN0b3MgZXN0YW4gb3JkZW5hZG9zIGRlIG1lbm9yIGEgbWF5b3IKIyMjIyMgTWVkaWFuYSBwb2JsYWNpb25hbApgYGB7cn0KbWVkaWFuYV9wb2JsYWNpb25hbCA8LSBtZWRpYW4ocG9ibGFjaW9uJFBhZ28pCm1lZGlhbmFfcG9ibGFjaW9uYWwKYGBgCgojIyMjIyBNZWRpYW5hIG11ZXN0cmFsCmBgYHtyfQogICBtZWRpYW5hX211ZXN0cmFsIDwtIG1lZGlhbihtdWVzdHJhJFBhZ28pCiAgIG1lZGlhbmFfbXVlc3RyYWwKICAjTW9kYTogVmFsb3IgcXVlIGFwYXJlY2UgY29uIG1heW9yIGZyZWN1ZW5jaWEgZW4gdW4gY29uanVudG8gZGUgZGF0b3MKICAgI0Z1bmNpb24gcGFyYSBjYWxjdWxhciBsYSBNb2RhCiAgIG1vZGUgPC0gIGZ1bmN0aW9uKHgpIHsKICAgICB1eCA8LSB1bmlxdWUoeCkgCiAgICAgdXhbd2hpY2gubWF4KHRhYnVsYXRlKG1hdGNoKHgsdXgpKSldIAogICAgIH0KYGBgCgojIyAqTm90YTogU2kgbmluZ3VuIGRhdG8gc2UgcmVwaXRlIGxhIGZ1bmNpw7NuIGNvbG9jYSBlbCBwcmltZXIgdmFsb3IgZW4gbHVnYXIgZGUgbWFyY2FyIGVycm9yKgojIyMjIyBNb2RhIHBvYmxhY2lvbmFsCmBgYHtyfQptb2RhX3BvYmxhY2lvbmFsIDwtIG1vZGUocG9ibGFjaW9uJFBhZ28pCiAgIG1vZGFfcG9ibGFjaW9uYWwKYGBgCiAgCiMjIyMjIE1vZGEgbXVlc3RyYWwKYGBge3J9Cm1vZGFfbXVlc3RyYWwgPC0gbW9kZShtdWVzdHJhJFBhZ28pCiAgIG1vZGFfbXVlc3RyYWwKYGBgCiAgIAojIyAqUmVsYWNpw7NuIGVudHJlIGxhIG1lZGlhLCBtZWRpYW5hIHkgbW9kYSoKIyMjIyMgU2kgbGEgbWVkaWEgPSBtZWRpYW5hID0gbW9kYSwgbG9zIGRhdG9zIHRpZW5lbiB1bmEgZGlzdHJpYnVjacOzbiBTSU1FVFJJQ0EKIyMjIyMgU2kgbGEgbWVkaWEgPCBtZWRpYW5hIDwgbW9kYSAsIGxvcyBkYXRvcyB0aWVuZW4gU0VTR08gTkVHQVRJVk8gICAKICAgCiMjICoqSGlzdG9ncmFtYSBkZSBQb2JsYWNpw7NuIChYLSBwb2JsYWNpw7NuKSAoWSAtIEZyZWN1ZW5jaWEpKioKIyMjIyMgTGEgcG9ibGFjacOzbiB0aWVuZSB1biBzZXNnbyBwb3NpdGl2bwpgYGB7cn0KaGlzdChwb2JsYWNpb24kUGFnbykKYGBgCgojIyAqKkhpc3RvZ3JhbWEgZGUgTXVlc3RyYSAoWC0gTXVlc3RyYSkgKFkgLSBGcmVjdWVuY2lhKSoqCmBgYHtyfQpoaXN0KG11ZXN0cmEkUGFnbykgCmBgYAoKIyMgKk1lZGlkYXMgZGUgZGlzcGVyc2nDs246IE1pZGVuIHF1ZSB0YW4gZXNwYXJjaWRvcyBzZSBlbmN1ZW50cmFuIGxvcyBkYXRvcyoKIyMjIyBSYW5nbzogaW50ZXJ2YWxvIG8gZGlmZXJlbmNpYSBlbnRyZSBlbCB2YWxvciBtYXhpbW8geSBlbCBtaW5pbW8gZGUgdW4gY29uanVudG8gZGUgZGF0b3MKIyMjIyMgUmFuZ28gcG9ibGFjaW9uYWwKYGBge3J9CnJhbmdvX3BvYmxhY2lvbmFsIDwtIG1heChwb2JsYWNpb24kUGFnbykgLSBtaW4ocG9ibGFjaW9uJFBhZ28pCiAgIHJhbmdvX3BvYmxhY2lvbmFsCiAgIAogICAgciA8LSAgICByYW5nZShwb2JsYWNpb24kUGFnbykKICAgIHIKYGBgCgojIyMjIyAqTGEgZnVuY2nDs24gcmFuZ2UoKSBkZXZ1ZWx2ZSBlbCB2YWxvciBtaW5pbW8geSBtYXhpbW8gcGVybyBubyBzdSBkaWZlcmVuY2lhLCBxdWUgZXMgZWwgdmFsb3IgcXVlIGJ1c2NhbW9zKgojIyMjIyBSYW5nbyBtdWVzdHJhbApgYGB7cn0KcmFuZ29fbXVlc3RyYWwgPC0gIG1heChtdWVzdHJhJFBhZ28pIC0gbWluKG11ZXN0cmEkUGFnbykKcmFuZ29fbXVlc3RyYWwKYGBgCgojIyMjIFZhcmlhbnphOiBFcyBlbCBwcm9tZWRpbyBlbGV2YWRvIGFsIGN1YWRyYWRvIGRlIGxhcyBkZXN2aWFjaW9uZXMgaW5kaXZpZHVhbGVzIGRlIGNhZGEgb2JzZXJ2YWNpb24gY29uIHJlc3BlY3RvIGEgbGEgbWVkaWEgZGUgdW5hIGRpc3RyaWJ1Y2nDs24uIAojIyMjICpOb3RhOiBTaSBlcyBwb2JsYWNpw7NuIHNlIGRpdmlkZSBlbnRyZSBOOiBTaSBlcyBtdWVzdHJhIHNlIGRpdmlkZSBlbnRyZSBuLTEuKgogICAgCiMjIyMjIFZhcmlhbnphIFBvYmxhY2lvbmFsIChTaWdtYSBjdWFkcmFkYSkKYGBge3J9CiB2YXJpYW56YV9wb2JsYWNpb25hbCA8LSB2YXIocG9ibGFjaW9uJFBhZ28pKihOLTEpL04KICAgIHZhcmlhbnphX3BvYmxhY2lvbmFsCmBgYAogICAKIyMjIyMgVmFyaWFuemEgTXVlc3RyYWwgKFNpZ21hIGN1YWRyYWRhKQpgYGB7cn0KdmFyaWFuemFfbXVlc3RyYWwgPC0gdmFyKG11ZXN0cmEkUGFnbykKICAgIHZhcmlhbnphX211ZXN0cmFsCmBgYAogIAogICAgCiMjIyMgRGVzdmlhY2nDs24gZXN0YW5kYXI6IFJhaXogY3VhZHJhZGEgZGUgbGEgdmFyaWFuemEKIyMjIyMgRGVzdmlhY2nDs24gZXN0YW5kYXIgcG9ibGFjaW9uYWwgKHNpZ21hKQpgYGB7cn0KIGRlc3ZpYWNpb25fZXN0YW5kYXJfcG9ibGFjaW9uYWwgPC0gc3FydCh2YXJpYW56YV9wb2JsYWNpb25hbCkKICAgIGRlc3ZpYWNpb25fZXN0YW5kYXJfcG9ibGFjaW9uYWwKYGBgCiAgIAojIyMjIyBEZXN2aWFjacOzbiBlc3RhbmRhciBtdWVzdHJhbCAoc2lnbWEpCmBgYHtyfQpkZXN2aWFjaW9uX2VzdGFuZGFyX211ZXN0cmFsIDwtIHNxcnQodmFyaWFuemFfbXVlc3RyYWwpCiAgICBkZXN2aWFjaW9uX2VzdGFuZGFyX211ZXN0cmFsCmBgYAogICAgCiMjICoqQ29uY2x1c2lvbmVzKioKClBhcmEgb2J0ZW5lciBpbmZvcm1hY2nDs24gbcOhcyBlc3BlY2lmaWNhIHkgQWNlcnRpdmEsIGJpZW4gc2Ugc2FiZSBxdWUgUnN0dWRpbyBlcyB1bmEgZXhjZWxlbnRlIGhlcnJhbWllbnRhIHF1ZSBub3MgYnJpbmRhIHVuIGVzdHVkaW8gcG9ibGFjaW9uYSBjb24gZXN0YWRpc3RpY2FzIGF1dG9tYXRpY2FzIHF1ZSBnZW5lcmEgdW5hIG1lam9yIHJlc3B1ZXN0YSBhIGxhIG9idGVuY2nDs24gZGUgaW5mb3JtYWNpw7NuLiBFcyBwb3IgZXN0byBxdWUgUnN0dWRpbyBlcyBtdXkgcmVjb21lbmRhZG8gcGFyYSBsYXMgZXN0YWRpc3RpY2FzIHBhcmEgZXZpdGFyIGVycm9yZXMgcXVlIHBvZHLDrWFuIG9jYXNpb25hciBtYWxhIGludGVycHJldGFjacOzbiBkZSBsb3MgZGF0b3MgeSBhZGVtw6FzIGRlIHF1ZSBSIHN0dWRpbyBnZW5lcmEgZGUgbWFuZXJhIGF1dG9tYXRpY2EgbGFzIGZvcm11bGFzIG8gbnVtZXJvIGNvbiBzaW1wbGVzIGNvbWFuZG9zIHNlbmNpbGxvcyBkZSB1dGlsaXphciB5IHBvbmVyIGVuIHByYWN0aWNhLiAKCiAgICAgICAg