compute_correlations: Compute the derived parameters correlations from a cubic...

In mcanouil/eggla: Early Growth Genetics Longitudinal Analysis

Compute the derived parameters correlations from a cubic splines mixed-effects model by time_model().

Description

Based on computed area under the curves (i.e., compute_aucs()) and slopes (i.e., compute_slopes()) for several intervals using a model fitted by time_model(), compute the correlations between each intervals derived parameters.

Usage

compute_correlations(
  fit,
  method,
  period = c(0, 0.5, 1.5, 3.5, 6.5, 10, 12, 17),
  knots = list(cubic_slope = NULL, linear_splines = c(0.75, 5.5, 11), cubic_splines =
    c(1, 8, 12))[[method]]
)

Arguments

fit	A model object from a statistical model such as from a call to time_model().
method	The type of model provided in fit, i.e., one of "cubic_slope", "linear_splines" or "cubic_splines".
period	The intervals knots on which AUCs are to be computed.
knots	The knots as defined fit and according to method.

Value

A list object with correlations between each intervals derived parameters.

Examples

data("bmigrowth")
ls_mod <- time_model(
  x = "age",
  y = "log(bmi)",
  cov = NULL,
  data = bmigrowth[bmigrowth[["sex"]] == 0, ],
  method = "linear_splines"
)
compute_correlations(
  fit = ls_mod,
  method = "linear_splines",
  period = c(0, 0.5, 1.5, 3.5, 6.5, 10, 12, 17)#,
  # knots = list(
  #   "cubic_slope" = NULL,
  #   "linear_splines" = c(0.75, 5.5, 11),
  #   "cubic_splines" = c(1, 8, 12)
  # )[[method]]
)

mcanouil/eggla documentation built on April 5, 2025, 3:06 a.m.