Limma moderated contrasts

Limma moderated contrasts

Limma moderated contrasts

See also

moderated_p_limma_long

Other modelling: AnovaExtractor, Contrasts, ContrastsDEqMSFacade, ContrastsDEqMSVoomFacade, ContrastsFirth, ContrastsFirthFacade, ContrastsLMFacade, ContrastsLMImputeFacade, ContrastsLMMissingFacade, ContrastsLimma, ContrastsLimmaFacade, ContrastsLimmaImputeFacade, ContrastsLimmaVoomFacade, ContrastsLimmaVoomImputeFacade, ContrastsLimpaFacade, ContrastsLmerFacade, ContrastsMissing, ContrastsModeratedDEqMS, ContrastsPlotter, ContrastsRLMFacade, ContrastsROPECA, ContrastsROPECAFacade, ContrastsTable, INTERNAL_FUNCTIONS_BY_FAMILY, LR_test(), Model, ModelFirth, ModelLimma, StrategyLM, StrategyLimma, StrategyLimpa, StrategyLmer, StrategyLogistf, StrategyRLM, build_contrast_analysis(), build_model(), build_model_glm_peptide(), build_model_glm_protein(), build_model_impute(), build_model_limma(), build_model_limma_impute(), build_model_limma_voom(), build_model_limma_voom_impute(), build_model_limpa(), build_model_logistf(), compute_borrowed_variance(), compute_borrowed_variance_limma(), compute_contrast(), compute_lmer_contrast(), contrasts_fisher_exact(), get_anova_df(), get_complete_model_fit(), get_p_values_pbeta(), group_label(), impute_refit_singular(), is_singular_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(), new_lm_imputed(), pivot_model_contrasts_to_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_limma(), strategy_limpa(), strategy_logistf(), summary_ROPECA_median_p.scaled()

Super class

prolfqua::ContrastsInterface -> ContrastsModerated

Public fields

Contrast

Class implementing the Contrast interface

model_name

name of model

subject_id

columns with subject_id (proteinID)

p.adjust

function to adjust p-values

Methods

Public methods

• ContrastsModerated$new()

• ContrastsModerated$get_contrast_sides()

• ContrastsModerated$get_linfct()

• ContrastsModerated$get_contrasts()

• ContrastsModerated$get_Plotter()

• ContrastsModerated$to_wide()

• ContrastsModerated$clone()

Inherited methods

• prolfqua::ContrastsInterface$column_description()
• prolfqua::ContrastsInterface$get_missing()

---

Method new()

initialize

Usage

ContrastsModerated$new(
  Contrast,
  model_name = paste0(Contrast$model_name, "_moderated"),
  p.adjust = prolfqua::adjust_p_values
)

Arguments

Contrast

class implementing the ContrastInterface

model_name

name of the model

p.adjust

function to adjust p-values - default BH

---

Method get_contrast_sides()

get both sides of contrasts

Usage

ContrastsModerated$get_contrast_sides()

---

Method get_linfct()

get linear functions from contrasts

Usage

ContrastsModerated$get_linfct(global = TRUE)

Arguments

global

logical TRUE - get the a linear functions for all models, FALSE - linear function for each model

---

Method get_contrasts()

applies limma moderation

Usage

ContrastsModerated$get_contrasts(all = FALSE)

Arguments

all

should all columns be returned (default FALSE)

global

use a global linear function (determined by get_linfct)

---

Method get_Plotter()

get ContrastsPlotter

Usage

ContrastsModerated$get_Plotter(fc_threshold = 1, fdr_threshold = 0.1)

Arguments

fc_threshold

fold change threshold to show in plots

fdr_threshold

FDR threshold to show in plots

---

Method to_wide()

convert to wide format

Usage

ContrastsModerated$to_wide(columns = c("p.value", "FDR", "statistic"))

Arguments

columns

value column default moderated.p.value

Returns

data.frame

---

Method clone()

The objects of this class are cloneable with this method.

Usage

ContrastsModerated$clone(deep = FALSE)

Arguments

deep

Whether to make a deep clone.

Examples

istar <- sim_lfq_data_protein_config(Nprot = 50)
#> creating sampleName from file_name column
#> completing cases
#> completing cases done
#> setup done
protIntensity <- istar$data
config <- istar$config


lProt <- LFQData$new(protIntensity, config)
lProt$rename_response("transformedIntensity")
modelFunction <-
  strategy_lm("transformedIntensity  ~ group_")
mod <- build_model(
 lProt,
 modelFunction)

Contr <- c("dil.b_vs_a" = "group_A - group_Ctrl")
contrast <- prolfqua::Contrasts$new(mod,
 Contr)
contrast <- ContrastsModerated$new(contrast)
bb <- contrast$get_contrasts()
#> determine linear functions:
#> get_contrasts -> contrasts_linfct
#> contrasts_linfct
#> Joining with `by = join_by(protein_Id, contrast)`
csi <- ContrastsMissing$new(lProt, contrasts = Contr)

contrast$get_contrasts() |> dim()
#> [1] 49 13
(xx <- csi$get_contrasts())   |> dim()
#> dil.b_vs_a=group_A - group_Ctrl
#> dil.b_vs_a=group_A - group_Ctrl
#> dil.b_vs_a=group_A - group_Ctrl
#> [1] 50 20
merged <- merge_contrasts_results(contrast, csi)
#> Joining with `by = join_by(protein_Id, contrast)`
#> Joining with `by = join_by(protein_Id, contrast)`

merged$more$get_contrasts() |> dim()
#> [1]  1 13
stopifnot(all(dim(merged$merged$get_contrasts() == c(50,13))))
stopifnot(all(dim(merged$same$get_contrasts()) == c(49,13)))

cs <- contrast$get_contrast_sides()
cslf <- contrast$get_linfct()
ctr <- contrast$get_contrasts()
ctrwide <- contrast$to_wide()
cp <- contrast$get_Plotter()

print(cp$histogram()$p.value, vp=NULL)

print(cp$histogram()$FDR, vp = NULL)


cp$volcano()
#> $FDR

#> 
cp$ma_plot()