I have identified differentially methylated positions from illumina EPIC methylation array data using a linear mixed model. Because I have related individuals in my sample, I needed to account for relatedness using a kinship matrix as a random effect. I would like to identify DMRs from the list of DMPs I have identified with my mixed model. But none of the approaches I have seen for identifying DMRs seem compatible with this.
Does anyone know of a function (preferably implemented in R) for identifying DMRs from a custom list of DMPs?