In Enlightengeo/GeoscienceBC_2019-008: Enlighten Geoscience BC 2019-008 Induced Seismicity Statistical Study Code

knitr::opts_chunk$set(echo = TRUE)
library(geosciencebc2019008)

HDBSCAN Plot

hdbscan = read_csv('../output/hdbscan_sensitivity_grid.csv')

seis_events = pivot_longer(hdbscan,-c(seismogenic_events,mean_seismogenic_mag))

ggplot(seis_events) +
  geom_point(aes(x=value, y=seismogenic_events, colour=mean_seismogenic_mag)) +
  scale_color_viridis_c(option='B', name='Mean Magnitude') +
  facet_wrap(. ~ name, scales = "free_x") +
  ylab('Seismogenic Events') +
  xlab('Value') +
  ggsave('../output/hdbscan_sensitivity_plot.jpeg', width=10, height=5, dpi=600)

Enlightengeo/GeoscienceBC_2019-008 documentation built on Feb. 4, 2021, 12:43 p.m.