View source: R/GeneExpression.R
GetCounts | R Documentation |
GetCounts it estimate the gene counts wrapping the Rsubread::featureCounts function
ft <- GetCounts("my.aligned_bam", feature = "exon")
sbjBAM |
the BAM file name |
feature |
string indicating if exons or genes are returned |
annotSource |
string possible values are "GENCODE19" or "inbuilt". If GENECODE19, it will use GENOCDE from Arriba software, if "inbuilt" it will use as in |
nthreads |
integer, total of treads to use, default parallel::detectCores() -1 |
an extended object as returned by featureCounts
, the extension include in the targets slots the feature extracted (exon or gene)
## Not run:
test.subject <- GetArribaRTest()
bam.subject <- RunSTAR(test.subject)
ft <- GetCounts(bam.subject)
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.