Each stratum is split into two halves; signs from a random Hadamard-like
matrix double one half and zero the other. For exact Hadamard ordering use
survey::as.svrepdesign(..., type = "BRR").
Usage
morie_weights_brr(weights, strata, n_replicates = NULL, seed = 42)
Arguments
- weights
Numeric vector of unit-level design weights.
- strata
Optional vector of stratum identifiers aligned with
weights.
- n_replicates
Integer; replicate count.
- seed
Integer RNG seed.
Examples
set.seed(1)
res <- morie_weights_brr(runif(16, 0.5, 2), rep(1:8, each = 2))
str(res, max.level = 1)
#> num [1:16, 1:8] 0 2.12 0 3.72 0 ...
