normalize {CGHcall}R Documentation

Normalization and cellularity adjustment for arrayCGH data.

Description

This function normalizes arrayCGH data using the global mode or median. It can also adjust for the cellularity of your data.

Usage

normalize(input, method = "median", cellularity = 1, smoothOutliers = TRUE, ...)

Arguments

input Object of class cghRaw.
method Normalization method, either 'median', 'mode', or 'none'.
cellularity A vector of cellularities ranging from 0 to 1 to define the contamination of your sample with healthy cells (1 = no contamination). See details for more information.
smoothOutliers Logical. Indicates whether outliers should be smoothed using the smooth.CNA function.
... Arguments for smooth.CNA.

Details

The cellularity parameter should be a vector of length n where n is the number of samples in your dataset. The vector is recycled if there are not enough values in it, or truncated if there are too many. For more information on the correction we refer to section 1.6 of the supplementary information for van de Wiel et al. 2006.

Value

This function returns a dataframe in the same format as the input with normalized and/or cellularity adjusted log2 ratios.

Author(s)

Sjoerd Vosse & Mark van de Wiel

Examples

  data(WiltingData)
  ## Convert to cghRaw object
  cgh <- cghRaw(WiltingData)
  ## First preprocess the data
  raw.data <- preprocess(cgh)
  ## Simple global median normalization for samples with 75% tumor cells
  perc.tumor <- rep(0.75, 3)
  normalized.data <- normalize(raw.data, cellularity=perc.tumor)

[Package CGHcall version 2.0.0 Index]