Expand description
Nested F-test classification of t-statistics (limma classifyTestsF).
Pure-Rust port of the classifyTestsF routine in limma’s decidetests.R.
Given a matrix of t-statistics (genes × contrasts) it computes, per gene, an
overall moderated F-statistic and — via a step-down procedure — classifies
each individual contrast as up (1), down (-1) or non-significant (0).
The optional correlation matrix is the estimated correlation of the coefficients (equivalently the covariance of the t-statistics). When given, the F-statistic whitens by it via the eigendecomposition; when omitted the contrasts are treated as independent and the F-statistic is the mean square.
Functions§
- classify_
tests_ f - Classify each contrast per gene as
1(up),-1(down) or0via the nested step-down F-test.dfis the denominator degrees of freedom, recycled when length 1 (f64::INFINITYselects the chi-squared limit). - classify_
tests_ fstat - Overall moderated F-statistic per gene (limma’s
fstat.only=TRUE). Returns the F-statistics together with the numerator degrees of freedomdf1 = r; the denominator df2 is the caller’sdf, unchanged.