Function to perform a One-way Anova statistical test on a MsnBase dataset
Source:R/anova_analysis.R
classic1wayAnova.RdFunction to perform a One-way Anova statistical test on a MsnBase dataset
Examples
# \donttest{
library(SummarizedExperiment)
data(subR25prot)
obj <- subR25prot
filter <- FunctionFilter('qMetacellOnConditions',
cmd = 'delete',
mode = 'AtLeastOneCond',
pattern = c("Missing POV", "Missing MEC"),
conds = design.qf(obj)$Condition,
percent = TRUE,
th = 0.8,
operator = '>')
obj <- filterFeaturesOneSE(obj, name = "Filtered", filters = list(filter))
qdata <- SummarizedExperiment::assay(obj[[3]])
conds <- design.qf(obj)$Condition
anova_tests <- apply(qdata, 1, classic1wayAnova, conditions = as.factor(conds))
anova_tests <- t(anova_tests)
# }