Write optional disk outputs for longitudinal mixed-effects analyses

Source: R/methylationGLMM_T1T2_steps.R

Write optional serialized outputs, summary tables, significant interaction tables, and annotated results from the longitudinal mixed-effects workflow.

writeMethylationGLMM_T1T2Outputs(
  modelResults,
  modelSummaries,
  annotatedResults,
  significantInteractions = NULL,
  outputRData,
  summaryTxtDir,
  significantInteractionDir,
  annotatedLMEOut,
  saveTxtSummaries = TRUE,
  saveSignificantInteractions = FALSE,
  verbose = FALSE,
  logs = FALSE,
  log_dir = NULL,
  log_file = "log_methylationGLMM_T1T2.txt"
)

Arguments

modelResults

Object returned by fitMethylationGLMM_T1T2Models().

modelSummaries

Object returned by summarizeMethylationGLMM_T1T2Models().

annotatedResults

Object returned by annotateMethylationGLMM_T1T2Summaries() or a compatible data frame.

significantInteractions

Object returned by collectSignificantInteractionsMethylationGLMM_T1T2() or NULL.

outputRData

Character. Directory used for serialized model and summary outputs.

summaryTxtDir

Character. Directory used for tab-delimited summary tables.

significantInteractionDir

Character. Directory used for significant interaction coefficient tables.

annotatedLMEOut

Character. Directory used for the annotated summary CSV file.

saveTxtSummaries

Logical. If TRUE, write tab-delimited summary tables.

saveSignificantInteractions

Logical. If TRUE, write significant interaction coefficient tables.

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_paths" containing the paths of the files written to disk.

Examples

ex <- dnaEPICO:::exampleMethylationGLMMStateDnaEpico()
annotation_data <- annotateMethylationGLMM_T1T2Summaries(
  modelSummaries = ex$modelSummaries,
  annotationObject = ex$annotationData,
  annotationCols = "Name,chr,pos",
  verbose = FALSE,
  logs = FALSE
)
significant_hits <- collectSignificantInteractionsMethylationGLMM_T1T2(
  modelResults = ex$modelResults,
  pvalThreshold = 1,
  verbose = FALSE,
  logs = FALSE
)
output_paths <- writeMethylationGLMM_T1T2Outputs(
  modelResults = ex$modelResults,
  modelSummaries = ex$modelSummaries,
  annotatedResults = annotation_data,
  significantInteractions = significant_hits,
  outputRData = file.path(ex$tempDir, "models"),
  summaryTxtDir = file.path(ex$tempDir, "summary"),
  significantInteractionDir = file.path(ex$tempDir, "significant"),
  annotatedLMEOut = file.path(ex$tempDir, "annotated"),
  saveTxtSummaries = TRUE,
  saveSignificantInteractions = TRUE,
  verbose = FALSE,
  logs = FALSE
)
names(output_paths)
#> [1] "modelFiles"                  "summaryFiles"               
#> [3] "summaryTxtFiles"             "significantInteractionFiles"
#> [5] "annotatedLME"