dmr_result <- champ.DMR(beta = dmr_matrix, pheno = group, arraytype = "EPIC" )

<< Find DMR with Bumphunter Method >> 3 cores will be used to do parallel Bumphunter computing. According to your data set, champ.DMR() detected 13322 clusters contains MORE THAN 7 probes within300 maxGap. These clusters will be used to find DMR.

Error in [<-.data.frame(*tmp*, list, value = 0.001) : new columns would leave holes after existing columns