proboscis {MAQCsubset}R Documentation

Produce a plot similar to Figure 2 of the Shippy MAQC paper (PMID 16964226).

Description

Produce a plot similar to Figure 2 of the Shippy MAQC paper (PMID 16964226).

Usage

proboscis(es, site=1, ABp=0.001, CDp=0.01, mmrad=100)

Arguments

es ExpressionSet instance with MAQC assay results
site numeric code – site to be assessed
ABp ABp – p-value threshold to declare concentration of gene in sample A to be different from ehe concentration in sample B
CDp CDp – p-value threshold to declare concentration of gene in sample C to be different from the concentration in sample D
mmrad numeric radius of the moving mean used to smooth the proportions differentially expressed

Details

Figure 2 of the Shippy paper consists of a collection of plots of estimated probabilities of self-consistent monotone titration – briefly, samples are such that A has 100% USRNA, B has 100% Ambion brain, C has 75% USRNA+25% brain, D has 25% USRNA, 75% brain. Self-consistent monotone titration holds for gene g if microarray measures for that gene satisfy A > C > D > B or B > C > D > A. The estimated probability functions look like a creature sticking its nose over a wall, thus the name of this function.

Value

an instance of proboStruct, for which a plot and lines method are available.

Author(s)

Vince Carey <stvjc@channing.harvard.edu>

References

PMID 16964226

Examples

data(afxsubRMAES)
NN2 = proboscis(afxsubRMAES, site=2)
plot(NN2)

[Package MAQCsubset version 1.0.13 Index]