compare_designs: Compare Survey Design Efficiencies

In stratifyR: Optimal Stratification of Univariate Populations

Compare Survey Design Efficiencies

Description

Given a fitted "strata" object (from strata.data or strata.distr), computes and compares the variance of the sample mean under three designs for the same total sample size n:

SRS

Simple random sampling without replacement (baseline).

Proportional

Stratified sampling with n_h \propto W_h (proportional allocation).

Neyman

Stratified sampling with n_h \propto W_h S_h (optimal/Neyman allocation), the allocation used by stratifyR.

Design effects (DEFF) follow Kish (1965): \mathrm{DEFF} = V_{\rm design} / V_{\rm SRS}. The equivalent SRS sample size is the number of observations an SRS design would need to match the precision of the Neyman design, and the cost saving is the corresponding percentage reduction.

Usage

compare_designs(object, ...)

## S3 method for class 'strata'
compare_designs(object, n = NULL, ...)

## S3 method for class 'compare_designs'
print(x, digits = 6, ...)

Arguments

object	An object of class "strata".
...	Currently unused.
n	Integer. Total sample size. Defaults to object$nhTot (the value used when the stratification was run).
x	A "compare_designs" object.
digits	Integer. Number of significant digits used when printing.

Value

An object of class "compare_designs" (an invisibly-printed list) with components:

n

Total sample size used.

H

Number of strata.

S2

Estimated population variance S^2.

V_within

Within-stratum variance component \sum W_h S_h^2.

V_srs, V_prop, V_opt

Variance of \bar{y} under SRS, proportional, and Neyman designs.

SE_srs, SE_prop, SE_opt

Corresponding standard errors.

deff_prop, deff_opt

Design effects relative to SRS.

n_srs_equiv

Equivalent SRS sample size for same precision as the Neyman design.

pct_saving

Percentage sample-size saving of Neyman over SRS.

WhShTot

\sum W_h S_h (objective function value).

See Also

strata.data, strata.distr

Examples

## Not run: 
res <- strata.data(data = anaemia$Iron, h = 3, n = 300)
cd  <- compare_designs(res)
cd

## End(Not run)

stratifyR documentation built on Sept. 10, 2026, 5:07 p.m.