U1A7.R

#Asignacion 7
setwd("~/Stat")

#cargar las bibliotecas con el paquete pacman
library(pacman)
library(gifski)
library(gganimate)

## Loading required package: ggplot2

p_load(readr, tidyverse, plotly)

#leer datos de la direccion general de epidemiologia (DGE)
datos <- read.csv("Casos_Diarios_Estado_Nacional_Confirmados_20200609.csv")

#dar formato a los diarios (absolutos) de Mexico
Mexico <- t(datos[datos$nombre=="MEXICO" ,])
Mexico <- as.vector(Mexico)
Mexico <- Mexico[4:153] #filtrar datos
Mexico <- as.numeric(Mexico) #para convertir a numerico
Mexico <- as.vector(Mexico)
cMexico <- cumsum(Mexico) #datos acumulados diarios

#Estadistica descriptiva
#figura1
plot(Mexico) #grafico de puntos de datos absoluto

summary(Mexico)

##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##     0.0     0.0    20.5   131.3   231.5   669.0

boxplot(Mexico)

hist(Mexico)

#dar formato a los 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 numerico
Jalisco <- as.vector(Jalisco)
cJalisco <- cumsum(Jalisco) #datos acumulados diarios

plot(Jalisco)

plot(cJalisco)

#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")
MexJal <- data.frame(Fecha, Mexico, Jalisco) #datos absolutos
cMexJal <- data.frame(Fecha, cMexico, cJalisco) #datos acumulados

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

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

#animacion grafica comparativa de frecuencias acumuladas
#por fechas para los estados de Mexico y Jalisco
#codificados por color 

ggplot(data=cMexJal)+
  geom_line(aes(Fecha, cJalisco, colour = 'Jalisco')) +
  geom_line(aes(Fecha, cMexico, colour = 'Mexico')) +
  xlab('Fecha') +
  ylab('Casos diarios acumulados') +
  labs(colour ='Estados') +
  transition_reveal(Fecha)

gMexJal <- ggplot(data=cMexJal)+
  ggtitle('COVID-19 en Mexico y Jalisco') +
  geom_line(aes(Fecha, cJalisco, colour = 'Jalisco')) +
  geom_line(aes(Fecha, cMexico, colour = 'Mexico')) +
  xlab('Fecha') +
  ylab('Casos diarios acumulados') +
  labs(colour ='Estados') 

gMexJal

ggplotly(gMexJal)