require(org.Mm.eg.db)
require(TxDb.Mmusculus.UCSC.mm9.knownGene)
txdb <- TxDb.Mmusculus.UCSC.mm9.knownGene
# loading H3K4me1 peaks as a GRanges object
# built based on the BED file from the GEO GSM1234488 sample
# limited to chr19:3200000-4000000
H3K4me1GR <- system.file("extdata", "H3K4me1GR.Rda", package="compEpiTools")
load(H3K4me1GR)
enh <- enhancers(H3K4me1GR, txdb)
m.enh <- matchEnhancers(enhGR=enh, minD=2e4, maxD=2e5,
txdb=txdb, EG2GS=org.Mm.eg.db)
Run the code above in your browser using DataLab