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