data.root

Used by the DP recurrences to evaluate the stratification objective at
a specific remaining distance 'd' and first-stratum width 'y', given
stratum cost 'c' and constants tucked in 'my_env.'

Determines Optimum Strata Boundaries (OSB) and Optimum Sample
Sizes (OSS) for univariate stratified sampling designs under Neyman
allocation. The stratification variable is described by a best-fitting
parametric distribution, selected automatically by AIC from a set of
continuous families (normal, log-normal, gamma, Weibull, exponential,
Cauchy, uniform, Pareto, triangular and right-triangular), and the
optimum boundaries are obtained by minimising the Neyman objective.
Version 2.0 keeps the original globally optimal Dynamic Programming (DP)
solver of Reddy and Khan (2020) as the default and adds two faster
derivative-free alternatives for interactive and large-scale use: a
multi-start 'COBYLA' solver and a two-phase 'global' solver that couples
'DIRECT-L' with 'COBYLA' refinement. It also provides cost-constrained
allocation with unequal per-stratum costs, a design-efficiency comparison
(compare_designs), two- and three-dimensional and interactive
visualisations, solution-quality diagnostics (a Cauchy-Schwarz optimality
gap and KKT first-order residuals for the derivative-free solvers) and a
self-contained 'shiny' application, while remaining backward compatible
with the strata.data() and strata.distr() interface of version 1.x. The
methodology follows Khan et al. (2008)
<https://www150.statcan.gc.ca/n1/pub/12-001-x/2008002/article/10761-eng.pdf>,
Reddy and Khan (2018) <doi:10.1111/anzs.12244> and Reddy and Khan (2020)
<doi:10.1111/anzs.12301>.

Karuna G. Reddy

stratifyR

Optimal Stratification of Univariate Populations

M. G. M. Khan

data.root function

<dl><dt>d</dt>
<dd>numeric: remaining distance/range on the (scaled) axis</dd>
<dt>y</dt>
<dd>numeric: first stratum width on the (scaled) axis</dd>
<dt>c</dt>
<dd>numeric: per stratum cost multiplier (Ch[k])</dd>
<dt>my_env</dt>
<dd>environment: holds distribution name, parameters and scaling</dd></dl>

Arguments

Calculate the objective function value for a given (d, y) — data.root

<dl>

<dt>d</dt>
<dd>numeric: remaining distance/range on the (scaled) axis</dd>


<dt>y</dt>
<dd>numeric: first stratum width on the (scaled) axis</dd>


<dt>c</dt>
<dd>numeric: per stratum cost multiplier (Ch[k])</dd>


<dt>my_env</dt>
<dd>environment: holds distribution name, parameters and scaling</dd>

</dl>

data.root: Calculate the objective function value for a given (d, y)

Description

Usage

Value

Arguments