In survey research, the design effect is a number that shows how well a sample of people may represent a larger group of people for a specific measure of interest (such as the mean). This is important when the sample comes from a sampling method that is different than just picking people using a simple random sample.
The design effect is a positive real number, represented by the symbol . If , then the sample was selected in a way that is just as good as if people were picked randomly. When , then inference from the data collected is not as accurate as it could have been if people were picked randomly.
When researchers use complicated methods to pick their sample, they use the design effect to check and adjust their results. It may also be used when planning a study in order to determine the sample size.
In survey methodology, the design effect (generally denoted as , , or ) is a measure of the expected impact of a sampling design on the variance of an estimator for some parameter of a population. It is calculated as the ratio of the variance of an estimator based on a sample from an (often) complex sampling design, to the variance of an alternative estimator based on a simple random sample (SRS) of the same number of elements. The (be it estimated, or known a priori) can be used to evaluate the variance of an estimator in cases where the sample is not drawn using simple random sampling. It may also be useful in sample size calculations and for quantifying the representativeness of samples collected with various sampling designs.
The design effect is a positive real number that indicates an inflation (), or deflation () in the variance of an estimator for some parameter, that is due to the study not using SRS (with , when the variances are identical). Intuitively we can get when we have some a-priori knowledge we can exploit during the sampling process (which is somewhat rare). And, in contrast, we often get when we need to compensate for some limitation in our ability to collect data (which is more common). Some sampling designs that could introduce generally greater than 1 include: cluster sampling (such as when there is correlation between observations), stratified sampling (with disproportionate allocation to the strata sizes), cluster randomized controlled trial, disproportional (unequal probability) sample (e.g. Poisson sampling), statistical adjustments of the data for non-coverage or non-response, and many others. Stratified sampling can yield that is smaller than 1 when using Proportionate allocation to strata sizes (when these are known a-priori, and correlated to the outcome of interest) or Optimum allocation (when the variance differs between strata and is known a-priori).
Many have been proposed in the literature for how a known sampling design influences the variance of estimators of interest, either increasing or decreasing it. Generally, the design effect varies among different statistics of interests, such as the total or . It also matters if the sampling design is correlated with the outcome of interest. For example, a possible sampling design might be such that each element in the sample may have a different probability to be selected. In such cases, the level of correlation between the probability of selection for an element and its measured outcome can have a direct influence on the subsequent design effect. Lastly, the design effect can be influenced by the distribution of the outcome itself. All of these factors should be considered when estimating and using design effect in practice.
The term "design effect" was coined by in his 1965 book "Survey Sampling." In it, Kish proposed the general definition for the design effect, as well as formulas for the (with intraclass correlation); and the famous . These are often known as "Kish's design effect", and were later combined into a single formula.
In a 1995 paper, Kish mentions that a similar concept, termed "Lexis ratio", was described at the end of the 19th century. The closely related was described by in 1950, while computations of ratios of variances were already published by Kish and others from the late 1940s to the 1950s. One of the precursors to Kish's definition was work done by Cornfield in 1951.
In his 1995 paper, Kish proposed that considering the design effect is necessary when averaging the same measured quantity from multiple surveys conducted over a period of time. He also suggested that the design effect should be considered when extrapolating from the error of simple statistics (e.g. the mean) to more complex ones (e.g. regression coefficients). However, when analyzing data (e.g., using survey data to fit models), values are less useful nowadays due to the availability of specialized software for analyzing survey data. Prior to the development of software that computes standard errors for many types of designs and estimates, analysts would adjust standard errors produced by software that assumed all records in a dataset were i.i.d by multiplying them by a (see Deft definition below).
The design effect, commonly denoted by (or , sometimes with additional subscripts), is the ratio of two theoretical variances for s of some ():
So that:
In other words, measures the extent to which the variance has increased (or, in some cases, decreased) because the sample was drawn and adjusted to a specific sampling design (e.g., using weights or other measures) compared to if the sample was from a (without replacement). Notice how the definition of is based on parameters of the population that are often unknown, and that are hard to estimate directly. Specifically, the definition involves the variances of estimators under two different sampling designs, even though only a single sampling design is used in practice.
For example, when estimating the population mean, the (for some sampling design p) is:
Where is the sample size, is the fraction of the sample from the population, is the (squared) (FPC), is the , and is some estimator of the variance of the mean under the sampling design. The issue with the above formula is that it is extremely rare to be able to directly estimate the variance of the estimated mean under two different sampling designs, since most studies rely on only a single sampling design.
There are many ways of calculation , depending on the parameter of interest (e.g. population total, population mean, quantiles, ratio of quantities etc.), the estimator used, and the sampling design (e.g. clustered sampling, stratified sampling, post-stratification, multi-stage sampling, etc.). The process of estimating for specific designs will be described in .
A related quantity to , proposed by Kish in 1995, is the Design Effect Factor, abbreviated as (or also ). It is defined as the square root of the variance ratios while also having the denominator use a simple random sample with replacement (SRSWR), instead of without replacement (SRSWOR):
In this later definition (proposed in 1995, vs 1965) Kish argued in favor of using over for several reasons. It was argued that SRS "without replacement" (with its positive effect on the variance) should be captured in the denominator part in the definition of the design effect, since it is part of the sampling design. Also, since often the use of the factor is in s), it was claimed that using will be simpler than writing . It is also said that for many cases when the population is very large, is (almost) the square root of (), hence it is easier to use than exactly calculating the (FPC).
Even so, in various cases a researcher might approximate the by calculating the variance in the numerator while assuming SRS with replacement (SRSWR) instead of SRS without replacement (SRSWOR), even if it is not precise. For example, consider a multistage design with primary sampling units (PSUs) selected systematically with probability proportional to some measure of size from a list sorted in a particular way (say, by number of households in each PSU). Also, let it be combined with an estimator that uses to match the totals for several demographic variables. In such a design, the joint selection probabilities for the PSUs, which are needed for a without replacement variance estimator, are 0 for some pairs of PSUs - implying that an exact design-based (i.e., repeated sampling) variance estimator does not exist. Another example is when a public use file issued by some government agency is used for analysis. In such a case the information on joint selection probabilities of first-stage units is almost never released. As a result, an analyst cannot estimate a with replacement variance for the numerator even if desired. The standard workaround is to compute a variance estimator as if the PSUs were selected with replacement. This is the default choice in software packages such as Stata, the R survey package, and the SAS survey procedures.
The effective sample size, defined by Kish in 1965, is calculated by dividing the original sample size by the design effect. Namely:
This quantity reflects what would be the sample size that is needed to achieve the current variance of the estimator (for some parameter) with the existing design, if the sample design (and its relevant parameter estimator) were based on a .
A related quantity is the effective sample size ratio (ESSR), which can be calculated by simply taking the inverse of (i.e., ).
For example, let the design effect, for estimating the population mean based on some sampling design, be 2. If the sample size is 1,000, then the effective sample size will be 500. It means that the variance of the based on 1,000 samples will be the same as that of a simple based on 500 samples obtained using a simple random sample.
Different sampling designs and statistical adjustments may have substantially different impact on the bias and variance of estimators (such as the mean).
An example of a design which can lead to estimation efficiency, compared to simple random sampling, is . This efficiency is gained by leveraging information about the composition of the population. For example, if it is known that gender is correlated with the outcome of interest, and also that the male-female ratio for some population is (say) 50%-50%, then sampling exactly half of the sample from each gender will reduce the variance of the outcome's estimator. Similarly, if a particular sub-population is of special interest, deliberately over-sampling from that sub-population will decrease the variance for estimations made about it.
Improvement in variance efficiency might sometimes be sacrificed for convenience or cost. For example, in the case the units may have equal or unequal selection probabilities, irrespective of their (and their negative effect of increasing the variance of the estimators). We might decide (for practical reasons) to collect responses from only 2 people of each household (i.e., a sampled cluster), which could lead to more complex post-sampling adjustment to deal with unequal selection probabilities. Also, such decisions could lead to less efficient estimators than just taking a fixed proportion of responses from a cluster.
When the sampling design isnâÂÂt set in advance and needs to be figured out from the data we have, this can lead to an increase of both the variance and bias of the weighted estimator. This might happen when making adjustments for issues like non-coverage, non-response, or an unexpected strata split of the population that wasnâÂÂt available during the initial sampling stage. In these cases, we might use statistical procedures such as post-stratification, raking, or inverse propensity score weighting (where the propensity scores are estimated), among other methods. Using these methods requires assumptions about the initial design model. For example, when we use post-stratification based on age and gender, it is assumed that these variables can explain a significant portion of the bias in the sample. The quality of these estimators is closely tied to the quality of the additional information and the assumptions used when making them. Either way, even when estimators (like propensity score models) do a good job capturing most of the sampling design, using the weights can make a small or a large difference, depending on the specific data-set.
Due to the large variety in sampling designs (with or without an effect on unequal selection probabilities), different formulas have been developed to capture the potential design effect, as well as to estimate the variance of estimators when accounting for the sampling designs. Sometimes, these different design effects can be compounded together (as in the case of unequal selection probability and cluster sampling, more details in the following sections). Whether or not to use these formulas, or just assume SRS, depends on the expected amount of bias reduction vs. the increase in estimator variance (and in the overhead of methodological and technical complexity).
There are various ways to sample units so that each unit would have the exact same probability of selection. Such methods are called (EPSEM) methods. Some of the more basic methods include (SRS, with or without replacement) and for getting a fixed sample size. There is also with a random sample size. More advanced techniques such as and can also be designed to be EPSEM. For example, in cluster sampling we can use a two stage sampling in which we sample each cluster (which may be of different sizes) with equal probability, and then sample from each cluster at the second stage using SRS with a fixed proportion (e.g. sample half of the cluster, the whole cluster, etc.). This method will yield EPSEM, but the specific number of elements we end up with is stochastic (i.e., non deterministic). Another strategy for cluster sampling that leads to EPSEM is to sample clusters in a way that is proportional to their sizes, and then sample a fixed number of elements inside each cluster.
In their works, and others highlight several known reasons that lead to unequal selection probabilities:
Adjusting for unequal probability selection through "individual case weights" (e.g. inverse probability weighting), yields various types of estimators for quantities of interest. Estimators such as yield unbiased estimators (if the selection probabilities are indeed known, or approximately known), for total and the mean of the population. Deville and Särndal (1992) coined the term "calibration estimator" for estimators using weights such that they satisfy some condition, such as having the sum of weights equal the population size. And more generally, that the weighted sum of weights is equal some quantity of an auxiliary variable: (e.g., that the sum of weighted ages of the respondents is equal to the population size in each age group).
The two primary ways to argue about the properties of calibration estimators are:
As we will see later, some proofs in the literature rely on the randomization-based framework, while others focus on the model-based perspective. When moving from the mean to the , more complexity is added. For example, in the context of , often the population size itself is considered an unknown quantity that is estimated. So in the calculation of the weighted mean is in fact based on a , with an estimator of the total at the numerator and an estimator of the population size in the denominator (making the variance calculation to be more complex).
There are many types (and subtypes) of weights, with different ways to use and interpret them. With some weights their absolute value has some important meaning, while with other weights the important part is the relative values of the weights to each other. This section introduces some of the more common types of weights so that they can be referenced in follow-up sections.
There are also indirect ways of applying "weighted" adjustments. For example, the existing cases may be duplicated to missing observations (e.g. from non-response), with variance estimated using methods such as . An alternative approach is to remove (assign a weight of 0 to) some cases. For example, when wanting to reduce the influence of over-sampled groups that are less essential for some analysis. Both cases are similar in nature to inverse probability weighting but the application in practice gives more/less rows of data (making the input potentially simpler to use in some software implementation), instead of applying an extra column of weights. Nevertheless, the consequences of such implementations are similar to just using weights. So while in the case of removing observations the data can easily be handled by common software implementations, the case of adding rows requires special adjustments for the uncertainty estimations. Not doing so may lead to erroneous conclusions(i.e., there is when using alternative representation of the underlying issues).
The term "Haphazard weights", coined by Kish, is used to refer to weights that correspond to , but ones that are not related to the expectancy or variance of the selected elements.
When taking an unrestricted sample of elements, we can then randomly split these elements into strata, each of them containing some size of elements so that . All elements in each stratum has some (known) non-negative weight assigned to them (). The weight can be produced by the inverse of some for elements in each stratum (i.e., following a procedure such as post-stratification). In this setting, Kish's design effect, for the increase in variance of the sample due to this design (reflected in the weights), versus of some outcome variable y (when there is no correlation between the weights and the outcome, i.e. haphazard weights) is:
By treating each item as coming from its own stratum , Kish (in 1992) simplified the above formula to the (well-known) following version:
This version of the formula is valid when one stratum had several observations taken from it (i.e., each having the same weight), or when there are just many strata were each one had one observation taken from it, but several of them had the same probability of selection. While the interpretation is slightly different, the calculation of the two scenarios comes out to be the same.
When using Kish's design effect for unequal weights, you may use the following simplified formula for "'s Effective Sample Size"
The above formula, by , gives the increase in the variance of the based on . This can also be written as the following formula where y are observations selected using (with no within-cluster correlation, and no relationship to the expectancy or variance of the outcome measurement), and y' are the observations we would have had if we got them from a :
It can be shown that the ratio of variances formula can be reduced to Kish's formula by using a . In it, Kish's formula will hold when all n observations () are (at least approximately) (), with the same () in the response variable of interest (y). It will also be required to assume the weights themselves are not a but rather some known constants (e.g. the inverse of probability of selection, for some pre-determined and known ).
The following is a simplified proof for when there are no clusters (i.e., no between element of the sample) and each stratum includes only one observation:
Transitions:
The conditions on y are trivially held if the y observations are with the same and . In such cases, , and we can estimate by using . If the y's are not all with the same expectations then we cannot use the estimated variance for calculation, since that estimation assumes that all s have the same expectation. Specifically, if there is a correlation between the weights and the outcome variable y, then it means that the expectation of y is not the same for all observations (but rather, dependent on the specific weight value for each observation). In such a case, while the design effect formula might still be correct (if the other conditions are met), it would require a different estimator for the variance of the weighted mean. For example, it might be better to use a .
If different s values have different variances, then while the weighted variance could capture the correct population-level variance, Kish's formula for the design effect may no longer be true.
A similar issue happens if there is some correlation structure in the samples (such as when using ).
Notice that Kish's definition of the design effect is closely tied to the (Kish also calls it relvariance or relvar for short) of the weights (when using the for ). This has several notations in the literature:
Where is the population variance of , and is the mean. When the weights are normalized to sample size (so that their sum is equal to n and their mean is equal to 1), then and the formula reduces to . While it is true we assume the weights are fixed, we can think of their variance as the variance of an defined by sampling (with equal probability) one weight from our set of weights (similar to how we would think about the correlation of x and y in a ).
Kish's original definition compared the variance under some sampling design to the variance achieved through a . Some literature provide the following alternative definition for Kish's design effect: "the ratio of the variance of the weighted survey mean under disproportionate stratified sampling to the variance under when all stratum unit variances are equal". Reflecting on this, Park and Lee (2006) stated that "The rationale behind [...][Kish's] derivation is that the loss in precision of [the weighted mean] due to haphazard unequal weighting can be approximated by the ratio of the variance under disproportionate stratified sampling to that under the proportionate stratified sampling".
Note that this alternative definition only approximated since if the denominator is based on "proportionate stratified sampling" (achieved via ) then such a selection will yield a reduced variance as compared with . This is since stratified sampling removes some of the variability in the specific number of elements per stratum, as occurs under SRS.
Relatedly, Cochran (1977) provides a formula for the proportional increase in variance due to deviation from optimum allocation (what, in Kish's formulas, would be called L).
Early papers used the term . As more definitions of the design effect appeared, was denoted (or ) or simply for short. Kish's design effect is also known as the "Unequal Weighting Effect" (or just UWE), termed by Liu et al. in 2002.
The estimator for the total is the "p-expanded with replacement" estimator (a.k.a.: pwr-estimator or ). It is based on a (with replacement, denoted SIR) of n items () from a population of size N. Each item has a probability of (k from 1 to N) to be drawn in a single draw (, i.e. it is a ). The probability that a specific will appear in the sample is . The "p-expanded with replacement" value is with the following expectancy: . Hence , the pwr-estimator, is an unbiased estimator for the sum total of y.
In 2000, Bruce D. Spencer proposed a formula for estimating the design effect for the variance of estimating the total (not the mean) of some quantity (), when there is correlation between the selection probabilities of the elements and the outcome variable of interest.
In this setup, a sample of size n is drawn (with replacement) from a population of size N. Each item is drawn with probability (where , i.e. ). The selection probabilities are used to define the : . Notice that for some random set of n items, the sum of weights will be equal to 1 only by expectation () with some variability of the sum around it (i.e., the sum of elements from a ). The relationship between and is defined by the following (population) :
Where is the outcome of element i, which linearly depends on with the intercept and slope . The residual from the fitted line is . We can also define the population variances of the outcome and the residuals as and . The correlation between and is .
Spencer's (approximate) design effect for estimating the total of y is:
Where:
This assumes that the regression model fits well so that the probability of selection and the residuals are , since it leads to the residuals, and the square residuals, to be uncorrelated with the weights, i.e., that and also .
When the population size (N) is very large, the formula can be written as:
(since , where )
This approximation assumes that the linear relationship between P and y holds. And also that the correlation of the weights with the errors, and the errors squared, are both zero. I.e., and .
We notice that if , then (i.e., the average of y). In such a case, the formula reduces to
Only if the variance of y is much larger than its mean, then the right-most term is close to 0 (i.e., ), which reduces Spencer's design effect (for the estimated total) to be equal to Kish's design effect (for the ratio means): . Otherwise, the two formulas will yield different results, which demonstrates the difference between the design effect of the total vs. the design effect of the mean.
In 2001, Park and Lee extended Spencer's formula to the case of the ratio-mean (i.e., estimating the mean by dividing the estimator of the total with the estimator of the population size). It is:
Where:
Park and Lee's formula is exactly equal to Kish's formula when . Both formulas relate to the design effect of the mean of y, while Spencer's relates to the estimation of the population total.
In general, the for the total () tends to be less efficient than the for the ratio mean () when is small. And in general, impacts the efficiency of both design effects.
For data collected using we assume the following structure:
When clusters are all of the same size , the design effect D<sub>eff</sub>, proposed by Kish in 1965 (and later re-visited by others), is given by:
It is sometimes also denoted as .
In various papers, when cluster sizes are not equal, the above formula is also used with as the average cluster size (which is also sometimes denoted as ). In such cases, Kish's formula (using the average cluster weight) serves as a conservative (upper bound) of the exact design effect.
Alternative formulas exists for unequal cluster sizes. Followup work had discussed the sensitivity of using the average cluster size with various assumptions.
In a 1987 paper, Kish proposed a combined design effect that incorporates both the effects due to weighting that accounts for unequal selection probabilities and cluster sampling:
The above uses notations similar to what is used in this article (the original 1987 publication used different notation). A justification for this formula was provided by Gabler et al.
In 2000, Liu and Aragon proposed a decomposition of unequal selection probabilities design effect for different strata in stratified sampling. In 2002, Liu et al. extended that work to account for stratified samples, where within each stratum is a set of unequal selection probability weights. The cluster sampling is either global or per stratum. Similar work was done also by Park et al. in 2003.
The Chen-Rust extends the model-based justification of KishâÂÂs 1987 formula for design effects proposed by Gabler, el. al., applying it to two-stage designs with stratification at the first stage and to three-stage designs without stratification. The modified formulae define the overall design effect using survey weights and population intracluster correlations. These formulae allow for insightful interpretations of design effects from various sources and can estimate intracluster correlations in completed surveys or predict design effects in future surveys.
Henry's proposes an extended model-assisted weighting design-effect measure for single-stage sampling and calibration weight adjustments for a case where , where is a vector of covariates, the model errors are independent, and the estimator of the population total is the general regression estimator (GREG) of Särndal, Swensson, and Wretman (1992). The new measure considers the combined effects of non-epsem sampling design, unequal weights from calibration adjustments, and the correlation between an analysis variable and the auxiliaries used in calibration.
Lohr's is for ordinary least squares (OLS) and generalized least squares (GLS) estimators in the context of cluster sampling, using a random coefficient regression model. Lohr presents conditions under which the GLS estimator of the regression slope has a design effect less than 1, indicating higher efficiency. However, the design effect of the GLS estimator is highly sensitive to model specification. If an underlying random coefficient model is incorrectly specified as a random intercept model, the design effect can be seriously understated. In contrast, the OLS estimator of the regression slope and the design effect calculated from a design-based perspective are robust to misspecification of the variance structure, making them more reliable in situations where the model specification may not be accurate.
may be used when planning a future data collection, as well as a diagnostic tool:
Considering the design effect is unnecessary when the source population is closely , or when the sample design of the data was drawn as a . It is also less useful when the sample size is relatively small (at least partially, for practical reasons).
While Kish originally hoped to have the design effect be as agnostic as possible to the underlying distribution of the data, sampling probabilities, their correlations, and the statistics of interest, followup research has shown that these do influence the design effect. Hence, these properties should be carefully considered when deciding which calculation to use, and how to use it.
The design effect is rarely applied when constructing confidence intervals. Ideally, one would be able to determine, for an estimator of a particular parameter, both the variance under Simple Random Sample (SRS) with replacement and the design effect (which accounts for all elements of the sampling design that change the variance). In such scenarios, the basic variance and the design effect could have been multiplied to compute the variance of the estimator for the specific design. This computed value can then be employed to form confidence intervals. However, in real-world applications, it is uncommon to estimate both values simultaneously. As a result, other methods are favored. For instance, Taylor linearization is utilized to construct confidence intervals based on the . More broadly, the bootstrap method, also known as , is applied for a range of weighted statistics.
Kish's design effect is implemented in various statistical software packages: