extendedMapSeqlevels {derfinder} | R Documentation |
If available, use the information from GenomeInfoDb for your species of
interest to map the sequence names from the style currently used to another
valid style. For example, for Homo sapiens map '2' (NCBI style) to 'chr2'
(UCSC style). If the information from GenomeInfoDb is not available, the
original sequence names will be returned. To disable this functionality
specify set chrsStyle
to NULL
.
extendedMapSeqlevels( seqnames, style = getOption("chrsStyle", "UCSC"), species = getOption("species", "homo_sapiens"), currentStyle = NULL, ... )
seqnames |
A character vector with the sequence names. |
style |
A single character vector specifying the naming style to use for renaming the sequence names. |
species |
A single character vector with the species of interest: it has
to match the valid species names supported in GenomeInfoDb. See
|
currentStyle |
A single character vector with the currently used
naming style. If |
... |
Arguments passed to other methods and/or advanced arguments. Advanced arguments:
|
This function is inspired from mapSeqlevels with the difference that it will return the original sequence names if the species, current naming style or target naming style are not supported in GenomeInfoDb.
If you want to disable this function, set chrsStyle
to NULL
.
From other functions in derfinder that pass the
...
argument to this function, use chrsStyle = NULL
.
This can be useful when working
with organisms that are absent from GenomeInfoDb as documented in
https://support.bioconductor.org/p/95521/.
A vector of sequence names using the specified naming style
.
L. Collado-Torres
## Without guessing any information extendedMapSeqlevels("2", "UCSC", "Homo sapiens", "NCBI") ## Guessing the current naming style extendedMapSeqlevels("2", "UCSC", "Homo sapiens") ## Guess species and current style extendedMapSeqlevels("2", "NCBI") ## Guess species while supplying the current style. ## Probably an uncommon use-case extendedMapSeqlevels("2", "NCBI", currentStyle = "TAIR10") ## Sequence names are unchanged when using an unsupported species extendedMapSeqlevels("seq2", "NCBI", "toyOrganism") ## Disable extendedMapSeqlevels. This can be useful when working with ## organisms that are not supported by GenomeInfoDb chrs <- c( "I", "II", "III", "IV", "IX", "V", "VI", "VII", "VIII", "X", "XI", "XII", "XIII", "XIV", "XV", "XVI", "XVII" ) extendedMapSeqlevels(chrs, chrsStyle = NULL) ## Not run: ## Set global species and style option options("chrsStyle" = "UCSC") options("species" = "homo_sapiens") ## Run using global options extendedMapSeqlevels("2") ## End(Not run)