Hello! This probably comes a lot later than you would like so. However, I would like to showcase my implementation. First of all, you will need extract the EPIC (or whichever array type you used) manifest for the mapping info.

manifest <- sesameDataGet("EPICv2.address")

From there you will create an annotation DF:

anno <- data.frame( CpG = manifest$ordering$Probe_ID, chr = manifest$hg38@seqnames@values, pos = manifest$hg38@ranges@start )

Inner join the two: anno_stat <- inner_join(dmp, anno, by ="CpG")

Since some of the variables are string type, we need to change them to int:

anno_stat$chr <- gsub("chr", "", anno_stat$chr)

anno_stat$chr <- recode(anno_stat$chr, X = "23", Y = "24", MT = "25")

anno_stat$chr <- as.numeric(anno_stat$chr)

Remove non-standard and NAs:

anno_stat <- anno_stat %>% filter(grepl("^chr[0-9XY]+$", paste0("chr", chr))) %>% filter(!is.na(chr), !is.na(pos), !is.na(P.Value))

And then we can finally plot:

library(qqman) manhattan(anno_stat, chr = "chr", bp = "pos", p = "P.Value", snp = "CpG")

Hope this helps, I tried sourcing the web, but I haven't stumbled on any CpG implementations yet.

Best, Dorde Racic

You can't really plot a manhattan plot with just 3 loci. They typically describe a distribution of values across an entire contig/chromosome/genome. Otherwise, a manhattan plot is just a scatter plot organised by genomic distance on the x and a given value on the y.

Here is a tutorial on how to do them in R with and without ggplot2. Halfway down the page are a few examples specifically of manhattan style plots.

EDIT: I realise there is the command head(sites), but it only printed 3 rows, when the default is 6 in R, hence my concern about the number of data points.