if(require("VanillaICE")){
library(oligoClasses)
library(GenomicRanges)
library(Biobase)
library(foreach)
registerDoSEQ()
data(oligoSetExample, package="oligoClasses")
oligoSet <- oligoSet[chromosome(oligoSet) == 1, ]
grl <- hmm(oligoSet,TAUP=1e10)
class(grl)## GRangesList
gr <- grl[[1]]
(frange <- makeFeatureGRanges(oligoSet))
## which features overlap with the second range in sample NA06993
subsetByOverlaps(frange, gr[2,])
}
Run the code above in your browser using DataLab