Home | About | Journals | Submit | Contact Us | Français |

**|**HHS Author Manuscripts**|**PMC2517157

Formats

Article sections

Authors

Related links

J Am Diet Assoc. Author manuscript; available in PMC 2008 August 18.

Published in final edited form as:

PMCID: PMC2517157

NIHMSID: NIHMS33120

Janet A. Tooze, PhD, MPH, Assistant Professor,^{} Douglas Midthune, MS, Mathematical Statistician, Kevin W. Dodd, PhD, Mathematical Statistician, Laurence S. Freedman, PhD, Director, Biostatistics Unit, Susan M. Krebs-Smith, PhD, MPH, RD, Nutritionist, Amy F. Subar, PhD, MPH, RD, Nutritionist, Patricia M. Guenther, PhD, RD, Nutritionist, Raymond J. Carroll, PhD, Distinguished Professor, Professor of Statistics, Nutrition and Toxicology, and Victor Kipnis, PhD, Mathematical Statistician

Janet A. Tooze, Section on Biostatistics, Department of Public Health Sciences, Wake Forest University School of Medicine, Medical Center Boulevard, Winston-Salem, NC 27157, Phone: (336) 716-3833, Fax: (336) 716-6427, Email: ude.cmbufw@ezootj;

The publisher's final edited version of this article is available at J Am Diet Assoc

See other articles in PMC that cite the published article.

We propose a new statistical method that uses information from two 24-hour recalls (24HRs) to estimate usual intake of episodically-consumed foods.

The method developed at the National Cancer Institute (NCI) accommodates the large number of non-consumption days that arise with foods by separating the probability of consumption from the consumption-day amount, using a two-part model. Covariates, such as sex, age, race, or information from a food frequency questionnaire (FFQ), may supplement the information from two or more 24HRs using correlated mixed model regression. The model allows for correlation between the probability of consuming a food on a single day and the consumption-day amount. Percentiles of the distribution of usual intake are computed from the estimated model parameters.

The Eating at America's Table Study (EATS) data are used to illustrate the method to estimate the distribution of usual intake for whole grains and dark green vegetables for men and women and the distribution of usual intakes of whole grains by educational level among men. A simulation study indicates that the NCI method leads to substantial improvement over existing methods for estimating the distribution of usual intake of foods.

The NCI method provides distinct advantages over previously proposed methods by accounting for the correlation between probability of consumption and amount consumed and by incorporating covariate information. Researchers interested in estimating the distribution of usual intakes of foods for a population or subpopulation are advised to work with a statistician and incorporate the NCI method in analyses.

When using dietary assessment among populations or individuals, investigators are often interested in capturing usual intakes – that is, long-term averages. The 24-hour dietary recall (24HR) provides rich details about dietary intake for a given day, but collecting more than two 24HRs per individual is impractical in large surveys such as the National Health and Nutrition Examination Survey (NHANES). Therefore, it is necessary to employ statistical methods to estimate usual dietary intake.

Researchers are interested in estimating the usual intake of foods to assess compliance with food-based dietary recommendations and to relate food intake to health parameters. Unlike most nutrients, which are consumed daily, estimating usual intake of episodically-consumed foods presents the following unique challenges for statistical modeling (see the Glossary for a definition of “statistical modeling” and related terms):

- accounting for days without consumption of a particular food or food group;
- allowing for consumption-day amount data that are generally positively skewed and have extreme values in the upper tail of the intake distribution;
- distinguishing within-person variability, which consists of day-to-day variation in intake and random reporting errors, from between-person variation;
- allowing for the correlation between the probability of consuming a food and the consumption-day amount; and,
- relating covariate information (for example, sex, age, race, ethnicity, or education level) to usual intake (1).

As discussed in Dodd et al. (1), two other methods have been used to estimate the distribution of usual intake of episodically-consumed foods with a few days of 24HRs: (i) the distribution of within-person means (WPM), and (ii) the method developed at Iowa State University for estimating the distribution of foods (ISUF) (2). The WPM method usually leads to biased estimates of the prevalence of either inadequate or excess food intake because it does not meet any of the challenges listed above. In particular, because the WPM method does not meet challenge (C), and thereby includes within-person variability, the variance of usual intake is inflated. The ISUF method meets the first three challenges, but it does not allow for correlation between probability and amount and cannot incorporate covariate information regarding usual intake.

This paper describes a new statistical method that was developed at the National Cancer Institute (NCI) to meet all five of the challenges noted above, using two 24HRs, and evaluates the new statistical method's application to estimating the distribution of usual intake of episodically-consumed foods.

In the NCI method, we assume that the 24HR is an unbiased instrument for usual intake of episodically consumed foods. This assumption has two components. First, we assume that the 24HR does not misclassify the respondent's food consumption, i.e., if a food actually is consumed on a surveyed day, the food will be reported on the 24HR; and if a food is not consumed, it will not be reported on the 24HR. Second, we assume that the 24HR is an unbiased measure for the amount of food consumed on the consumption day. This does not mean that the 24HR captures the amount of food consumed by an individual exactly on each recall – at a given time an individual may report more or less than was actually consumed – but over many days it produces the correct average intake. This assumption will be discussed in detail later in the paper. In addition to the assumptions made about the 24HR, we make the usual assumptions for parametric regression analysis in our models. In particular, we assume that, after an appropriate transformation, the amount of food consumed on a consumption day is approximately normally distributed.

The NCI method for estimating the usual intake of foods has two steps. The first step consists of fitting a two-part statistical model that describes the relationship between usual intake and covariates and estimates the variability of intake both within and between individuals. We adapted a two-part model with correlated person-specific random effects, developed by Tooze et al. (3), for this purpose. Similar to the ISUF method, the statistical model represents usual intake as the product of the probability to consume a food on a given day and the usual consumption-day amount. The amount data are transformed to approximate normality, using the Box-Cox transformation (4), as part of the model-fitting process. To account for the correlation between the probability and amount that exists for most foods, as described below, the two parts of the model are linked.

The second step of the method involves additional statistical procedures that, depending on the application of interest, are used to obtain the final “product” of the analysis. Examples of those products include estimates of the distribution of usual intake in a population or subpopulation of interest, or estimates of individual intakes. The latter may be used to assess diet-health relationships. Because these are clearly varying endpoints, different statistical procedures are required. However, because the data share a common structure in each case, the same statistical model is used to obtain parameter estimates, which are the inputs for the final step. This paper focuses on the statistical model and its application for estimating the distribution of usual intake for populations and subpopulations. Predicting individual usual intake and relating it to health outcomes are beyond the scope of this paper.

The first part of the model estimates the *probability* of consuming a food (positive intake reported on the 24HR) using logistic regression with a person-specific random effect (mixed model). The logistic regression model incorporates covariates to represent the effect of personal characteristics, such as age, sex, or body mass index, on the probability of food consumption. The person-specific effect is a factor that allows an individual's consumption probability to differ from the population level. It may be thought of as the individual's personal tendency to consume a food. The probability of consumption is estimated from two or more 24HRs, accounting for covariates. Symbolically, Part I may be represented as:

$${\text{Logit}(\text{24HR Probability})=\text{Intercept}}_{I}+{\text{Slope}}_{I}\times {\text{Covariate}+\text{Person-Specific Effect}}_{\text{I}}$$

[1]

where, for probability *p*,
$\text{logit}\phantom{\rule{0.2em}{0ex}}(\text{p})=\text{log}\phantom{\rule{0.2em}{0ex}}\left(\stackrel{\text{p}}{\phantom{\rule{0ex}{.6}}}/\underset{1-\text{p}}{\phantom{\rule[--.8]{0ex}{0ex}}}\right)$. The intercept, slope, and variance of the person-specific effect are the model parameters, and subscript I indicates their association with Part I. Although one covariate is shown in equation [1], the model allows for multiple covariates or no covariates.

The second part of the model specifies the consumption-day *amount* of a food using the 24HR data on a transformed scale. Similar to Part I, Part II may incorporate covariate information to estimate amount. As before, the covariates, which need not be the same covariates as in Part I, represent the effect of personal characteristics on the consumption-day amount. This part of the model also includes a person-specific effect as well as within-person variability due to day-to-day variation in an individual's intake and other sources of random error. The model for Part II is:

$${\text{Transformed 24HR Amount}=\text{Intercept}}_{\mathrm{II}}+{\text{Slope}}_{\mathrm{II}}\times {\text{Covariate}+\text{Person-specific Effect}}_{\mathrm{II}}+{\text{Within-person Variability}}_{\mathrm{II}}$$

[2]

where subscript II indicates that these parameters are associated with Part II, and differ from those in Part I. Two or more 24HRs on a number of individuals with reports of the food of interest are required to distinguish between- and within-person variation. The model is specified on the transformed scale where the person-specific effect and within-person random variability are normally distributed.

Unlike the ISUF method, in which the two parts of the model are assumed to be independent and are estimated separately, the NCI method fits both parts simultaneously, which associates probability to amount in two different ways. First, the two person-specific effects are modeled as correlated random variables. Second, some covariates may be the same in both parts of the model, inducing correlation between them. By linking the two parts of the statistical model, the relationship between probability and amount is accounted for, meeting challenge (D) stated previously.

The model is fit by the maximum likelihood method, using a SAS software (version 8.2, SAS Institute, Cary, NC) macro. To account for the correlation, all of the model parameters are estimated at the same time using a non-linear mixed effects model. In addition, the Box-Cox transformation parameter is also estimated as part of the likelihood maximization procedure. The advantage of using the normality transformation within the modeling step, not before modeling as in other methods, is that the 24HR-reported amount is transformed to normality conditionally on the covariates in the model. Estimates are obtained for the model parameters presented in equations [1] and [2], and for the correlation between the person-specific effects.

As described by Subar et al. (5), frequencies from an FFQ, NCI's Diet History Questionnaire (DHQ), are generally positively related to the proportion of 24HRs with reported consumption of those foods. This demonstrates that food frequency information could be a useful covariate in estimating the probability of consumption. Additionally, because of the correlation between the probability of consumption and the amount consumed, frequency responses can contribute to estimating not only probability to consume, but amount as well (5). The sum of the frequencies of several individual FFQ items may be used to represent the frequency of consumption of a food group containing them. The relationship between the FFQ and the 24HR is often non-linear. Consequently, a polynomial model may be used to model the relationship between the FFQ and 24HR.

To estimate the distribution of usual intake, the estimated model parameters are used to simulate a population that has the same characteristics (as described by the values of the covariates) and between-person variability as the sample on which the model was fit. The within-person variation in model [2] is not included since, by definition, it does not contribute to long-term intake.

First, the estimated intercept and slope(s) for the covariate(s) are used to obtain for each individual in the sample an estimate based on the covariate values used to fit the model. To these, it is necessary to add an estimate of the person-specific random effects. Because these effects are unobservable, we generate an estimate from their statistical distribution (bivariate normal with mean zero and variance parameters estimated from fitting the statistical model). To improve the precision of the estimated usual intake distribution, we generate 100 pseudo-persons for each individual in the sample, each with the same covariate values but different simulated person-specific effects.

Because we transform the consumption-day amount data using the Box-Cox transformation during the model fitting process, it is necessary to back-transform the amount data to the original scale before we can obtain estimates of the distribution of usual intake. The back-transformation is similar to the approach used in the method developed at Iowa State University for estimating usual nutrient intake distributions (6). It adds an adjustment term to make the mean of the back-transformed variable match the mean on the original scale as described by Dodd et al (1). Lastly, the mean, standard deviation, and percentiles are estimated empirically from this simulated population.

Estimates of the distribution of usual intake for a subpopulation are made in the same way as those for the total population, except that covariates that define the subgroup are included in the NCI model. When making the estimates for a subpopulation, only the covariate values differ; all other variance components remain the same. This specification leads to smaller standard errors of estimated parameters than stratifying by subpopulation.

One of the challenges in developing this statistical model was to account for the phenomenon that occurs for some food groups: that those individuals who eat a food most frequently tend to eat more of it (Challenge D). To determine how often this happens, we used the Eating at America's Table Study (EATS) data to assess the proportion of food groups that exhibit a positive correlation between the probability of consumption and consumption-day amount on 24HRs. Details of the EATS, conducted in 1997-1998, are published elsewhere (7). We used data from 965 men and women, 20 to 70 years of age, from a national sample, who completed four 24HRs three months apart, followed by NCI's DHQ. The study was approved by the NCI Special Studies Institutional Review Board.

To determine what proportion of 27 Pyramid food groups (described in (5)) exhibited a significant correlation between probability of consumption and consumption-day amount, we calculated the proportion of respondents who reported consumption on none, one, two, three or four of the four 24HRs and the median portion size consumed for each category by sex. Spearman correlation coefficients between the number of recalls with consumption and the consumption-day amount were then computed.

We also used the EATS data to (i) illustrate the application of the NCI method, in comparison with the WPM and ISUF methods, for estimating the distribution of usual intake for two food groups, whole grains and dark green vegetables, for men and women and (ii) estimate the distribution of usual intake by education level for whole grain consumption by men. In these analyses, only DHQ reported frequencies were used; the average portion size information (small, medium, large) was not used. This variable of reported frequencies is very similar to the Food Propensity Questionnaire (FPQ) used in the 2003-2006 NHANES (5). Daily frequencies were determined for each respondent for each food group by summing the line items belonging to each group. The covariate for whole grains included the frequency data for the following DHQ items: breads, crackers, hot cereals, popcorn, ready-to-eat cereal, rice or other cooked grains, and potato/tortilla/corn chips. The covariate for dark green vegetables included frequencies of broccoli, raw greens, cooked greens, and lettuce. The proportions of the line items in each category (e.g., the proportion of bread that was whole grain) were considered when deriving the frequency variables. (See (5) for further details.) To accommodate the non-linear relationship between the 24HR and the food-frequency variables, a polynomial model was used, that is, square root, linear, and quadratic functions of the food frequency variable were included as covariates in the statistical model.

We also estimated the distribution of usual intake of whole grain consumption for three subpopulations of men: those with a high school education or less, those with some college, and college graduates. This was achieved by including the indicator variables for education level as covariates in both parts of the statistical model.

Simulation studies provide a means of comparing statistical methods to each other, as well as to a measure of true usual intake. We conducted a simulation study to evaluate the NCI method and to compare it to other currently available methods for estimating usual food intake distributions. Because we were interested in the performance of the methods for the minimum number of 24HR days, two days were chosen to compare six different methods of estimating usual intake: (i) the two-day WPM, (ii) the ISUF method, (iii) the standard NCI method with correlated person-specific effects and the FPQ as a covariate (as described above), (iv) the NCI method with correlated person-specific effects but without the FPQ (i.e., the covariate(s) generated by the FFQ were removed from the model fittings of equations [1] and [2]), (v) the NCI method with uncorrelated person-specific effects (i.e., the person-specific effects in equations [1] and [2] were assumed to be uncorrelated in the model fittings) with the FPQ, and (vi) the NCI method with uncorrelated person-specific effects without the FPQ. By comparing the NCI method both with and without the correlated person-specific effects and with and without the FPQ, we were able to isolate the effects of these components of the model in this simulation study.

We simulated 200 data sets, each with 2000 pseudo-persons, based on whole grain consumption by women in EATS. For each pseudo-person in a simulated data set, an FFQ frequency value was selected from the actual values in EATS. Next, a probability was generated for that individual, using the average probability from the EATS data stratified by FFQ frequency group (Figure 1a). For example, a pseudo-person who fell into FFQ group 5 (corresponding to a value between the 20^{th} to 25^{th} percentile, or approximately four times a week), would have about a 52% chance of consuming whole grains, based on the recall data for everyone in that group. Using this probability, 365 pseudo-days were generated for this person, with each one having the underlying chance of 52% of consuming a whole grain on that day. Next, for the days that were simulated to be consumption days, an average amount was generated for each consumption day using the mean consumption-day amounts from EATS represented in Figure 1b. For example, for a woman in group 5 (Figure 1b), the mean value of about 0.8 servings, plus or minus a randomly-generated value reflecting the variability about the mean, was used. Finally, two correlated person-specific effects (from a bivariate normal distribution), corresponding to probability and amount, were generated for each pseudo-person with a correlation equivalent to that found in the EATS data.

Figure 1a. Average probability of whole grains consumption by Diet History Questionnaire (DHQ) whole grains frequency group for women in the Eating at America's Table Study.

After combining the simulated probability, consumption-day amount, and person-specific effects, each person in the data set has 365 days of pseudo-data, the mean of which is used to estimate true intake. The average percentiles estimated from all six methods above were compared to the percentiles of this true intake.

Table 1 presents the proportion of respondents in EATS who consumed a food group on none, one, two, three, or four of the recall days, and the corresponding median amount, expressed as the number of Pyramid servings (8), consumed per day. Forty-five of the 54 (83%) food group by sex combinations (27 food groups for each sex) exhibited a positive correlation between the probability of consuming the food and the average consumption-day amount.

The smoothed distribution curve for whole grain consumption by women from the EATS data is represented in Figure 2. This figure illustrates the differences between the WPM, ISUF, and NCI methods using four days of 24HRs. First, even with four days of recall, it is still difficult to estimate empirically the lower tail of the distribution using WPMs, resulting in a large spike at zero. Second, because the WPM method does not distinguish within-person from between-person variability, its distribution has a longer tail. Although the ISUF software produced a warning that it should not be used for whole grains because significant correlation existed between the probability of consuming and the amount consumed, its results are presented here for illustration purposes. The curve for the NCI method is to the left of the ISUF curve below about 1.7 servings, at which point it shifts to the right of the ISUF curve for larger amounts. The difference between these curves is due to the positive correlation between probability and amount. Because those women who are less likely to eat whole grains eat smaller amounts when they do eat them, the area under the NCI method curve is larger in the lower part of the distribution than the area below the curve produced using the ISUF method, which assumes no relationship between probability and amount. For the same reason, this relationship reverses in the upper part of the area under the NCI method curve, reflecting that women who are more likely to eat whole grains eat larger amounts when they do eat them.

Table 2 illustrates the differences among the three methods. For both dark green vegetables and whole grains, the four-day WPM method produces estimates that are higher than the ISUF and NCI methods in both of the tails, due to the large spike at zero and the longer tail in the upper end of the distribution.

Percentage of persons who consume below or above a cutoff number of servings using different methods.

Other differences between the NCI method and the ISUF method are determined by the strength of the correlation between probability of consumption and the amount consumed. This correlation coefficient is not significantly different from zero for dark green vegetable consumption by men in EATS, and, as can be seen in Table 2, the estimates from the NCI and ISUF methods are similar for men. In contrast, due to substantially correlated probability and amount for whole grains (r=0.29 for women and r=0.34 for men), the estimates of the number of people consuming less than 0.5 serving and greater than three servings of this food group are considerably lower in the ISUF method than the NCI method for both women and men.

The percentages of men falling above or below specified cutoff values for servings of whole grains by education levels are given in Table 3. It is clear from this table that the proportion of men in EATS who consume whole grains differs by level of educational attainment.

Figure 3 illustrates the average bias, defined as the difference between the estimate and simulated truth (the 365-day mean), of each of the methods. Except at the mean of the distribution, the two-day WPM estimate has a much greater bias than any of the other methods. The NCI method with correlated person-specific effects either with or without the FPQ produced estimates that are very close to the estimate from the 365-day mean, with essentially no bias for all percentiles.

Bias of percentile estimates from simulations based on whole grains for women (from the Eating at America's Table Study). The dashed line at zero represents no bias. (2 day mean: within-person mean (WPM) of 2 days of simulated 24-hour recalls; ISUF: Iowa **...**

Figure 4 illustrates that the ISUF method and the NCI method with uncorrelated person-specific effects and without the FPQ produced very similar curves that are shifted from the 365-day mean curve. By ignoring correlation between probability and amount, these two methods tend to overestimate the amounts consumed by those with a low probability of consumption, and underestimate the amounts consumed by those with a high probability of consumption, leading to biased estimates. The NCI method with uncorrelated person-specific effects but with the FPQ led to the best results when compared to the same method without the FPQ. Including the shared covariate (the FPQ information) in both parts of the model captured some of the correlation between probability and consumption-day amount, although not as much as incorporating the correlated person-specific effects. Figures 3 and and44 illustrate that modeling the correlation between probability of consumption and consumption-day amount leads to the largest gains over the other methods.

A statistical model used for estimating the usual intake of episodically consumed foods using two or more 24HRs per subject needs to appropriately account for the characteristics of such data. First, it must account for the spike at zero due to non-consumption of a food on the recalled days. This is achieved in the NCI method by representing usual intake as the product of the probability of consuming a food and the amount consumed on a consumption day. Additionally, the model must transform skewed distributions of the consumption-day amounts to approximate normality. It also must have the ability to distinguish within-person variability that results from day-to-day differences in intake and random reporting error from the variability among individuals. Like the ISUF method, the NCI method addresses all of these challenges. The NCI method overcomes the limitations of the ISUF method, however, by incorporating covariate information and by accounting for the correlation between the probability of consuming a food and the consumption-day amount. Because probability and amount are substantially correlated in a majority of food groups, as our findings from EATS illustrate, the ISUF method should not be used on most foods according to its own criterion.

The importance of using covariates in the model is highly dependent upon the application of interest. When interest is on estimating the distribution of usual intake in a subpopulation or population, it is not important how the between-person variation is partitioned between the part that is explained by covariates and the unexplained component captured by the person-specific effect, but how well person-specific effect and within-person random error can be transformed to normality. Including covariates in the model may make normality more realistic, with consequent improvement in estimating the distribution, especially its tails. This was not the case with our example. After allowing for correlated person-specific effects, our simulation study results indicated little difference between the NCI method with and without the FPQ. This may be because the 24HR reported amounts for whole grains (the basis for the simulated data) can be transformed to normality almost equally well unconditionally or conditionally on the FPQ. For foods without this characteristic, including the FPQ in the model may lead to some improvement in estimating the tails of the distribution. Because the FPQ may not be necessary to estimate the usual intake distribution of foods, it should be possible to estimate the usual intake distribution of foods from previous survey data with at least 2 24HRs per participant.

When the distribution of usual intake is estimated in subpopulations, incorporating covariates that characterize the subpopulation, such as age, sex, race, income, or education, may provide a substantial improvement by leading to more efficient estimation than does stratification. This advantage in efficiency is expected to increase as size of the subpopulation decreases. In addition, covariates may be used to adjust for temporal effects, such as seasonality and day-of-week effects, and the reduction in mean levels of intake that can occur with repeat 24HRs, the time-in-sample effect. This is done simply by including an indicator variable into the statistical model as a covariate to indicate that the second recall is being modeled, allowing its mean to be adjusted for repeat application of the recall, if necessary.

For predicting individual usual intake and relating it to health outcomes, the goal is to reduce unexplained between-person variation of intake, which may be achieved by including appropriate covariates in the model. Different types of covariates may be incorporated in the estimation of usual intake for this purpose, including personal characteristics associated with intake, as well as FPQ information. Because the FPQ frequencies are associated both with the probability of consuming a food and the amount consumed on the consumption day (5), incorporating food frequency information as a covariate in both parts of the model can explain at least part of between-person variation in the 24HR, therefore providing a better estimate of diet-health relations. Although the NCI statistical model is the first step in estimating individual usual intake, the development of additional statistical methodology for the second step is necessary; this work is currently underway.

Another benefit of including covariates in the statistical model is the improved ability to make inferences regarding the effects of these covariates on food consumption, similar to the purpose of the statistical models described by Haines et al (9) and Guenther et al. (10). Like these approaches, the NCI model allows the analyst to separate the effect of covariates on the decision to consume a food from the effect of covariates on the amount that is consumed when the food is eaten. The NCI model may be used to approximate the effects of individual covariates to determine which variables are associated with the probability of consumption and the consumption-day amount by making inferences about the strength of the relationships of the covariates with individual intake.

The NCI method was developed to meet the special challenges for estimating intake of episodically-consumed foods. Although some of those challenges (e.g., (A) and (D) listed in the introduction) are unique to episodically-consumed foods, others apply as well to nutrients or foods that are consumed on a daily basis. Part II of the NCI model (equation [2]) alone could be used to estimate usual intake for foods or nutrients that are consumed on a daily basis by nearly everyone. The ability to easily incorporate covariates in the estimation of usual intake makes the NCI method an attractive alternative to the ISU nutrient method.

Although the NCI method appears to present a substantial improvement over existing methods for estimating the distribution of usual intake for foods, it does have some limitations. First, the model never produces a true zero intake because the logistic regression that is used to model the probability of consumption does not predict a zero value. Furthermore, the model requires that a sufficient number of people consume a given food on at least two recalled days. For foods that are consumed very episodically in the population, such as organ meats, this condition may not be satisfied.

Most importantly, the model is based on the assumption that the 24HR is an unbiased instrument for measuring usual food intake. Many recent studies with doubly labeled water have found misreporting of energy intake on both the 24HR and FFQ, almost always in the direction of underreporting (11-14). This suggests that at least some foods are underreported as well. A few studies have investigated the extent and type of underreporting by food, and it appears that underreporting may be differential by food (15); although, unfortunately, it is impossible to know which foods, and by how much, are misreported on the 24HR. Due to this uncertainty, we follow the only available practical convention and assume that the 24HR is unbiased. For those foods that are reported with bias on the 24HR, the estimated intake will be biased as well. When the FFQ is used as a covariate in the model, it is allowed to involve systematic bias as well as random measurement error. It is important to note that, in the NCI method, the FFQ does not replace information from the 24HR but is being calibrated using the 24HR as a reference instrument. The method, therefore, includes useful information from the FFQ without subjecting the final estimates to its measurement error.

The NCI method has been developed and illustrated here using data from what we implicitly assumed to be a random sample of independent individuals, not a complex survey sample. We are presently working on its extension for complex surveys, such as NHANES, so that estimates of the distribution of usual intake in the US population may be made using the NCI method.

The NCI method to estimate usual intake of even episodically consumed foods using two 24-hour recalls represents an advance in dietary assessment. It provides distinct advantages over previously proposed methods by accounting for the correlation between probability of consumption and amount consumed and by incorporating covariate information. The software for implementing the method and the guidance for using it will be publicly available in the future. Researchers interested in estimating the distribution of usual intakes of food groups for a population or subpopulation or the percent of persons above or below a given standard, are advised to work with a statistician to incorporate this method in their analyses. A further application of the NCI method for estimating individual usual intake and relating it to health outcomes is forthcoming.

The authors would like to thank Phillip S. Kott, Joseph D. Goldman, Richard P. Troiano, and Amy Millen for their thoughtful reviews and Anne Brown Rodgers for her expert editing assistance.

**Publisher's Disclaimer: **This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

Janet A. Tooze, Section on Biostatistics, Department of Public Health Sciences, Wake Forest University School of Medicine, Medical Center Boulevard, Winston-Salem, NC 27157, Phone: (336) 716-3833, Fax: (336) 716-6427, Email: ude.cmbufw@ezootj.

Douglas Midthune, National Cancer Institute, 6130 Executive Blvd, EPN 3131, MSC 7354, Bethesda, MD 20892, Phone: (301) 496-7463, Fax: (301)-402-0816, Email: vog.hin.liam@dnuhtdim.

Kevin W. Dodd, National Cancer Institute, 6116 Executive Boulevard, MSC 8317, Bethesda, MD 20892, Phone: (301) 435 - 1834, Fax:, Email: vog.hin.liam@kddod.

Laurence S. Freedman, Gertner Institute for Epidemiology and Health Policy Research, Sheba Medical Center, Tel Hashomer 52161, Israel, Telephone: 972-3-5305390 (W), 972-2-6793472 (H), Fax: 972-3-5349607, Email: li.oc.moctca@fsl.

Susan M. Krebs-Smith, National Cancer Institute, 6130 Executive Blvd, EPN 4005, MSC 7344, Bethesda, MD 20892, Phone:, Fax:, Email: vog.hin.liam@smssberk.

Amy F. Subar, National Cancer Institute, 6130 Executive Blvd, EPN 4005, MSC 7344, Bethesda, MD 20892, Phone: (301) 594-0831, Fax: (301) 435-3710, Email: vog.hin.liam@arabus.

Patricia M. Guenther, US Department of Agriculture, Center for Nutrition Policy and Promotion, Address: 3101 Park Center Dr., Ste. 1034, Phone: (703) 605-0253, Fax: (703) 305-3300, Email: vog.adsu.ppnc@rehtneuG.aicirtaP.

Raymond J. Carroll, Department of Statistics, Texas A&M University, TAMU 3143, College Station, TX 77843-3143, Phone: (979)845-3141, Fax: (979)845-3144, Email: ude.umat.tats@llorrac..

Victor Kipnis, National Cancer Institute, 6130 Executive Blvd, EPN 3131, MSC 7354, Bethesda, MD 20892, Phone: (301) 496-7464, Fax: (301) 402-0816, Email: vog.hin.liam@vsinpik..

1. Dodd K, Guenther PM, Freedman LS, Subar AF, Kipnis V, Midthune D, Tooze JA, Krebs-Smith SM. Statistical methods for estimating usual intake of nutrients and foods: A review of the theory. under review. [PubMed]

2. Nusser SM, Fuller WA, Guenther PM. Estimation of usual dietary intake distributions: adjusting for measurement error and nonnormality in 24-hour food intake data. In: Trewin D, editor. Survey Measurement and Process Quality. New York, NY: Wiley; 1996. pp. 689–709.

3. Tooze JA, Grunwald GK, Jones RH. Analysis of repeated measures data with clumping at zero. Stat Methods Med Res. 2002;11(4):341–355. [PubMed]

4. Box GEP, Cox DR. An analysis of transformations. J R Stat Soc Ser B. 1964;26:211–252.

5. Subar AF, Dodd K, Guenther PM, Kipnis V, Midthune D, McDowell M, Tooze JA, Freedman LS, Krebs-Smith SM. The Food Propensity Questionnaire (FPQ): concept and development. Journal of the American Dietetic Association. in press. [PubMed]

6. Nusser SM, Carriquiry AL, Dodd KW, Fuller WA. A semi-parametric transformation approach to estimating usual nutrient intake distributions. J Am Stat Assoc. 1996;91:1440–1449.

7. Subar AF, Thompson FE, Kipnis V, Midthune D, Hurwitz P, McNutt S, McIntosh A, Rosenfeld S. Comparative validation of the Block, Willett, and National Cancer Institute food frequency questionnaires: the Eating at America's Table Study. Am J Epidemiol. 2001;154:1089–1099. [PubMed]

8. [4/14/2006]. http://www.ba.ars.usda.gov/cnrg/services/foodlink.html.

9. Haines PS, Guilkey DK, Popkin BM. Modeling food consumption decisions as a two-step process. Am J Agric Econ. 1988;70:543–552.

10. Guenther PM, Jensen HH, Batres-Marquez SP, Chen CF. Sociodemographic, knowledge, and attitudinal factors related to meat consumption in the United States. J Am Diet Assoc. 2005;105:1266–1274. [PubMed]

11. Hill RJ, Davies PS. The validity of self-reported energy intake as determined using the doubly labelled water technique. Br J Nutr. 2001;85:415–30. [PubMed]

12. Trabulsi J, Schoeller DA. Evaluation of dietary assessment instruments against doubly labeled water, a biomarker of habitual energy intake. Am J Physiol Endocrinol Metab. 2001;281:E891–9. [PubMed]

13. Macdiarmid J, Blundell J. Assessing dietary intake: who, what and why of under-reporting. Nutr Res Rev. 1998;11:231–253. [PubMed]

14. Subar AF, Kipnis V, Troiano RP, Midthune D, Schoeller DA, Bingham S, Sharbaugh CO, Trabulsi J, Runswick S, Ballard-Barbash R, Sunshine J, Schatzkin A. Using intake biomarkers to evaluate the extent of dietary misreporting in a large sample of adults: the OPEN Study. Am J Epidemiol. 2003;158:1–13. [PubMed]

15. Krebs-Smith SM, Graubard BI, Kahle LL, Subar AF, Cleveland LE, Ballard-Barbash R. Low energy reporters vs others: a comparison of reported food intakes. Eur J Clin Nutr. 2000;54:281–7. [PubMed]

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |