Exposure measurement error in time-series studies of air pollution: concepts and consequences.

Misclassification of exposure is a well-recognized inherent limitation of epidemiologic studies of disease and the environment. For many agents of interest, exposures take place over time and in multiple locations; accurately estimating the relevant exposures for an individual participant in epidemiologic studies is often daunting, particularly within the limits set by feasibility, participant burden, and cost. Researchers have taken steps to deal with the consequences of measurement error by limiting the degree of error through a study's design, estimating the degree of error using a nested validation study, and by adjusting for measurement error in statistical analyses. In this paper, we address measurement error in observational studies of air pollution and health. Because measurement error may have substantial implications for interpreting epidemiologic studies on air pollution, particularly the time-series analyses, we developed a systematic conceptual formulation of the problem of measurement error in epidemiologic studies of air pollution and then considered the consequences within this formulation. When possible, we used available relevant data to make simple estimates of measurement error effects. This paper provides an overview of measurement errors in linear regression, distinguishing two extremes of a continuum-Berkson from classical type errors, and the univariate from the multivariate predictor case. We then propose one conceptual framework for the evaluation of measurement errors in the log-linear regression used for time-series studies of particulate air pollution and mortality and identify three main components of error. We present new simple analyses of data on exposures of particulate matter < 10 microm in aerodynamic diameter from the Particle Total Exposure Assessment Methodology Study. Finally, we summarize open questions regarding measurement error and suggest the kind of additional data necessary to address them. ImagesFigure 1Figure 2Figure 3

Misclassification of exposure has long been recognized as an inherent limitation of epidemiologic studies of the environment and disease (1). For many agents of interest, exposures take place over time and in multiple locations so that it is difficult to accurately estimate the relevant exposures for individual study participants, particularly within the limits set by feasibility, participant burden, and cost. In general, exposure measurement error tends to blunt the sensitivity of epidemiologic studies for detecting the effects of environmental agents, although the specific impact of exposure error on effect estimates depends on several factors induding the study design, the types of error, and the relationships between the outcome and the independent variables (1,2). As the problem of exposure error has become well recognized, researchers have taken steps to control its consequences by limiting the degree of error through careful study design and data collection, by estimating the degree of error using a nested validation study, and by making adjustments for measurement error in statistical analyses.
In this paper, we address the problem of exposure error in observational ecologic time-series studies of air pollution and health.
The pollution of outdoor air is a public health concern throughout the world. For decades, epidemiologic studies have been a cornerstone of our approach to investigating the health effects of air pollution and have been a principal basis for setting regulations to protect the public against adverse health effects. Two broad types of observational study designs have been used in research on air pollution: ecologic or aggregate-level studies, either crosssectional or time-series in design, and individual-level studies, primarily of the cross-sectional or cohort designs. In ecologic studies, population-level indicators of exposure are typically drawn from centrally sited air pollution monitors. In individual-level cross-sectional and cohort studies, exposure estimates for individual participants may be based on centrally located monitors, on the combination of central monitors with personal records of environments where participants spend time, or on personal exposure monitoring (3).
Regardless of study design, any pollution exposure assessment strategy introduces some degree of exposure measurement error. For example, in the Six Cities Study (4,5), a prospective cohort study of air pollution and respiratory health and mortality, exposure estimates for persons from each of the six cities were based on centrally sited monitors.
Exposures were further characterized for samples of participants using personal monitors and monitors placed in their homes; the resulting data provide an understanding of the components of error associated with using the central site data for all participants.
The problem of measurement errors in predictor variables in regression analysis has been carefully studied in the statistics and epidemiologic literature for several decades. Fuller (6) summarized early research on linear regression with so-called "errors-in-x" variables. Carroll et al. (2) extended this literature to generalized linear models including Poisson, logistic, and survival regression analyses. Thomas et al. (2) presented an overview of the exposure error or misdassification problem from the general epidemiologic perspective. Spiegelman et al. (8), Willett (95, and Pierce et al. (10) provided recent illustrations of statistical approaches to measurement error in epidemiologic research.
In one of the early papers on the topic of exposure error in studies of air pollution, Shy et al. (11) described the problem and addressed its consequences in an epidemiologic framework. Goldstein and Landovitz (12,13) recognized that a single monitoring station may not adequately represent a geographic area and conducted an analysis of correlations among concentration data from several monitors in New York City. In the ensuing decades, there has been deepening understanding of measurement error in general and of its potential implications for the study of air pollution (14,15). Articles * Zeger et al.
During the 1 990s, substantial new evidence, largely from ecologic time-series analyses of air pollution and mortality, showed that daily variation in ambient measures of particulate air pollution within the current standards of the U.S. Environmental Protection Agency was associated with daily mortality levels (16). There are strong concerns about interpreting these associations in view of potential errors in the exposure measurements. In a series of papers, Lipfert and Wyzga (17) and Lipfert (18,19) suggested that the central monitoring data used in the time-series analyses have uncertain relationships with the exposures of individuals in the study communities; they further argued that those errors vary among pollutants, complicating interpretation of any multipollutant models. Lipfert and Wyzga (17) referred specifically to an analysis by Schwartz et al. (20) that attributed effects on mortality to fine rather than coarse partides, based in part on the results of multivariable models which included variables for both particulate measures.
A number of exposure assessment studies found sizable differences between actual personal exposures to particles and estimates based on central monitor values (21). Some investigators have questioned whether the observed associations are plausible given these findings. However, Schwartz et al. (20) responded that as the number of deaths per day is calculated over the population, the relevant exposure measure is the mean of personal exposures on that day, which is probably more tightly correlated with central station monitoring than individual exposures. Janssen et al. (22) reported that much of the variation in particulate matter < 10 ,um in aerodynamic diameter (PMl) measurements is between people and that the longitudinal correlation between average and ambient PMIO measures is relatively high. The debate over measurement error and its consequences has taken place, however, without the development of a more comprehensive formulation of the problem.
Because exposure measurement error may have substantial implications for interpreting epidemiologic studies on air pollution, particularly the time-series analyses, we developed one systematic conceptual formulation of the problem of exposure error in epidemiologic time-series studies of air pollution and considered the possible consequences for relative risk estimation. We used available and relevant data to obtain rough estimates of the magnitudes of the effects of measurement error for one city.

Overview of Measurement Error Effects in Regression Models
The fundamental concepts of how exposure error can affect an epidemiologic study of pollution and health can be shown by considering the effects of exposure measurement error in a standard linear Gaussian regression model. The effects in Gaussian models have been discussed in full detail elsewhere (2,6,7,23,24). For simplicity, consider a regression of the health response (e.g., log mortality rate on day t) and predictors (e.g., PM10, O31 and weather): yt= a + xxt+£ [1] where a and x are regression coefficients to be estimated, and et represents residual error that is assumed to be independent of xt. Here Ox is the expected change in mortality per unit change in true exposure. Given observations (xp y), t = 1, ... T and appropriate assumptions about the distribution of the residuals, ordinary least-squares estimation provides optimal (unbiased and minimum varianced) estimates of the regression coefficients. Now we assume that instead of the true exposure levels x,, we have only an imperfect measure of exposure, denoted z. The overall difference between xt and zt comprises multiple components of error, induding differences between individual-and population-average exposures; between population-average exposures and ambient levels at central sites; and between actual ambient levels and the measurements of those levels. Suppose we regress the health outcome yt on the imperfect zt rather than xp which is unavailable: yt= a*+ Zt + Et. [2] How will ,z differ from xv To answer this question, we first assume that zt is a surrogate for xt, which means that, given xp there is no additional information in zt about yt. We then can distinguish two fundamentally distinct types of relationships between the true and measured exposures, which represent poles of a measurement error continuum. The first type is referred to as the classical error model (7), in which we assume that z is an imperfect measure of x, so that the average z within each x stratum equals x [E(zlx) = x]. Then it follows that the measurement error zx is uncorrelated with the true value x. This classical model is a reasonable one for the difference between measured ambient levels of pollution and the true values for a measuring device that is unbiased. That is, when the true level of pollution is x, an unbiased instrument will measure x on average, even if individual measurements z differ from x. The second type ofmodel for measurement error is the Berkson error model (2). In this model, we assume that the average value of the true exposure x within each stratum of measured level z equals z [E(Alz) = z]. This Berkson model is appropriate when z represents a measurable environmental factor that is shared by a group of participants whose individual exposures x might vary because of time-activity patterns. For example, z might be the spatially averaged ambient level of a pollutant without major indoor sources and x might be the personal exposures that, when averaged across people, match the ambient level.
Classical and Berkson models for exposure measurement errors represent two extremes of a continuum. Most exposure errors combine elements of each, but because the consequences on risk assessment of classical and Berkson errors differ, it is useful to consider each in turn. In the case of the Berkson error, if we regress yt on zp rather than on xp the estimate z is an unbiased estimate of the coefficient x that would be obtained by regressing yt on the actual exposure x,. Having zt rather than xt does not lead to bias in the regression coefficients under the surrogacy assumption. The exposure measurement error increases the variance of the regression coefficient, however, because having zt rather than xt is obviously not as informative about the coefficient P, Bias is not introduced, however. The same is true if the average x at each value of z differs from z by a fixed amount a, i.e., E(xlz) = z -a.
In contrast, under the classical error model Z, obtained by regressing yt on the imperfect measure exposure zp, is a biased estimate of Rx. In the simple linear regression with one explanatory variable, 0Z is expected to be smaller than P,, or attenuated. The degree of attenuation increases as the variance of the exposure error increases. Again, a constant difference in the expected values of the two measures does not change this result.
It is useful to establish these results on the effects of exposure error on simple linear regression coefficients and helpful to do so in advance of considering a multiple regression case. The model of interest is Equation 1, but because xt is unobserved we instead might regress yt on zt (Equation 2).
The question is how will Z from Equation 2 estimate x in Equation 1. Under the Berkson model, xt is assumed to vary about zt, so that by Equation 1, E(ytizt) = a + 0 E(xtIz) = a + xz1, [3] Comparing Equation 3 and Equation 2 shows that P = Px in the Berkson error case; that is, z is an unbiased estimate of fx. Thus, the estimated regression coefficient is biased toward zero. In one pertinent case, x = 0, the naive estimate Z is unbiased with E(Iz)= X = 0; that is, under the classical error model, measurement error does not lead to spurious associations if there is truly no association. Random variation, of course, can produce such associations by chance, as it can when there is no measurement error. How-ever, the probability of such false positive associations (the type 1 error rate) remains the same.
Realistic models for estimating the effects of air pollution on mortality have elements of both classical and Berkson error models. In general, the effect of such exposure errors is intermediate between the two extreme models. The effect of measurement error, therefore, likely depends on the direction and magnitude of the correlation of measurement errors with the measured exposures and not just on the variance of the measurement errors.
More complex multipollutant models are often applied in an attempt to estimate the independent effect of a pollutant present in a mixture with other pollutants. For example, in an analysis of air pollution and mortality in Philadelphia, Kelsall et al. (25) regress mortality on as many as five pollutants. Because little empirical evidence about the simultaneous errors in multiple pollutants is currently available, we only lay a foundation that can inform the design of future studies, as discussed in "Summary and Research Recommendations." Confining our attention to the classical and the Berkson error cases, we again assume a linear regression model of the form given by Equation 1, where x, now represents a vector of exposure variables, with a corresponding vector of regression coefficients Px and zt denotes a vector of measurements of each exposure variable. In the Berkson error case, the assumption that xt is an imprecise version of zt or E(xtIz) = z, still assures that the estimates of the regression coefficients are unbiased, as in the univariate instance. Under the classical error model, however, the multiple regression extension is not so straightforward. We again assume that Zt is an imprecise measure of xp i.e., E(ztlx) = x,. To compute E(xtlz), the average x, at each zp let V denote the covariance matrix of xt and let Tdenote the covariance matrix of the difference 8t = Zt -Xt, and, as before, we assume that 6 and x are independent. The matrix generalization of the earlier result is that 0 , I = 0ixC, where C = T (T + W1). Now it is no longer true that Ozj < , j for each component (j) and estimates of regression coefficients can be biased toward or away from the null; that is, positive associations can be produced when the component is correlated with at least one component having a nonzero effect, even though the true coefficient for a particular component is zero. Table 1 illustrates the magnitude of bias that can result from regressing yt on two predictors zlt and z2t instead of on xlt and x2, This example might refer to estimating the effects of PM 0 and 03 on mortality when ambient values (z values) instead of personal exposure (x values) are available. We assume zlt = xlt + 6It and z2t = Z2t + 62? V 1 = var(xl) = V22 = var(x2) = 1. Table 1 Table 1 is intended to illustrate the consequences of measurement error in the two-predictor model.
The first line of Table 1 refers to an example in which there is no correlation between xlt and x2t and there is equal variability of the two exposure errors 6 t and 62 and these errors are not correlated; that is, the error in one predictor does not predict the error in the other. Here, there is an equal degree of attenuation in the coefficients for the two variables. With unequal variances but no correlation, i.e., the sixth row, the degree of attenuation is greater for the variable with greater variance. If the exposures are correlated but the errors are uncorrelated (the second and third rows), the two effect estimates are similarly altered with the direction of the effect depending on the sign of the correlation. Introducing correlation between the errors, i.e., the fourth and fifth rows, has an effect that depends on the pattern of correlation. The bottom half of Table   1 shows more complex patterns with differing patterns of correlation and variation of the two errors. Some of the scenarios introduce substantially different effects of the two variables, but none yield effect estimates above the true value of one, even with more extreme differences in error variances or the two correlations. Table 2 also addresses the consequences of measurement error in a two-variable model, but in this example only one variable (x2) has a true effect; the other exposure (xl) has no effect on the health outcome (y). Either correlation between xltand x2, or their errors can introduce an apparent effect of xl on y. Some scenarios of variance and correlation even bring the apparent effects of the two variables quite close (e.g., the tenth and eleventh rows), but in every case, including more extreme situations than shown, the estimate for the true predictor ([2) is always larger than for the null predictor (PI).
Some general conclusions can be offered concerning multipollutant models under this simple classical error model. Conclusion 1. There is a general tendency for the coefficient from the regression on zt to be smaller than the corresponding coefficient from the regression on xp i.e., kjz> I^j if all 4>0. Conclusion 2. The degree of attenuation of each coefficient depends, in large part, on its measurement error variance relative to the variance of the true exposure-i.e., TJy V. .. Thus, the coefficients for variables that are measured with considerable error will be attenuated more than those of variables with less error. Conclusion 3. Depending on the correlation structure of the attenuation matrix C, some of the effect of one variable, _ may be transferred to the estimate of another variable's effect, f.,k Such transfers of effect are generally from a more poorly measured variable to a better measured variable. However, for such transfers to be large, the true exposure Environmental Health Perspectives * VOLUME 108 1 NUMBER  variables or their measurement errors need to be substantially correlated. Conclusion 4. As a consequence of conclusion 3, the estimate of a parameter can be biased away from the true value. However, this type of bias generally arises only with a very strong negative correlation between the measurement errors (e.g., rows 9-11 ofTable 2).
Conclusion 5. Also as a consequence of conclusion 3, there will generally be spurious associations for a variable x; that, in fact, has no effect only if xj is substantially correlated with one or more variables which actually have an effect. Generally, the correlation among the errors has a larger influence on the bias than the correlation among the true pollutant levels.
These condusions are obtained from and therefore pertain to the dassical linear regression model with two predictors, assuming that zt is a surrogate for xt (nondifferential errors). The actual exposure measurement situation in the air pollution-mortality context is obviously more complex. First, log-linear, not linear, models are used, although the degree of nonlinearity is usually small in mortality studies. Second, the measurement errors are not purely of the classical nondifferential type. For example, the degree of error for gaseous pollutants may depend on temperature or other covariates. Finally, errors may be multiplicative rather than additive. Nonetheless, the linear regression with dassical measurement error is a leading case that provides insight into the major possible consequences of exposure errors.

Framework for Assessing Measurement Error Effects in Pollution-Mortality Studies
We now build on the fundamental concepts underlying statistical models of exposure measurement error and focus on the specific log-linear regressions used for assessing the pollutant-mortality association, controlling for weather variables. Our discussion is based on the premise that the ideal investigation of the health effects of air pollution would be conducted at the individual level with measurements of personal exposure to pollutants. However, exposure and mortality data are generally only available after aggregation to a municipal level; little or no data from indoor air monitoring are available. Finally, air pollutant measurements are imprecise and this imprecision has consequences for estimates of pollutant effects on mortality.
To investigate the effects of exposure error in the log-linear regressions widely used to assess the pollutant-mortality association, consider the following model for an individual's risk ofmortality: lit= Xoeitexp(x2tp) [4] where Xt is the risk of death for person i on day i Xoit is that individual's baseline risk in the absence of exposure, i.e., xit = 0, and exp(xit P) is the relative risk of death associated with the explanatory variables xit. Let yit = 1 if person i dies on day t and 0 if that person does not. We typically observe the total number ofdeaths for a population nt yt = y=;yi.' where nt n is the population size on day t.

By Equation 4
, the expected total number of deaths It in a community is Xt = Eyt= Xit oit exp(xij). [5] l l In analyzing population-level data on mortality and air pollution, log-linear regressions of the following form have been fit it= exp[(t) + zpJ + ufj [6] where s(t) is an arbitrary but smooth function of time introduced to control for the confounding of longer-term trends and seasonality, zt is the average of multiple monitor measurements of ambient pollution measurement for day t, and ut are other possible confounders such as temperature and dew point temperature on the same and previous days.
If the regression coefficient jx for a pollutant in the personal risk model Equation 4 is  x,t for person i on day t and the available ambient values zt measured with error by monitors. Assuming, for simplicity, a high degree of spatial homogeneity in ambient levels, personal exposure is contributed to by Zt' the true outdoor level, and wi, the indoor level, which is also influenced by Zt from penetration of the pollutant in outdoor air into indoor spaces. For example, personal exposure to PMIO is determined by the time spent outdoors, the concentration during that time, and by the concentrations in indoor environments that are determined by indoor sources such as cigarette smoking and the penetration of particles indoors because air is exchanged between the outdoors and the indoor environments. Figure 1 further shows that the personal risk of dying is influenced by an individual's baseline risk in addition to the unobserved personal exposure to pollutant x1.. Only the measured ambient pollution data are observed and are therefore shown in a rectangular box.
In considering the consequences for z as an estimate of x with an imprecise measure of ambient pollution zt, rather than actual personal exposure xip it is useful to begin by decomposing the pollution measurement difference between xit and zt into three components: X= Zt+ (Xit-x) +(X t + (ZtZt)' [7] Here, (xit --) is the error due to having aggregated rather than individual exposure data; (xt -Z) is the difference between the average personal exposure and the true ambient pollutant level; and (Ztz) represents the difference between the true and the measured ambient concentration.

The first term (x,tx) is an example of
Berksonian error so that, in a simple linear model, having aggregate rather than individual exposure does not itself lead to bias into the regression coefficient. The second term (X -Z*) is not Berksonian and is likely to be a source of bias. The final term (ztz) is largely of the Berkson type if the average of the available monitors zt is an unbiased estimate of the true spatially averaged ambient level Zt* We can now further study the effects of these three terms on risk estimation by substituting the decomposition in Equation 7 into Equation 5. After some straightforward calculations detailed in the "Appendix," the expected number ofdeaths on day tcan be written Eyt = exp[log(nt 0) + zJ3x + {(5(w) -+ (Xt-z) + (zt -Zf}j. [ Articles * Measurement error in time-series studies Here P, is the personal log-relative risk of interest from Equation 5. Note the approximation Equation 8 retains only linear terms in the expansion of an exponential function. The second-order terms are an order of magnitude smaller and are ignored to simplify the exposition. For studies of particulate pollution effects on mortality, the effect sizes are on the order of 1 or 2% so that ignoring secondorder terms should not qualitatively affect the results. In studies of morbidity, higher order terms may be more important.
The total baseline risk (n,X0) almost certainly varies smoothly over time because it is an average risk over a large population. Hence, it will be appropriately controlled for in log-linear regressions by inclusion of the smooth s(4 in Equation 6. We now consider Z1I3x and the three components of error in turn. The first error term 3t(w) -xt is the difference between the baseline risk-weighted average personal exposure and the unweighted average personal exposure. It derives from the Berkson error (xix,) and produces no bias in the linear unaggregated model. This difference due to risk weighting in our log-linear model with person-specific baseline risks is likely to be small and to vary slowly over time. Hence, it can be adequately controlled by inclusion of the smooth function s(t) in the log-linear regression of yt on z, One scenario in which this difference would vary from day to day and therefore not be adequately controlled would occur if the more frail individuals were to follow pollution reports (or a correlate such as weather) and reduce their exposures to ambient air on high pollution days by, for example, staying indoors. Current warning systems for air pollution alerts are intended to reduce exposures of susceptible persons in this fashion.
The second error term Xt-Zt is non-Berksonian and has the greatest potential to introduce bias in the estimate z when Z is correlated with X -Zt Even if the terms are  uncorrelated so that r will be a roughly unbiased estimate of 0,,, it will reduce efficiency relative to a study in which x, is available because z, andzt share the same coefficient in Equation 8. The difference Xt-Z between average personal exposures and the true ambient value can be analyzed further by considering an individual personal exposure xi, Because individual i's exposure on day tderives either from indoor or ambient sources, we can write xit= Ct z* + (1 -ao.)la where Ii is the concentration of pollutant generated by indoor sources such as tobacco smoke and pets and ait is his or her fraction ofexposure from ambient sources that take place either outdoors or result from the penetration of ambient pollution indoors. It follows that x-=a* + It where It = (1al)IJn1. That is, the average personal exposure is proportional to the ambient level offset by the effects of the population average of the non-ambient indoor sources.
Wllson and Suh (26) argued that the daily population average concentrations of fine particles derived from indoor sources It are approximately independent of ambient levels zt across time. When this is true, failure to measure indoor sources will not introduce further bias in the estimation of [x because the deviations due to indoor air exposure are a second example of Berkson error, and the errors will tend to cancel one another out when averaged over the population. Never-theless, Z is only proportional to xtso that even ifa varied little over time (at a), the coefficient tZ from a regression of yt on z would estimate a,, not fx. Hence, if 20% of daily exposure results from indoor sources independent of the ambient levels, the regression on ambient levels will yield coefficients that are roughly 20% smaller than would have occurred with actual personal exposures. However, this may be the appropriate coefficient for policymakers seeking an estimate of the effect of an inarguable measure of ambient levels. This, however, assumes that particles from indoor sources and outdoor sources are identical; that is, they are similar in composition and toxicity. If this is not the case, then the two types of partides are more appropriately treated as separate pollutants, and the personal exposure measure desired would be aotz, the personal exposure to particles from outdoor sources. Studies that use sulfates as a tracer for particles from outdoor sources indicate that indoor/outdoor ratios are < 1. Because people spend most of their time indoors, this suggests that axt will be < 1 and that the second term in Equation 8 will be negatively correlated with z,, and will bias the estimated coefficient downward. This also illustrates that the model is not restricted to cases where E(x) = E(z).
The final of the three error terms in Equation 8, zt -z1 represents the instrument measurement error in the ambient levels; like xit-xt, it is dose to the Berkson type. This term would tend to be cancelled out by spatial averaging across multiple unbiased ambient monitors in a region. For example, Kelsall et al. (25) averaged daily total suspended particulate data from up to nine monitors in their analysis of the effects of partides on mortality in Philadelphia. However, in many cities there is only one monitor or a few monitors operating concurrently. Even with a small number of monitors, longer term drift in instruments will not substantially affect estimates of Px because the time-series models control for such trends by indusion of s(4 in Equation 6.
For this final error term to cause substantial bias in 0, the error it -Zt must be strongly correlated with zt at shorter time scales. Further investigations of this correlation in cities with many monitors are warranted.
We have discussed three components of measurement error: a) an individual's deviation from the risk-weighted average personal exposure; b) the difference between the average personal exposure and the true ambient level; and c) the difference between the measured and the true ambient levels, which includes spatial variation and instrument error. We argue that the first and third components are of the Berkson type and therefore are likely to have smaller effects on the relative risk estimates. However, the second component can be a source of substantial bias if, for example, there are short-term associations of the contributions of indoor sources with ambient concentrations. We present one simple analysis of the Particle Total Assessment Methodology (PTEAM) data (27) that illustrates how we can further study the effects of the most important second component.

Evaluating Potential Measurement Error Bias in Pollutant-Mortality Relative Risk Estimates
The "Framework for Assessing Measurement Error Effects in Pollution-Mortality Studies" can be used, in combination with data on the components of error, to quantify the consequences of exposure measurement error. One of the few available data sets with ambient and personal measurements will be used to illustrate one approach. We used daily measurements of personal exposure for persons followed in the PTEAM Study (27) to quantify the difference between concentration measured by an ambient monitor and the average of personal exposures. We studied one approach for estimating the size of bias in estimated PM10-mortality regression coefficients Z as an estimate of the true relative risk for personal exposure P., with data Environmental Health Perspectives * VOLUME 1081 NUMBER  November 1990. In addition, a daily average PM10 value from an ambient monitor positioned near the homes was also collected; Pellizzari and Spengler (29) provided details on the methods used to collect these data. We used the PTEAM Study data to estimate the correlation between the daily PM1O concentration for the ambient monitor zt and the difference between the average personal exposure and concentration measured by the ambient monitor xt -z( These estimates correctly account for the varying number of observations on a given day; however, the average personal exposure value is based on relatively few measurements and is therefore more variable across time than the actual mean exposure. Equation 8 indudes a weighted average of personal exposures, with weights determined by the baseline personal risk for each individual. Those weights were unavailable in the PTEAM Study and hence, we used an unweighted average. Figure 2 displays a time-series plot of the daily ambient values and the average personal exposures. The correlation across time of these two series is estimated to be 0.58 [95% confidence interval (CI), 0.35-0.74]. This correlation is much greater than the more widely cited cross-sectional correlation from the same study. It would likely be even greater if the mean personal exposure was calculated on a larger number of persons each day. The corresponding correlation across time between the ambient monitor concentrations and the daily differences between the personal and ambient values is -0.63 (CI, -0.77 to -0.42). Hence, the hypothesis that the measurement errorzt is uncorrelated with zt is not consistent with the PTEAM Study data. Some bias in the regression coefficient is therefore expected. Because the correlation of xtz, and zt is negative, the coefficient i in the regression on zt wiul tend to underestimate the co-efficient in the regression on xt in a single-pollutant analysis. We now assess the size of the bias that will result from this measurement error.
Addressing the Bias in PM1-Mortality Regression

Coefficients
The PTEAM Study results or other, perhaps more appropriate, data sets on the difference between average risk-weighted personal exposure and ambient monitor concentrations, can be used to estimate bias in the results of loglinear regression models.
If available, we would have used the average personal exposure series, X-for at-risk residents of each city in the standard log-linear regression model rather than zp as was used in the original analyses. We would then have compared the regression coefficients obtained when xt is the predictor with those using zt to assess the bias.
Obviously, xt is not available except in special circumstances. However, from the PTEAM Study data (shown in Figure 2) or similar data, we can estimate the relationship of xtand zp for example, by assuming: xt=0o + 0zt+et [9] where 00 and 01 are the intercept and slope to be estimated from the available data. We can then use the fitted Equation 9 to predict the unobserved xt from the available zt and then use the predicted value xt as the desired exposure values when estimating the pollution-mortality relative risk ., In fact the estimate of Px has the simple form = 11 This well-known approach to adjust for exposure measurement error is called regression calibration (7). As an illustration, we applied this strategy to a regression of daily mortality Calibration is easy to implement and apply. Its limitations are that confidence intervals for x depend on large sample theory and that it does not extend easily to situations where multiple sources of information about the xp, zt relationship are available.
It is simple to overcome the possible limitations of calibration by using a simulated *~~~~~~~~~~~A value`* rather than the predicted value xt from Equation 9. That is, we use Equation 9 to simulate the average personal exposure, ', from the ambient exposure, zp for a city or period of interest when xt is not available, under the assumption that the estimated Os and var(e) are applicable. This simulated series -* is then used instead of zt in the loglinear regression. The result is one estimate of _i3all it Ox, If we then repeatedly simulate xts and fit the log-linear regression for each to obtain X , , we obtain a distribution of 0s. The difference between the mean of the simulated xs and the 0_, derived from the log-linear regression of mortality on Zt is a measure of the bias resulting from having zt rather than the average personal exposure for that city. By simulating xts rather than using a fixed predicted value xt, we properly account for nonlinearities and sources of variation in 0 and can extend the analysis to more complicated situations. Figure 3 shows the distribution of the s for Riverside (solid curve). Also shown is the normal approximation of the likelihood func-tion for the coefficient . from the log-linear regression of mortality directly on z, (dotted curve). Solid and dotted lines are at the centers of these distributions. We find that the 0,ls have a mean 1.42% 1.76). Hence, measurement error has biased the result toward the null. Second, the distribution of the jxs is more dispersed than the distribution of fiZ. This is because we have taken into account the variability due to having zp, not x, i.e., arising from var(s,) in Equation 9. The results are very similar to what we obtain from calibration. This calculation assumes that the estimated relationship between xt and zt for the PTEAM Study is the true one, and hence, we ignored a second component of uncertainty due to estimation of the relationship between k and it from the finite sample size of the PTEAM Study data taken at one site and at a particular time period. That is, even if we assume that the relationship between and z, is known, estimating the association of mortality with is less precise than with zp given only zt in that particular city. Of course, the relationship of and zt is not precisely known and needs to be quantified further. Dominici et al. (30) provided a more complete analysis of the bias in . as an estimate of x using the PTEAM Study and four other data sets and a more complete statistical model. Their findings were qualitatively similar to those presented here. Finally, our assessment of bias assumed that the health effects of personal exposure to partides originating outdoors and indoors are the same. To assume otherwise would require substantially more detailed data and modeling.

Recommendations
The differences between true personal exposure for every individual (xi) and measured ambient concentrations, averaged over a few fixed imprecise monitors (z), are inherently complex, as is the effect of this exposure measurement error on estimates ofpollution-mortality relative risks. Nonetheless, it is useful and imperative to analyze these effects in light of our current understanding of the measurement process. This paper presented one framework for doing so. We distinguished two extremes of a continuum of types of measurement errors: Berkson and dassical errors. The former is likely to create little bias in mortality-relative risk estimates; the latter has more serious consequences.
We posited a relative risk model in which an individual's hazard of death on a given day is expressed as a function of his or her personal exposure, which is decomposed to highlight three types of exposure errors. We then aggregated the model to produce the model for the expected total deaths in a population used in most time-series analyses. This model showed that a risk-weighted average personal exposure measure is the desired exposure measure. The likely consequence of using ambient concentrations instead is to underestimate the pollution effects. In contrast, differences between individual exposures on a given day and the risk-weighted average of personal exposures are examples of Berkson error and are not likely to cause substantial bias in coefficients from time-series morbidity studies. Our analysis suggested that the largest biases in inferences about the mortality-personal exposure relative risk will occur because of the more complex errors between ambient and average personal exposure measures. If indoor sources produce partides of similar composition and toxicity as outdoor source particles, indoor sources may be a major component ofthis error. Finally, we used the best available data (from the PTEAM Study in Riverside) with both personal exposure and ambient time-series data to quantify the size of this error. Our analysis indicated that the coefficient obtained from regressing mortality on measured ambient levels (z) is smaller than what we expect ifwe regress mortality on average personal exposure (-x).
For tractability and clarity, we conducted a first-order analysis of exposure errors and ignored possible second-and higher order effects in which daily fluctuations in the variance of personal exposures across a population or in the covariations among the measurement errors could introduce additional biases. Second-order terms will be insignificant in studies of particulate effects on mortality where the first-order terms are on the order of percent. Such higher order analyses for other studies of, for example, morbidity, are beyond the scope of this paper and will require substantially more detailed models and data. It is, however, possible that higher order effects are important; further investigation is necessary.
Epidemiologic research is necessarily limited by the quality of the health outcome and risk factor measurements (31). Time-series studies of the acute effects of air quality on mortality are subject to the limitations posed by the available measurements of pollution levels. The generic criticism-that measurement errors render the results of such time-series models uninterpretable-is inaccurate. The consequences of measurement error can be quantified, although only a few informative data sets are presently available. Differences between the average personal exposure and ambient measurements are the most likely source of substantial bias. Data should be collected for the comparison of risk-weighted average personal exposure with ambient levels in several cities with varying Relative rate (%1O0 ig/rm3) Figure 3. The solid line is the distribution of the relative rate , obtained when the simulated series it of the total personal exposure is the predictor in the log-linear regression. The dotted line is a normal approximation of the distribution of the relative rate , obtained when the ambient concentration zt in Riverside, California, is the predictor in the log-linear regression.
degrees of spatial heterogeneity in ambient levels, population composition, and indoor pollution sources. Given such data, models like those summarized by Dominici et al. (32) can be used to quantify more precisely the biases due to pollutant measurement errors. This paper focuses on the effects on relative risk estimates of using zt (measured ambient particle levels) rather than xit (actual personal exposures in log-linear regressions). Such effects are important from a scientific perspective to quantify the health risks of exposure to particulate pollution. From a regulatory perspective, the effect of having the imprecise zt rather than the true ambient value zt may be of greater interest because it is ambient levels that may or may not be regulated further. A more detailed error analysis of the zt-zt difference would investigate the spatial variation in particulate levels and how the number of monitors used to calculate zt reduced this source of measurement error. In measurement errors in a single pollutant measure, PM1O, simultaneous errors in several pollutants can complicate the analysis. However, qualitative biases-that is, changes in the sign of a coefficient-can occur only when the measurement errors for different pollutants are highly correlated with one another. This level of correlation might arise if two or more pollutants are measured by the same instrument (e.g., different fractions of particulate matter) or if multiple instruments are housed in the same location, which is subject to atypical exposure patterns. The possibility nevertheless requires detailed investigation, because in this case the findings of epidemiologic studies could be misleading. Personal exposure studies that collect Environmental Health Perspectives * VOLUME 108 1 NUMBER  We considered the effects of exposure measurement error on regression coefficients from log-linear models in which serial correlation is accounted for using flexible smoothing splines. An alternate analytic strategy is to fit a linear regression with time-series errors [ARIMA model (33)]. In certain specific time-series models, the degree of attentuation due to dassical error might be reduced because to account for the autocorrelated errors, the ARIMA filters or smooths both the responses and the predictors that might reduce the degree of measurement error. Further research on this possibility is warranted.
The measurement error framework and the illustrative calculations discussed here make apparent several open questions and opportunities for additional data collection. These opportunities would enable more accurate quantification of the effects of measurement error in assessing the air pollution-mortality relationship. In relation to single-pollutant models, the two most important questions are a) Is the average personal exposure to pollutants from indoor sources correlated over time with ambient levels? and b) Does the difference between baseline riskweighted average exposure and population average exposure vary slowly over time?
For models with multiple pollutants, the additional key question follows: How do the components of error identified in Equation 5 covary across pollutants? For example, how do the differences between actual ambient levels and the measured levels correlate across the different pollutants and how do these differences depend on the true values of other pollutants or covariates? Wldson and Suh (26) conducted a metaanalysis of data from multiple sites and conduded, in answer to the first question above, that concentrations of fine partides originating from indoor sources are independent of ambient levels over time. To