exportBigWigs: Export bigwig files from values contained in a boost mode...

Description Usage Arguments Value Methods (by class) Author(s) Examples

Description

In boost mode, binding values of experiments are stored in the form of GRanges tables in the object. exportBigWigs output these values giving a proper bigwig file name corresponding to the transformations already performed.

Usage

1
2
3
4
5
6
7
8
9
exportBigWigs(theObject, verbose = TRUE)
            
            
## S4 method for signature 'ChIPSeqSpikeDatasetBoost'
exportBigWigs(theObject, verbose = TRUE)

## S4 method for signature 'ChIPSeqSpikeDatasetListBoost'
exportBigWigs(theObject, 
verbose = TRUE)

Arguments

theObject

ChIPSeqSpike dataset in boost mode (see ?spikeDataset)

verbose

If FALSE, do not output processing messages. Default is TRUE

Value

Output bigwig files of binding values.

The suffix of the bigwig file reflects the transformation steps performed on the object. If all steps were performed, the file name will be of the form: 'expName-RPM-BGSub-reverse-spiked.bw'. The suffixes 'RPM', 'BGSub', 'reverse' and 'spiked' stands for RPM scaling, input subtraction, RPM scaling reversal and exogenous scaling respectively.

Methods (by class)

Author(s)

Nicolas Descostes

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
## Mock example on a restricted number of reads
info_file_csv <- system.file("extdata/info.csv", package="ChIPSeqSpike")
bam_path <- system.file(c("extdata/bam_files"), package="ChIPSeqSpike")
bigwig_path <- system.file(c("extdata/bigwig_files"), package="ChIPSeqSpike")

if(.Platform$OS.type != 'windows') {
    csds <- spikeDataset(infoFile = info_file_csv, bamPath = bam_path, 
                         bigWigPath = bigwig_path, boost = TRUE)

    ## Creating test folder
    dir.create("./test_chipseqspike")

    csds <- estimateScalingFactors(csds)

    ## Apply RPM scaling
    csds <- scaling(csds, outputFolder = "test_chipseqspike")

    ## output RPM scaled files
    exportBigWigs(csds)

    ## Apply input subtraction
    csds <- inputSubtraction(csds)

    ## output input subtracted files
    exportBigWigs(csds)

    ## Delete all files generated in this example
    unlink("test_chipseqspike/", recursive = TRUE)
}

descostesn/ChIPSeqSpike documentation built on Aug. 6, 2019, 3:51 p.m.