man/examples/example_abYlm.Mlm.R

In abima: Adaptive Bootstrap Inference for Mediation Analysis with Enhanced Statistical Power

\donttest{
  ## Load libraries
  library(abima)

  set.seed(2)


  ## Set up parameters
  M.family <- gaussian()
  Y.family <- gaussian()

  simulation <- function(alpha_S = 0, beta_M = 0) {
    data <- generate_all_data(
      n = 200,
      alpha_S = alpha_S,
      beta_M = beta_M,
      M.family = M.family,
      Y.family = Y.family
    )
    S <- data$S
    M <- data$M
    Y <- data$Y
    X <- data$X

    out <- abYlm.Mlm(S, M, Y, X, lambda = 2, B = 199)
    out
  }


  simulation(1 / 8, 1 / 8)

  simulation(0, 0)



}