Create limma modelling strategy

Analogous to strategy_lm but for limma's matrix-based pipeline. Returns a list consumed by build_model_limma.

Usage

strategy_limma(modelstr, model_name = "limma", trend = FALSE, robust = FALSE)

Arguments

modelstr

model formula as string (e.g. "abundance ~ group_")

model_name

name of model

trend

logical, passed to eBayes

robust

logical, passed to eBayes

See also

Other modelling: Contrasts, ContrastsFirth, ContrastsLimma, ContrastsMissing, ContrastsModerated, ContrastsModeratedDEqMS, ContrastsPlotter, ContrastsProDA, ContrastsROPECA, ContrastsTable, INTERNAL_FUNCTIONS_BY_FAMILY, LR_test(), Model, ModelFirth, ModelLimma, build_model(), build_model_limma(), build_model_logistf(), contrasts_fisher_exact(), get_anova_df(), get_complete_model_fit(), get_p_values_pbeta(), group_label(), isSingular_lm(), linfct_all_possible_contrasts(), linfct_factors_contrasts(), linfct_from_model(), linfct_matrix_contrasts(), merge_contrasts_results(), model_analyse(), model_summary(), moderated_p_deqms(), moderated_p_deqms_long(), moderated_p_limma(), moderated_p_limma_long(), my_contest(), my_contrast(), my_contrast_V1(), my_contrast_V2(), my_glht(), pivot_model_contrasts_2_Wide(), plot_lmer_peptide_predictions(), sim_build_models_lm(), sim_build_models_lmer(), sim_build_models_logistf(), sim_make_model_lm(), sim_make_model_lmer(), strategy_logistf(), summary_ROPECA_median_p.scaled()

Examples

strat <- strategy_limma("abundance ~ group_")
strat$formula
#> abundance ~ group_
#> <environment: 0x562ab64f3ee8>
strat$model_name
#> [1] "limma"