Object
~~|
~~+--ParametersInterface
~~~~~~~|
~~~~~~~+--MultiSourceCopyNumberNormalization
Directly known subclasses:
public static class MultiSourceCopyNumberNormalization
extends ParametersInterface
The multi-source copy-number normalization (MSCN) method [1] is a
normalization method that normalizes copy-number estimates measured
by multiple sites and/or platforms for common samples. It normalizes the
estimates toward a common scale such that for any copy-number level
the mean level of the normalized data are the same.MultiSourceCopyNumberNormalization(dsList=NULL, fitUgp=NULL, subsetToFit=NULL, targetDimension=1, align=c("byChromosome", "none"), tags="*", ...)list of K AromaUnitTotalCnBinarySet:s.AromaUgpFile that specifies the
common set of loci used to normalize the data sets at.fitUgp
object) to be used to fit the normalization functions.
If NULL, loci on chromosomes 1-22 are used, but not on ChrX and ChrY.character specifying type of alignment applied, if any.
If "none", no alignment is done.
If "byChromosome", the signals are shifted chromosome
by chromosome such tgetAllNames Gets the names of all unique samples across all sources.
getAsteriskTags -
getInputDataSets Gets the list of data sets to be normalized.
getOutputDataSets -
getTags -
nbrOfDataSets -
process Normalizes all samples.
}Methods inherited from ParametersInterface: getParameterSets, getParameters, getParametersAsString
Methods inherited from Object: $, $<-, [[, [[<-, as.character, attach, attachLocally, clearCache, clearLookupCache, clone, detach, equals, extend, finalize, gc, getEnvironment, getFieldModifier, getFieldModifiers, getFields, getInstantiationTime, getStaticInstance, hasField, hashCode, ll, load, objectSize, print, registerFinalizer, save, asThis
This means that if the transformation applied by a particular preprocessing method is not the same for the sex chromosomes as the autosomal chromosomes, the normalization applied on the sex chromosomes is not optimal one. This is why multi-source normalization sometimes fails to bring sex-chromosome signals to the same scale across sources. Unfortunately, there is no automatic way to handle this. The only way would be to fit a specific normalization function to each of the sex chromosomes, but that would require that there exist copy-number abberations on those chromosomes, which could be a too strong assumption.
A more conservative approach is to normalize the signals such that
afterward the median of the smoothed copy-number levels are the same
across sources for any particular chromosome.
This is done by setting argument align="byChromosome".
However, the current implementation is such that it first generates smoothed data for all samples/arrays. Then, it normalizes the sample one by one.