propexpr: Estimate Proportion of Expressed Probes

In richierocks/limma2: Linear Models for Microarray Data

Description

Estimate the proportion of microarray probes which are expressed in each array.

Usage

propexpr(x, neg.x=NULL, status=x$genes$Status, labels=c("negative","regular"))

Arguments

x	matrix or similar object containing raw intensities for a set of arrays.
neg.x	matrix or similar object containing raw intensities for negative control probes for the same arrays. If NULL, then negative controls must be provided in x.
status	character vector giving probe types.
labels	character vector giving probe type identifiers.

Details

This function estimates the proportion of expressed in a microarray by utilizing the negative control probes. Illumina BeadChip arrays contain 750~1600 negative control probes. The expression profile of these control probes can be saved to a separate file by the Illumina BeadStudio software when using it to output the expression profile for regular probes. The control probe profile could be re-generated if it was not generated when the regular probe profile was created by BeadStudio. Other microarray platforms can also use this function to estimate the proportion of expressed probes in each array, provided that they have a set of negative control probes.

labels can include one or two probe type identifiers. Its first element should be the identifier for negative control probes (negative by default). If labels only contains one identifier, then it will be assumed to contain the identifier for negative control probes. By default, regular is the identifier for regular probes.

Value

Numeric vector giving the proportions of expressed probes in each array.

Author(s)

Wei Shi and Gordon Smyth

References

Shi, W, de Graaf, C, Kinkel, S, Achtman, A, Baldwin, T, Schofield, L, Scott, H, Hilton, D, Smyth, GK (2010). Estimating the proportion of microarray probes expressed in an RNA sample. Nucleic Acids Research 38, 2168-2176.

See Also

Description to the control probes in Illumina BeadChips can be found in read.ilmn.

Examples

## Not run: 
x <- read.ilmn(files="sample probe profile.txt",ctrlfiles="control probe profile.txt")
propexpr(x, )

## End(Not run)

richierocks/limma2 documentation built on May 27, 2019, 8:47 a.m.