BamBigwig_to_chipProfile {profileplyr} | R Documentation |
Generate a soGGi ChIPprofile object with multiple BAM/bigWig files or multiple BED files as the input
BamBigwig_to_chipProfile(signalFiles, testRanges, format, style = "percentOfRegion", nOfWindows = 100, bin_size = 20, ...)
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 |
A profileplyr object
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", )