edmr
edmr copied to clipboard
get.dmr.genes() is wrong
I believe the correct code should be: get.dmr.genes() <- function (myDMR, subject, id.type = "gene.symbol") { ind = findOverlaps(myDMR, subject) unique(values(subject)[unique(ind@to), id.type]) }
The code package offer is get.dmr.genes() <- function (myDMR, subject, id.type = "gene.symbol") { ind = findOverlaps(subject, myDMR) unique(values(subject)[unique(ind@to), id.type]) } And this would not return correct result.