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.
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.