BamBigwig_to_chipProfile {profileplyr}R Documentation

BamBigwig_to_chipProfile

Description

Generate a soGGi ChIPprofile object with multiple BAM/bigWig files or multiple BED files as the input

Usage

BamBigwig_to_chipProfile(signalFiles, testRanges, format,
  style = "percentOfRegion", nOfWindows = 100, bin_size = 20, ...)

Arguments

signalFiles

paths to either BAM files or bigwig files. More than one path can be in this character vector, but all paths in one function call must point to be either all BAM files or all bigWig files, not a combination of the two.

testRanges

Either a character vector with paths to BED files.

format

character string of "bam", "bigwig", "RleList" or "PWM"

style

a character string, "percentOfRegion" (default) for normalised length divided into bins set by the 'nOfWindows' argument, "point" for per base pair plot where the number of base pairs per bin is set by the 'bin_size' argument, and "region" for combined plot

nOfWindows

The number of windows/bins the normalised ranges will be divided into if 'style' is set to 'percentOfRegion'. Default is 100.

bin_size

If 'style' is set to 'point' then this will determine the size of each bin over which signal is quantified. The default is 20 base pairs.

...

pass to regionPlot() within the soGGi package

Value

A profileplyr object

Examples

signalFiles <- c(system.file("extdata",
                             "Sorted_Hindbrain_day_12_1_filtered.bam",
                              package = "profileplyr"))
require(Rsamtools)
for (i in seq_along(signalFiles)){
 indexBam(signalFiles[i])
}
testRanges <- system.file("extdata", 
                          "newranges_small.bed", 
                          package = "profileplyr")
BamBigwig_to_chipProfile(signalFiles, 
                         testRanges, 
                         format = "bam",
                         paired=FALSE,
                         style="percentOfRegion",
                         )

[Package profileplyr version 1.0.1 Index]