Skip to contents

Sweep Envirome Size and Quality (Power Analysis)

Source: R/env_prediction.R
sim_W_grid.Rd

Accessory to sim_W. Repeatedly simulates enviromes across a grid of noise, n_var and collinearity, returning the mean recovery and effective rank at each design point.

Usage

sim_W_grid(
  C_env,
  noise = seq(0, 1, 0.25),
  n_var = c(10, 50, 200),
  collinearity = c(0, 0.5, 0.9),
  n_rep = 10,
  calibrate = TRUE,
  cal.nrep = 15L,
  seed = NULL,
  verbose = TRUE
)

Arguments

C_env

q x q correlation AMONG ENVIRONMENTS (see sim_W).

noise, n_var, collinearity

numeric vectors defining the grid.

n_rep

integer. Replicate simulations per design point.

calibrate

logical. Passed to sim_W.

cal.nrep

integer. Calibration replicates per sim_W call.

seed

integer. RNG seed; each design point gets a reproducible sub-seed.

verbose

logical. Print per-design-point progress. Default TRUE.

Value

A data.frame with one row per design point: noise, n_var, collinearity, mean recovery r and its SD r_sd, r2, and mean eff_rank.

See also

sim_W, sim_met_C

Examples

# \donttest{
C <- sim_met_C(q = 10, min.cor = 0.2, max.cor = 0.8, seed = 1)
#> Warning: Rejection sampling failed in 200 tries; C was repaired by eigenvalue clamping.
#>   Requested off-diagonal range [0.200, 0.800]; realised [0.235, 0.712].
grid <- sim_W_grid(C, noise = c(0, 0.3, 0.6), n_var = c(20, 100),
                   collinearity = c(0, 0.8), n_rep = 5)
#>   [1/12] noise=0.00  n_var=20  collinearity=0.00
#>   [2/12] noise=0.30  n_var=20  collinearity=0.00
#>   [3/12] noise=0.60  n_var=20  collinearity=0.00
#>   [4/12] noise=0.00  n_var=100  collinearity=0.00
#>   [5/12] noise=0.30  n_var=100  collinearity=0.00
#>   [6/12] noise=0.60  n_var=100  collinearity=0.00
#>   [7/12] noise=0.00  n_var=20  collinearity=0.80
#>   [8/12] noise=0.30  n_var=20  collinearity=0.80
#> Warning: Calibration failed; using noise as the raw variance share.
#> Warning: Calibration failed; using noise as the raw variance share.
#> Warning: Calibration failed; using noise as the raw variance share.
#> Warning: Calibration failed; using noise as the raw variance share.
#> Warning: Calibration failed; using noise as the raw variance share.
#>   [9/12] noise=0.60  n_var=20  collinearity=0.80
#> Warning: Calibration failed; using noise as the raw variance share.
#> Warning: Calibration failed; using noise as the raw variance share.
#> Warning: Calibration failed; using noise as the raw variance share.
#> Warning: Calibration failed; using noise as the raw variance share.
#> Warning: Calibration failed; using noise as the raw variance share.
#>   [10/12] noise=0.00  n_var=100  collinearity=0.80
#>   [11/12] noise=0.30  n_var=100  collinearity=0.80
#> Warning: Calibration failed; using noise as the raw variance share.
#> Warning: Calibration failed; using noise as the raw variance share.
#> Warning: Calibration failed; using noise as the raw variance share.
#> Warning: Calibration failed; using noise as the raw variance share.
#> Warning: Calibration failed; using noise as the raw variance share.
#>   [12/12] noise=0.60  n_var=100  collinearity=0.80
#> Warning: Calibration failed; using noise as the raw variance share.
#> Warning: Calibration failed; using noise as the raw variance share.
#> Warning: Calibration failed; using noise as the raw variance share.
#> Warning: Calibration failed; using noise as the raw variance share.
#> Warning: Calibration failed; using noise as the raw variance share.
# }