In most epidemiological studies on diet and disease risk, food-frequency questionnaire (FFQs) are used to assess dietary habits (Hu & Willett, Reference Hu and Willett2002). There is an on-going debate about the accuracy of dietary measurements in large-scale cohort studies, since the observations on the association between diet and cancer were sometimes inconsistent (Kipnis et al. Reference Kipnis, Midthune, Freedman, Bingham, Day, Riboli, Ferrari and Carroll2002; Bingham et al. Reference Bingham, Luben, Welch, Wareham, Khaw and Day2003; Schatzkin & Kipnis, Reference Schatzkin and Kipnis2004). An adequate balance between accuracy and effort has to be achieved.
Due to the long latency period for the development of most chronic diseases, the validity and reproducibility of dietary intake data are important issues in nutritional epidemiology. The validation against standard reference methods is used to appraise the ability to discriminate between subjects with true exposure differences (Willett, Reference Willett1998). The reproducibility describes the consistency of FFQ measurements at different times for the same individual. The European Prospective Investigation into Cancer and Nutrition (EPIC) is a prospective multicentre study to further investigate the association between diet, lifestyle and chronic diseases with emphasis on cancer (Riboli & Kaaks, Reference Riboli and Kaaks1997). At the time of the baseline assessment the problem of measurement error was addressed in the Heidelberg component of EPIC by performing validity and reproducibility studies (Boeing et al. Reference Boeing, Bohlscheid-Thomas, Voss, Schneeweiss and Wahrendorf1997; Bohlscheid-Thomas et al. Reference Bohlscheid-Thomas, Hoting, Boeing and Wahrendorf1997a, Reference Bohlscheid-Thomas, Hoting, Boeing and Wahrendorfb).
Diet–disease associations can be explored for food groups and nutrients often by applying a ranking system (e.g. quantiles). The ability to rank individuals along the distribution of intake and the stability of ranking over time is a prerequisite of FFQs. Measurement error in dietary exposure assessment and instability in individual dietary intake may hide true diet–disease associations.
The aims of the present analyses were (1) to assess the long-term reproducibility of the FFQ used in EPIC-Heidelberg and (2) to appraise the dietary changes in the EPIC-Heidelberg cohort.
Material and methods
Study population
The EPIC-Heidelberg cohort comprises 25 540 subjects, of whom women were aged 35–65 years and men 40–65 years at recruitment. Participants were recruited between 1994 and 1998 from the general population in the Heidelberg region (Boeing et al. Reference Boeing, Korfmann and Bergmann1999a, Reference Boeing, Wahrendorf and Beckerb).
A self-administered FFQ was used to attain data on the habitual diet at baseline. The FFQ comprised questions on the frequency and portion size of 148 food items eaten during the year preceding enrolment (FFQ-1; Bohlscheid-Thomas et al. Reference Bohlscheid-Thomas, Hoting, Boeing and Wahrendorf1997a). The second dietary assessment was performed consecutively along with the date of study entry. After a mean follow-up time of 68·9 (sd 4·1) months in men and 68·7 (sd 4·0) months in women the second follow-up (2001–3) was performed. The repeated dietary measurement was performed with the baseline instrument (FFQ-2), with a supplement of forty-five questions on food preparation methods (Rohrmann & Becker, Reference Rohrmann and Becker2002).
From 25 540 participants with complete FFQ-1, 2505 (10 %) subjects omitted FFQ-2, 1003 (4 %) were lost to follow-up or withdrew, and 570 (2 %) had died. In total 21 462 subjects completed FFQ-2, of whom 729 subjects were excluded based on extreme energy intake ( < 1st or >99th percentile of total energy intake), leaving 20 733 subjects (9530 men and 11 203 women) for the analysis.
At recruitment and during the follow-up, data on lifestyle factors such as smoking habits, physical activity and anthropometric measures were collected (Klipstein-Grobusch et al. Reference Klipstein-Grobusch, Georg and Boeing1997; Bergmann et al. Reference Bergmann, Bussas and Boeing1999).
Variable definition
Dietary intake estimates for food groups were based on the results of the self-administered, quantitative FFQ. Nutrient intake data were calculated based on the food intake data of the Bundesanstalt für gesundheitlichen Verbraucherschutz und Veterinärmedizin (1999). The validity of the questionnaire against 24 h dietary recalls as standard method was determined (Kroke et al. Reference Kroke, Klipstein-Grobusch, Voss, Möseneder, Thielecke, Noack and Boeing1999). Foods were summarized into seventeen food groups based on the EPIC-SOFT classification (Slimani et al. Reference Slimani, Deharveng and Charrondiere1999). Energy adjustment of food and nutrient intake was performed by means of the residual method regression analysis for both dietary assessments (Willett & Stampfer, Reference Willett and Stampfer1986). The quintile cut-points were based on the individual consumption levels of the corresponding percentiles.
Educational level was categorized as low, middle and high. BMI was calculated as weight (kg)/height2 (m2) and classified in < 18·5, 18·5–24·9, 25–29·9, 30–34·9 and ≥ 35 kg/m2. Smoking status was summarized as non-smoker and current smoker at the time of the assessment.
Despite the visual (FFQ-1) and programmed (FFQ-2) data check at data-entry, missing values in the FFQ occurred (Bohlscheid-Thomas et al. Reference Bohlscheid-Thomas, Hoting, Boeing and Wahrendorf1997a). Missing values were replaced according to the following manner in both assessments: (1) food items with no marks were regarded as not being consumed; (2) if a food item was indicated as consumed, but the information on portion size or food frequency was missing, then the most frequently indicated portion size or food frequency given in the entire cohort was used; (3) if a food item was marked as not consumed, but either portion size or food frequency was indicated then this item was set as not consumed (Bohlscheid-Thomas et al. Reference Bohlscheid-Thomas, Hoting, Boeing and Wahrendorf1997a).
Statistical method
Comparative evaluation of absolute and energy-adjusted values of food and nutrient intake by gender was performed. Median, range and percentiles were calculated using the SAS procedure UNIVARIATE. The relative changes between both assessments were calculated and shown as percentage of difference to total intake. The distribution of the food group intake was not normally distributed, thus non-parametric tests were carried out. Spearman correlation coefficients between FFQ-1 and FFQ-2 for food groups and selected nutrients were calculated using the SAS procedure CORR (Masson et al. Reference Masson, McNeill, Tomany, Simpson, Peace, Wei, Grubb and Bolton-Smith2003).
Quintiles were calculated based on gender-specific percentile cut-points. The agreement of individual quintile classification by food groups between baseline and follow-up was evaluated and weighted κ values were calculated by means of the SAS procedure PROC FREQ (Fleiss & Cohen, Reference Fleiss and Cohen1973). The results of the κ statistic indicate excellent agreement for values over 0·80, good agreement for values between 0·61 and 0·80, moderate agreement for values between 0·41 and 0·60, fair agreement for values between 0·21 and 0·40, and poor agreement for values equal or less than 0·20 (Altman, Reference Altman1991). The statistical software package SAS release 9·1 (SAS Institute, Cary, NC, USA) was used for all analyses. All statistical tests and corresponding P values were two-sided, and P < 0·05 was considered statistically significant.
Results
Table 1 shows the cohort characteristics by gender of the 20 733 subjects included in the analyses. During the follow-up of the cohort, the age distribution of the cohort shifted to older age groups and the rate of retired persons increased from 9·8 to 20·3 % in women and from 16·5 to 36·5 % in men.
The absolute and the energy-adjusted distribution of daily intake for food groups and selected nutrients for men are given in Table 2 and for women in Table 3. In both genders, the distribution of the BMI did not change substantially. The increase in median daily consumption of non-alcoholic beverages in men (+11 % relative increase) and in women (+14 % relative increase) was most striking. The median consumption of alcoholic beverages declined in men by 1 % and in women by 2 %. Among foods, the median daily intake of fruits (+6 % in men, +6 % in women), vegetables (+3 % in men, +4 % in women) and cereals/cereal products (+5 % in men, +5 % in women) was higher in the FFQ-2 than in the baseline data for both genders. In contrast, the daily median consumption of potatoes ( − 5 % in men, − 3 % in women), cakes ( − 9 % in men, − 14 % in women), added fats and oils ( − 6 % in men, − 11 % in women) and sugar/confectionary ( − 7 % in men, − 8 % in women) was lower in FFQ-2 than in the FFQ-1. The median intake of dairy products was lower in FFQ-2 than FFQ-1 ( − 2 %) in women only. The observed shifts in consumption were more prominent in men than in women for intake of potatoes, cereals/cereal products and alcoholic beverages. A marginal reduction in the median meat/meat products consumption ( − 2 %) was observed only in men, whereas the shift for vegetable consumption was more pronounced in women. The median daily energy intake declined by 1 % in both men and women. Among the macronutrients the reduction in the daily fat intake was most prominent ( − 3 % in men, − 5 % in women) whereas among the micronutrients, the average intakes of β-carotene (+2 % in men, +4 % in women), vitamin C (+6 % in men, +7 % in women) and calcium (+3 % in men, +3 % in women) were slightly higher in the second assessment.
* All correlations reached significance level P < 0·0001.
* All correlations reached significance level P < 0·0001.
Analyses by 10-year age-strata revealed no clear differences for most dietary factors, except for the meat intake among men, which increased in younger age groups by 1 %, but decreased by 4 % among men in the age of 50 years or more, and the consumption of sugar and confectionery, for which a stronger decrease ( − 11 %) in men aged 40–49 years than in older men ( − 5 % and − 1 %) was observed (data not shown). Concerning the consumption of sugar and confectionery also stronger reduction was observed among younger (age 35–39 years) women ( − 13 %) than in women of older age ( − 9 %, − 3 % and − 9 %).
Spearman's correlation coefficients between FFQ-1 and FFQ-2 ranged for food groups and selected nutrients from 0·41 (vegetables) to 0·73 (alcoholic beverages consumption) in men (Table 2). The values in women ranged from 0·44 (cereals/cereal products and fat) to 0·77 (alcohol beverages consumption) (Table 3). In both genders, the correlations reached statistical significance. After adjustment for energy intake, most of the correlation coefficients remained similar or were slightly attenuated. An exception was total fat intake, with an attenuation of the correlation coefficient by energy adjustment from 0·53 to 0·43 in men and from 0·53 to 0·41 in women. Overall, the long-term reproducibility of the FFQ appears best for alcoholic beverages and moderate to good for food groups and selected nutrients. Notably, after energy adjustment the correlation coefficients for the macronutrients protein (in both genders) and carbohydrates (in women) did not improve.
In order to appraise the test–retest reproducibility, the ranking of individuals according to quintile cut-points in both assessment periods was compared (Table 4). The proportions of subjects classified into the comparable quintile ranged from 31·2 % for vegetables to 50·3 % for the consumption of alcoholic beverages in men and from 32·1 % for added fats and oils to 50·3 % for alcoholic beverages in women. Concerning nutrients, reclassification to the same quintile was lowest for vitamin C (31·7 % in men, 32·4 % in women) and highest for ethanol intake (50·1 % in men, 50·5 % in women). The percentages for the re-classification into an opposite quintile ranged from 1·0 % (alcoholic beverages) to 3·2 % (miscellaneous) in men and from 0·7 % (alcoholic beverages) to 2·6 % (vegetables, added fats and oils) in women. Overall, 60–70 % of the participants in both genders were re-classified to the same or adjacent quintile, indicating that the agreement for all food groups and selected nutrients is high.
In order to evaluate the test–retest agreement, the weighted κ statistic was calculated. The lowest weighted κ value with 0·28 was observed for vegetable intake in men and 0·29 for cereals/cereal products intake in women. The highest weighted κ values in both genders (0·57 in men and 0·60 in women) were calculated for alcoholic beverage or ethanol intake, respectively. For most food groups or nutrients, the weighted κ ranged around 0·35 in both genders, indicating fair concordance. Moderate agreement was observed for the intake of meat/meat products (0·44), sugar/confectionery (0·42) and cakes (0·41) in women and for the intake of dairy products (0·44), sugar/confectionery (0·43) and cakes (0·42) in men. Overall, the present results indicate fair test–retest agreement for most food groups and the selected nutrients.
Discussion
In both genders, the median food and nutrient intake levels remained relatively stable over time. An increase in consumption was found for non-alcoholic beverages, cereals/cereal products, fruits and vegetables, which was more pronounced in women than in men. The median daily consumption of potatoes, dairy products (in women), alcohol and meat (in men) was lower according to the FFQ-2 data. To some extent the observed dietary changes may be attributed to real changes of dietary habits. On the ecological level, the observed trends are consistent with other reports on dietary trends based on food consumption data (Mensink, Reference Mensink2002; Deutsche Gesellschaft für Ernährung, 2004; Jungjohann et al. Reference Jungjohann, Luhrmann, Bender, Blettner and Neuhauser-Berthold2005), indicating a slight increase in fruit, vegetable, fish and cereals/cereal products consumption, and a decrease in the consumption of potatoes, meat/meat products and alcoholic beverages between 1995 and 2002 in western Germany (Deutsche Gesellschaft für Ernährung, 2004). The trend for increasing dairy product intake in Germany was in contrast to the present observations in women. The present result may be influenced by different selection patterns of dairy foods such as cheese or yoghurt. Data on dietary consumption of middle-aged men in 1984–5 and 1994–5 from the MONICA-Project in Augsburg are available (Döring et al. Reference Döring, Honig-Blum, Winkler, Kammerlohr, Fischer and Keil1998). In the MONICA data, the mean intake of fat and alcohol were lower in the second survey and the mean carbohydrate intake was higher than in the first survey (Döring et al. Reference Döring, Honig-Blum, Winkler, Kammerlohr, Fischer and Keil1998). The present observations on the intake of carbohydrates, protein, calcium, β-carotene and vitamin C are in line with the national trends for these nutrients (Deutsche Gesellschaft für Ernährung, 2004). Regarding total fat intake, the present data suggest a reduction of the median daily intake, which is also observed in some other reports (Döring et al. Reference Döring, Honig-Blum, Winkler, Kammerlohr, Fischer and Keil1998; Mensink, Reference Mensink2002; Deutsche Gesellschaft für Ernährung, 2004; Jungjohann et al. Reference Jungjohann, Luhrmann, Bender, Blettner and Neuhauser-Berthold2005). The increased consumption of non-alcoholic beverages may reflect a trend to more health-conscious choices (Jungjohann et al. Reference Jungjohann, Luhrmann, Bender, Blettner and Neuhauser-Berthold2005; Roddam et al. Reference Roddam, Spencer, Banks, Beral, Reeves, Appleby, Barnes, Whiteman and Key2005) but may also be affected by the past hot summers. Changes in dietary habits depend on the age group, sociodemographic factors and the cultural setting (Colditz et al. Reference Colditz, Willett, Stampfer, Sampson, Rosner, Hennekens and Speizer1987). This limits a direct comparison with the results from other surveys. Looking at the changes in median intakes it has to be considered that if the measurement error in both measurements was the same, then the effect of the flattened slope would underestimate the true changes in median intakes (Kipnis et al. Reference Kipnis, Subar, Midthune, Freedman, Ballard-Barbash, Troiano, Bingham, Schoeller, Schatzkin and Carroll2003).
In the present data, correlation coefficients for food and nutrient intakes ranged from 0·5 to 0·7, indicating a moderate to high correlation for the food groups and nutrients. This is in line with other studies on long-term reproducibility measured for periods from 1 to 10 years (Pietinen et al. Reference Pietinen, Hartman, Haapa, Rasanen, Haapakoski, Palmgren, Albanes, Virtamo and Huttunen1988; Jarvinen et al. Reference Jarvinen, Seppanen and Knekt1993; Goldbohm et al. Reference Goldbohm, van't Veer, van den Brandt, van't Hof, Brants, Sturmans and Hermus1995; Willett, Reference Willett1998; Masson et al. Reference Masson, McNeill, Tomany, Simpson, Peace, Wei, Grubb and Bolton-Smith2003). In a former reproducibility study over 1 year, using a precursor of the current FFQ, the correlation coefficients of food group intake ranged from 0·49 for bread to 0·89 for alcohol with a median coefficient of 0·70 (Bohlscheid-Thomas et al. Reference Bohlscheid-Thomas, Hoting, Boeing and Wahrendorf1997b). Goldbohm et al. (Reference Goldbohm, van't Veer, van den Brandt, van't Hof, Brants, Sturmans and Hermus1995) collected FFQ data annually over 5 years. The results of the regression analysis indicated a slight decline in the median consumption of most nutrients and food groups over time. Concerning macronutrients, the correlation coefficients for the unadjusted FFQ data were mostly between 0·60 and 0·70 (Bohlscheid-Thomas et al. Reference Bohlscheid-Thomas, Hoting, Boeing and Wahrendorf1997a). Low correlation can be caused by the method of dietary measurement or real changes of dietary habits. In the validation study, the FFQ-derived measurements were compared with data collected with the 24 h diet recall method used twelve times over 1 year as standard measurement. The correlation coefficients for food group intake between FFQ and 24 h recall data had a median of 0·45 (range 0·19 for legumes and 0·90 for alcoholic beverages) (Bohlscheid-Thomas et al. Reference Bohlscheid-Thomas, Hoting, Boeing and Wahrendorf1997a, Reference Bohlscheid-Thomas, Hoting, Boeing and Wahrendorfb). The ratio of within- to between-person variations was larger than 1 for almost all food groups, however, the authors noted a lack of precision of the reference measurements (Bohlscheid-Thomas et al. Reference Bohlscheid-Thomas, Hoting, Boeing and Wahrendorf1997a). Kipnis et al. (Reference Kipnis, Subar, Midthune, Freedman, Ballard-Barbash, Troiano, Bingham, Schoeller, Schatzkin and Carroll2003) proposed as a structure of measurement error besides within-person random error, two biases (intake-related and person-specific). The intake-related bias may lead to the flattening slope phenomenon which was found to be more pronounced in energy-adjusted data. The person-specific bias reflects the difference between total within-person bias and its intake-related component (Kipnis et al. Reference Kipnis, Subar, Midthune, Freedman, Ballard-Barbash, Troiano, Bingham, Schoeller, Schatzkin and Carroll2003). Measurement error of dietary intake often attenuates the disease-related risk estimates and may consequently miss diet–disease associations. The development of refined biological markers such as doubly labelled water for the assessment for the energy expenditure may help to be better able to distinguish between real changes of diet and measurement errors (Kaaks et al. Reference Kaaks, Ferrari, Ciampi, Plummer and Riboli2002; Subar et al. Reference Subar, Kipnis and Troiano2003).
On the food group level, correlation coefficients remained similar or were slightly attenuated after adjustment for energy intake. Differing from the observation in food groups and micronutrients, in macronutrients the correlation coefficients for total fat and protein intake in men and fat, protein and carbohydrates intake in women declined after energy adjustment. Adjustment for energy was suggested to correct partially for measurement error (Willett, Reference Willett1998). Energy adjustment had the strongest impact on the energy-contributing macronutrients such as total fat intake in both genders, whereas in women the correlation coefficient for carbohydrates was worse in the unadjusted data (0·51 and 0·49). The present observation suggests that the macronutrients were differently prone to reporting errors (Voss et al. Reference Voss, Kroke, Klipstein-Grobusch and Boeing1998).
The agreement between the same or adjacent ranks of subjects according to food or nutrient intake levels reached generally about 70 %. Extreme misclassification (opposite quintiles) occurred at maximum in 3 % of the participants. The present results indicate a favourable performance of the FFQ and give confidence in the ranking of subjects according to dietary exposure. Comparing the agreement of classification according to food group intake levels obtained by 24 h recalls and FFQ, Bohlscheid-Thomas (et al. Reference Bohlscheid-Thomas, Hoting, Boeing and Wahrendorf1997a) found values in a similar range as the present values. Self-selection of subjects whose motivation is above average is likely for the whole cohort. However, the participation rate during follow-up was high, indicating that no further selection process during follow-up may have happened.
When reporting on changes of food group and nutrient intake, several limitations need to be considered. Though identical dietary questionnaires were used, some time-dependent differences in over- or underreporting related to the conceptions of healthy diet may have occurred. In the baseline assessment this systematic error was explored by comparison with a standard method. About 8 % of the initial study sample provided additionally 24 h dietary recall data for calibration purpose of the EPIC subcohorts (Slimani et al. Reference Slimani, Deharveng and Charrondiere1999), which were analysed on the food group level.
Measuring the long-term reproducibility may be influenced by the decreasing ability to complete the FFQ due to increasing age. However, the strength of the present study is the high participation rate, the application of an identical FFQ, and the calculation of food groups and nutrients according to an identical procedure.
In conclusion, the present results suggest a moderate to good long-term reproducibility of the FFQ used in the EPIC-Heidelberg cohort. After a mean follow-up time of 68·8 months, we observed indications for some modifications in the consumption level of food groups and selected nutrients. The observed deviations may reflect true changes in dietary behaviour, but may also be related to measurement error. For several nutrients the dietary changes were consistent with the results of independent national surveys over time, suggesting that the FFQ results may partly reflect real dietary changes over time. In addition, the observed agreement of classification according to food and nutrient intake levels after several years gives confidence that the dietary data can be used as long-term exposure variables to investigate diet and disease associations. Future analyses on diet and disease risk can additionally benefit from more accurate classification of long-term habitual diet using individual exposure levels based on two dietary measurements, especially for those participants with instable FFQ results.
Acknowledgements
The authors thank the volunteers who participated in EPIC-Heidelberg. EPIC-Heidelberg was supported by the German Cancer Research Centre, grants from the European Community (SI2.326938) and Deutsche Krebshilfe e.V. (70-2488).