orQA provides methods for quality assessment microarray titration data. These include permutation based trend and shape tests useful for the assessment of accuracy and cross platform consistency as well as methods for the estimation of variance components under order restrictions.

orQA provides methods for quality assessment of microarray titration
data. Exploiting the monotonic nature of such measurements accuracy,
precision and cross-platform agreement can be derived. For the assessment
of accuracy using shape tests see `pttest`

. To evaluate
precision using variance component estimates see
`est.lme`

. Order restricted inference of monotonic trends
for the purpose of cross-platform comparison is provided by
`e2test`

. For an example see below.

Florian Klinglmueller

Maintainer: Florian Klinglmueller <float_at_lefant.net>

Klinglmueller, F., Tuechler, T., Posch, M. (2010) "Cross Platform Comparison Of Microarray Data Using Order Restricted Inference" Under Review

Pinheiro, J., Bates, D., DebRoy, S., Sarkar, D.,R Development Core Team (2010) "nlme: Linear and Nonlinear Mixed Effects Models"

Barlow, R. E., Bartholomew, D. J., Bremner, J. M., and Brunk, H. D. (1972) "Statistical inference under order restrictions"; Wiley, London.

Robertson, T., Wright,F. T. and Dykstra, R. L. (1988) "Order Restricted Statistical Inference"; Wiley, New York.

Guo W., Sarkar SK., Peddada SD. (2010) "Controlling False Discoveries in Multidimensional Directional Decisions, with Applications to Gene Expression Data on Ordered Categories" Biometrics

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 | ```
## Let's assume a titration study with 10 genes 4 titration groups and
## 10 replicates which come from either of 2 biologically different
## samples (i.e. two times 4 technical replicates)
## WARNING: examples are run with a very limited number of permutations
## Titration groups, random factor
g <- rep(1:4,each=8)
r <- rep(rep(1:2,each=4),4)
## No differences (global null)
nulldata <- matrix(rnorm(320),nc=32)
## Differences between titration levels in each gene
altdata <- t(t(nulldata)+g)
## Accuracy are there any significantly non monotonous trends
res <- pttest(nulldata,g,1000,r) # apply shape test
## no significant trends at all
sigdirPttest(res)
## with alternatives
res <- pttest(altdata,g,1000,r) # apply shape test
## some significant trends
sigdirPttest(res)
## type of monotonicity (no significant trend, up, down,
## "anti-monotonous")
table(monotonicity(res))
## Precision - estimate variance components
res <- est.lme(nulldata,g,r)
round(apply(res,2,summary),2)
res <- est.lme(altdata,g,r)
round(apply(res,2,summary),2)
## Trend test and agreement
## some other dataset with some alternatives in both directions
tdir <- sample((1:3)-2,10,rep=TRUE)
altdata2 <- matrix(rnorm(320),nc=32)+ (tdir %*% t(g))
res1 <- e2test(altdata,g,1000,r)
res2 <- e2test(altdata2,g,1000,r)
## trinomial coding for directions (-1 down, 0 non sig., 1 up)
sigdir1 <- sigdirE2test(res1)
sigdir2 <- sigdirE2test(res2)
## contingency table of directional decisions
table(sigdir1,sigdir2)
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.