Growth in early yyears: statistical and clinical insights

Similar documents
STA Module 6 The Normal Distribution

STA Module 6 The Normal Distribution. Learning Objectives. Examples of Normal Curves

The multivariate piecewise linear growth model for ZHeight and zbmi can be expressed as:

SENIOR VCAL NUMERACY INVESTIGATION SENIOR VCAL NUMERACY INVESTIGATION Only A Little Bit Over. Name:

wine 1 wine 2 wine 3 person person person person person

Dietary Diversity in Urban and Rural China: An Endogenous Variety Approach

BORDEAUX WINE VINTAGE QUALITY AND THE WEATHER ECONOMETRIC ANALYSIS

Activity 10. Coffee Break. Introduction. Equipment Required. Collecting the Data

CHAPTER 2 ANNUAL RETAIL FOOD PRICE MOVEMENTS

Using Growing Degree Hours Accumulated Thirty Days after Bloom to Help Growers Predict Difficult Fruit Sizing Years

Predicting Wine Quality

Final Exam Financial Data Analysis (6 Credit points/imp Students) March 2, 2006

From VOC to IPA: This Beer s For You!

IT 403 Project Beer Advocate Analysis

Lollapalooza Did Not Attend (n = 800) Attended (n = 438)

ICC September 2018 Original: English. Emerging coffee markets: South and East Asia

Business Statistics /82 Spring 2011 Booth School of Business The University of Chicago Final Exam

Chapter 3. Labor Productivity and Comparative Advantage: The Ricardian Model. Pearson Education Limited All rights reserved.

Emerging Local Food Systems in the Caribbean and Southern USA July 6, 2014

Comparing R print-outs from LM, GLM, LMM and GLMM

Online Appendix to The Effect of Liquidity on Governance

Panel A: Treated firm matched to one control firm. t + 1 t + 2 t + 3 Total CFO Compensation 5.03% 0.84% 10.27% [0.384] [0.892] [0.

1.3 Box & Whisker Plots

Chapter 3. Labor Productivity and Comparative Advantage: The Ricardian Model

Preview. Chapter 3. Labor Productivity and Comparative Advantage: The Ricardian Model

Wine-Tasting by Numbers: Using Binary Logistic Regression to Reveal the Preferences of Experts

The Financing and Growth of Firms in China and India: Evidence from Capital Markets

REPORT OF THE WORKING GROUP ON TEA BREW BY T C CHAUDHURI N MURALEEDHARAN ANOOP KUMAR BAROOAH

Transportation demand management in a deprived territory: A case study in the North of France

Analysis of Things (AoT)

EXECUTIVE SUMMARY. 1. When do Asian clams reproduce in Lake George? 2. How fast do Asian clams grow in Lake George?

DOES BEER PLAY A SOLE ROLE IN ALCOHOL AND HEALTH SYMPHONY?

Biologist at Work! Experiment: Width across knuckles of: left hand. cm... right hand. cm. Analysis: Decision: /13 cm. Name

Fairtrade Buying Behaviour: We Know What They Think, But Do We Know What They Do?

*p <.05. **p <.01. ***p <.001.

Haddock. Seafood Industry Factsheet. Market overview: haddock

Napa County Planning Commission Board Agenda Letter

Handling Missing Data. Ashley Parker EDU 7312

Gail E. Potter, Timo Smieszek, and Kerstin Sailer. April 24, 2015

November K. J. Martijn Cremers Lubomir P. Litov Simone M. Sepe

STAT 5302 Applied Regression Analysis. Hawkins

Gasoline Empirical Analysis: Competition Bureau March 2005

HW 5 SOLUTIONS Inference for Two Population Means

Preview. Introduction. Chapter 3. Labor Productivity and Comparative Advantage: The Ricardian Model

Multiple Imputation for Missing Data in KLoSA

Pitfalls for the Construction of a Welfare Indicator: An Experimental Analysis of the Better Life Index

A latent class approach for estimating energy demands and efficiency in transport:

Imputation of multivariate continuous data with non-ignorable missingness

Missing value imputation in SAS: an intro to Proc MI and MIANALYZE

Can You Tell the Difference? A Study on the Preference of Bottled Water. [Anonymous Name 1], [Anonymous Name 2]

Which of the following are resistant statistical measures? 1. Mean 2. Median 3. Mode 4. Range 5. Standard Deviation

Structural Reforms and Agricultural Export Performance An Empirical Analysis

Which of your fingernails comes closest to 1 cm in width? What is the length between your thumb tip and extended index finger tip? If no, why not?

Problem Set #3 Key. Forecasting

Preview. Introduction (cont.) Introduction. Comparative Advantage and Opportunity Cost (cont.) Comparative Advantage and Opportunity Cost

PEEL RIVER HEALTH ASSESSMENT

Lack of Credibility, Inflation Persistence and Disinflation in Colombia

Experimental characterization of short channel bulk MOSFET devices with different gate LER

Preview. Introduction. Chapter 3. Labor Productivity and Comparative Advantage: The Ricardian Model

EFFECT OF TOMATO GENETIC VARIATION ON LYE PEELING EFFICACY TOMATO SOLUTIONS JIM AND ADAM DICK SUMMARY

Recent U.S. Trade Patterns (2000-9) PP542. World Trade 1929 versus U.S. Top Trading Partners (Nov 2009) Why Do Countries Trade?

WINE GRAPE TRIAL REPORT

Grape Growers of Ontario Developing key measures to critically look at the grape and wine industry

Chapter 3: Labor Productivity and Comparative Advantage: The Ricardian Model

Structures of Life. Investigation 1: Origin of Seeds. Big Question: 3 rd Science Notebook. Name:

The Sources of Risk Spillovers among REITs: Asset Similarities and Regional Proximity

Update to A Comprehensive Look at the Empirical Performance of Equity Premium Prediction

Specialty Coffee Market Research 2013

ARE THERE SKILLS PAYOFFS IN LOW AND MIDDLE-INCOME COUNTRIES?

OF THE VARIOUS DECIDUOUS and

An application of cumulative prospect theory to travel time variability

Online Appendix to Voluntary Disclosure and Information Asymmetry: Evidence from the 2005 Securities Offering Reform

The Elasticity of Substitution between Land and Capital: Evidence from Chicago, Berlin, and Pittsburgh

Price & Time Symmetry DENNIS W. WILBORN, SR.

Application & Method. doughlab. Torque. 10 min. Time. Dough Rheometer with Variable Temperature & Mixing Energy. Standard Method: AACCI

FOOD SAFETY RISK ASSESSMENT FOR CCS CANTEENS AND KITCHENS. General purpose Catering kitchen, Hot & Cold Drinks. Created on 22/11/2010

RESULTS OF THE MARKETING SURVEY ON DRINKING BEER

Internet Appendix. For. Birds of a feather: Value implications of political alignment between top management and directors

EAT TOGETHER EAT BETTER BEAN MEASURING ACTIVITY

Credit Supply and Monetary Policy: Identifying the Bank Balance-Sheet Channel with Loan Applications. Web Appendix

To: Professor Roger Bohn & Hyeonsu Kang Subject: Big Data, Assignment April 13th. From: xxxx (anonymized) Date: 4/11/2016

Duration of resource: 17 Minutes. Year of Production: Stock code: VEA12062

Food Allergies on the Rise in American Children

Lesson 23: Newton s Law of Cooling

INFLUENCE OF ENVIRONMENT - Wine evaporation from barrels By Richard M. Blazer, Enologist Sterling Vineyards Calistoga, CA

Evaluation of the Weltech PW-2050 dry matter assessment system

The People of Perth Past, Present and Future

Labor Supply of Married Couples in the Formal and Informal Sectors in Thailand

DRY BEANS PRICES AND INCOME DISTRIBUTION IN RWANDA: A NON- PARAMETRIC ANALYSIS

Decision making with incomplete information Some new developments. Rudolf Vetschera University of Vienna. Tamkang University May 15, 2017

Missouri State University

1/17/manufacturing-jobs-used-to-pay-really-well-notanymore-e/

Rail Haverhill Viability Study

This appendix tabulates results summarized in Section IV of our paper, and also reports the results of additional tests.

COMPARISON OF EMPLOYMENT PROBLEMS OF URBANIZATION IN DISTRICT HEADQUARTERS OF HYDERABAD KARNATAKA REGION A CROSS SECTIONAL STUDY

Determining the Optimum Time to Pick Gwen

Missing Data Treatments

Religion and Life - Year 8 ISBL

A Comparison of X, Y, and Boomer Generation Wine Consumers in California

Table A.1: Use of funds by frequency of ROSCA meetings in 9 research sites (Note multiple answers are allowed per respondent)

Transcription:

Growth in early yyears: statistical and clinical insights Tim Cole Population, Policy and Practice Programme UCL Great Ormond Street Institute of Child Health London WC1N 1EH UK

Child growth Growth is the unique paediatric indicator of well-being, which can monitor a child for endocrine, nutritional, emotional and physical health Ian Jefferson

How to assess growth Measure child Plot on growth chart Read growth chart Take action

Centiles and centile crossing high centile low centile centile tracking 99.6 98 91 75 50 25 9 2 0.4 centiles UK-WHO chart

What is a centile? Centile - percentage point of the frequency distribution Cut-off identifies that percentage of children with measurements below it Examples 50% of children lie below 50th centile (median) 25% of children lie below 25th centile 99.6% of children lie below 99.6th centile (0.4% above) The cut-offs vary by age, so the centiles appear as curves on the growth chart

Growth distance and velocity Growth chart is road to health Current size (i.e. centile) indicates distance travelled Centile crossing indicates velocity of travel Growth chart quantifies size/distance Centile Growth chart does not quantify growth velocity Centile crossing is uncalibrated Ironic - growth chart does not measure growth

A new concept: growth acceleration Growth distance One measurement Centile Growth velocity Two measurements Centile crossing Growth acceleration Three measurements Change in centile crossing

Distance - one measurement

Velocity two measurements

Acceleration three measurements

Growth pattern many measurements

Growth pattern many measurements Modelling growth curves with SITAR Provides simple summary of individual growth patterns

Aims To show how statistics helps chart assessment for: 1. Growth distance 2. Growth velocity 3. Growth acceleration 4. Growth pattern

Growth distance One measurement

Constructing growth charts Growth charts designed to assess single measurements Compare measurement to distribution of reference measurements for age and sex LMS method popular way to construct growth charts Worked example for weight in girls Cole TJ, Green PJ. Smoothing reference centile curves: the LMS method and penalized likelihood. Stat Med 1992;11:1305-19.

Constructing growth charts Weight in 4000 girls Age 1-21 years Aim: to define weight distribution at each age 95% below 95 th centile 50% below 50 th centile 5% below 5 th centile etc

Constructing growth charts Weight in 4000 girls Age 1-21 years Aim: to define weight distribution at each age Construct smooth centile curves 50% below 50 th centile 3% below 3 rd centile 97% below 97 th centile etc

LMS method Cole, JRSS A (1988) Split into narrow age groups Summarise distribution in each group Need to adjust for skewness Raise weight to Box-Cox power λ Calculate mean μ and coefficient of variation σ So λ μ and σ vary by age

LMS method Plot λ μ and σ against age and fit smooth curves L curve for Box-Cox power λ M curve for median μ S curve for coeff of variation σ Hence LMS method

LMS method Centile curves are functions of L M and S curves So if L M and S curves are smooth, centiles are too

Cole-Green LMS method Peter Green (1988) proposed using maximum penalized likelihood to improve LMS method Elegantly avoids arbitrary age groupings See Cole & Green, Stat Med (1992) Now the standard method Peter Green FRS

40 countries use LMS method WHO growth standard

Growth velocity Two measurements

Growth velocity Velocity appears as centile crossing on chart Two problems with chart centiles They assess distance not velocity Light babies grow faster, heavy babies slower Regression to the mean So velocity depends on starting weight Only experience can tell if centile crossing is abnormal Need a way to flag abnormal centile crossing on chart

Velocity and centile crossing Show line on chart whose slope corresponds to 5 th velocity centile over 4 weeks Depends on age and initial centile

Centile crossing over 4 weeks 95 5 95 50 5 95 5

Statistics of centile crossing Two weights 4 weeks apart Convert to z-scores z 1 and z 2 Expected mean of z 2 is r.z 1 where r is correlation between z-scores SD of z 2 is 1-r 2 So 5 th centile for z 2 is z 2 = Mean 1.64 SD = r.z 1 1.64 1-r 2 So z 2 depends on z 1 and r

Thrive lines for growth velocity For ages 0-4, 4-8, 8-12 weeks calculate correlations r 0-4 4, r 4-8 8, r 8-12 etc Choose baseline value z 0 Then using formula z 0 > z 4 > z 8 > z 12 defines a curve Call the curve a thrive line as it defines failure to thrive Cole TJ. Presenting information on growth distance and conditional velocity in one chart: practical issues of chart design. Stat Med 1998;17:2697-707.

Thrive line overlay - 5th centile weight gain Boys weight 5th centile for 4-week intervals

Thrive lines Thrive lines assess weight velocity 5th velocity centile Over a 4-week period A child s plot that tracks along the thrive lines for 4 weeks is growing on the 5th velocity centile Tracking for longer is worse: e.g. for 8 weeks, growth < 1st velocity centile Thrive lines presented as plastic overlay to superimpose on chart

Mild centile crossing - 1 channel width over 8 weeks

Weight gain above 5th centile

Weight gain above 5th centile

Moderate centile crossing - 2 channel widths over 8 weeks

Weight gain below 5th centile

Weight gain below 5th centile

Thrive 95 lines Rapid infant weight gain also a concern Useful to identify rapid weight gain Hence Thrive 95 lines Define 95th centile for weight gain

Thrive 95 lines

Thrive 5 and Thrive 95 lines 5th centile weight gain 95th centile weight gain

Benefit of thrive lines Plastic overlay designed to fit on British 1990 chart format Distance and velocity both assessed yet data plotted just once No need for separate distance and velocity charts Useful addition to weight chart

Thrive lines and electronic charts Now easy to add thrive lines to electronic charts Thrive lines can be drawn for any required velocity centile, e.g. 1 st or 99 th Switch between distance and velocity centiles

Distance centiles

Velocity centiles

Growth acceleration Three measurements

Question You observe an infant grow over 4 weeks They show upward or downward centile crossing Ask yourself: How will they grow over the next 4 weeks? Will they stay on the same centile? Will they continue to cross centiles the same way? Or will they cross centiles the other way?

Centile crossing over 4 weeks

Statistics of centile crossing and deviation As before, convert weights to z-scores The change in z-score over 4 weeks is deviation e.g. from birth to 4 weeks: Deviation = z 4 z 0 = d 04 Deviation the same as centile crossing Research question: What is the correlation between successive deviations? e.g. correlation between d 04 and d 48 Possible answers: zero, positive or negative

Two growth studies of Cambridge infants Widdowson Study (1959-65) 1094 infants measured monthly from 0-12 months Representative of Cambridge infants in ~1960 Weights obtained from child welfare clinics Cambridge Infant Growth Study (1984-87) 255 infants measured every 4 weeks from 0-52 weeks Families more selected and of higher social class Infants weighed by experienced research nurse In brief, monthly weights in infancy

Correlation between successive deviations Widdowson Study C I G Study Age 3-4 months positive correlation

Correlation between successive deviations Widdowson Study C I G Study Age 10-11 months negative correlation

Surprise - deviations are correlated! At 3-4 months there is a positive correlation Infants crossing centiles one month are likely to cross centiles in the same direction the next month At 10-11 months there is a negative correlation Infants crossing centiles are likely to cross centiles in the opposite direction the next month Cole TJ, Singhal A, Fewtrell MS, et al. Weight centile crossing in infancy: correlations between successive months show evidence of growth feedback and an infant-child growth transition. Am J Clin Nutr 2016;104:1101-9.

How does the correlation change with age?

Deviation and feedback Before 6 months infants crossing centiles tend to continue to cross centiles After 6 months they tend to cross back again Examples of feedback Positive feedback before 6 months Negative feedback after 6 months

Positive feedback Before 6 months, some young infants want to shift to a different centile Mismatch between fetal growth and target size? So need to cross centiles in same direction for a time But eventually reach their target Example of positive feedback

Negative feedback Older infants depart from growth trajectory due to some exposure e.g. infection leads to downward centile crossing Response is to compensate the following month e.g. catch-up following infection Example of negative feedback

Implications for chart assessment Centile crossing predicts centile crossing But depends on age Early centile crossing (before 6 months) Expect more centile crossing Late centile crossing (after 6 months) Expect reverse centile crossing Mid-age centile crossing Expect centile tracking Easiest to see on weight z-score scale

Growth velocity and growth assessment

Growth acceleration and feedback Assessment of acceleration a novel idea Highlights change from positive to negative feedback Reflects how and why centile crossing becomes less common with increasing age in infancy

Growth pattern Many measurements

Variation in growth pattern Interesting to look at individual growth curves To see how they differ, and how they are the same Here are a sample of growth curves from the Cambridge Infant Growth Study

Summarising growth pattern Curves largely the same shape But differing in position Some high, some low Some steep, some shallow SITAR is a growth curve model that adjusts each curve for being high/low h/l (size) early/late (timing) steep/shallow (intensity)

All growth curves, colour-coded measured every 4 weeks

SITAR SITAR adjustment makes all curves like the mean curve High curves shifted down, low curves up (size) Steep curves made shallower, shallow steeper (intensity) Early curves shifted later, late curves earlier (timing) Size, timing and intensity estimated as random effects Net effect is to superimpose curves Then fit mean curve through superimposed curves

All growth curves, colour-coded measured every 4 weeks

All growth curves, colour-coded after SITAR adjustment

All growth curves, colour-coded with SITAR mean curve

SITAR growth patterns SITAR converts growth curves: to a mean curve: and a growth pattern for each individual: size, timing, intensity Summary like growth distance or growth velocity SITAR - a useful instrument for growth curve analysis. Cole TJ, Donaldson MD, Ben-Shlomo Y. Int J Epidemiol 2010;39:1558-66.

SITAR growth patterns SITAR explains over 95% of variance Very ygood fit So random effects define individual growth pattern Can be used as individual growth summary To relate to earlier exposures or later life course BUT note that SITAR not useful clinically It needs whole growth curve Comes too late to make clinical decisions

Conclusions Growth summary for one, two, three and many measurements Distance, velocity, acceleration and pattern Useful to assess growth in individuals Improving decision making Shows how statistics can help in the assessment of growth