## ----Load a regionDB-------------------------------------------------------
library("LOLA")
dbPath = system.file("extdata", "hg19", package="LOLA")
regionDB = loadRegionDB(dbPath)

## --------------------------------------------------------------------------
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")