M-HPCA Simulation

Function for simulating data as described in the supplementary material of Campos et al. (202?)

MHPCA_simulation(
  sig_eps,
  n_d,
  J,
  D,
  num_reg,
  num_time,
  missing_level = FALSE,
  K = 2,
  L = 2,
  P = 2,
  M = 2
)

Arguments

sig_eps

error standard deviation (scalar)
n_d

group sample size (scalar)
J

number of repetitions per subject (scalar)
D

number of groups (scalar)
num_reg

number of regions (scalar)
num_time

number of functional points (scalar)
missing_level

allow for missing repetitions (logical)
K

number of marginal level 1 eigenvectors (scalar)
L

number of marginal level 1 eigenfunctions (scalar)
P

number of marginal level 2 eigenvectors (scalar)
M

number of marginal level 2 eigenfunctions (scalar)

Value

A list with the simulated components.

  • mu - The simulated mean function

  • eta - The simulated group-region-repetition shift from the mean

  • v_k - The simulated level 1 eigenvectors

  • phi_l - the simulated level 1 eigenfunctions

  • lambda_kl - The simulated level 1 eigenvalues

  • v_p - The simulated level 2 eigenvectors

  • phi_m - The simulated level 2 eigenfunctions

  • lambda_pm - The simulated level 2 eigenvalues

  • xi - The simulated subject-specific scores

  • zeta - The simulated subject/repetition-specific scores

  • data - The simulated dataframe with columns for

  • sig_eps - The simulated measurement error variance