#Ejercicio 5
setwd("~/Estadistica2020")

#cargar las bibliotecas con el paquete pacman
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
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 datos diarios (absolutos) de Campeche
Campeche <- t(datos[datos$nombre=="CAMPECHE" ,])
Campeche <- as.vector(Campeche)
Campeche <- Campeche[4:153] #filtrar datos
Campeche <- as.numeric(Campeche) #para convertir a numérico
Campeche <- as.vector(Campeche)
cCampeche <- cumsum(Campeche) #Datos acumulados diarios

#estadÃdstica descriptiva
#figura 1 
plot(Campeche)  #gráfico de puntos de datos absolutos

summary(Campeche)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   0.000   0.000   0.000   5.427   7.000  38.000
boxplot(Campeche)

hist(Campeche)

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

plot(Chiapas)

plot(cChiapas)

#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")
CamChi <- data.frame(Fecha, Campeche, Chiapas) #datos absolutos
cCamChi <- data.frame(Fecha, cCampeche, cChiapas) #datos acumulados

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

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

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

ggplot(data=cCamChi)+
  geom_line(aes(Fecha, cChiapas, colour = 'Chiapas')) +
  geom_line(aes(Fecha, cCampeche, colour = 'Campeche')) +
  xlab('Fecha') +
  ylab('Casos diarios acumulados') +
  labs(colour ='Estados') +
  transition_reveal(Fecha)

gCamChi <- ggplot(data=cCamChi)+
  ggtitle('COVID-19 en Campeche y Chiapas') +
  geom_line(aes(Fecha, cChiapas, colour = 'Chiapas')) +
  geom_line(aes(Fecha, cCampeche, colour = 'Campeche')) +
  xlab('Fecha') +
  ylab('Casos diarios acumulados') +
  labs(colour ='Estados') 

gCamChi

ggplotly(gCamChi)