Add numerical integration points for unanchored comparison

Generate quasi-Monte Carlo integration points using Sobol sequences and a Gaussian copula to account for correlations between covariates in the AgD.

Usage

add_integration_unanchored(
  data,
  n_int = 64,
  cor = NULL,
  cor_adjust = NULL,
  ...
)

Arguments

data

An unanchored_data object from combine_unanchored

n_int

Number of integration points (default 64, powers of 2 recommended)

cor

Correlation matrix for covariates. If NULL, computed from IPD.

cor_adjust

Adjustment method: "spearman", "pearson", or "none"

...

Distribution specifications for each covariate using distr

Value

An unanchored_data object with integration points added

Examples

if (FALSE) { # \dontrun{
network <- add_integration_unanchored(
  network,
  n_int = 64,
  x1 = distr(qnorm, mean = x1_mean, sd = x1_sd),
  x2 = distr(qbern, prob = x2_mean)
)
} # }