COVID-19 Cases

The following is the distribution by department of the number of COVID-19 cases as of September 24, 2020.

## install.packages("leaflet")
library(leaflet)

df <- data.frame(Department= c("Bogotá D.C.", "Antioquia", "Atlántico", "Valle del Cauca", "Cundinamarca", "Santander", "Bolívar", "Córdoba", "Cesar", "Nariño", "Norte de Santander", "Meta", "Magdalena", "Sucre", "Tolima", "Risaralda", "Huila", "Cauca", "Caquetá", "La Guajira", "Boyacá", "Caldas", "Chocó", "Putumayo", "Quindío", "Amazonas", "Casanare", "Arauca", "San Andrés", "Vaupés", "Guainía", "Guaviare", "Vichada"),
                          Confirmed = c(257679, 107331, 66862, 58588, 31575, 28751, 28304, 23340, 17878, 17742, 14932, 14850, 14814, 13785, 11422, 10239, 10082, 8652, 7843, 7609, 6410, 5117, 3967, 3679, 2941, 2739, 1972, 1519, 1055, 731, 720, 677, 463),
                          deaths = c(6568, 2266, 3037, 2154, 929, 1209, 761, 1505, 526, 663, 856, 381, 806, 570, 293, 220, 310, 238, 282, 286, 129, 114, 155, 164, 87, 116, 41, 39, 10, 6, 10, 11, 4),
                 lat= c( 4.60971, 6.25184,  10.96854, 3.43722, 4.57937, 7.12539, 10.39972, 8.74798, 10.46314, 1.21361, 7.89391,  4.142, 11.24079, 9.30472, 4.43889, 4.81333, 2.9273, 2.43823, 1.61438, 11.54444, 5.53528, 5.06889, 5.69188, 1.15284, 4.53389, -4.21528, 5.33775, 7.08471, 12.58317, 1.25509, 3.867, 2.57286, 6.18493),
                          lng= c( -74.08175, -75.56359, -74.78132, -76.5225, -74.21682, -73.1198, -75.51444, -75.88143, -73.25322, -77.28111, -72.50782, -73.62664, -74.19904, -75.39778, -75.23222, -75.69611, -75.28189, -76.61316, -75.60623, -72.90722, -73.36778, -75.51738, -76.65835, -76.65208, -75.68111, -69.94056, -72.39586, -70.75908, -81.70636, -70.235, -67.917, -72.64591, -67.4894))

df %>% 
    leaflet () %>%
    addTiles () %>%
    addCircles(weight = 2, radius = sqrt(df$Confirmed)*300,  popup = c(df$Department, df$Confirmed))
## Assuming "lng" and "lat" are longitude and latitude, respectively

COVID-19 Deaths

Distribution of deaths from COVID-19 in Colombia by September 24, 2020.

df %>% 
    leaflet () %>%
    addTiles () %>%
    addCircles(weight = 2, color= "red", radius = sqrt(df$deaths)*400,  popup = c(df$Department, df$deaths))
## Assuming "lng" and "lat" are longitude and latitude, respectively