SamplingStrata (version 1.5-4)

bethel: Multivariate optimal allocation

Description

Multivariate optimal allocation for different domains of interest in stratified sample design under a given stratification of the sampling frame

Usage

bethel  (
          stratif,
          errors, 
          minnumstrat=2,          
          maxiter=200,
          maxiter1=25,
          printa=FALSE,
		  realAllocation=FALSE,
          epsilon=1e-11
          )

Value

A vector containing the computed optimal allocation

Arguments

errors

Data frame of coefficients of variation for each domain

stratif

Data frame of survey strata

minnumstrat

Minimum number of units per strata (default=2)

maxiter

Maximum number of iterations of the algorithm (default=200)

maxiter1

Maximum number of iterations (default=25) of the general procedure. This kind of iteration may be required by the fact that when in a stratum the number of allocated units is greater or equal to its population, that stratum is set as "census stratum", and the whole procedure is re-initialised

printa

If TRUE then two attributes are added to the resulting vector. The first ('confr') is a comparison between results obtained with 3 different allocation methods: Bethel, proportional and equal. The second ('outcv') is a table reporting planned and actual CV, together with a sensitivity analysis

realAllocation

If FALSE, the allocation is based on INTEGER values; if TRUE, the allocation is based on REAL values

epsilon

Epsilon (default=1e-11)): this value is used to compare the difference in results from one iteration to the other; if it it is lower than "epsilon", then the procedure stops

Author

Daniela Pagliuca with contributions from Teresa Buglielli and Giulio Barcaroli

Examples

if (FALSE) {
library(SamplingStrata)
data(strata)
data(errors)
n <- bethel(strata, errors, printa=TRUE)
sum(n)
attributes(n)$confr
attributes(n)$outcv
}