Estimate Proportion of Expressed Probes
Estimate the proportion of microarray probes which are expressed in each array.
propexpr(x, neg.x=NULL, status=x$genes$Status, labels=c("negative","regular"))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 |
status |
character vector specifying control type of each probe. Only used if |
labels |
character vector giving the |
This function estimates the overall proportion of probes on each microarray that are correspond to expressed genes using the method of Shi et al (2010). The function is especially useful for Illumina BeadChips arrays, although it can in principle be applied to any platform with good quality negative controls.
The negative controls can be supplied either as rows of x or as a separate matrix.
If supplied as rows of x, then the negative controls are identified by the status vector.
x might also include other types of control probes, but these will be ignored in the calculation.
Illumina BeadChip arrays contain 750~1600 negative control probes.
If read.idat is used to read Illumina expression IDAT files, then the control probes will be populated as rows of the output EListRaw object, and the vector x$genes$Status will be set to identify control probes.
Alternatively, expression values can be exported from Illumina's GenomeStudio software as tab-delimited text files. In this case, the control probes are usually written to a separate file from the regular probes.
Numeric vector giving the proportions of expressed probes in each array.
Wei Shi and Gordon Smyth
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(7), 2168-2176. https://www.ncbi.nlm.nih.gov/pubmed/20056656
Description to the control probes in Illumina BeadChips can be found in read.ilmn.
## Not run:
# Read Illumina binary IDAT files
x <- read.idat(idat, bgx)
propexpr(x)
# Read text files exported from GenomeStudio
x <- read.ilmn(files = "sample probe profile.txt",
ctrlfiles = "control probe profile.txt")
propexpr(x)
## End(Not run)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.