AEUCI_9

---

library(pacman)
p_load("fdth","base64enc", "htmltools", "mime", "xfun", "prettydoc","readr", "knitr","DT","dplyr", "ggplot2","modeest")

library(normtest) 
library(nortest) 
library(moments)

## 
## Attaching package: 'moments'

## The following object is masked from 'package:modeest':
## 
##     skewness

Introducción

En esta investigacion se plantea buscar cual pais entre México y Rusia emite más CO2 a la atmosfera anualmente en un rango de años de 1891 a 2019.

Analisis de la contaminacion por CO2 en mexico con respecto al pais ruso en diferentes años. Para este analisis se tomaron los datos del siguiente link:

https://drive.google.com/file/d/1u_Ad4PgLrM3LZ_npBp34nJFw8igwwVpr/view?usp=sharing. Para los cuales se tomaron los datos de Mexico y rusia (año y percapita de la contaminacion con respecto a este)

Preguntas y respuestas ¿QUÉ PAÍS EMITE MÁS? ¿POR QUÉ? ¿QUÉ IMPLICA ESTO?

Ademas para ello se analizaran pruebas para compararlos entre si utilizando la analogia visual y textual por medio del programa R estudio.

PRUEBA DE HIPOTESIS ESTADISTICA DESCRIPTIVA REGRESION LINEAL CONTRASTE CON GRAFICO DE CAJA Y BIGOTE

#Importar datos Año y Co2 (percapita) Para comenzar este estudio lo primero que se hara, es importar los datos que se analizaran.

library(readr)
mexico <- read_csv("mexico.csv")

## Rows: 129 Columns: 4

## -- Column specification --------------------------------------------------------
## Delimiter: ","
## chr (2): Entity, Code
## dbl (2): Year, Annual CO2 emissions (per capita)

## 
## i Use `spec()` to retrieve the full column specification for this data.
## i Specify the column types or set `show_col_types = FALSE` to quiet this message.

View(mexico)

library(readr)
rusia <- read_csv("rusia.csv")

## Rows: 129 Columns: 4

## -- Column specification --------------------------------------------------------
## Delimiter: ","
## chr (2): Entity, Code
## dbl (2): Year, Annual CO2 emissions (per capita)

## 
## i Use `spec()` to retrieve the full column specification for this data.
## i Specify the column types or set `show_col_types = FALSE` to quiet this message.

View(rusia)

Emisión de CO2 percapita entre México y Rusia

Prueba de hipotesis

¿Cuál pais emite mayor concentración de CO2?

H0: Rusia emite mayor concentración de CO2 que México. H1: México emite mayor concentración de CO2 que Rusia.

Tabla de Mexico

datatable(mexico)

Tabla de Rusia

datatable(rusia)

Tablas de distribucion de frecuencia

Distribucion de frecuenacia Mexico

tablaMEX = fdt(mexico$"Annual CO2 emissions (per capita)" )
tablaMEX

##     Class limits  f   rf rf(%)  cf  cf(%)
##  [0.0392,0.7156) 23 0.18 17.83  23  17.83
##   [0.7156,1.392) 31 0.24 24.03  54  41.86
##    [1.392,2.068) 16 0.12 12.40  70  54.26
##    [2.068,2.745)  6 0.05  4.65  76  58.91
##    [2.745,3.421)  5 0.04  3.88  81  62.79
##    [3.421,4.097) 31 0.24 24.03 112  86.82
##    [4.097,4.774) 14 0.11 10.85 126  97.67
##     [4.774,5.45)  1 0.01  0.78 127  98.45
##     [5.45,6.126)  2 0.02  1.55 129 100.00

Distribucion de frecuenacia Rusia

tablaRUS = fdt(rusia$"Annual CO2 emissions (per capita)" )
tablaRUS

##      Class limits  f   rf rf(%)  cf  cf(%)
##  [0.24176,2.1358) 46 0.36 35.66  46  35.66
##   [2.1358,4.0298) 13 0.10 10.08  59  45.74
##   [4.0298,5.9239)  6 0.05  4.65  65  50.39
##   [5.9239,7.8179)  7 0.05  5.43  72  55.81
##    [7.8179,9.712)  4 0.03  3.10  76  58.91
##    [9.712,11.606) 29 0.22 22.48 105  81.40
##     [11.606,13.5)  7 0.05  5.43 112  86.82
##     [13.5,15.394)  9 0.07  6.98 121  93.80
##   [15.394,17.288)  8 0.06  6.20 129 100.00

Histograma

hist(mexico$"Annual CO2 emissions (per capita)", main = "Mexico", xlab = "CO2")

hist(rusia$"Annual CO2 emissions (per capita)", main = "Rusia", xlab = "CO2")

plot(tablaMEX,type="fh", main = "Histograma de frecuencia absoluta Mexico")   

plot(tablaRUS,type="fh", main = "Histograma de frecuencia absoluta Rusia")   

plot(tablaMEX,type="rfh", main = "Histograma de frecuencia relativa Mexico")

plot(tablaRUS,type="rfh", main = "Histograma de frecuencia relativa Rusia")

Medidas de tendencia

Media

• Mexico

mean(mexico$"Annual CO2 emissions (per capita)")

## [1] 2.252332

*Rusia

mean(rusia$"Annual CO2 emissions (per capita)")

## [1] 6.576326

Mediana

• Mexico

median(mexico$"Annual CO2 emissions (per capita)")

## [1] 1.776

• Rusia

mean(rusia$"Annual CO2 emissions (per capita)")

## [1] 6.576326

Moda

• Mexico

mlv(mexico$"Annual CO2 emissions (per capita)")

## Warning: argument 'method' is missing. Data are supposed to be continuous. 
##             Default method 'shorth' is used

## [1] 0.8591754

• Rusia

mlv(rusia$"Annual CO2 emissions (per capita)")

## Warning: argument 'method' is missing. Data are supposed to be continuous. 
##             Default method 'shorth' is used

## [1] 1.517985

Valor minimo

• Mexico

min(mexico$"Annual CO2 emissions (per capita)")

## [1] 0.0396

• Rusia

min(rusia$"Annual CO2 emissions (per capita)")

## [1] 0.2442

Valor Maximo

• Mexico

max(mexico$"Annual CO2 emissions (per capita)")

## [1] 6.0657

• Rusia

max(rusia$"Annual CO2 emissions (per capita)")

## [1] 17.117

Cuartiles

• Mexico

summary(mexico$"Annual CO2 emissions (per capita)")

##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##  0.0396  1.0456  1.7760  2.2523  3.8218  6.0657

• Rusia

summary(rusia$"Annual CO2 emissions (per capita)")

##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##  0.2442  0.8390  5.6756  6.5763 11.2014 17.1170

Caja y bigote

• Mexico

boxplot(mexico$"Annual CO2 emissions (per capita)",main="Mexico co2 per Capita")

* Rusia

boxplot(rusia$"Annual CO2 emissions (per capita)",main="Mexico co2 per Capita")

Varianza

• Mexico

var(mexico$"Annual CO2 emissions (per capita)")

## [1] 2.52578

• Rusia

var(rusia$"Annual CO2 emissions (per capita)")

## [1] 30.86872

Desviacion extandar

• Mexico

sd(mexico$"Annual CO2 emissions (per capita)")

## [1] 1.58927

• Rusia

sd(rusia$"Annual CO2 emissions (per capita)")

## [1] 5.555962

Diagrama de dispersion

dd= read_csv("MXvsRS.csv")

## Rows: 129 Columns: 2

## -- Column specification --------------------------------------------------------
## Delimiter: ","
## dbl (2): mexico, rusia

## 
## i Use `spec()` to retrieve the full column specification for this data.
## i Specify the column types or set `show_col_types = FALSE` to quiet this message.

pairs(mexico$Year~mexico$"Annual CO2 emissions (per capita)")

pairs(rusia$Year~rusia$"Annual CO2 emissions (per capita)")

pairs(dd)

Pruebas de normalidad.

Kolmogorov-Smirnov

• Mexico

ks.test(dd$mexico, "pnorm", mean = mean(dd$mexico), sd = sd((dd$mexico)))

## Warning in ks.test(dd$mexico, "pnorm", mean = mean(dd$mexico), sd =
## sd((dd$mexico))): ties should not be present for the Kolmogorov-Smirnov test

## 
##  One-sample Kolmogorov-Smirnov test
## 
## data:  dd$mexico
## D = 0.15388, p-value = 0.004443
## alternative hypothesis: two-sided

• Rusia

ks.test(dd$rusia, "pnorm", mean = mean(dd$rusia), sd = sd((dd$rusia)))

## 
##  One-sample Kolmogorov-Smirnov test
## 
## data:  dd$rusia
## D = 0.17978, p-value = 0.0004781
## alternative hypothesis: two-sided

Prueba de Lilliefors

• Mexico

lillie.test(dd$mexico)

## 
##  Lilliefors (Kolmogorov-Smirnov) normality test
## 
## data:  dd$mexico
## D = 0.15388, p-value = 6.064e-08

• Rusia

lillie.test(dd$rusia)

## 
##  Lilliefors (Kolmogorov-Smirnov) normality test
## 
## data:  dd$rusia
## D = 0.17978, p-value = 4.993e-11

Analisis de correlacion

regMEX = lm(mexico$"Annual CO2 emissions (per capita)" ~ mexico$Year,  data =  mexico )
regRUS = lm(rusia$"Annual CO2 emissions (per capita)" ~ rusia$Year,  data =  rusia )

• Mexico

# Mexico
summary(regMEX)

## 
## Call:
## lm(formula = mexico$"Annual CO2 emissions (per capita)" ~ mexico$Year, 
##     data = mexico)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -1.0645 -0.5990 -0.2598  0.3741  4.9118 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -60.906478   4.794832  -12.70   <2e-16 ***
## mexico$Year   0.032306   0.002452   13.18   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.037 on 127 degrees of freedom
## Multiple R-squared:  0.5775, Adjusted R-squared:  0.5741 
## F-statistic: 173.6 on 1 and 127 DF,  p-value: < 2.2e-16

• Rusia

summary(regRUS)

## 
## Call:
## lm(formula = rusia$"Annual CO2 emissions (per capita)" ~ rusia$Year, 
##     data = rusia)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.4721 -2.1445 -0.6154  1.5325  5.9864 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -2.478e+02  1.246e+01  -19.89   <2e-16 ***
## rusia$Year   1.301e-01  6.372e-03   20.42   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 2.695 on 127 degrees of freedom
## Multiple R-squared:  0.7666, Adjusted R-squared:  0.7647 
## F-statistic:   417 on 1 and 127 DF,  p-value: < 2.2e-16

Recta de minimos cuadrados

plot(mexico$Year, mexico$`Annual CO2 emissions (per capita)`, xlab="Año", ylab = "CO2 per capita",main="Mexico")
abline(regMEX)

plot(rusia$Year, rusia$`Annual CO2 emissions (per capita)`, xlab="Año", ylab = "CO2 per capita",main="Rusia")
abline(regRUS)

Distribucion normal

• Mexico

pnorm(dd$mexico,mean = mean(dd$mexico),sd =sd(dd$mexico))

##   [1] 0.08235559 0.08276751 0.08231735 0.08231735 0.08191665 0.08249911
##   [7] 0.08325807 0.08353798 0.08359599 0.08531088 0.08751141 0.08824418
##  [13] 0.08911330 0.09018306 0.08997856 0.09139661 0.09377566 0.10503689
##  [19] 0.10421523 0.10768136 0.14505871 0.14722124 0.18008762 0.16277597
##  [25] 0.21072418 0.24816696 0.38378136 0.45974638 0.63435426 0.96655359
##  [31] 0.99179005 0.98455621 0.93958008 0.90990251 0.80762676 0.63582050
##  [37] 0.42883239 0.31104595 0.28107293 0.25516961 0.21088776 0.20228581
##  [43] 0.20493882 0.22383675 0.24707282 0.24880474 0.27642330 0.26492569
##  [49] 0.24343124 0.24339182 0.23122374 0.19488421 0.19049259 0.19363803
##  [55] 0.20709044 0.22391202 0.24572363 0.24693375 0.22847349 0.23260576
##  [61] 0.25668785 0.26568870 0.25488672 0.24505053 0.25920829 0.26684564
##  [67] 0.30295612 0.32290095 0.33889025 0.35686357 0.35803768 0.33746378
##  [73] 0.34048050 0.37234793 0.36494464 0.38219639 0.41639557 0.42173022
##  [79] 0.44830147 0.49010276 0.53138756 0.54168599 0.57623446 0.60355257
##  [85] 0.62389340 0.68116210 0.69454223 0.77345139 0.80774717 0.85762649
##  [91] 0.87606570 0.89880626 0.83786363 0.82249959 0.83240002 0.83225798
##  [97] 0.84504590 0.83269963 0.90940039 0.83129323 0.84421997 0.83831107
## [103] 0.83719871 0.84994066 0.80484452 0.82080884 0.84954409 0.87007032
## [109] 0.86506882 0.86490487 0.87683654 0.87118469 0.89524168 0.88969014
## [115] 0.90937985 0.91476140 0.91118417 0.91651385 0.89350660 0.87297481
## [121] 0.88799744 0.89355283 0.88062778 0.86393204 0.85802275 0.85434319
## [127] 0.81824549 0.79729472 0.77202804

• Rusia

pnorm(dd$rusia,mean = mean(dd$rusia),sd =sd(dd$rusia))

##   [1] 0.1295232 0.1300899 0.1320735 0.1322083 0.1342257 0.1343698 0.1368353
##   [8] 0.1390203 0.1423280 0.1454870 0.1459933 0.1447236 0.1438568 0.1466785
##  [15] 0.1431805 0.1466702 0.1500528 0.1499521 0.1515352 0.1505192 0.1514930
##  [22] 0.1559544 0.1591320 0.1544264 0.1508855 0.1528785 0.1495874 0.1300442
##  [29] 0.1281873 0.1272058 0.1278897 0.1299604 0.1311059 0.1335186 0.1331699
##  [36] 0.1400436 0.1454582 0.1484051 0.1505907 0.1630734 0.1727137 0.1799531
##  [43] 0.1876079 0.2038919 0.2163589 0.2326305 0.2344660 0.2409726 0.2414215
##  [50] 0.2705442 0.2465759 0.1800570 0.1824715 0.2163273 0.2220921 0.2364871
##  [57] 0.2651151 0.2751025 0.3003457 0.3242463 0.3422651 0.3578244 0.3703288
##  [64] 0.3944308 0.4356061 0.4675425 0.5016065 0.5260876 0.5442676 0.5582242
##  [71] 0.5677681 0.5874262 0.6203038 0.6500632 0.6796624 0.7077156 0.7306725
##  [78] 0.7440047 0.7650225 0.7913941 0.8218472 0.8462336 0.8645148 0.8798952
##  [85] 0.8998774 0.9131808 0.9238105 0.9337280 0.9343466 0.9447929 0.9363554
##  [92] 0.9402720 0.9425509 0.9431021 0.9622296 0.9610727 0.9581955 0.9631006
##  [99] 0.9558551 0.9710981 0.9582067 0.8835269 0.8581393 0.7903954 0.7807878
## [106] 0.7696702 0.7302622 0.7248643 0.7379215 0.7339552 0.7507734 0.7485469
## [113] 0.7626599 0.7667611 0.7748281 0.7974232 0.7972606 0.8086105 0.7691777
## [120] 0.7994381 0.8163901 0.8201954 0.7983106 0.7980877 0.7970116 0.7942706
## [127] 0.8029725 0.8173306 0.8125204

Preguntas y respuestas

¿Qué país emite más?

Rusia, confirmando la hipotersis nula.

¿Por qué?

Rusia aún a ser un país más grande que México pensamos que por tener una población algo cercana (128.9 millones en México y 144.1 millones en Rusia) y cantidad de habitantes por kilometro cuadrado, este emitiria menor concentración que mexico en cuanto a la concentración de poblacio, sin embargo, debido a que Rusia es un país grande y frío los servicos de electricidad, combustible y calefacción estan a mayor uso que México, sin mencionar que tambien es un país armamentistico y con servicios espaciales, esto dio como resultado una mayor cantidad de concentración de CO2 en comparacion con México que es un país mayormente en producción de materia prima y manufactura.

¿Qué implica esto?

México es un país más pequeño pero con la población más concentrada a comparación con Rusia, sin embargo, Rusia utiliza más recursos por condiciones de clima y publación que son 10 millones de habitantes mayor a México, contribuyendo más a la contaminación de CO2 en comparación a México.

Conclusión

En conclusión, el país que emite más contaminación atmosférica debido al gas CO2 es el Russo, los datos registrados que se analizaron son del año 1891 al 2019. Como se puede apreciar en los histogramas de frecuencias y en los modelos durante el estudio, se puede notar que los valores de contaminación en Mexico, el maximo fue en el año 1921 con un valor de 6.06 de percapita en la población mexicana( Se establece una Secretaría de Estado, que se denominará Secretaría de Educación Pública, se inaugura universiadades de bellas arte y se abre la posibilidad de dar pensiones en el extranjero), por otro lado el país Ruso en el año de 1990 registró 17.11 de percapita en su población(un año despues empieza a disminuir constantemente debido a la caida de la Union Sovietica), en cuanto a los registos más actuales en 2019 se registro que México tenia 3.4372 y Rusia 11.5057 teniando casi el triple en comparación con México, donde estos paises an aplicado leyes sobre la contaminación de CO2.

Entonces podemos decir que los valores de Mexico en la contaminación son mínimos pero frecuentes a comparación del Ruso que son menos frecuentes pero en mayor cantidad en el análisis, estos niveles altos de contaminación del CO2 probablemente se deben a que el país Russo al ser un país frio, utiliza una mayor cantidad de recursos energeticos que tienen a aparte de sus pruebas armamentisticas que un alto impacto en el ambiente.