makeCgh: Constructs a 'cghRaw', 'cghSeg', or 'cghCall' object

Description Usage Arguments Value Author(s) Examples

Description

Constructs a 'cghRaw', 'cghSeg', or 'cghCall' object.

Usage

1
makeCgh(object, filter=TRUE, chromosomeReplacements=c(X = 23, Y = 24, MT = 25), ...)

Arguments

object

A QDNAseqCopyNumbers object.

filter

If TRUE, bins are filtered, otherwise not.

chromosomeReplacements

A named integer vector of chromosome name replacements to be done. QDNAseq stores chromosome names as characters, but CGHcall expects them to be integers. Defaults to c(X=23, Y=24, MT=25) for human. Value of "auto" will use running numbers in order of appearance in the bin annotations.

...

Not used.

Value

Returns a cghRaw if the object has not been segmented, a cghSeg if it has been segmented but not called, or cghCall if it has been called as well.

Author(s)

Ilari Scheinin

Examples

1
2
3
4
5
6
7
8
data(LGG150)
readCounts <- LGG150
readCountsFiltered <- applyFilters(readCounts)
readCountsFiltered <- estimateCorrection(readCountsFiltered)
copyNumbers <- correctBins(readCountsFiltered)
copyNumbersNormalized <- normalizeBins(copyNumbers)
copyNumbersSmooth <- smoothOutlierBins(copyNumbersNormalized)
cgh <- makeCgh(copyNumbersSmooth)

Example output

38,819	total bins
38,819	of which in selected chromosomes
36,722	of which with reference sequence
33,347	final bins
Calculating correction for GC content and mappability
    Calculating fit for sample LGG150 (1 of 1) ...
Done.
Applying median normalization ...
Smoothing outliers ...

QDNAseq documentation built on Nov. 8, 2020, 6:57 p.m.