## ---- echo=FALSE----------------------------------------------------------- # These settings make the vignette prettier knitr::opts_chunk$set(results="hold", message=FALSE) ## ----Load a regionDB------------------------------------------------------- library("LOLA") dbPath = system.file("extdata", "hg19", package="LOLA") regionDB = loadRegionDB(dbPath) ## ----Look at the elements of a regionDB------------------------------------ names(regionDB) ## ----Load sample user sets and universe------------------------------------ data("sample_input", package="LOLA") # load userSets data("sample_universe", package="LOLA") # load userUniverse ## ----Run the calculation--------------------------------------------------- locResults = runLOLA(userSets, userUniverse, regionDB, cores=1) ## -------------------------------------------------------------------------- colnames(locResults) head(locResults) ## -------------------------------------------------------------------------- locResults[order(support, decreasing=TRUE),] ## -------------------------------------------------------------------------- locResults[order(maxRnk, decreasing=TRUE),] ## ----Write results--------------------------------------------------------- writeCombinedEnrichment(locResults, outFolder= "lolaResults") ## ----Write split results--------------------------------------------------- writeCombinedEnrichment(locResults, outFolder= "lolaResults", includeSplits=TRUE) ## ----Extracting overlaps--------------------------------------------------- oneResult = locResults[2,] extractEnrichmentOverlaps(oneResult, userSets, regionDB) ## ----Grabbing individual region sets--------------------------------------- getRegionSet(regionDB, collections="ucsc_example", filenames="vistaEnhancers.bed") ## ----Grabbing individual region sets from disk----------------------------- getRegionSet(dbPath, collections="ucsc_example", filenames="vistaEnhancers.bed")