sigdirPttest: Summarize results from 'pttest'

In orQA: Order Restricted Assessment Of Microarray Titration Experiments

Description

Translate results from pttest into trinomial coding where 1 signifies an upward trend, 0 a none significant result, and -1 a downward trend.

Usage

sigdirPttest(o,alpha)

Arguments

o	result object from pttest
alpha	desired alpha level for the test

Details

Translate results from e2test into trinomial coding where 1 signifies an upward trend, 0 a none significant result, and -1 a downward trend.

Value

a matrix with rows corresponding to the number of tests and columns corresponding to the number treatment groups coding each result as either 1,0, or -1.

Author(s)

Florian Klinglmueller <float_at_lefant.net>

Examples

groups <- rep(1:4,each=18)
ind <- rep(rep(1:3,each=6),4)
tdir <- sample((1:3)-2,100,rep=TRUE)
data <- matrix(rnorm(7200),nc=72)+(tdir %*% t(groups))
out <- pttest(data,groups,B=1000,rep=ind)
sigdir <- sigdirPttest(out)
head(sigdir)

orQA documentation built on May 1, 2019, 10:31 p.m.