The generalized exponential model for sampling weight calibration for extreme values, nonresponse, and poststratification
The generalized exponential model for sampling weight calibration for extreme values, nonresponse, and poststratification was written by R.E. Folsom and A.C. Singh in 2000. It was part of the proceedings of the American Statistical Association Section on Survey Research Methods.
The authors build on Deville and Särndal (1992) (notated DS in this article). Establishing some notation:
individuals are indexed by k
the sample of individuals is s
yk is the study variable
xk is a vector of p predictor variables
Tx is a vector of p (i.e., same dimension as xk) post-stratification controls
dk is a design weight
wk is an adjusted weight, i.e. dkak
λ be a vector of p (i.e., same dimension as xk) parameters
The logit-type SD model for the aforementioned weight adjustments is:
where l and u are set as lower/upper bounds, for which the only requirement is that l < 1 < u since 1 is implicitly the 'center'; and where A = (u - l)/(u - 1)(1 - l). The λ parameters are then estimated iteratively through satisfying the equation Σs xkdkak(λ) - Tx = 0 and minimizing Δ(wk,dk):
The authors extend this model for centers other than 1, variable bounds, and post-stratification controls designed to address problems other than sampling error. They call this the generalized exponential model (GEM).
bounds and center are set such that lk < ck < uk
it follows that Ak = (uk - lk)/(uk - ck)(ck - lk)
Σs xkdkak(λ) - T̃k = 0 (note the tilde over T)
The λ parameters can be estimated through Newton-Raphson iterative steps.
Let v denote the iteration.
Set initial values as φk(0) = 1 and λ(0) = 0.
Within each iteration, let Γφ,v = diag(dkφk(v)) for convenience.
Update the parameters as λ(v) = λ(v-1) + (X'Γφ,v-1X)-1 (Tx - T̂x(v-1)).
For subsequent iterations, φk(v) = (uk - ak(v))(ak(v) - lk)/(uk - ck)(ck - lk).
Iterate until ||Tx - T̂x(v)|| increases, instead of decreasing.
Reading Notes
Section 4 gets into asymptotic properties, will need to come back to this too.
