analyzeDuprates: Read in a BAM file and count the tags falling on the features...

Description Usage Arguments Details Value Examples

View source: R/analyzeDuprates.R

Description

analyzeDuprates returns a data.frame with tag counts

Usage

1
2
analyzeDuprates(bam, gtf, stranded = 0, paired = FALSE, threads = 1,
  verbose = FALSE, ...)

Arguments

bam

The bam file containing the duplicate-marked reads

gtf

The gtf file describing the features

stranded

Whether the reads are strand specific

paired

Paired end experiment?

threads

The number of threads to be used for counting

verbose

Whether to output Rsubread messages into the console

...

Other params sent to featureCounts

Details

This function makes use of the Rsubread package to count tags on the GTF features in different scenarios. The scenarios are the 4 possible combinations of allowing multimappers (yes/no) and duplicate reads (yes/no).

Value

A data.frame with counts on features, with and without taking into account multimappers/duplicated reads

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
bam <- system.file("extdata",
                   "wgEncodeCaltechRnaSeqGm12878R1x75dAlignsRep2V2_duprm.bam",
                   package="dupRadar")
gtf <- system.file("extdata","genes.gtf",package="dupRadar")
stranded <- 2    # '0' (unstranded), '1' (stranded) and '2' (reverse)
paired   <- FALSE
threads  <- 4

# call the duplicate marker and analyze the reads
dm <- analyzeDuprates(bam,gtf,stranded,paired,threads)

Bioconductor-mirror/dupRadar documentation built on June 1, 2017, 7:37 a.m.