Collect significant longitudinal model terms from fitted mixed-effects models — collectSignificantInteractionsMethylationGLMM

Collect raw coefficient tables for CpGs whose phenotype main effect or requested interaction p-value passes the chosen threshold.

collectSignificantInteractionsMethylationGLMM_T1T2(
  modelResults,
  pvalThreshold = 0.05,
  interactionTerm = NULL,
  verbose = FALSE,
  logs = FALSE,
  log_dir = NULL,
  log_file = "log_methylationGLMM_T1T2.txt"
)

Arguments

modelResults: Object returned by fitMethylationGLMM_T1T2Models().
pvalThreshold: Numeric. Threshold applied to the extracted phenotype or interaction p-values.
interactionTerm: Character or NULL. Optional interaction term. When NULL, phenotype main effects are used.
verbose: Logical. If TRUE, emit progress messages with message().
logs: Logical. If TRUE, write the same messages to a log file.
log_dir: Character or NULL. Directory used for the log file when logs = TRUE.
log_file: Character. File name used when logs = TRUE.

Value

A list with class "dnaEPICO_methylationGLMM_T1T2_significant" containing the retained coefficient tables for each phenotype.

Examples

ex <- dnaEPICO:::exampleMethylationGLMMStateDnaEpico()
significant_hits <- collectSignificantInteractionsMethylationGLMM_T1T2(
  modelResults = ex$modelResults,
  pvalThreshold = 1,
  verbose = FALSE,
  logs = FALSE
)
names(significant_hits)
#> [1] "score"