I am conducting a methylation analysis and attempting to identify Differentially Methylated Regions (DMRs). To achieve this, I've used the following code:
grs <- makeGenomicRatioSetFromMatrix(mat2, rownames = NULL, pData = dades.df,
array = "IlluminaHumanMethylationEPIC",
annotation="ilm10b2.hg19",
mergeManifest = FALSE, what = c("M"))
phenotype <- pData(grs)$SMOKER
designMatrix <- model.matrix(~ phenotype)
dmrs <- bumphunter(grs, design = designMatrix,
cutoff = 0.2, B = 0, type = "M")
My issue lies in the output obtained, which consists of variables such as chr, start, end, value, area, cluster, indexStart, indexEnd, L, and clusterL. However, I do not receive a column containing p-values. This is crucial for my analysis, and I need to ascertain the p-values for each DMR.
Does anyone know how I can obtain this p-value? Is it possible to obtain it using the bumphunter function? Alternatively, is there another method for obtaining DMRs that provides p-values in the output?
Any help is appreciated
Irene