Variance estimation from replicate estimates.

method selects the rescaling: "JK1", "JKn", "BRR", "Fay", "bootstrap", "SDR".

Usage

morie_weights_replicate_variance(
  full_estimate,
  replicate_estimates,
  method = c("JK1", "JKn", "BRR", "Fay", "bootstrap", "SDR"),
  fay_coefficient = 0,
  strata = NULL
)

Arguments

full_estimate

Numeric scalar; full-sample point estimate.

replicate_estimates

Numeric vector of replicate point estimates.

method

Character; calibration / smoothing / variance method. Allowed values depend on the caller.

fay_coefficient

Fay's coefficient (\(k\)) for Fay's-BRR.

strata

Optional vector of stratum identifiers aligned with weights.

Value

A named list (see Details).

Examples

set.seed(1)
reps <- rnorm(50, 1.2, 0.05)
res <- morie_weights_replicate_variance(1.2, reps, method = "bootstrap")
str(res, max.level = 1)
#> List of 4
#>  $ variance: num 0.00173
#>  $ se      : num 0.0416
#>  $ ci_lower: num 1.12
#>  $ ci_upper: num 1.28