Entering edit mode
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?
Thanks!
Thank you. That's great to hear and I look forward to it.
Do you know of any other package that allows to do that?