bam2rpm: Compute the normalized number of reads (rpm) fall into...

Description Usage Arguments Value Examples

View source: R/bam2rpm.R

Description

bam2rpm computes the normalized number of reads (rpm) fall into specific genomic region such as promoter, enhancer, genebody

Usage

1
bam2rpm(bamFile, region, fragLength = 180)

Arguments

bamFile

Aligned bam file as input.

region

The GRanges object defined by user to calculate the number of reads fall into this specific region. For ChIP-Seq of histone modifications they are usually promoter, enhancer and genebody regions.

fragLength

Extend reads toward the 3'-end to the average DNA fragment size obtained after DNA size selection

Value

a vector of numbers

Examples

1
2
3
data("promoter")
file.bam <- system.file("extdata", "SRR925640.bam", package = "intePareto")
bam2rpm(bamFile = file.bam, region = promoter, fragLength = 180)

intePareto documentation built on June 18, 2021, 5:08 p.m.