Description Usage Arguments Details Value Author(s) See Also Examples
Generic functions to extract genomic features from an object. This page documents the methods for OrganismDb objects only.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 | ## S4 method for signature 'MultiDb'
transcripts(x, columns=c("TXID", "TXNAME"), filter=NULL)
## S4 method for signature 'MultiDb'
exons(x, columns="EXONID", filter=NULL)
## S4 method for signature 'MultiDb'
cds(x, columns="CDSID", filter=NULL)
## S4 method for signature 'MultiDb'
genes(x, columns="GENEID", filter=NULL)
## S4 method for signature 'MultiDb'
transcriptsBy(x, by, columns, use.names=FALSE,
outerMcols=FALSE)
## S4 method for signature 'MultiDb'
exonsBy(x, by, columns, use.names=FALSE, outerMcols=FALSE)
## S4 method for signature 'MultiDb'
cdsBy(x, by, columns, use.names=FALSE, outerMcols=FALSE)
## S4 method for signature 'MultiDb'
getTxDbIfAvailable(x, ...)
## S4 method for signature 'MultiDb'
asBED(x)
## S4 method for signature 'MultiDb'
asGFF(x)
## S4 method for signature 'MultiDb'
disjointExons(x, aggregateGenes=FALSE,
includeTranscripts=TRUE, ...)
## S4 method for signature 'MultiDb'
microRNAs(x)
## S4 method for signature 'MultiDb'
tRNAs(x)
## S4 method for signature 'MultiDb'
promoters(x, upstream=2000, downstream=200, use.names=TRUE, ...)
## S4 method for signature 'GenomicRanges,MultiDb'
distance(x, y, ignore.strand=FALSE,
..., id, type=c("gene", "tx", "exon", "cds"))
## S4 method for signature 'BSgenome'
extractTranscriptSeqs(x, transcripts, strand = "+")
## S4 method for signature 'MultiDb'
extractUpstreamSeqs(x, genes, width=1000, exclude.seqlevels=NULL)
## S4 method for signature 'MultiDb'
intronsByTranscript(x, use.names=FALSE)
## S4 method for signature 'MultiDb'
fiveUTRsByTranscript(x, use.names=FALSE)
## S4 method for signature 'MultiDb'
threeUTRsByTranscript(x, use.names=FALSE)
## S4 method for signature 'MultiDb'
isActiveSeq(x)
|
x |
A MultiDb object. Except for the extractTranscriptSeqs method. In that case it's a BSgenome object and the second argument is an MultiDb object. |
... |
Arguments to be passed to or from methods. |
by |
One of |
columns |
The columns or kinds of metadata that can be retrieved from the
database. All possible columns are returned by using the |
filter |
Either |
use.names |
Controls how to set the names of the returned
GRangesList object.
These functions return all the features of a given type (e.g.
all the exons) grouped by another feature type (e.g. grouped by
transcript) in a GRangesList object.
By default (i.e. if Finally, |
upstream |
For |
downstream |
For |
aggregateGenes |
For |
includeTranscripts |
For |
y |
For |
id |
A |
type |
A |
ignore.strand |
A |
outerMcols |
A |
transcripts |
An object representing the exon ranges of each transcript to extract.
It must be a GRangesList or MultiDb
object while the |
strand |
Only supported when Can be an atomic vector, a factor, or an Rle object,
in which case it indicates the strand of each transcript (i.e. all the
exons in a transcript are considered to be on the same strand).
More precisely: it's turned into a factor (or factor-Rle)
that has the "standard strand levels" (this is done by calling the
|
genes |
An object containing the locations (i.e. chromosome name, start, end, and
strand) of the genes or transcripts with respect to the reference genome.
Only GenomicRanges and MultiDb objects
are supported at the moment. If the latter, the gene locations are obtained
by calling the |
width |
How many bases to extract upstream of each TSS (transcription start site). |
exclude.seqlevels |
A character vector containing the chromosome names (a.k.a. sequence levels) to exclude when the genes are obtained from a MultiDb object. |
These are the range based functions for extracting transcript information from a MultiDb object.
a GRanges or GRangesList object
M. Carlson
MultiDb-class for how to use the
simple "select" interface to extract information from a
MultiDb
object.
transcripts for the original
transcripts
method and related methods.
transcriptsBy for the original
transcriptsBy
method and related methods.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 | ## extracting all transcripts from Homo.sapiens with some extra metadata
library(Homo.sapiens)
cols = c("TXNAME","SYMBOL")
res <- transcripts(Homo.sapiens, columns=cols)
## extracting all transcripts from Homo.sapiens, grouped by gene and
## with extra metadata
res <- transcriptsBy(Homo.sapiens, by="gene", columns=cols)
## list possible values for columns argument:
columns(Homo.sapiens)
## Get the TxDb from an MultiDb object (if it's available)
getTxDbIfAvailable(Homo.sapiens)
## Other functions listed above should work in way similar to their TxDb
## counterparts. So for example:
promoters(Homo.sapiens)
## Should give the same value as:
promoters(getTxDbIfAvailable(Homo.sapiens))
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.