vignettes/Earthquake_package.R

In flyeye/earthquakes: Cleaning and visualizing data from NOAA Significant Earthquake Database

## ---- echo=TRUE, message=FALSE, warning=FALSE----------------------------
  require(readr)
  filename <- system.file("extdata", "earthquakes.tsv.gz", package = "earthquakes")
  raw_data <- readr::read_delim(filename, delim = "\t")

## ---- echo=TRUE, message=FALSE, warning=FALSE----------------------------
 require(earthquakes)
 head(earthquakes)
 str(earthquakes)

## ---- fig.height=5, fig.width=7------------------------------------------
 require(earthquakes)
 ggplot2::ggplot() +
   ggplot2::theme_void() +
   geom_timeline(data = earthquakes, ggplot2::aes(event_date = DATE,
                                   xmindate = 2000,
                                   xmaxdate = 2003,
                                   richter = EQ_PRIMARY,
                                   death = DEATHS,
                                   transparency = 0.3,
                                   colour = "red",
                                   scale = 0.5))+
   ggplot2::theme(legend.position="none")

## ---- fig.height=5, fig.width=7------------------------------------------
 require(earthquakes)
 ggplot2::ggplot() +
   ggplot2::theme_void() +
   geom_timeline(data = earthquakes, ggplot2::aes(event_date = DATE,
                                   xmindate = 2000,
                                   xmaxdate = 2003,
                                   richter = EQ_PRIMARY,
                                   death = DEATHS,
                                   transparency = 0.3,
                                   colour = "red",
                                   scale = 0.5, 
                                   layers = COUNTRY))+
   ggplot2::theme(legend.position="none")

## ---- fig.height=6, fig.width=8------------------------------------------
 require(earthquakes)
   ggplot2::ggplot() +
   ggplot2::theme_void() +
   geom_timeline(data = earthquakes, ggplot2::aes(event_date = DATE,
                                   xmindate = 2000,
                                   xmaxdate = 2005,
                                   richter = EQ_PRIMARY,
                                   death = DEATHS,
                                   transparency = 0.3,
                                   colour = "red",
                                   scale = 0.5,
                                   layers = COUNTRY))+
   ggplot2::theme(legend.position="none") +
   geom_timeline_label(data = earthquakes, ggplot2::aes(event_date = DATE,
                                         xmindate = 2000,
                                         xmaxdate = 2005,
                                         richter = EQ_PRIMARY,
                                         layers = COUNTRY,
                                         labels = LOCATION_NAME))


## ---- eval=FALSE, message=FALSE, warning=FALSE, include=FALSE------------
#   require(earthquakes)
#   filename <- system.file("extdata", "earthquakes.tsv.gz", package = "earthquakes")
#   raw_data <- read_delim(filename, delim = "\t")
#   data <- eq_clean_data(raw_data)
#   head(data)

## ---- message=FALSE, warning=FALSE---------------------------------------
 require(earthquakes)
 filename <- system.file("extdata", "earthquakes.tsv.gz", package = "earthquakes")
 raw_data <- read_delim(filename, delim = "\t")
 data <- eq_location_clean(raw_data)
 head(data$LOCATION_NAME)

## ---- message=FALSE, warning=FALSE---------------------------------------
require(earthquakes)
labels <- eq_create_label(earthquakes)
head(labels)

## ---- fig.height=5, fig.width=7, message=FALSE, warning=FALSE------------
 require(earthquakes)

 filename <- system.file("extdata", "earthquakes.tsv.gz", package = "earthquakes")
 readr::read_delim(filename, delim = "\t") %>%
   eq_clean_data() %>%
   dplyr::filter(COUNTRY == "MEXICO" & as.integer(YEAR) >= 2000) %>%
   dplyr::mutate(popup_text = eq_create_label(.)) %>%
   eq_map(annot_col = "popup_text")