I want to detect differentially methylated regions (DMRs) on a whole-genome bisulfite sequencing (BS-seq) dataset.
However, instead of having two groups to compare, I'm working with a clinical parameter that is a continuous variable and would like to identify changes for unit of increase of my clinical parameter.
Software packages for RNA-seq analysis such as DESeq2 allow that by simply including the continuous variable in the model formula.
How can I perform a differential methylation analysys of BS-seq data using a continuous variable? What Bioconductor package(s) can I use?