## -----------------------------------------------------------------------------------------------------------------------------#库(“ Lola”)###QueillregionDB(“ Lolacore/HG19”)###geirionseta = readbed(“ lola_vignette_data/seta_100.bed)”)#egionionSetb = readbed(“ lola_vignette_data/setb_100.bed”)#egiemionsEtc = readbed(“ lola_vignette_data/setc_100.bed”)#activedhs = readbed(lola_vignette_data/activeed_data/activeed_data/activeed_data/activeed_data/activeedhs_universe.universe.bed”-----------------------------------------------------------------------------------------------#userets = grangeslist(egionionseta,gearnionsetb,eginionsetc)#locresults = runlola(userest,activedhs,actigiveDhs,regiondb,cores = 1)## ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------#locresults [1:10,] ## ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- # locResults[order(meanRnk, decreasing=FALSE),][1:20,]#locresults [order(maxrnk,decteasing = false),] [1:20,]#locresults [collection ==“ sheffield_dnase”,],] [order(maxrnk,decrosing = false),]##########,] [订单(最大rnk,降低= false),],] [1:10,]#locresults [usereset == 3,] [order(maxrnk,decroding = false),] [1:10,]#####locresults [userset == 1,][collection ==“ sheffield_dnase”,] [1:15,]#locresults [userset == 2,] [collection ==“ sheffield_dnase”,] [1:15,]#locresults [usereset == 3,]==“ sheffield_dnase”,] [1:15,] ## ---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------#WriteCombinedEnrichment(Locresults,Outfolder =“ LolareSults”,includeplits = true)