## ----eval = FALSE-------------------------------------------------------------
# activeDHS = unlist(regionDB$regionGRL[which(regionDB$regionAnno$collection == "sheffield_dnase")])
# activeDHS = disjoin(activeDHS)
# activeDHS

## ----eval = FALSE-------------------------------------------------------------
# restrictedUniverse = unlist(userSets)

## ----eval = FALSE-------------------------------------------------------------
# locResults = runLOLA(userSets, activeDHS, regionDB, cores=1)
# locResultsRestricted = runLOLA(userSets, restrictedUniverse, regionDB, cores=1)

## ----eval = FALSE-------------------------------------------------------------
# locResults[userSet==2,][order(maxRnk, decreasing=FALSE),][1:10,]
# locResultsRestricted[userSet==2,][order(maxRnk, decreasing=FALSE),][1:10,]
# 
# locResults[userSet==3,][order(maxRnk, decreasing=FALSE),][1:10,]
# locResultsRestricted[userSet==3,][order(maxRnk, decreasing=FALSE),][1:10,]