estimateScalingFactors {ChIPSeqSpike} | R Documentation |
Compute scaling factors for endogenous and exogenous experiment from 'Experiment', 'ExperimentLoaded', 'ChIPSeqSpikeDataset', 'ChIPSeqSpikeDatasetBoost', 'ChIPSeqSpikeDatasetList', and 'ChIPSeqSpikeDatasetListBoost'
estimateScalingFactors(theObject, paired = FALSE, verbose = TRUE) ## S4 method for signature 'Experiment' estimateScalingFactors(theObject, paired = FALSE, verbose = TRUE) ## S4 method for signature 'ExperimentLoaded' estimateScalingFactors(theObject, paired = FALSE, verbose = TRUE) ## S4 method for signature 'ChIPSeqSpikeDataset' estimateScalingFactors(theObject, paired = FALSE, verbose = TRUE) ## S4 method for signature 'ChIPSeqSpikeDatasetBoost' estimateScalingFactors(theObject, paired = FALSE, verbose = TRUE) ## S4 method for signature 'ChIPSeqSpikeDatasetList' estimateScalingFactors(theObject, paired = FALSE, verbose = TRUE) ## S4 method for signature 'ChIPSeqSpikeDatasetListBoost' estimateScalingFactors(theObject, paired = FALSE, verbose = TRUE)
theObject |
|
paired |
Indicate if sequences are single- or paired-ended. Default is FALSE |
verbose |
If False, do not output processing messages. Default is TRUE |
Estimating scaling factors is the first step to perform on a dataset. A scaling factor is defined as:
1/(bam_count/1000000)
bam_count being the number of reads aligned to the genome. The count is determined for the endogenous and exogenous experiments.
Scaling factors will the be applied to the bigwig files in the following steps of the procedure. After estimating scaling factors, RPM normalization should be performed.
Return an object of the same class of the input object containing computed scaling factors.
Experiment
: Method for signature theObject='Experiment'
ExperimentLoaded
: Method for signature theObject='ExperimentLoaded
'
ChIPSeqSpikeDataset
: Method for signature theObject=
'ChIPSeqSpikeDataset'
ChIPSeqSpikeDatasetBoost
: Method for signature theObject=
'ChIPSeqSpikeDatasetBoost'
ChIPSeqSpikeDatasetList
: Method for signature theObject=
'ChIPSeqSpikeDatasetList'
ChIPSeqSpikeDatasetListBoost
: Method for signature theObject=
'ChIPSeqSpikeDatasetListBoost'
Nicolas Descostes
spikeSummary
scaling
spikePipe
## 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") csds <- spikeDataset(infoFile = info_file_csv, bamPath = bam_path, bigWigPath = bigwig_path) csds <- estimateScalingFactors(csds)