Plot the amount of sequencing data from SRA

plot_sra_stat.R

library(ggplot2)
library(tidyverse)
library(readr)
sra_stat <- read_csv("https://www.ncbi.nlm.nih.gov/Traces/sra/sra_stat.cgi")
sra_stat$date <- as.Date(sra_stat$date, format = "%m/%d/%Y")

ggplot(sra_stat, aes(x = date, group = 1)) +
  geom_area(aes(y = bases / 10^15), fill = "orange", alpha = 0.8) + 
  geom_area(aes(y = bytes / 10^15), fill = "blue", alpha = 0.8) + 
  labs(title = "Amount of Sequencing Data in SRA", 
       x = "Year", y = "Petabases (orange) or petabytes (blue)") +
  scale_x_date(limits = as.Date(c("2013-01-01", Sys.Date()))) +
  theme_bw(base_size = 12)