Introduction. The Journal of Nutrition Methodology and Mathematical Modeling

Similar documents
Statistics 11 Lecture 18 Sampling Distributions (Chapter 6-2, 6-3) 1. Definitions again

Statistics Lecture 13 Sampling Distributions (Chapter 18) fe1. Definitions again

Estimation and Confidence Intervals

Objectives. Sampling Distributions. Overview. Learning Objectives. Statistical Inference. Distribution of Sample Mean. Central Limit Theorem

Chapter 21. Recall from previous chapters: Statistical Thinking. Chapter What Is a Confidence Interval? Review: empirical rule

Review for Chapter 9

5.1 Description of characteristics of population Bivariate analysis Stratified analysis

Maximum Likelihood Estimation of Dietary Intake Distributions

Chapter 8 Descriptive Statistics

Caribbean Examinations Council Secondary Education Certificate School Based Assessment Additional Math Project

Supplemental Material can be found at: 9.DC1.html

Measures of Spread: Standard Deviation

Lecture Outline. BIOST 514/517 Biostatistics I / Applied Biostatistics I. Paradigm of Statistics. Inferential Statistic.

5/7/2014. Standard Error. The Sampling Distribution of the Sample Mean. Example: How Much Do Mean Sales Vary From Week to Week?

Ovarian Cancer Survival

STATISTICAL ANALYSIS & ASTHMATIC PATIENTS IN SULAIMANIYAH GOVERNORATE IN THE TUBER-CLOSES CENTER

Should We Care How Long to Publish? Investigating the Correlation between Publishing Delay and Journal Impact Factor 1

CHAPTER 8 ANSWERS. Copyright 2012 Pearson Education, Inc. Publishing as Addison-Wesley

Methodology National Sports Survey SUMMARY

Measuring Dispersion

APPROVAL REQUIRED. By approving this proof you are confirming that the contact information is correct.

Chapter - 8 BLOOD PRESSURE CONTROL AND DYSLIPIDAEMIA IN PATIENTS ON DIALYSIS

How is the President Doing? Sampling Distribution for the Mean. Now we move toward inference. Bush Approval Ratings, Week of July 7, 2003

Appendix C: Concepts in Statistics

Statistical Analysis and Graphing

Reporting Checklist for Nature Neuroscience

Bayesian Sequential Estimation of Proportion of Orthopedic Surgery of Type 2 Diabetic Patients Among Different Age Groups A Case Study of Government

23.3 Sampling Distributions

Sampling Distributions and Confidence Intervals

DISTRIBUTION AND PROPERTIES OF SPERMATOZOA IN DIFFERENT FRACTIONS OF SPLIT EJACULATES*

Concepts Module 7: Comparing Datasets and Comparing a Dataset with a Standard

Basic Requirements. of meeting cow herd production and profitability goals for the beef cattle enterprise.

Practical Basics of Statistical Analysis

Statistics for Managers Using Microsoft Excel Chapter 7 Confidence Interval Estimation

A Supplement to Improved Likelihood Inferences for Weibull Regression Model by Yan Shen and Zhenlin Yang

Technical Assistance Document Algebra I Standard of Learning A.9

The Effect of Question Order on Reporting Physical Activity and Walking Behavior

Drug use in Ireland and Northern Ireland

Guidance on the use of the Title Consultant Psychologist

Children and adults with Attention-Deficit/Hyperactivity Disorder cannot move to the beat

Standard deviation The formula for the best estimate of the population standard deviation from a sample is:

Sec 7.6 Inferences & Conclusions From Data Central Limit Theorem

Chapter 8 Student Lecture Notes 8-1

Nutritional Models. Food-Based Dietary Guidelines Can Be Developed and Tested Using Linear Programming Analysis

EDEXCEL NATIONAL CERTIFICATE UNIT 28 FURTHER MATHEMATICS FOR TECHNICIANS OUTCOME 1- ALGEBRAIC TECHNIQUES TUTORIAL 3 - STATISTICAL TECHNIQUES

What are minimal important changes for asthma measures in a clinical trial?

GSK Medicine Study Number: Title: Rationale: Study Period: Objectives: Primary Secondary Indication: Study Investigators/Centers: Research Methods

Study No.: Title: Rationale: Phase: Study Period: Study Design: Centres: Indication: Treatment: Objectives: Primary Outcome/Efficacy Variable:

GOALS. Describing Data: Numerical Measures. Why a Numeric Approach? Concepts & Goals. Characteristics of the Mean. Graphic of the Arithmetic Mean

Pilot and Exploratory Project Support Grant

Ida Leida M.Thaha, Mega Marindrawati Rochka 1, Muh. Syafar 2

Automatic reasoning evaluation in diet management based on an Italian cookbook

Plantar Pressure Difference: Decision Criteria of Motor Relearning Feedback Insole for Hemiplegic Patients

Biopharmaceutics Classification System (BCS) Biowaiver Assessment Report

Intro to Scientific Analysis (BIO 100) THE t-test. Plant Height (m)

Risk factors for repetition and suicide following self-harm in older adults: multicentre cohort study {

The Efficiency of the Denver Developmental Screening Test with Rural Disadvantaged Preschool Children 1

Dietary glycemic load, added sugars, and carbohydrates as risk factors for pancreatic cancer: the Multiethnic Cohort Study 1 4

Primary: To assess the change on the subject s quality of life between diagnosis and the first 3 months of treatment.

Pilot and Exploratory Project Support Grant

Evaluation of C-14 Based Radiation Doses from Standard Food Ingestion in Korea

GSK Medicine: Study Number: Title: Rationale: Study Period: Objectives: Indication: Study Investigators/Centers: Research Methods:

Repeatability of the Glaucoma Hemifield Test in Automated Perimetry

Vanuatu NCD Risk Factors

Relationship between Established Breast Cancer Risk Factors and Risk of Seven Different Histologic Types of Invasive Breast Cancer

Measurement Variability in Duplex Scan Assessment of Carotid Atherosclerosis

Modified Early Warning Score Effect in the ICU Patient Population

The Nutritional Density Ratio Dilemma: Developing a Scale for Nutritional Value Paul D. Q. Campbell

Hypertension in patients with diabetes is a well recognized

The relationship between hypercholesterolemia as a risk factor for stroke and blood viscosity measured using Digital Microcapillary

Finite Element Simulation of a Doubled Process of Tube Extrusion and Wall Thickness Reduction

Solomon Islands NCD Risk Factors STEPS REPORT

JUST THE MATHS UNIT NUMBER STATISTICS 3 (Measures of dispersion (or scatter)) A.J.Hobson

Measuring the Ability to Identify One s Own Emotions: The Development and Initial Psychometric Evaluation of a Maximum-Performance Test

Clinical Usefulness of Very High and Very Low Levels of C-Reactive Protein Across the Full Range of Framingham Risk Scores

Body Mass Index and Disability Pension in Middle-Aged Men Non-Linear Relations

Randomised controlled trial of a brief alcohol intervention in a general hospital setting

2008 STEPwise Approach to Chronic Disease Risk Factor Survey Report

MONGOLIAN STEPS SURVEY ON THE PREVALENCE OF NONCOMMUNICABLE DISEASE AND INJURY RISK FACTORS

Associations of self-efficacy, social support, and knowledge with fruit and vegetable consumption in Japanese workers

Visual Acuity Screening of Children 6 Months to 3 Years of Age

ETHIOPIA STEPS REPORT ON RISK FACTORS FOR NON-COMMUNICABLE DISEAES AND PREVALENCE OF SELECTED NCDs

1 Barnes D and Lombardo C (2006) A Profile of Older People s Mental Health Services: Report of Service Mapping 2006, Durham University.

A longitudinal study of self-assessment accuracy

Federated States of Micronesia (Chuuk) NCD Risk Factors STEPS REPORT

Objectives. Types of Statistical Inference. Statistical Inference. Chapter 19 Confidence intervals: Estimating with confidence

An algorithm for prioritizing the maintenance of power transformers

Minimum skills required by children to complete healthrelated quality of life instruments for asthma: comparison of measurement properties

IMPAIRED THEOPHYLLINE CLEARANCE IN PATIENTS WITH COR PULMONALE

International Journal of Mathematical Archive-4(3), 2013, Available online through ISSN

Self-Reported Reasons Men Decide Not to Participate in Free Prostate Cancer Screening

Family Mealtimes, Dietary Quality, and Body Mass Index in Children

Research on the effects of aerobics on promoting the psychological development of students based on SPSS statistical analysis

Long-term vitamin C supplement use and prevalence of early age-related lens opacities1

CURRENT ALCOHOL USE IS ASSOCIATED WITH A REDUCED RISK OF HOT FLASHES IN MIDLIFE WOMEN

Autism Awareness Education. April 2018

Kingdom of Tonga NCD Risk Factors STEPS REPORT

S3: Ultrasensitization is Preserved for Transient Stimuli

The incidence of tumors of the digestive tract is high

NIH Public Access Author Manuscript J Low Genit Tract Dis. Author manuscript; available in PMC 2014 April 01.

Transcription:

The Joural of Nutritio Methodology ad Mathematical Modelig The Populatio Distributio of Ratios of Usual Itakes of Dietary Compoets That Are Cosumed Every Day Ca Be Estimated from Repeated 24-Hour Recalls 1 3 Laurece S. Freedma, 4 * Patricia M. Guether, 5 Kevi W. Dodd, 6 Susa M. Krebs-Smith, 6 ad Douglas Midthue 6 4 Gerter Istitute for Epidemiology ad Health Policy Research, Tel Hashomer 52161, Israel; 5 Ceter for Nutritio Policy ad Promotio, USDA, Alexadria, VA 22302; ad 6 Natioal Cacer Istitute, Bethesda, MD 20892 Abstract Estimatig the populatio distributio of the usual itake of a utriet relative to that of aother utriet requires determiatio of idividual-level ratios. If itake data are available o a per-day basis, as with 24-h dietary recalls, those ratios ca be determied i 1 of 2 ways: as the usual ratio of itakes or the ratio of usual itakes. Each of these ratios has its ow meaig ad determiatio; the ratio of usual itakes is coceptually cosistet with determiatios obtaied from FFQ data. We preset a method for estimatig the ratio of usual itakes that uses bivariate modelig of the 2 utriet itakes i questio. Applicatio of the method to the NHANES data for the years 2001 2004 yielded estimated distributios for percet of usual eergy itake from total fat, percet of usual eergy itake from saturated fat, ad usual sodium itake per 1000 kcal (4184 kj) of usual eergy itake. Distributios for both the total populatio ad for age-geder subgroups were estimated. Approximately 60% of adults (.19 y) had a usual total fat itake that was withi the recommeded rage of 20 35% of total eergy, but oly ~34% had a usual saturated fat itake,10% of total eergy. The results chaged oly miimally whe the other defiitio of usual itake, the usual ratio of itakes, was adopted. J. Nutr. 140: 111 116, 2010. Itroductio 1 L.S.F. was supported through a cotract held by Iformatio Maagemet Systems Ic with the US Natioal Cacer Istitute; the remaiig authors were supported by their respective istitutios: USDA ad the Natioal Cacer Istitute. 2 Author disclosures: L. S. Freedma, P. M. Guether, K. W. Dodd, S. M. Krebs-Smith, ad D. Midthue, o coflicts of iterest. 3 A Supplemetal Appedix describig the statistical method used is available with the olie postig of this paper at j.utritio.org. * To whom correspodece should be addressed. E-mail: lsf@actcom.co.il. A major purpose of dietary surveillace or moitorig is to evaluate dietary itakes relative to some stadard. Stadards may be averages, aroud which the populatio s itakes should be distributed, or thresholds, above or below which the populatio s itake should fall, but they are all established with regard to usual itake, geerally defied as the log-ru average daily itake over a period of time (1,2). This is importat because diets vary cosiderably from day to day. Noetheless, the primary assessmet method used i dietary surveillace is the 24-h dietary recall (3), a method that iheretly captures itake oe day at a time ad thus yields a large amout of withi-perso variatio. I previous work o estimatig the proportios of the populatio meetig recommeded dietary itake levels, ivestigators have either preseted the distributios of the reported values (4) or have adjusted the variace of the distributio to exclude withi-perso variatio (5 8). Oly the latter methods successfully remove the withi-perso variatio ad produce distributios of usual itake that have approximately the correct spread or variace. However, these methods have heretofore bee available oly for itakes of sigle dietary compoets that are cosumed early every day by early everyoe [see Dodd et al. (9) for a review]. I dietary surveillace, it is ofte of iterest to examie the populatio distributio of usual itake for a ratio of 2 dietary compoets. Examples iclude the distributio of itakes of saturated fat, total fat, or sodium, each expressed as a ratio of total eergy itake. Estimatig such distributios from atioal survey data that are based upo dietary itake reports from 24-h recalls is ot trivial. It has bee oted previously (10,11) that there are 2 ways of describig the distributio of ratios of itakes. Likewise, 2 defiitios of a idividual s usual itake ratio are available. First, oe could take the log-term mea of the daily ratio of itakes, which we will call the usual ratio (of itakes). Secod, oe could take the ratio of the log-term meas of the daily itakes, which we will call the ratio of usual itakes. A simple (but rather extreme) example shows that these are differet. 0022-3166/08 $8.00 ã 2010 America Society for Nutritio. Mauscript received May 18, 2009. Iitial review completed August 3, 2009. Revisio accepted November 1, 2009. 111 First published olie November 18, 2009; doi:10.3945/j.109.110254.

TABLE 1 Summary of the distributio of percet of usual eergy itake from total fat i the total US populatio ad selected age-geder subgroups, 2001 2004 1,2 Percetile Mea 6 SE 5 10 25 50 75 90 95 ad females % % 1 3 y 1515 32.57 6 0.31 25.49 26.94 29.44 32.38 35.50 38.42 40.24 4 8 y 1701 32.30 6 0.31 25.68 27.03 29.37 32.13 35.06 37.80 39.54 9 13 y 1061 33.27 6 0.23 24.41 26.26 29.40 33.03 36.91 40.57 42.91 14 18 y 1424 32.98 6 0.33 24.74 26.46 29.37 32.76 36.36 39.76 41.84 19 30 y 1100 31.98 6 0.49 23.97 25.69 28.55 31.81 35.27 38.48 40.59 31 50 y 1466 33.41 6 0.42 25.18 26.91 29.82 33.20 36.77 40.18 42.32 51 70 y 1252 35.01 6 0.40 25.98 27.85 31.06 34.76 38.69 42.51 44.89 $71 y 832 34.06 6 0.34 24.62 26.57 29.90 33.78 37.92 41.89 44.41 $19 y 4650 33.55 6 0.24 24.95 26.73 29.76 33.29 37.05 40.72 43.03 9 13 y 1112 33.38 6 0.36 24.94 26.69 29.68 33.19 36.84 40.34 42.54 14 18 y 1362 33.44 6 0.51 25.01 26.74 29.72 33.23 36.93 40.44 42.62 19 30 y 1325 32.53 6 0.43 24.39 26.08 28.95 32.32 35.90 39.23 41.37 31 50 y 1595 34.33 6 0.42 25.69 27.48 30.58 34.11 37.86 41.45 43.66 51 70 y 1284 34.93 6 0.35 25.94 27.78 30.95 34.66 38.65 42.44 44.74 $71 y 860 33.94 6 0.29 24.69 26.63 29.90 33.71 37.66 41.57 43.95 $19 y 5064 34.03 6 0.28 25.28 27.05 30.17 33.79 37.62 41.33 43.63 $1 y 17,889 33.56 6 0.17 25.13 26.85 29.82 33.29 37.01 40.64 42.93 2 Percet of usual eergy itake from fat = 100 3 [usual itake of fat (kj)] / [usual itake of eergy (kj)]. Suppose that a perso cosumes o alterate days 10,000 kj, of which 4000 kj are from fat, ad 6000 kj, of which oly 600 kj are from fat. Usig the usual ratio, the usual percet eergy from fat is the average of 40 ad 10%, i.e. 25%. Usig the ratio of usual itakes, the usual fat itake is 2300 kj, the usual eergy itake is 8000 kj, ad their ratio is 29%. The questio of which defiitio is the more appropriate for purposes of dietary surveillace is ot a easy matter to resolve, TABLE 2 Summary of the distributio of percet of usual eergy itake from saturated fat i the total US populatio ad selected age-geder subgroups, 2001 2004 1,2 Percetile Mea 6 SE 5 10 25 50 75 90 95 ad females % % 1 3 y 1515 12.71 6 0.15 9.05 9.78 11.04 12.57 14.23 15.82 16.85 4 8 y 1701 11.55 6 0.18 8.37 8.99 10.10 11.43 12.88 14.26 15.14 9 13 y 1061 11.78 6 0.12 8.31 9.01 10.23 11.67 13.23 14.68 15.62 14 18 y 1424 11.36 6 0.12 8.16 8.82 9.94 11.26 12.67 14.03 14.88 19 30 y 1100 10.77 6 0.22 7.71 8.34 9.44 10.70 12.04 13.30 14.13 31 50 y 1466 10.89 6 0.16 7.79 8.43 9.52 10.80 12.16 13.48 14.30 51 70 y 1252 11.18 6 0.14 7.83 8.52 9.69 11.07 12.55 13.99 14.86 $71 y 832 11.03 6 0.20 7.53 8.24 9.46 10.91 12.46 13.97 14.90 $19 y 4650 10.95 6 0.09 7.76 8.41 9.54 10.85 12.26 13.62 14.49 9 13 y 1112 11.76 6 0.15 8.22 8.93 10.16 11.64 13.21 14.75 15.74 14 18 y 1362 11.26 6 0.17 7.81 8.50 9.70 11.13 12.69 14.19 15.12 19 30 y 1325 10.79 6 0.17 7.51 8.16 9.30 10.67 12.15 13.56 14.48 31 50 y 1595 11.15 6 0.17 7.70 8.39 9.61 11.03 12.57 14.07 15.00 51 70 y 1284 10.94 6 0.16 7.43 8.13 9.35 10.80 12.40 13.94 14.90 $71 y 860 10.73 6 0.15 7.12 7.84 9.10 10.60 12.19 13.79 14.79 $19 y 5064 10.96 6 0.12 7.50 8.20 9.41 10.84 12.39 13.89 14.85 $1 y 17,889 11.16 6 0.07 7.76 8.45 9.63 11.03 12.55 14.03 14.98 2 Percet of usual eergy itake from fat = 100 3 [usual itake of fat (kj)] / [usual itake of eergy (kj)]. 112 Freedma et al.

because the decisio may deped o biological kowledge that is ofte lackig. The usual ratio is simpler to calculate, but the ratio of usual itakes is coceptually cosistet with the ratio that would be determied from FFQ, because they query logterm itake. For this reaso aloe, it is desirable to have a statistical method for estimatig the ratio of usual itakes ad derivig its distributio i the populatio from 24-h recall data. Therefore, our purposes i this article are to describe the methodology for estimatig the populatio distributio of the ratio of usual itakes ad to apply the method to estimate the distributios of several ratios. A feature of the method is the ability to estimate distributios i subpopulatios. Accordigly, we also preset distributios accordig to selected age ad geder groups. Fially, we compare the results of estimatig the populatio distributio of the ratio of usual itakes with those for the usual ratio ad shows that the differece was small i the examples. Methods Data. The data for this study were obtaied from 17,889 participats i the 2001 2004 NHANES. All recalls obtaied from idividuals aged $1 y were icluded except those deemed ureliable by survey staff (297 idividuals), those reportig breast milk cosumptio (39 idividuals) because breast milk itake was ot quatified, ad 1 idividual who reported fastig o the survey day. Dietary itakes were assessed via computer-assisted, itervieweradmiistered 24-h recalls. I 2001, a sigle recall was requested from all idividuals. I 2003 2004, 2 recalls were requested of all participats. The first dietary iterview was coducted i a mobile examiatio ceter ad the secod was coducted 3 10 d later via telephoe. The NHANES has a complex, multistage, probability desig. For this study, oe of the authors (K.W.D.) created replicate weight sets suitable for the balaced repeated replicatio method of variace estimatio. The replicate weights were based o the dietary survey weights accompayig the NHANES data files, which adjust for survey desig, orespose, day of week, ad ay sequece effect i the secod recall. The NHANES 2001 2004 surveys were approved by the NHANES Istitutioal Review Board uder protocol o. 98 12. Further iformatio about the NHANES is available elsewhere (12). Estimatig the populatio distributio of the ratio of usual itakes. The method, described briefly here, is based o a joit bivariate model for the 2 utriets (the umerator ad the deomiator of the ratio) i questio. A more precise ad techical descriptio is provided i the Supplemetal Appedix. Each reported utriet itake was first mathematically trasformed to approximate ormality usig a Box-Cox (power) trasformatio. The trasformatio was chose to miimize the mea squared error aroud a straight lie fit to a weighted QQ plot usig the samplig survey weights of each participat. The data were the aalyzed usig a bivariate liear mixed effects model, icludig the followig terms for each utriet: a itercept (fixed effect), a idicator for whether the reported day was a weekday or a weeked (fixed effect), a idicator for the sequece umber (first vs. secod) of the report (fixed effect), idicators for age group (fixed effects), a subject-specific term (radom effect), ad a withi-subject error term (radom error). For the aalysis of childre s itakes, a extra covariate for geder (fixed effect) was icluded. For adults, me ad wome were aalyzed separately. The covariate weekday/weeked was icluded to accommodate the differece i itake that ofte occurs betwee the weeked (Friday to Suday) ad the rest of the week. The covariate sequece umber was icluded to accommodate the possibility that participats report less fully o the secod recall tha o the first. The first recall was take to be the ubiased report. The covariate age group, which is a factor with several levels, was icluded to allow estimatio of the distributio withi specific age groups. TABLE 3 Summary of the distributio of usual sodium itake relative to usual eergy itake i the total U.S. populatio ad selected age-geder subgroups, 2001 2004 1,2 Percetile Mea 6 SE 5 10 25 50 75 90 95 ad females mg/1000 kcal 3 mg/1000 kcal 1 3 y 1515 1392 6 10 1041 1112 1234 1380 1537 1685 1779 4 8 y 1701 1509 6 12 1166 1235 1355 1498 1652 1798 1893 9 13 y 1061 1551 6 23 1208 1277 1396 1538 1688 1842 1937 14 18 y 1424 1513 6 19 1198 1261 1372 1502 1644 1781 1865 19 30 y 1100 1505 6 21 1192 1255 1366 1494 1633 1765 1855 31 50 y 1466 1551 6 15 1233 1297 1408 1540 1682 1822 1911 51 70 y 1252 1608 6 14 1261 1330 1452 1594 1750 1903 2003 $71 y 832 1609 6 20 1245 1316 1443 1594 1759 1921 2025 $19 y 4650 1560 6 8 1228 1294 1410 1547 1696 1844 1939 9 13 y 1112 1520 6 18 1148 1220 1348 1505 1673 1838 1946 14 18 y 1362 1475 6 21 1109 1182 1307 1460 1626 1789 1894 19 30 y 1325 1558 6 22 1186 1258 1386 1541 1711 1880 1986 31 50 y 1595 1571 6 20 1186 1263 1396 1554 1728 1901 2010 51 70 y 1284 1624 6 14 1218 1297 1437 1604 1791 1977 2094 $71 y 860 1618 6 22 1195 1278 1423 1599 1790 1982 2107 $19 y 5064 1588 6 12 1195 1272 1407 1569 1749 1928 2043 $1 y 17,889 1554 6 7 1187 1260 1388 1538 1702 1867 1973 2 Usual itake of sodium (mg) per 1000 kcal of usual eergy itake = 1000 3 [usual itake of sodium (mg)] / [usual itake of eergy (kcal)]. 3 1 kcal = 4.184 kj. Distributio of ratios of usual itakes 113

TABLE 4 Proportios i the total U.S. populatio ad selected age-geder subgroups with usual itake of total fat equal to 20 35% of usual eergy itake ad proportios with usual itake of saturated fat less tha 10% of usual eergy itake, 2001 2004 1 Proportio with total fat = 20 35% of eergy 6 SE Proportio with saturated fat,10% of eergy 6 SE ad females 1 3 y 1515 0.72 6 0.03 0.12 6 0.02 4 8 y 1701 0.74 6 0.03 0.23 6 0.03 9 13 y 1061 0.63 6 0.02 0.22 6 0.02 14 18 y 1424 0.66 6 0.02 0.26 6 0.03 19 30 y 1100 0.73 6 0.03 0.36 6 0.05 31 50 y 1466 0.63 6 0.03 0.34 6 0.03 51 70 y 1252 0.52 6 0.03 0.30 6 0.02 $71 y 832 0.58 6 0.02 0.34 6 0.03 $19 y 4650 0.62 6 0.02 0.33 6 0.02 9 13 y 1112 0.63 6 0.02 0.23 6 0.03 14 18 y 1362 0.63 6 0.03 0.30 6 0.03 19 30 y 1325 0.69 6 0.03 0.38 6 0.03 31 50 y 1595 0.56 6 0.03 0.31 6 0.03 51 70 y 1284 0.52 6 0.02 0.36 6 0.03 $71 y 860 0.58 6 0.02 0.40 6 0.03 $19 y 5064 0.58 6 0.02 0.35 6 0.02 $1 y 17,889 0.62 6 0.01 0.31 6 0.02 Icludig a covariate i the model allows estimatio of distributios i subpopulatios ad results i much greater precisio tha aalyzig the model o oly the subset of data from that subpopulatio (13). The subject-specific terms of the 2 utriets were assumed to have a bivariate ormal distributio, with 2 variaces ad a covariace estimated from the data. Likewise, the withi-subject terms were assumed to be bivariate ormal ad idepedet of the subject-specific terms ad required aother 2 variaces ad a covariace to be estimated. The model parameters could be estimated, because may participats completed more tha oe 24-h recall. If there are o repeated determiatios, the withi-subject variace parameters could ot be estimated. I the data from NHANES 2001 2004, 69% of the participats completed 2 24-h recalls. The model was fitted to the data usig the NLMIXED procedure i the statistical software package SAS ad the samplig weights of each participat were icorporated ito the aalysis. The output yielded estimates of the fixed-effect ad radom-effect model parameters, icludig the withi-subject covariaces. Three rus of the NLMIXED procedure for each type of ratio were performed o childre aged 1 8 y, males aged $9 y, ad females $9 y. Mote Carlo simulatios were the ru usig the parameter values estimated from the model. These simulatios geerated usual itakes of each utriet for a large umber of pseudo-idividuals (the details of these simulatios are described i the Supplemetal Appedix). The ratio of usual itakes was calculated for each pseudo-idividual ad, usig the samplig weights, the percetiles of the distributio of ratios of usual itakes were estimated. SE were estimated usig balaced repeated replicatio. These procedures were used to estimate distributios for 10 age-geder subgroups. The distributios for me ad wome $19 y were estimated by combiig the distributios of 4 fier age groupigs accordig to the populatio proportio i each subgroup. Results are show for the followig ratios: total fat:eergy, saturated fat:eergy, ad sodium:eergy. Values i the text are meas or proportios 6 1 SE, or estimated percetiles. The SAS code for performig these aalyses is available at http://riskfactor. cacer.gov/diet/usualitakes/. Estimatig the populatio distributio of the usual ratio. The method for estimatig the populatio distributio of the usual ratio is simpler tha for the ratio of usual itakes, because it is based o the uivariate distributio of the ratio of daily itakes rather tha o the bivariate distributio of the umerator ad deomiator. Several methods ca be applied directly to idividuals 24-h recall reported daily ratios (5,6,13) ad would be expected to give very similar results. For this report, the oe-part versio of a method developed at the Natioal Cacer Istitute (13) was used. Results The estimated mea percetage of usual eergy itake from total fat i the total populatio was 33.6 6 0.2% ad the media 33.3% (Table 1). The 5th ad 95th percetiles were 25.1 ad 42.9%, respectively. There was modest variatio with age, with older persos tedig to have a higher percetage itake ad very little differece betwee me ad wome. The estimated mea percetage of usual eergy itake from saturated fat i the total populatio was 11.2 6 0.1% ad the media 11.0% (Table 2). The 5th ad 95th percetiles were 7.8 ad 15.0%, respectively. The yougest childre (1 3 y) had the highest percetage; the percetage was lower i childre aged 4 18 y ad lowest i adults. 19 y, with little variatio amog adults by age group or betwee me ad wome. The estimated mea usual itake of sodium per 1000 kcal (4184 kj) of usual eergy itake i the total populatio was 1554 6 7 mg ad the media was 1538 mg (Table 3). Itakes were lower amog youg childre (1 3 y). Itakes were also slightly lower amog youger adults (31 50 y) tha older adults (.50 y). Differeces betwee me ad wome were mior. Supplemetal Tables providig SE for the idividual percetile TABLE 5 Compariso of the mea ratio of usual itakes to the mea usual ratio for percet eergy from total fat i the total U.S. populatio ad selected age-geder subgroups, 2001 2004 1 Ratio of usual itakes Usual ratio Differece ad females % 1 3 y 1515 32.57 6 0.31 32.40 6 0.32 20.17 6 0.09 4 8 y 1701 32.30 6 0.31 32.09 6 0.32 20.21 6 0.08 9 13 y 1061 33.27 6 0.23 32.49 6 0.23 20.78 6 0.15 14 18 y 1424 32.98 6 0.33 32.50 6 0.32 20.48 6 0.12 19 30 y 1100 31.98 6 0.49 31.84 6 0.44 20.14 6 0.14 31 50 y 1466 33.41 6 0.42 33.18 6 0.38 20.22 6 0.11 51 70 y 1252 35.01 6 0.40 34.60 6 0.40 20.41 6 0.16 $71 y 832 34.06 6 0.34 33.43 6 0.36 20.64 6 0.14 $19 y 4650 33.55 6 0.24 33.26 6 0.22 20.29 6 0.09 9 13 y 1112 33.38 6 0.36 32.71 6 0.36 20.67 6 0.09 14 18 y 1362 33.44 6 0.51 32.87 6 0.46 20.58 6 0.13 19 30 y 1325 32.53 6 0.43 32.17 6 0.41 20.36 6 0.10 31 50 y 1595 34.33 6 0.42 34.02 6 0.44 20.30 6 0.10 51 70 y 1284 34.93 6 0.35 34.42 6 0.32 20.52 6 0.12 $71 y 860 33.94 6 0.29 33.13 6 0.28 20.81 6 0.13 $19 y 5064 34.03 6 0.28 33.60 6 0.26 20.43 6 0.09 $1 y 17,889 33.56 6 0.17 33.18 6 0.17 20.39 6 0.05 1 Values are meas 6 SE, estimated from data o 17,889 participats i NHANES. 114 Freedma et al.

TABLE 6 Compariso of the mea ratio of usual itakes to the mea usual ratio for percet eergy from saturated fat i the total U.S. populatio ad selected age-geder subgroups, 2001 2004 1 Ratio of usual itakes Usual ratio Differece ad females % 1 3 y 1515 12.71 6 0.15 12.67 6 0.15 20.05 6 0.05 4 8 y 1701 11.55 6 0.18 11.50 6 0.17 20.05 6 0.03 9 13 y 1061 11.78 6 0.12 11.53 6 0.13 20.26 6 0.06 14 18 y 1424 11.36 6 0.12 11.21 6 0.13 20.15 6 0.04 19 30 y 1100 10.77 6 0.22 10.73 6 0.21 20.04 6 0.05 31 50 y 1466 10.89 6 0.16 10.77 6 0.15 20.12 6 0.04 51 70 y 1252 11.18 6 0.14 11.00 6 0.15 20.18 6 0.05 $71 y 832 11.03 6 0.20 10.78 6 0.21 20.25 6 0.06 $19 y 4650 10.95 6 0.09 10.82 6 0.08 20.13 6 0.04 9 13 y 1112 11.76 6 0.15 11.56 6 0.15 20.20 6 0.03 14 18 y 1362 11.26 6 0.17 11.08 6 0.17 20.18 6 0.05 19 30 y 1325 10.79 6 0.17 10.65 6 0.16 20.14 6 0.04 31 50 y 1595 11.15 6 0.17 11.05 6 0.19 20.11 6 0.04 51 70 y 1284 10.94 6 0.16 10.77 6 0.16 20.17 6 0.05 $71 y 860 10.73 6 0.15 10.48 6 0.15 20.24 6 0.05 $19 y 5064 10.96 6 0.12 10.82 6 0.12 20.15 6 0.04 $1 y 17,889 11.16 6 0.07 11.02 6 0.07 20.14 6 0.02 1 Values are meas 6 SE, estimated from data o 17,889 participats i NHANES. For the 3 examples give, differeces betwee me ad wome were mior. Thus, although absolute itakes of these utriets differed substatially betwee me ad wome, after adjustmet for eergy itake, the itakes were very similar. Because there are 2 differet defiitios of usual itake of a ratio, ivestigators may ask which to use ad whether it matters which is used. The followig formula provides the ivestigator with a method of judgig whether the 2 defiitios will lead to similar results. Deote a idividual s ratio o a give day by x/y. If CV y is the withi-perso CV of the deomiator (expressed as a proportio, ot as a percetage) ad CV x is that of the umerator, ad if r xy is the withi-perso correlatio betwee the umerator ad deomiator, the: idividual 0 s usual ratio ffi ð1 þ CV 2 2r y xycv x CV y Þ 3 idividual 0 s ratio of usual itakes: This approximate result is obtaied usig Taylor s expasio ad should provide reasoable accuracy for CV,0.75. It is similar to ad mathematically cosistet with the expressio give by Krebs-Smith et al. (10), although those authors were cocered with ratios at the populatio level rather tha at the idividual level. The closer the bracketed expressio is to 1.0, the closer the usual ratio will be to the ratio of usual itakes ad, cosequetly, the closer the 2 methods of estimatig the populatio distributio of usual itake will be for the ratio of utriets i questio. Data from repeat 24-h recalls i the Eatig at America s Table Study (15) idicate that the withi-perso CV for eergy is ~0.32, for saturated fat is 0.37, ad for sodium 0.43 [Supplemetal Table 1 i (11)]. The withi-perso correlatios with eergy were 0.78 for saturated fat ad 0.71 for estimates for these ratios ad the estimated distributios for other ratios are available. 7 The methodology described also allows estimatio of the proportio of the populatio whose ratio of usual itakes falls withi a chose iterval. A estimated 62 6 1% of the populatio had a usual total fat cosumptio that was withi the acceptable macroutriet distributio rage for total fat (20 35% of eergy) as set by the Istitute of Medicie (14) (Table 4). The proportio meetig this recommedatio was lowest i the 51- to 70-y age group, where oly 52 6 2% did so. Approximately 34% of the adult populatio s usual saturated fat cosumptio was,10% of eergy, as recommeded by the Dietary Guidelies for Americas (2). The proportio of childre meetig that recommedatio was smaller. The differece betwee the estimated distributio of the ratio of usual itakes ad the distributio of the usual ratio of itakes for the same 3 ratios was small (Tables 5 7). The differeces were i both directios but very small compared with the quatity estimated. All differeces were,3% of the estimated value. Discussio The method described i this article has provided for the first time, to our kowledge, estimates of the distributio of ratios of usual itakes of utriets to usual eergy itake for the U. S. populatio ad age-geder subgroups. Ratios give more direct iformatio o the quality of the compositio of the diet cosumed tha is provided by absolute itakes of utriets ad distributios of such ratios have heretofore bee uavailable. TABLE 7 Compariso of the mea ratio of usual itakes to the mea usual ratio for sodium relative to eergy itake i the total U.S. populatio ad selected age-geder subgroups, 2001 2004 1 Ratio of usual itakes Usual ratio Differece ad females mg/1000 kcal 2 1 3 y 1515 1391.7 6 10.4 1393.6 6 11.4 1.9 6 5.3 4 8 y 1701 1509.3 6 12.2 1521.8 6 12.3 12.5 6 4.7 9 13 y 1061 1550.5 6 23.0 1558.5 6 24.0 8.0 6 8.0 14 18 y 1424 1513.3 6 18.5 1525.6 6 18.1 12.3 6 5.3 19 30 y 1100 1504.6 6 21.0 1519.0 6 20.6 14.4 6 4.9 31 50 y 1466 1551.4 6 14.6 1571.9 6 15.1 20.5 6 6.1 51 70 y 1252 1607.9 6 14.1 1624.7 6 15.0 16.8 6 6.0 $71 y 832 1609.2 6 20.1 1617.3 6 21.2 8.1 6 6.9 $19 y 4650 1560.4 6 8.4 1577.4 6 7.5 17.0 6 5.2 9 13 y 1112 1519.7 6 18.4 1539.3 6 18.8 19.6 6 5.3 14 18 y 1362 1475.4 6 20.6 1491.4 6 20.2 16.0 6 6.4 19 30 y 1325 1557.8 6 21.5 1576.9 6 22.9 19.2 6 5.2 31 50 y 1595 1570.8 6 20.0 1590.9 6 19.9 20.1 6 4.7 51 70 y 1284 1624.4 6 14.0 1641.9 6 14.4 17.4 6 5.0 $71 y 860 1618.3 6 21.7 1632.1 6 20.8 13.8 6 7.0 $19 y 5064 1587.8 6 12.0 1606.3 6 12.1 18.4 6 3.8 $1 y 17,889 1553.5 6 7.0 1569.6 6 6.0 16.1 6 2.6 1 Values are meas 6 SE, estimated from data o 17,889 participats i NHANES. 2 1 kcal = 4.184 kj. Distributio of ratios of usual itakes 115

sodium [Supplemetal Table 1 i (11)]. Thus, the bracketed expressio above for the saturated fat:eergy ratio is equal to 1.010 ad for the sodium:eergy ratio is 1.004, values that are close eough to 1 that oe would expect very little differece betwee the methods. That the values are so close to 1 is due to the relatively high correlatio betwee umerator ad deomiator, together with the similarity i their CV. As predicted by the above formula, the results preseted i Tables 6 ad 7 show very little differece betwee the 2 methods for the saturated fat:eergy ratio ad the sodium: eergy ratio. Table 5 shows that the differece was also small for the total fat:eergy ratio. The formula predicts, however, that differeces could be more importat i other situatios. For example, i a case where the deomiator has a high CV ad the umerator has little or o correlatio with the deomiator, oe might expect larger differeces. Thus, ratios of microutriets, such as vitami C, to eergy may display larger differeces. I view of the possibility of larger differeces beig foud betwee the 2 methods, it would seem appropriate to idetify if ad whe oe of the defiitios is more appropriate. As metioed i the Itroductio, this choice is difficult i the absece of relevat biological iformatio. For example, high itakes of saturated fat are cosidered deleterious ad it is recommeded that Americas cosume,10% of eergy from such fatty acids. What is less clear is whether this ratio of 10% of eergy should ot be exceeded o a diural basis or over time. If each day s ratio has some relevace, the the usual ratio may be more appropriate. If, however, the relevat ratio is the proportio of all eergy cosumed over a log period of time derived from saturated fat, the the ratio of usual itakes may be more appropriate. For most ratios of iterest i utritio, biological iformatio to guide choosig the most relevat ratio is ot yet available. It may be that oe or the other method will be foud to be more appropriate for some ratios tha others. Most utritioists we have questioed have stated a preferece for the ratio of usual itakes, perhaps reflectig a ituitio that this term better reflects the steady state of a idividual. As metioed earlier, the ratio of usual itakes has the added advatage of beig coceptually cosistet with determiatios based o food frequecy data. The method described ca be applied to data gathered by repeated 24-h recalls o utriets or other dietary costituets that are cosumed early every day by early everyoe i the populatio. The method is iteded for use o quite large datasets with sample sizes of at least 1000 or more, especially if distributios i populatio subgroups are to be estimated. If oe is iterested i the distributio of a total populatio oly, the a sample of several hudred may suffice. O occasio there is a eed to examie the populatio distributio of the ratio of a episodically cosumed food or utriet (oe that is ot typically cosumed every day) to aother food or eergy. For example, oe may be iterested i cholesterol itake from eggs as a proportio of total cholesterol itake. I such cases, the required statistical modelig becomes more complex. Methods usig a extesio of the 2-part model described by Kipis et al. (16) are ow beig developed to allow the estimatio of these types of distributios. Ackowledgmets L.S.F., P.M.G., ad S.M.K-S. desiged research; D.M. ad K.W.D. performed statistical aalysis; ad L.S.F., P.M.G., ad S.M.K-S. wrote the paper. L.S.F. had primary resposibility for fial cotet. All authors read ad approved the fial mauscript. Literature Cited 1. Istitute of Medicie, Food ad Nutritio Board. Dietary referece itakes: applicatios i dietary assessmet. Washigto, DC: Natioal Academy Press; 2000. 2. U.S. Departmet of Health ad Huma Services ad USDA. Dietary guidelies for Americas 2005. 6th ed. Washigto, DC: US Govermet Pritig Office; 2005. 3. Moshfegh AJ, Raper N, Igwerse I, Clevelad L, Aad J, Goldma J, LaComb R. A improved approach to 24-hour dietary recall methodology. A Nutr Metab. 2001;45 Suppl 1:156. 4. Moshfegh AJ. Food ad utriet itakes by idividuals i the Uited States, by sex ad age, 1994 96. Natiowide Food Surveys Report No. 96 2. Washigto, DC: USDA; 1998. 5. NRC, Subcommittee o Criteria for Dietary Evaluatio, Coordiatig Committee o Evaluatio of Food Cosumptio Surveys, Food ad Nutritio Board. Nutriet adequacy: assessmet usig food cosumptio surveys. Washigto, DC: Natioal Academy Press; 1986. 6. Nusser SM, Carriquiry AL, Dodd KW, Fuller WA. A semiparametric trasformatio approach to estimatig usual daily itake distributios. J Am Stat Assoc. 1996;91:1440 9. 7. Guether PM, Kott PK, Carriquiry AL. Developmet of a approach for estimatig usual utriet itake distributios at the populatio level. J Nutr. 1997;127:1106 12. 8. Moshfegh A, Goldma J, Clevelad L. 2005. What We Eat i America, NHANES 2001 2002: usual utriet itakes from food compared to dietary referece itakes. USDA, Agricultural Research Service [cited 2009 Apr 27]. Available from: http://www.ars.usda.gov/sp2userfiles/ Place/12355000/pdf/usualitaketable s2001 02.pdf. 9. Dodd KW, Guether PM, Freedma LS, Subar AF, Kipis V, Midthue D, Tooze JA, Krebs-Smith SM. Statistical methods for estimatig usual itake of utriets ad foods: a review of the theory. J Am Diet Assoc. 2006;106:1640 50. 10. Krebs-Smith SM, Kott PS, Guether PM. Mea proportio ad populatio proportio: two aswers to the same questio? J Am Diet Assoc. 1989;89:671 6. 11. Freedma LS, Guether PM, Krebs-Smith SM, Kott PS. A populatio s mea Healthy Eatig Idex-2005 scores are best estimated by the populatio ratio whe oe 24-hour recall is available. J Nutr. 2008;138: 1725 9. 12. NHANES. [Website o the Iteret]. [cited 2009 28 Oct]. Available from: http://www.cdc.gov/chs/haes.htm. 13. Tooze JA, Midthue D, Dodd KW, Freedma LS, Krebs-Smith SM, Subar AF, Guether PM, Carroll RJ, Kipis V. A ew statistical method for estimatig the usual itake of episodically cosumed foods with applicatio to their distributio. J Am Diet Assoc. 2006;106:1575 87. 14. Istitute of Medicie, Food ad Nutritio Board. Dietary referece itakes for eergy, carbohydrate, fiber, fat, fatty acids, cholesterol, protei, ad amio acids. Washigto, DC: Natioal Academies Press; 2005. 15. Subar AF, Thompso FE, Kipis V, Midthue D, Hurwitz P, McNutt S, McItosh A, Rosefeld S. Comparative validatio of the Block, Willett, ad Natioal Cacer Istitute food frequecy questioaires: The Eatig at America s Table Study. Am J Epidemiol. 2001;154:1089 99. 16. Kipis V, Midthue D, Buckma DW, Dodd KW, Guether PM, Krebs- Smith SM, Subar AF, Tooze JA, Carroll RJ, et al. Modelig data with excess zeros ad measuremet error: applicatio to evaluatig relatioships betwee episodically cosumed foods ad health outcomes. Biometrics. Epub 2009 Feb 26. 116 Freedma et al.