Helper functions for scads

Helper functions for scads

compute_lfc_stats2(
  X,
  F,
  L,
  f0,
  D = matrix(rnorm(1000 * k), 1000, ncol(F)),
  U = matrix(runif(1000 * k), 1000, ncol(F)),
  M = matrix(sample(k, 1000 * k, replace = TRUE), 1000, k) - 1,
  lfc.stat = "le",
  conf.level = 0.68,
  rw = 0.3,
  e = 1e-15,
  verbose = TRUE
)