paCalls(object, method, ..., verbose=TRUE)
## S3 method for class 'ExonFeatureSet':
paCalls(object, method, verbose = TRUE)
## S3 method for class 'GeneFeatureSet':
paCalls(object, method, verbose = TRUE)
## S3 method for class 'ExpressionFeatureSet':
paCalls(object, method, ..., verbose = TRUE)
method
are 'DABG' (p-values for each probe) and 'PSDABG'
(p-values for each probeset). For Expression arrays, the only option
currently available for method
is 'MAS5'.
ABOUT MAS5 CALLS:
The additional arguments that can be passed to MAS5 are:
alpha1
: a significance threshold in (0, alpha2);alpha2
: a significance threshold in (alpha1, 0.5);tau
: a small positive constant;ignore.saturated
: if TRUE, do the saturation correction described in
the paper, with a saturation level of 46000;This function performs the hypothesis test:
H0: median(Ri) = tau, corresponding to absence of transcript H1: median(Ri) > tau, corresponding to presence of transcript
where Ri = (PMi - MMi) / (PMi + MMi) for each i a probe-pair in the probe-set represented by data.
The p-value that is returned estimates the usual quantity:
Pr(observing a more "present looking" probe-set than data | data is absent)
So that small p-values imply presence while large ones imply absence of transcript. The detection call is computed by thresholding the p-value as in:
call "P" if p-value < alpha1 call "M" if alpha1
Liu, W. M. and Mei, R. and Di, X. and Ryder, T. B. and Hubbell, E. and Dee,
S. and Webster, T. A. and Harrington, C. A. and Ho, M. H. and Baid, J. and
Smeekens, S. P. (2002) Analysis of high density expression microarrays with
signed-rank call algorithms, Bioinformatics, 18(12), pp. 1593--1599.
Liu, W. and Mei, R. and Bartell, D. M. and Di, X. and Webster, T. A. and
Ryder, T. (2001) Rank-based algorithms for analysis of microarrays,
Proceedings of SPIE, Microarrays: Optical Technologies and Informatics, 4266.
Affymetrix (2002) Statistical Algorithms Description Document, Affymetrix
Inc., Santa Clara, CA, whitepaper.
if (require(oligoData) & require(pd.huex.1.0.st.v2)){
data(affyExonFS)
## Get only 2 samples for example
dabgP = paCalls(affyExonFS[, 1:2])
dabgPS = paCalls(affyExonFS[, 1:2], "PSDABG")
head(dabgP) ## for probe
head(dabgPS) ## for probeset
}
Run the code above in your browser using DataLab