FigS2B_GlobalMethylation.R

library(ggplot2)

global <- read.csv("/Users/victoriasugrue/Desktop/globalmethaverage.csv") #contains file of average global methylation across all probes in each individual sample

ggplot(global, aes(x=Age, y=average, colour=Sex)) +
  geom_jitter(width=0.1, size=4, alpha=1/1.1) + 
  geom_smooth(method=lm, se=FALSE, size=1) + 
  theme_classic() + 
  ggtitle("Global methylation") + 
  ylab("Average methylation (%)") +
  xlab("Age (years)")