Package: LOLA 1.43.0

Nathan Sheffield

LOLA: Locus overlap analysis for enrichment of genomic ranges

Provides functions for testing overlap of sets of genomic regions with public and custom region set (genomic ranges) databases. This makes it possible to do automated enrichment analysis for genomic region sets, thus facilitating interpretation of functional genomics and epigenomics data.

Authors:Nathan Sheffield <http://www.databio.org> [aut, cre], Christoph Bock [ctb]

LOLA_1.43.0.tar.gz
LOLA_1.43.0.zip(r-4.7)LOLA_1.43.0.zip(r-4.6)LOLA_1.43.0.zip(r-4.5)
LOLA_1.43.0.tgz(r-4.6-any)LOLA_1.43.0.tgz(r-4.5-any)
LOLA_1.43.0.tar.gz(r-4.7-any)LOLA_1.43.0.tar.gz(r-4.6-any)
LOLA_1.43.0.tgz(r-4.6-emscripten)
manual.pdf |manual.html✨
DESCRIPTION |NEWS
card.svg |card.png
LOLA/json (API)

# Install 'LOLA' in R:
install.packages('LOLA', repos = c('https://bioc.r-universe.dev', 'https://cloud.r-project.org'))

Bug tracker:https://github.com/nsheff/lola/issues

Datasets:
  • userSets - An example set of regions, sampled from the example database.
  • userUniverse - A reduced GRanges object from the example regionDB database

On BioConductor:LOLA-1.43.0(bioc 3.24)LOLA-1.42.0(bioc 3.23)

genesetenrichmentgeneregulationgenomeannotationsystemsbiologyfunctionalgenomicschipseqmethylseqsequencing

9.73 score 80 stars 249 scripts 56 mentions 20 exports 18 dependencies

Last updated from:a12914e6d5. Checks:1 NOTE, 9 OK. Indexed: yes.

TargetResultTimeFilesSyslog
bioc-checksNOTE172
linux-devel-x86_64OK239
source / vignettesOK189
linux-release-x86_64OK235
macos-release-arm64OK117
macos-oldrel-arm64OK123
windows-develOK148
windows-releaseOK145
windows-oldrelOK125
wasm-releaseOK148

Exports:buildRestrictedUniversecheckUniverseAppropriatenessextractEnrichmentOverlapsgetRegionFilegetRegionSetlistRegionSetsloadRegionDBmergeRegionDBsplotTopLOLAEnrichmentsreadBedreadCollectionreadCollectionAnnotationreadCollectionFilesreadRegionGRLreadRegionSetAnnotationredefineUserSetsrunLOLAsetSharedDataDirsplitFileIntoCollectionwriteCombinedEnrichment

Dependencies:BiocGenericsclidata.tablegenericsGenomicRangesglueIRangeslifecyclemagrittrplyrRcppreshape2rlangS4VectorsSeqinfostringistringrvctrs

Getting Started with LOLA
LOLA bioconductor package | Preparing analysis | Run the analysis | Exploring LOLA Results | Extracting certain region sets from a database
Rendered fromgettingStarted.Rmdusingknitr::rmarkdown

Last update: 2019-06-14
Started: 2015-04-09

Choosing a LOLA Universe
What's the universe?
Rendered fromchoosingUniverse.Rmdusingknitr::rmarkdown

Last update: 2018-02-06
Started: 2015-04-20

Using LOLA Core
Loading full-scale example data | Running the enrichment | Exploring results | Mystery set identity
Rendered fromusingLOLACore.Rmdusingknitr::rmarkdown

Last update: 2018-02-06
Started: 2015-04-16

Readme and manuals

Help Manual

Help pageTopics
If you want to test for differential enrichment within your usersets, you can restrict the universe to only regions that are covered in at least one of your sets. This function helps you build just such a restricted universebuildRestrictedUniverse
Check universe appropriatenesscheckUniverseAppropriateness
cleanws takes multi-line, code formatted strings and just formats them as simple stringscleanws
Just a reverser. Reverses the order of arguments and passes them untouched to countOverlapsAny - so you can use it with lapply.countOverlapsAnyRev
Given a single row from an enrichment table calculation, finds the set of overlaps between the user set and the test set. You can then use these, for example, to get sequences for those regions.extractEnrichmentOverlaps
Grab the filename for a a single region set from a database specified by filename.getRegionFile
Grab a single region set from a database, specified by filename.getRegionSet
Function to run lapply or mclapply, depending on the option set in getOption("mc.cores"), which can be set with setLapplyAlias().lapplyAlias
Lists the region sets for given collection(s) in a region database on disk.listRegionSets
converts a list of GRanges into a GRangesList; strips all metadata.listToGRangesList
Helper function to annotate and load a regionDB, a folder with subfolder collections of regions.loadRegionDB
Genome locus overlap analysis.LOLA
Given two regionDBs, (lists returned from loadRegionDB()), This function will combine them into a single regionDB. This will enable you to combine, for example, LOLA Core databases with custom databases into a single analysis.mergeRegionDBs
Named list function.nlist
Given some results (you grab the top ones on your own), this plots a barplot visualizing their odds ratios.plotTopLOLAEnrichments
Imports bed files and creates GRanges objects, using the fread() function from data.table.readBed
Given a bunch of region set files, read in all those flat (bed) files and create a GRangesList object holding all the region sets. This function is used by readRegionGRL to process annotation objects.readCollection
Read collection annotationreadCollectionAnnotation
Given a database and a collection, this will create the region annotation data.table; either giving a generic table based on file names, or by reading in the annotation data.readCollectionFiles
This function takes a region annotation object and reads in the regions, returning a GRangesList object of the regions.readRegionGRL
Given a folder containing region collections in subfolders, this function will either read the annotation file if one exists, or create a generic annotation file.readRegionSetAnnotation
This function will take the user sets, overlap with the universe, and redefine the user sets as the set of regions in the user universe that overlap at least one region in user sets. this makes for a more appropriate statistical enrichment comparison, as the user sets are actually exactly the same regions found in the universe otherwise, you can get some weird artifacts from the many-to-many relationship between user set regions and universe regions.redefineUserSets
This will change the string in filename to have a new extensionreplaceFileExtension
Enrichment CalculationrunLOLA
Function to sample regions from a GRangesList object, in specified proportionsampleGRL
To make parallel processing a possibility but not required, I use an lapply alias which can point at either the base lapply (for no multicore), or it can point to mclapply, and set the options for the number of cores (what mclapply uses). With no argument given, returns intead the number of cpus currently selected.setLapplyAlias
setSharedDataDir Sets global variable specifying the default data directory.setSharedDataDir
Efficiently split a data.table by a column in the tablesplitDataTable
This function will take a single large bed file that is annotated with a column grouping different sets of similar regions, and split it into separate files for use with the LOLA collection format.splitFileIntoCollection
An example set of regions, sampled from the example database.userSets
A reduced GRanges object from the example regionDB databaseuserUniverse
Wrapper of write.table that provides defaults to write a simple .tsv file. Passes additional arguments to write.tablewrite.tsv
Function for writing output all at once: combinedResults is an table generated by "locationEnrichment()" or by rbinding category/location results. Writes all enrichments to a single file, and also spits out the same data divided into groups based on userSets, and Databases, just for convenience. disable this with an option.writeCombinedEnrichment
Given a data table and a factor variable to split on, efficiently divides the table and then writes the different splits to separate files, named with filePrepend and numbered according to split.writeDataTableSplitByColumn