#U1A7
setwd("~/PyE")
#El estado de Jalisco es el estado con mayor tradición de México; muchos de sus elementos culturales 
#se han convertido en sinónimo de mexicanidad, como el mariachi, la charrería y el tequila. jalisco, es uno de los estados que mas se ha visto afectados 
#ante la epidemia del COVID-19, además de que ha tenido una administracion muy polémica durante la pandemia.
#En esta ocasión vamos a realizar el analisis y comparacion de Confirmados por día de COVID-19 ante el estado de Nayarit. 
#Otro estado muy reconocido en México por su diversidad cultural.

#cargar las bibliotecas 
library(pacman)
library(gifski)
library(gganimate)
## Loading required package: ggplot2
library(plotly)
## 
## Attaching package: 'plotly'
## The following object is masked from 'package:ggplot2':
## 
##     last_plot
## The following object is masked from 'package:stats':
## 
##     filter
## The following object is masked from 'package:graphics':
## 
##     layout
#Leer datos de la dirección general de epidemiología (DGE)

datos <- read.csv("Casos_Diarios_Estado_Nacional_Confirmados_20200609.csv")

#dar formato a los datos diarios (absolutos) de Jalisco
Jalisco <- t(datos[datos$nombre=="JALISCO" ,])
Jalisco <- as.vector(Jalisco)
Jalisco <- Jalisco[4:153] #filtrar datos
Jalisco <- as.numeric(Jalisco) #para convertir a numérico
Jalisco <- as.vector(Jalisco) 
cJalisco <- cumsum(Jalisco) #Datos acumulados diarios

#estadística descriptiva
#figura 1 
plot(Jalisco)  #gráfico de puntos de datos absolutos

summary(Jalisco)#Resumen de datos
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##    0.00    0.00    5.00   21.45   21.50  207.00
boxplot(Jalisco)#gráfico de Caja y Bigote

hist(Jalisco)#Histograma de los datos de Jalisco

#dar formato a los datos diarios (absolutos) de Nayarit
Nayarit <- t(datos[datos$nombre=="NAYARIT" ,])
Nayarit<- as.vector(Nayarit)
Nayarit <- Nayarit [4:153] #filtrar datos
Nayarit <- as.numeric(Nayarit) #para convertir a numérico
Nayarit<- as.vector(Nayarit)
cNayarit <- cumsum(Nayarit) #Datos acumulados diarios

plot(Nayarit)

plot(cNayarit)#En esta ocasión, podemos observar que es una gráfica que va creciendo de manera exponencial

#Construir un vector de Fechas
Fecha <- seq(from = as.Date("2020-01-12"), to = as.Date("2020-06-09"), by= "day")

#crear un arreglo matricial (marco de datos "data frame", "tibble")
NarJal <- data.frame(Fecha, Jalisco, Nayarit) #datos absolutos
cNarJal <- data.frame(Fecha, cJalisco, cNayarit) #datos acumulados

#grafica para Jalisco de datos absolutos con ggplot
ggplot(data = NarJal) +
  ggtitle("Casos diarios COVID-19 en Jalisco (DGE)") +
  geom_line (mapping = aes(x= Fecha, y= Jalisco))

#grafica para Nayarit de datos absolutos con ggplot
ggplot(data = NarJal) +
  ggtitle("Casos diarios COVID-19 en Nayarit (DGE)") +
  geom_line (mapping = aes(x= Fecha, y= Nayarit))

#animación gráfica comparativa de frecuencias acumuladas
#por fechas para los estados de Jalisco y Nayarit
#codificados por color 

ggplot(data=cNarJal)+
  geom_line(aes(Fecha, cNayarit, colour = 'Nayarit')) +
  geom_line(aes(Fecha, cJalisco, colour = 'Jalisco')) +
  xlab('Fecha') +
  ylab('Casos diarios acumulados') +
  labs(colour ='Estados') +
  transition_reveal(Fecha)

#Generamos la gráfica interactiva para mayor comodidad del usuario 
gNarJal <- ggplot(data=cNarJal)+
  ggtitle('COVID-19 en Jalisco y Nayarit') +
  geom_line(aes(Fecha, cNayarit, colour = 'Nayarit')) +
  geom_line(aes(Fecha, cJalisco, colour = 'Jalisco')) +
  xlab('Fecha') +
  ylab('Casos diarios acumulados') +
  labs(colour ='Estados') 

gNarJal

ggplotly(gNarJal)
#El día 3 de Marzo del presente año, el estado de Jalisco registró el primer caso de COVID-19,en cambio 7 días despues el estado de Nayarit
#registró el suyo, desde entonces, el estado de Jalisco ha tomado la delantera respecto a confirmados por día
#por otro lado, esto puede tener una respuesta ya que Jalisco cuenta con 8,256 millones de habitantes aproximadamente con una extensión territorial de 
#78.588 km² y Nayarit 1,181 millones de habitantes y una extensión territorial de 27.857 km²
#el contagio por día en el estado de Jalisco es muy alto
#en comparacion al estado de Nayarit, más de 2,000 confirmados de diferencia aproximadamente. 





#DATOS RECOPILADOS DE: https://coronavirus.gob.mx/datos/