## ---- 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,]