Rank Invariant Normalization

Share:

Description

This function basically adjusts the samples to the same background level and then optionally scales to the same foreground level.

Usage

1
rankinvariant(x.lumi, targetArray = NULL, rrc = .05, lowRank = seq(.5, .25, -.05), highRank = .9, minSize = .02, maxit=200)

Arguments

x.lumi

an ExpressionSet inherited object or a data matrix with columns as samples and rows as genes

targetArray

A target chip is the model for other chips to normalize. It can be a column index, a vector or a LumiBatch object with one sample.

rrc

The relative rank change allowed for a gene to be selected as rank invariant

lowRank

A vector with, in decreasing order, the minimum ranks where candidate genes can be selected as rank invariant

highRank

The maximum rank where candidate genes can be selected as rank invariant

minSize

Fraction of genes required to be selected as rank invariant

maxit

Maximum number of iterations for rlm to reach convergence

Details

Rank invariant normalization uses a set of genes that are rank invariant between a given sample and a target sample. The target sample can be predefined by setting the targetArray argument. If targetArray is NULL the average expression of all samples will be the target. Rank invariant genes are found for each sample seperately by calculation the relative rank change for each gene. Furthermore, only genes with ranks between the lowRank and highRank are considered. If the number of probes is less than minSize multiplies by the number of genes the next lowRank value tried. If no rank invariant set can be found an error is thrown.

The default settings of this function are the same as used Genomstudio (Illumina). The results produced by this method are similar, but not identical to Genomestudio.

Value

Return an object with expression values normalized. The class of the return object is the same as the input object x.lumi.

Author(s)

Arno Velds (contact: a.velds (at) nki.nl)

See Also

lumiN

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.