Indique su nombre y cédula

Nombre: David Pizarro

Cédula: 31.187.123

Estadística I, Tarea 01

A continuación hay una serie de ejercicios que debe completar. Algunos de los códigos ya están vaciados para que se facilite realizar estar tarea.

Se recomienda para cada uno de los conjuntos de datos a trabajar, revise los datos que contienen, así como también usar las funciones de ayuda ? obtener información sobre los mismos y de esta manera lograr facilitar la comprensión sobre los datos a trabajar.

Recomendaciones

Usar funciones apropiadas para obtener los valores solicitados. Revise los ejercicios de codificación anteriormente presentados en el curso, las láminas, así como el libro Open Intro Stats.

Si realiza algún cálculo manualmente deberá añadir las celdas de código con tales cálculo, tantas como sea necesario.

En cada una de las celdas de código, añada comentarios con el símbolo numeral (#) al inicio de la celda para indicar los procedimientos realizados.

Recuerde cuidar la redacción y ortografía en la información que presente.

En caso de tener alguna duda sobre los procedimientos a seguir, usar el grupo de whatsapp.

Se recuerda que la terea debe ser realizada de forma individual.

# instalar y cargar paquetes con los que se trabajará. No modificar esta celda
install.packages ('openintro')
install.packages ('gtsummary')
library ('openintro')
library ('gtsummary')
library ('dplyr')
Installing package into ‘/usr/local/lib/R/site-library’
(as ‘lib’ is unspecified)

also installing the dependencies ‘airports’, ‘cherryblossom’, ‘usdata’


Installing package into ‘/usr/local/lib/R/site-library’
(as ‘lib’ is unspecified)

also installing the dependencies ‘V8’, ‘litedown’, ‘reactR’, ‘bigD’, ‘bitops’, ‘juicyjuice’, ‘markdown’, ‘reactable’, ‘cards’, ‘gt’


Loading required package: airports

Loading required package: cherryblossom

Loading required package: usdata


Attaching package: ‘dplyr’


The following objects are masked from ‘package:stats’:

    filter, lag


The following objects are masked from ‘package:base’:

    intersect, setdiff, setequal, union

1) Tratamiento de tumores con etanol

Nombre del conjunto de datos

ethanol

Descripción

Experimento en el que se probaron 3 tratamientos diferentes de etanol para el tratamiento de tumores de cáncer oral en hámsters.

La columna “regress” especifica si el tratamiento tuvo efecto con la categoría “yes” o “no” en caso de no tener algún efecto sobre el tumor.

La columna “treatment” indica el tratamiento que recibió el hámster.

Detalles

Los tratamientos de etil_celulosa y etanol_puro consistían en aproximadamente una cuarta parte del volumen de los tumores, mientras que el tratamiento de etanol_puro_16x era 16 veces mayor, es decir, unas 4 veces el tamaño de los tumores.

Fuente

Morhard R, et al. 2017. Desarrollo de la ablación mejorada con etanol como alternativa a la cirugía en el tratamiento de tumores sólidos superficiales.

Asignación

  1. ¿Cómo se presentan los datos?
  2. ¿Cuántas observaciones tiene el conjunto de datos?
  3. ¿Cuántas variables tiene el conjunto de datos y cómo se llaman?
  4. Crear una tabla de contigencia que represente los totales por tratamiento. Calcular la frecuencia relativa de cada uno de estos experimentos y su éxito (puede usar los códigos revisados o hacerlo manualmente).
  5. Realice una explicación escrita sobre los resultados obtenidos (tipo de estudio, siginficado de los resultados, generalizable, etc.). Puede llegar a alguna conclusión con los datos presentados? Sea lo más explícito posible.
ethanol
A tibble: 24 × 2
treatment regress
<fct> <fct>
pure_ethanol_16x yes
pure_ethanol_16x yes
pure_ethanol_16x yes
pure_ethanol_16x yes
pure_ethanol_16x no
pure_ethanol_16x no
pure_ethanol_16x no
pure_ethanol_16x no
pure_ethanol_16x no
pure_ethanol_16x no
pure_ethanol_16x no
pure_ethanol_16x no
pure_ethanol no
pure_ethanol no
pure_ethanol no
pure_ethanol no
pure_ethanol no
ethyl_cellulose yes
ethyl_cellulose yes
ethyl_cellulose yes
ethyl_cellulose yes
ethyl_cellulose yes
ethyl_cellulose yes
ethyl_cellulose no

Respuesta 1.1:

Los datos se presentan de forma Categórica/Nominales, tanto la columna de “treatment” como la de “regress”.

Observación: los datos se presentan en dos columnas donde cada una corresponde a una variable

Respuesta 1.2:

El conjunto de datos contiene 24 observaciones.

Respuesta 1.3: El conjunto de datos contiene 3 variables:

  1. ethyl_cellulose
  2. pure_ethanol
  3. pure_ethanol_16x

Observación: no, sólo contiene dos variables que son las que se encuentran en cada columna. Lo que indica anteriormente son el tipo de tratamiento. Es equivalente a que una variable fuese el nombre de una país y luego indique que existen 185 variables, una por cada país, siendo lo correcto que la variable país puede asumir 1 de 185 posibles valores.

#Respuesta 1.4:
table(ethanol)
                  regress
treatment          no yes
  ethyl_cellulose   1   6
  pure_ethanol      5   0
  pure_ethanol_16x  8   4
tbl_summary(ethanol)%>%
            as_kable()


|**Characteristic** | **N = 24** |
|:------------------|:----------:|
|treatment          |            |
|ethyl_cellulose    |  7 (29%)   |
|pure_ethanol       |  5 (21%)   |
|pure_ethanol_16x   |  12 (50%)  |
|regress            |  10 (42%)  |
gtsummary::tbl_summary(ethanol,
                       include = 'treatment',
                       by='regress',
                       percent='row')%>%
                          as_kable()


|**Characteristic** | **no**  N = 14 | **yes**  N = 10 |
|:------------------|:--------------:|:---------------:|
|treatment          |                |                 |
|ethyl_cellulose    |    1 (14%)     |     6 (86%)     |
|pure_ethanol       |    5 (100%)    |     0 (0%)      |
|pure_ethanol_16x   |    8 (67%)     |     4 (33%)     |

observacion: si se quiere verificar para cada tratamiento, en términos relativos su efectividad, este sería el código

gtsummary::tbl_summary(ethanol,
                       include = 'treatment',
                       by='regress',
                       percent='column')%>%
                          as_kable()

Respuesta 1.5:

Tipo de estudio: Experimental

Los resultados nos muestran que hubo una mayor cantidad negativa de exito en las pruebas realizadas con las tres variables.

observación: es absolutamente incoherente lo que está indicando. Sino se obtuvo un resultado positivo, no se puede decir “una mayor cantidad negativa de éxito”. No entiendo. Adicionalmente, el análisis debe ser por tipo de tratamiento

2) Indicadores Económicos y de Sostenibilidad para Sudáfrica.

Nombre del Conjunto de Datos

sa_gdp_elec

Descripción

Incluye datos anuales sobre pib, gni, emisiones de co2, costes de puesta en marcha.

Detalles

año: Año de recogida de los datos.

access_elec: Acceso a la electricidad en porcentaje de la población.

puesta en marcha: Coste de los procedimientos de creación de empresas en porcentaje de la GNI.

co2: Emisión de CO2 en kt (kilotoneladas).

gdp: PIB per cápita, PPA en dólares internacionales constantes de 2017.

gni: GNI per cápita, PPA en dólares internacionales constantes de 2017. GNI (gross national income) es la cantidad total de dinero que ganan los ciudadanos y las empresas de un país.

co2_kg_ppp: kg por PIB en dólares PPA de 2017.

Fuente

World Bank I

Carbon Dioxide Information Analysis Center, Environmental Sciences Division

Asignación

  1. ¿Cuántas observaciones tiene el conjunto de datos?

  2. Indique el tipo de cada una de ls variables.

  3. Realice un histograma para la variable co2 e indique cuántas barras recomienda aplicar para visualizar este gráfico.

  4. Explique si tiene sentido crear un histograma para la variable year

  5. Muestre dos gráficos de dispersión entre variables que pueden estar asociadas, el tipo de relación que presenten e intente identificar cuáles pueden ser las variables explicativas y de respuesta

  6. ¿Puede indentificar algunas variables que resulten independientes?. Justifique

  7. Identifique el tipo de estudio que se revisó en esta sección

sa_gdp_elec
A tibble: 16 × 7
year access_elec startup co2 gdp gni co2_kg_ppp
<dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
2003 78.8 9.4 353110 10646.83 10388.18 0.7098952
2004 80.9 9.1 379970 10997.01 10807.17 0.7306183
2005 80.8 8.6 377640 11434.92 11233.35 0.6897401
2006 80.7 6.9 379560 11924.07 11711.35 0.6564602
2007 82.0 7.1 396900 12402.05 12010.29 0.6515254
2008 81.9 6.0 426560 12628.20 12251.73 0.6785602
2009 82.6 5.9 404020 12262.14 11999.96 0.6527439
2010 82.9 6.0 425110 12452.34 12195.42 0.6665559
2011 83.6 0.3 409120 12666.71 12346.55 0.6210867
2012 85.3 0.3 426710 12743.94 12404.11 0.6337627
2013 85.2 0.3 436870 12852.78 12523.90 0.6331184
2014 86.0 0.3 447980 12884.48 12549.90 0.6374456
2015 85.3 0.3 424880 12840.04 12528.09 0.5974440
2016 84.2 0.2 425180 12702.92 12355.72 0.5954893
2017 84.4 0.2 435140 12701.35 12320.43 0.6009385
2018 84.7 0.2 433250 12627.93 12229.11 0.5936559

Respuesta 2.1:

El conjunto de datos “sa_gdp_elec” presenta 16 observaciones.

Respuesta 2.2: Los tipos de variables presentados son los siguientes:

year: Numerica/Discreta

access_elec: Numerica/Continua

startup: Numerica/Continua

co2: Numerica/Discreta

gdp: Numerica/Continua

gni: Numerica/Continua

co2_kg_ppp: Numerica/Continua

#Respuesta 2.3: Recomiendo aplicar una cantidad de 7 barras para tener una vision mas clara para hallar valores concretos.
hist(co2, breaks=7)

observación: si sólo hay 16 observaciones, usar una cantidad de 7 barras es un valor muy alto, que termina generando poca agrupación, es decir, que tenemos mucho detalle y no una representación con menos detalle de los datos

Respuesta 2.4:

Hacer un histograma para la variable “year” no es recomendable puesto que es una variable temporal. El histograma nos muestra frecuencias de años, pero en los datos datos solo aparecen una sola vez.

observación: no es apropiado que digas “variable temporal” sino de tiempo. Adicionalmente, no es recomendable es porque solo hay una observación por año y no se tendría mucho que agrupar

Respuesta 2.6:

Las variables independientes son: Gni y Co2_Kg_ppp

Respuesta 2.7:

Tipo de estudio: Observacional longitudinal

3) Representación de Histogramas

A continuación se presenta el código para generar 3 histogramas donde cada uno corresponde a un conjunto de datos distinto. Para cada uno, realizar las siguientes actividades:

  1. Generar el gráfico y determinar si variando la cantidad de barras mediante el parametro breaks puede mejorar la representación.

  2. Identificar el conjunto de dato usado para generar los gráficos así como el nombre de la variable, especificando lo que representa y el tipo de variable.

  3. Identifique la modalidad de la distribución, así como posibles sesgos, según su selección de cantidad de barras

hist(london_murders$age, breaks = 15)

respuestas:

usar tantas celdas de código y texto como sea necesario

#Con una cantidad de barras de 8 nos permite observar un sesgo a la derecha y la distribucion es unimodal.
hist(london_murders$age, breaks = 8)

hist(LAhomes$sqft)

respuestas:

usar tantas celdas de código y texto como sea necesario

observación: incompleta esta sección

hist(life_exp$income)

respuestas:

usar tantas celdas de código y texto como sea necesario

hist(lego_population$pieces, breaks=5)

respuestas:

usar tantas celdas de código y texto como sea necesario

#En este histograma al agregarle mas barras nos nos da una mejor vision de el, aunque con 20 barras podemos ya decir que esta sesgada a la derecha y unimodal.
hist(lego_population$pieces, breaks=20)

observación: eso mismo lo podías apreciar con menos barras. Igualmente, no entraste en análisis sobre lo que quiere decir el comportamiento de la variable representada