a volcano plot example

[raynamharris]

A volcano plot example example

# set working directory to `training-rstudio-binder/GTEx/`
# install.packages("ggplot2")
library(ggplot2)

results <- read.table("./data/GTEx_Heart_20-29_vs_50-59.tsv")

ggplot(results, aes(x = logFC, y = -log10(adj.P.Val))) +
  geom_point(aes(color = ifelse( adj.P.Val < 0.05, "p < 0.05", "NS"))) +
  geom_hline(yintercept = -log10(0.05))  +
  theme(legend.position = "bottom") +
  labs(color = "20-29 vs 50-59 year olds", 
       subtitle = "Heart Tissue Gene Expression")