countReads: quantify reads in consecutive bins

Description Usage Arguments

View source: R/countReads.R

Description

quantify reads in consecutive bins

Usage

1
2
3
4
countReads(samplenames, gtf, fragmentLength = 150, bamFolder,
  outputDir = NA, modification = "m6A", binSize = 50,
  strandToKeep = "opposite", paired = FALSE, threads = 1,
  saveOutput = FALSE)

Arguments

samplenames

The names of each sample (prefix for bam files)

gtf

The gtf format gene annotation file

fragmentLength

The RNA fragment length (insert size of the library).

bamFolder

Path to the folder where bam file locates

outputDir

The directory to save output files

modification

The middle file name used to differentiate IP sample from Input sample. This should correspond to the file name of BAM files. For example, modification="m6A" for BAM files named as "samplenames.m6A.bam".

binSize

The size of consecutive bins to slice the transcripts

strandToKeep

According to library preparation protocol, choose which strand to count. Stranded RNA library usually seq the "ooposite" strand. Small RNA library seq the "same" strand.

paired

Logical indicating whether the input bam files are from paired end sequencing. Default is FALSE. If using paired end data, the read length will be estimated from the data and only good mate are counted.

threads

The number of threads to use for hyperthreading


scottzijiezhang/RADAR documentation built on Feb. 11, 2021, 8:35 p.m.