setwd("~/ONLINEDS")
#Analisis de datos abiertos de COVI19 EN MEXICO 

library(pacman)
## Warning: package 'pacman' was built under R version 3.6.3
p_load(markdown, knitr, dplyr, tidyr,hashmap, tidyverse, lubridate,
       summarytools, ggpubr, kableExtra, reshape2,
       sf, tmap, readr, devtools, plotly, gganimate, gifski)
## Installing package into 'C:/Users/Daniel/Documents/R/win-library/3.6'
## (as 'lib' is unspecified)
## Warning: package 'hashmap' is not available (for R version 3.6.2)
## Warning: unable to access index for repository http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/3.6:
##   no fue posible abrir la URL 'http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/3.6/PACKAGES'
## Warning: 'BiocManager' not available.  Could not check Bioconductor.
## 
## Please use `install.packages('BiocManager')` and then retry.
## Warning in p_install(package, character.only = TRUE, ...):
## Warning in library(package, lib.loc = lib.loc, character.only = TRUE,
## logical.return = TRUE, : there is no package called 'hashmap'
## Installing package into 'C:/Users/Daniel/Documents/R/win-library/3.6'
## (as 'lib' is unspecified)
## also installing the dependencies 'sf', 'classInt'
## Warning: unable to access index for repository http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/3.6:
##   no fue posible abrir la URL 'http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/3.6/PACKAGES'
## Warning: package 'sf' is in use and will not be installed
## package 'classInt' successfully unpacked and MD5 sums checked
## Warning: cannot remove prior installation of package 'classInt'
## Warning in file.copy(savedcopy, lib, recursive = TRUE): problema al copiar C:
## \Users\Daniel\Documents\R\win-library\3.6\00LOCK\classInt\libs\x64\classInt.dll
## a C:\Users\Daniel\Documents\R\win-library\3.6\classInt\libs\x64\classInt.dll:
## Permission denied
## Warning: restored 'classInt'
## package 'tmap' successfully unpacked and MD5 sums checked
## 
## The downloaded binary packages are in
##  C:\Users\Daniel\AppData\Local\Temp\RtmpQngdC4\downloaded_packages
## 
## tmap installed
## Warning: package 'tmap' was built under R version 3.6.3
## Warning in p_load(markdown, knitr, dplyr, tidyr, hashmap, tidyverse, lubridate, : Failed to install/load:
## hashmap, tmap
##devtools::install_github("nathan-russell/hashmap")

library(ggplot2)
#Primera parte datos diarios y acumulados

#Leer datos 
datos<- read.csv("Casos_Diarios_Estado_Nacional_Confirmados.csv")
#Casos diarios Baja California
BC <- t(datos[datos$nombre=="BAJA CALIFORNIA" ,])
BC <- as.vector(BC)
BC <- BC[4:133]
BC <- as.numeric(BC)
BC <- as.vector(BC)
cBC <- cumsum(BC)

#Casos diarios Chiapas
CHIAPAS <- t(datos[datos$nombre=="CHIAPAS" ,])
CHIAPAS <- as.vector(CHIAPAS)
CHIAPAS <- CHIAPAS[4:133]
CHIAPAS <- as.numeric(CHIAPAS)
CHIAPAS <- as.vector(CHIAPAS)
cCHIAPAS <- cumsum(CHIAPAS)

#Fecha
Fecha = seq(from = as.Date("2020-01-05"), to = as.Date("2020-05-13"), by = 'day')

BCCH <-data.frame(Fecha, BC, CHIAPAS)
cBCCH <-data.frame(Fecha, cBC, cCHIAPAS)

#Graficos de casos diarios confirmados en Baja california
plot(BC)

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

# Grafica comparativa Baja california y Chiapas
ggplot(data=cBCCH) +
  geom_line(aes(Fecha, BC, colour = 'Baja California')) + 
  geom_line(aes(Fecha, CHIAPAS, colour ='Chiapas') )+ 
  xlab('Fecha')+
  ylab('Casos Diarios') +
  labs(colour = "Estados")+
  ggtitle("Comparativa Baja California y Chiapas Covid-19")+
  transition_reveal(Fecha)

# Grafica interactiva
p <- ggplot(data=cBCCH) +
  geom_line(aes(Fecha, cBC, colour = 'Baja California')) + 
  geom_line(aes(Fecha, cCHIAPAS, colour ='Chiapas') )+ 
  xlab('Fecha')+
  ylab('Casos Diarios') +
  labs(colour = "Estados")+
  ggtitle("Comparativa Baja California y Chiapas Covid-19")
ggplotly(p)
BC1 <-data.frame(Fecha, BC, cBC)

#Diarios y acumulados para Baja California
g2 <- ggplot(data=BC1)+
  geom_col(aes(Fecha, cBC))+
  xlab('Fecha')+
  ylab('Casos diarios acumulados')+
  ggtitle("Confirmados de COVID-19 en Baja California 13 de Mayo 2020")
 

g3 <- ggplot(data=BC1)+
  geom_line(aes(Fecha, cBC))+
  xlab('Fecha')+
  ylab('Casos diarios')+
  ggtitle("Confirmados de COVID-19 en Baja California 13 de Mayo 2020")

p_load(gridExtra)
grid.arrange(g2, g3)

CH1 <-data.frame(Fecha, CHIAPAS, cCHIAPAS)

#Diarios y acumulados para Chiapas
g4 <- ggplot(data=CH1)+
  geom_col(aes(Fecha, cCHIAPAS))+
  xlab('Fecha')+
  ylab('Casos diarios acumulados')+
  ggtitle("Confirmados de COVID-19 en Chiapas 13 de Mayo 2020")


g5 <- ggplot(data=CH1)+
  geom_line(aes(Fecha, cCHIAPAS))+
  xlab('Fecha')+
  ylab('Casos diarios')+
  ggtitle("Confirmados de COVID-19 en Chiapas 13 de Mayo 2020")

p_load(gridExtra)
grid.arrange(g4, g5)