U1A7.R

setwd("~/PYE VERANO")

#cargar las bibliotecas con el paquete pacman
library(pacman)
p_load(readr, tidyverse, plotly, gganimate,gifski)

#Leer datos de la dirección general de epidemologia (DGE)
datos <- read.csv("Casos_Diarios_Estado_Nacional_Confirmados_20200609.csv")

#dar formato a los datos diarios de baja california
BC <- t(datos[datos$nombre=="BAJA CALIFORNIA" , ])
BC <- as.vector(BC)
BC <- BC[4:153]
BC <- as.numeric(BC)
BC<- as.vector(BC)
cBC <- cumsum(BC)

plot(BC)

summary(BC)

##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##    0.00    0.00    7.00   41.98   85.75  183.00

hist(BC)

#dar formato a los datos diarios de yucatan
Yucatan <- t(datos[datos$nombre=="YUCATAN" , ])
Yucatan <- as.vector(Yucatan)
Yucatan <- Yucatan[4:153]
Yucatan <- as.numeric(Yucatan)
Yucatan <- as.vector(Yucatan)
cYucatan <- cumsum(Yucatan)

Fecha <- seq(from= as.Date("2020-01-12"), to= as.Date("2020-06-09"), by="day")
  
  BCYu <- data.frame(Fecha, BC, Yucatan)
  cBCYu <- data.frame(Fecha, cBC, cYucatan)

ggplot(data=BCYu) +
  ggtitle("Casos diarios COVID-19 en Baja California (DGE)") +
  geom_line(mapping=aes(x=Fecha, y=BC))

ggplot(data=BCYu) +
  ggtitle("Casos diarios COVID-19 en Yucatan (DGE)") +
  geom_line(mapping=aes(x=Fecha, y=Yucatan))

ggplot(data = cBCYu)+
  geom_line(aes(Fecha, cYucatan, colour='Yucatan'))+
  geom_line(aes(Fecha, cBC, colour='Baja California'))+
  xlab('Fecha')+
  ylab('Casos diarios acumulados')+
  labs(colour="Estados")+
  transition_reveal(Fecha)

gBCYu <- ggplot(data = cBCYu)+
  ggtitle("COVID-19 en Baja California y Yucatan") +
  geom_line(aes(Fecha, cYucatan, colour='Yucatan'))+
  geom_line(aes(Fecha, cBC, colour='Baja California'))+
  xlab('Fecha')+
  ylab('Casos diarios acumulados')+
  labs(colour="Estados")

ggplotly(gBCYu)