Introduction to Statistical Modeling with SAS/STAT Software


Inference Principles for Survey Data

Design-based and model-assisted statistical inference for survey data requires that the randomness due to the selection mechanism be taken into account. This can require special estimation principles and techniques.

The SURVEYMEANS, SURVEYFREQ, SURVEYREG, SURVEYLOGISTIC, and SURVEYPHREG procedures support design-based and/or model-assisted inference for sample surveys. Suppose $\pi _ i$ is the selection probability for unit i in sample S. The inverse of the inclusion probability is known as sampling weight and is denoted by $w_ i$. Briefly, the idea is to apply a relationship that exists in the population to the sample and to take into account the sampling weights. For example, to estimate the finite population total $T_ N = \sum _{i \in U_ N y_ i}$ based on the sample S, you can accumulate the sampled values while properly weighting: $\widehat{T}_{\pi } = \sum _{i \in S} w_ i y_ i$. It is easy to verify that $\widehat{T}_{\pi }$ is design-unbiased in the sense that $\mr{E}[\widehat{T}_{\pi } | {\mathcal F}_ N] = T_ N$ (see Cochran 1977).

When a statistical model is present, similar ideas apply. For example, if $\beta _{N0}$ and $\beta _{N1}$ are finite population quantities for a simple linear regression working model that minimize the sum of squares

\[ \sum _{i \in U_ N} \left(y_ i - \beta _{0N} - \beta _{1N}x_ i\right)^2 \]

in the population, then the sample-based estimators $\widehat{\beta }_{0S}$ and $\widehat{\beta }_{1S}$ are obtained by minimizing the weighted sum of squares

\[ \sum _{i \in S} w_ i (y_ i - \hat\beta _{0S} - \hat\beta _{1S}x_ i)^2 \]

in the sample, taking into account the inclusion probabilities.

In model-assisted inference, weighted least squares or pseudo-maximum likelihood estimators are commonly used to solve such estimation problems. Maximum pseudo-likelihood or weighted maximum likelihood estimators for survey data maximize a sample-based estimator of the population likelihood. Assume a working model with uncorrelated responses such that the finite population log-likelihood is

\[ \sum _{i \in U_ N} l(\theta _{1N}, \ldots , \theta _{pN}; y_ i), \]

where $\theta _{1N}, \ldots , \theta _{pN}$ are finite population quantities. For independent sampling, one possible sample-based estimator of the population log likelihood is

\[ \sum _{i \in S} w_ i l(\theta _{1N}, \ldots , \theta _{pN}; y_ i) \]

Sample-based estimators $\widehat{\theta }_{1S}, \ldots , \widehat{\theta }_{pS}$ are obtained by maximizing this expression.

Design-based and model-based statistical analysis might employ the same statistical model (for example, a linear regression) and the same estimation principle (for example, weighted least squares), and arrive at the same estimates. The design-based estimation of the precision of the estimators differs from the model-based estimation, however. For complex surveys, design-based variance estimates are in general different from their model-based counterpart. The SAS/STAT procedures for survey data (SURVEYMEANS, SURVEYFREQ, SURVEYREG, SURVEYLOGISTIC, and SURVEYPHREG procedures) compute design-based variance estimates for complex survey data. See the section Variance Estimation, in ChapterĀ 14: Introduction to Survey Sampling and Analysis Procedures, for details about design-based variance estimation.