get_japan_regional_cases: Japan regional cases, daily

In epiforecasts/NCoVUtils: Data Extraction Tools for the Covid-19 Outbreak

Description

Extracts regional case counts for Japan. Source.

Usage

get_japan_regional_cases()

Examples

get_japan_regional_cases()

## Not run: 

library(ggplot2)
library(dplyr)

 # Mapping: set up
 data <- get_japan_regional_cases()

 regions <- rnaturalearth::ne_states("Japan", returnclass = "sf")

 regions_with_data <- regions %>%
   dplyr::left_join(data, by = c("name" = "region")) %>%
   dplyr::mutate(cases =  ifelse(is.na(cases), 0, cases)

   # Map inset: Okinawa
   jp_okinawa <- dplyr::filter(regions_with_data, name == "Okinawa")

   jp_okinawa <- ggplot2::ggplot(jp_okinawa) +
   ggplot2::geom_sf(ggplot2::aes(fill = cases)) +  ggplot2::coord_sf(datum = NA) +
   ggplot2::xlab(jp_okinawa$region) +
   ggplot2::theme_bw() +
   ggplot2::theme(legend.position = "none")

   ## Map: mainland with insets
   jp_main <- dplyr::filter(regions_with_data, name != "Okinawa") %>%
   ggplot2::ggplot() +
   ggplot2::geom_sf(ggplot2::aes(fill = cases)) +
   ggplot2::coord_sf(crs = sf::st_crs(4326), xlim = c(127, 146), ylim = c(29, 46)) +
   ggplot2::theme_bw()

   jp_main +
   ggplot2::annotation_custom(
   grob = ggplot2::ggplotGrob(jp_okinawa), xmin = 140,  xmax = 146,  ymin = 24,  ymax = 37)


## End(Not run)

epiforecasts/NCoVUtils documentation built on June 26, 2020, 1:20 a.m.