Calculate Absolute Copy Number from Sequencing Counts
Description
Taking into account mappability and GC content biases, the absolute copy number is calculated, by assuming that the median read depth is a copy number of 1.
Usage
1 2 3 4 5 6  ## S4 method for signature 'data.frame,matrix,GCAdjustParams'
GCadjustCopy(input.windows, input.counts,
gc.params, ...)
## S4 method for signature 'GRanges,matrix,GCAdjustParams'
GCadjustCopy(input.windows, input.counts,
gc.params, verbose = TRUE)

Arguments
input.windows 
A 
input.counts 
A matrix of counts. Rows are genomic windows and columns are samples. 
gc.params 
A 
... 

verbose 
Whether to print the progess of processing. 
Details
First, the mappability of all counting windows is calculated, and windows that have mappability less than the cutoff specified by in the parameters object are ignored in further steps. The remaining windows have their counts scaled by multiplying their counts by 100 / percentage mappability.
The range of GC content of the counting windows is broken into a number of bins, as specified by the user in the parameters object. A probability density function is fitted to the counts in each bin, so the mode can be found. The mode is taken to be the counts of the copy neutral windows, for that GC content bin.
A polynomial function is fitted to the modes of GC content bins. Each count is divided by its expected counts from the polynomial function to give an absolute copy number estimate. If the ploidy has been provided in the parameters object, then all counts within a sample are multiplied by the ploidy for that sample. If the sample ploidys were omitted, then no scaling for ploidy is done.
Value
A AdjustedCopyEstimate
object describing the input windows and their
estimates.
Author(s)
Dario Strbenac
Examples
1 2 3 4 5 6 7 8 9 10 11 12 13 14  ## Not run:
library(BSgenome.Hsapiens.UCSC.hg18)
library(BSgenome.Hsapiens36bp.UCSC.hg18mappability)
load("inputsReads.RData")
windows < genomeBlocks(Hsapiens, chrs = paste("chr", c(1:22, 'X', 'Y'), sep = ''),
width = 20000)
counts < annotationBlocksCounts(inputsReads, anno = windows, seq.len = 300)
gc.par < GCAdjustParams(genome = Hsapiens, mappability = Hsapiens36bp,
min.mappability = 50, n.bins = 10, min.bin.size = 10,
poly.degree = 4, ploidy = c(2, 4))
abs.cn < GCadjustCopy(input.windows = windows, input.counts = counts, gc.params = gc.par)
## End(Not run)
