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,
      distanceAround = 100,
      ...,
      quant_params = NULL
    )

### 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` | A character vector with paths to BED files.
`format` | character string of "bam", "bigwig", "RleList" or "PWM"
`style` | a character string, "percentOfRegion" (default) for normalized 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.
`distanceAround` | If 'style' is 'percentOfRegion', then this controls the distance around the region that is included. Default is 100, meaning that a distance equal to 100 percent of that particular region on either side of the region will be included in the heatmap.
`...` | pass to regionPlot() within the soGGi package
`quant_params` | An optional `BiocParallelParam` instance determining the parallel back-end to be used during evaluation. When this argument is set to NULL (default) SerialParam() will be used. For parallelization, MulticoreParam() can be used.

### 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.4.3 Index]