A function for quickly calculating f statistics for use in sva
This function does simple linear algebra to calculate f-statistics for each row of a data matrix comparing the nested models defined by the design matrices for the alternative (mod) and and null (mod0) cases. The columns of mod0 must be a subset of the columns of mod.
fstats(dat, mod, mod0)
dat |
The transformed data matrix with the variables in rows and samples in columns |
mod |
The model matrix being used to fit the data |
mod0 |
The null model being compared when fitting the data |
fstats A vector of F-statistics one for each row of dat.
library(bladderbatch) data(bladderdata) dat <- bladderEset[1:50,] pheno = pData(dat) edata = exprs(dat) mod = model.matrix(~as.factor(cancer), data=pheno) mod0 = model.matrix(~1,data=pheno) fs <- fstats(edata, mod, mod0)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.