Description Usage Arguments Value Author(s) Examples
The function quantify a aligned Bam file from RNA-Seq with user-specified genome regions. The reads from bam file will be filtered by the regions. The output can be counts or maxdepths for each region.
1 |
bam |
The aligned bam file from RNA-Seq experiment. |
GR |
|
method |
The method to quantify reads. It can be count, maxdepth or both. |
Reduce |
The methods to reduce the overlapped regions. The regions from different transcripts but the same gene can be overlapped. If it is "union", the overlapped regions from the same gene will be united together. If it is "intersection", the returned object will be the most shared regions among the transcripts in the same gene. If it is NULL, all the regions will be quantified separately. |
ovtype |
The overlap type that reads map to genes. Only the reads mapping with the type can be used in the quantification. By default, the reads that aligned within the input regions is used. The other reads are filtered out. See more details in "type" section of help(findOverlaps). |
revstrand |
For paired reads. If TRUE, the paired reads that are in the same direction will be filtered out. If FALSE, the direction of reads will be ignored. |
ingene |
For paired reads. If TRUE, the paired reads that are not in the same gene will be filtered out. If FALSE, all the paired reads will be used. |
... |
More options when read bam file. See more details in help(readBamGappedAlignments) |
The results return a GRanges object. If reduce option is used, the return GRanges object is reduced ranges by gene_id. Columns "counts" or/and "maxdepths" will be added to the elementMetadata of the GRanges object.
Qiang Hu
1 2 3 | ##---- Should be DIRECTLY executable !! ----
##-- ==> Define data, use random,
##-- or do help(data=index) for the standard data sets.
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.