classifyTests {limma} | R Documentation |
Classify a series of related t-statistics as up, down or not significant.
classifyTests(object, cor.matrix=NULL, design=NULL, contrasts=diag(ncol(design)), df=Inf, p.value=0.01, fstat.only=FALSE)
object |
numeric matrix of t-statistics or an MArrayLM object from which the t-statistics may be extracted. |
cor.matrix |
covariance matrix of each row of t-statistics. Defaults to the identity matrix. |
design |
full rank numeric design matrix. Not used if cor.matrix is specified. |
contrasts |
numeric matrix with columns specifying contrasts of the coefficients of interest. Not used if cor.matrix is specified. |
df |
numeric vector giving the degrees of freedom for the t-statistics.
May have length 1 or length equal to the number of rows of tstat . |
p.value |
numeric value between 0 and 1 giving the desired size of the test |
fstat.only |
logical, if TRUE then return the overall F-statistic as for FStat instead of classifying the test results |
FStat
computes the F-statistic for testing all the contrasts equal to zero.
It is equivalent to classifyTests
with fstat.only=TRUE
.
classifyTests
classifies using a nested F-test approach giving particular attention to correctly classifying genes which have two or more significant t-statistics, i.e., are differential expressed under two or more conditions.
The overall F-statistic used is computed by FStat
.
At least one constrast will be classified as significant if and only if the overall F-statistic is significant.
classifyTestsT
and classifyTestsP
implement simpler classification schemes based on threshold or critical values for the individual t-statistics in the case of classifyTestsT
or p-values obtained from the t-statistics in the case of classifyTestsP
.
Rows of tstat
correspond to genes and columns to coefficients or contrasts.
For each row of tstat
, F-statistics are constructed from the t-statistics.
If the overall F-statistic is significant, then the function makes a best choice as to which t-statistics contributed to this result.
The methodology is based on the principle that any t-statistic should be called significant if the F-test is still significant for that row when all the larger t-statistics are set to the same absolute size as the t-statistic in question.
If tstat
is an MArrayLM
object, then all arguments except for p.value
are extracted from it.
cor.matrix
is the same as the correlation matrix of the coefficients from which the t-statistics are calculated.
If cor.matrix
is not specified, then it is calculated from design
and contrasts
if at least design
is specified or else defaults to the identity matrix.
In terms of design
and contrasts
, cor.matrix
is obtained by standardizing the matrix
t(contrasts) %*% solve(t(design) %*% design) %*% contrasts
to a correlation matrix.
An object of class TestResults
.
This is essentially a numeric matrix with elements -1
, 0
or 1
depending on whether each t-statistic is classified as significantly negative, not significant or significantly positive respectively.
Gordon Smyth
An overview of linear model functions in limma is given by 5.LinearModels.