WPS8529 Policy Research Working Paper 8529 Joint Effects of Parenting and Nutrition Status on Child Development Evidence from Rural Cambodia Jan Berkes Abbie Raikes Adrien Bouguen Deon Filmer Development Economics Development Research Group July 2018 Policy Research Working Paper 8529 Abstract Substantial work has demonstrated that early nutrition these inequalities can be explained by differences in paren- and home environments, including the degree to which tal stimulation and early nutrition status. Better educated children receive cognitive stimulation and emotional sup- parents engage in better parental activities that stimulate port from parents, play a profound role in influencing early children’s development. However, the positive association childhood development. Yet, less work has documented between parental activities and child outcomes is particu- the joint influences of parenting and nutrition status on larly strong for non-stunted children, and parental activities child development among children in the preschool years can only explain about 8–14 percent of the cognitive gap living in low-income countries. Using panel data on par- between the lowest and highest wealth quintiles. The results enting, nutrition status, and early developmental outcomes highlight the need for integrated interventions that address of about 7,000 Cambodian preschool-age children, this both parenting and early nutrition, also suggesting that paper demonstrates that inequities in early development parenting interventions for the most disadvantaged fam- associated with family wealth are evident at the start of ilies should be carefully designed and evaluated to ensure preschool and increase over time. A significant share of maximum effectiveness. This paper is a product of the Development Research Group, Development Economics. It is part of a larger effort by the World Bank to provide open access to its research and make a contribution to development policy discussions around the world. Policy Research Working Papers are also posted on the Web at http://www.worldbank.org/research. The authors may be contacted via jberkes@diw.de or dfilmer@worldbank.org. The Policy Research Working Paper Series disseminates the findings of work in progress to encourage the exchange of ideas about development issues. An objective of the series is to get the findings out quickly, even if the presentations are less than fully polished. The papers carry the names of the authors and should be cited accordingly. The findings, interpretations, and conclusions expressed in this paper are entirely those of the authors. They do not necessarily represent the views of the International Bank for Reconstruction and Development/World Bank and its affiliated organizations, or those of the Executive Directors of the World Bank or the governments they represent. Produced by the Research Support Team Joint Effects of Parenting and Nutrition Status on Child Development: Evidence from Rural Cambodia Jan Berkes, Abbie Raikes, Adrien Bouguen, Deon Filmer* Keywords: parental stimulation, parenting, nutrition, early development, wealth gap, Cambodia JEL codes: J13; I14; I24 * Jan Berkes, Department of Education and Family, German Institute for Economic Research (DIW Berlin). Abbie Raikes, University of Nebraska. Adrien Bouguen, University of California Berkeley. Deon Filmer, World Bank. The authors do not have any interests that might be interpreted as influencing the research. Correspondence: Jan Berkes, Department of Education and Family, DIW Berlin, Mohrenstraße 58, 10117 Berlin, Germany; email: jberkes@diw.de The authors thank Tsuyoshi Fukao and members of the Early Child Development Department of the Royal Government of Cambodia’s Ministry of Education for input and assistance. Useful comments from Sophie Naudeau and Magdalena Bendini on an earlier version of this paper are gratefully acknowledged. This study benefited from funding from the World Bank’s Strategic Impact Evaluation Fund (TF0A0326) and the Ministry of Education, Youth and Sport of Cambodia. The authors are responsible for any errors. The findings, interpretations, and conclusions expressed in this paper are those of the authors and do not necessarily represent the views of the World Bank, its executive directors, or the governments it represents. Joint Effects of Parental Stimulation and Nutrition Status on Child Development: Evidence from Rural Cambodia 1) Introduction and Background Early childhood development sets the course for lifelong learning, health and well-being (Black et al., 2017). Despite the potential social capital that can be built through investments in early childhood development (e.g. Heckman, 2006), recent estimates suggest that 250 million children living in low- and middle-income countries are at risk of not fulfilling their developmental potential (Black et al., 2017). Family poverty is associated with several conditions that negatively impact child development. Children who experience early stunting, poor health or lack of educational opportunities are less likely to achieve their developmental potential than other children (Britto et al., 2017), with a common underlying influence of family poverty. Family poverty has been shown to have a wide range of negative effects on child development in many countries, with gradients in child development due to family income evident as early as 4 months of age and growing throughout the course of childhood (Paxson and Schady, 2007; Naudeau et al., 2011, Fernald, Kariger, Hidrobo & Gertler, 2012; Schady et al., 2015; World Bank 2018). Yet, the mechanisms accounting for the correlation between family wealth and child development remain imperfectly described and likely include a range of factors, including health and nutrition as well as the quality of home environments (e.g., McLoyd, 1998; Black et al., 2017). While perinatal (in-vitro and during the first months after birth) health and nutrition conditions are known to impact later cognition (Veena, Gale, Krishnaveni, Kehoe, Srinivasan, & 2 Fall, 2016), little is known about how parenting may help compensate for disparities in development due to poverty. Impacts of parenting and family wealth on child development Defining the mechanisms that link family assets to child development is critical for developing effective interventions. One such mechanism is parenting. Parenting quality has been shown to have strong and long-term impacts on child development (Bradley & Corwyn, 2002; see Bornstein, 2005, for a complete review). Parenting has been posited to be a central mediator of the relationship between family poverty and negative child outcomes (Blair & Raver, 2010). Parenting quality may be lower among families living in poverty due to the stress of living without adequate economic resources and the concomitant characteristics of poverty, such as low levels of parental education, which in turn contributes to the negative impacts of poverty on children. Thus, a key hypothesis is that parenting quality mediates the association between family poverty and child development. Several studies have examined the hypothesis that undernutrition arises in part from poor- quality parenting, particularly mediated through lack of maternal responsiveness. Results have been mixed. Maternal depression and other common mental disorders in the first two years after a child’s birth have been associated with child underweight and stunting (Nyugen et al., 2013; Surkan, Kennedy, Hurley & Black, 2011). This effect may be mediated through lack of engaged and stimulating parenting which in turn influences child feeding and growth, but sources of this association are not yet clear. A recent meta-analysis of interventions focused on improving “responsive feeding” yielded few studies documenting positive impacts of responsive feeding, which suggests that responsive parenting is not a primary determinant of undernutrition (Bentley, 3 Wasser & Creed-Kanashiro, 2011). In sum, while existing literature suggests that maternal mental disorders and undernutrition may be associated, the evidence supporting a general association between parenting quality and undernutrition is mixed, suggesting that parenting and undernutrition may influence child development through mostly distinct pathways. Parenting quality clearly has important impacts on child development overall. Both child development and economic theories suggest two sources of variation in parenting on child development: first, that different types of parenting may affect different domains of child development and, second, that parenting may have variable effects on child development depending on child characteristics and environmental conditions. Role of different types of parenting on child development. In interactions with young children, parents engage in multiple forms of teaching and support. Parenting is the primary mechanism by which children’s health and nutrition, as well as overall growth and development are facilitated (Phillips & Shonkoff, 2000). Emotionally responsive parenting refers to parents’ abilities to respond sensitively to children both verbally and physically, while cognitively stimulating parenting refers to parents’ stimulation of learning through use of language, introduction of new games and ideas, as well as exposure to high-quality learning environments (Bornstein & Putnick, 2012). Existing research suggests some degree of differentiation between domains of parenting: parents who are emotionally responsive are not necessarily cognitively stimulating; in addition, there is some degree of differentiation demonstrating that emotionally-responsive parenting is especially important for social/emotional development, while cognitively-stimulating parenting promotes young children’s learning and development. “Negative parenting,” or reliance on harsh, punitive physical and verbal parenting, relates to disruptions in social/emotional 4 development and also delayed cognitive development (Chang, Schwartz, Dodge & McBride- Chang, 2003; Hughes et al., 2006). Maternal education and, specifically, executive function (Cuevas et al., 2014) show relations with the development of children’s executive function. An overall tendency for parents to engage in more emotionally-supportive parenting than cognitively stimulating parenting has been documented across countries (Bornstein & Putnick, 2012). Defining differential impacts of parenting on child development. Recent work has documented that environmental conditions may influence children’s development differently. Describing the interaction between the child and their environment has emerged as a seminal element of understanding child development, especially in relation to the interaction between biological and environmental factors (Ellis, Boyce, Belsky, Bakersmans-Kranenburg, & van Ijzendoorn, 2011). Recent research suggests that children’s responsiveness to environmental conditions, and thus reactivity to both positive and negative stimuli, is not uniform across children (e.g., Belsky, Bakersmans-Kranenburg, & van Ijzendoorn, 2007). While most of this work has been conducted in high-income countries, the concept of studying child–environment interactions may also apply to children who have experienced extreme environmental deprivation resulting in negative health outcomes such as stunting. Economic analyses have also suggested that children with very low cognitive and emotional resources may respond less to parenting, based on the theory that low vectors of skill development persist over time, even when environmental resources are enhanced (Heckman, Stixrud & Urzua, 2006). Little work to date has examined this interaction in the context of environmental conditions that characterize child development in low- and middle-income countries. Nutrition and child development 5 A large body of work demonstrates strong associations between early nutrition status and child development (Prado & Dewey, 2014). Children with poor early nutrition fail to reach expected linear growth milestones, leading to stunting, long-term cognitive and motor deficits, and eventually worse economic outcomes (e.g., Galasso, Wagstaff, Naudeau & Shekar, 2017; Sudfield et al., 2015). Some work similarly demonstrates that cognitive stimulation can ameliorate some of the effects of poor early nutrition on child development, especially when coupled with parenting interventions (Paxson & Schady, 2007; Grantham-McGregor, Powell, Walker & Hines, 1993). Most evidence to date has been produced through intervention studies, with less work conducted using large-scale samples, where effects may be more representative of typical interactions between parenting and nutrition status. The present study examines the emergence of wealth gradients in child development and the joint roles of nutrition status and parenting in explaining the gaps between high- and low- wealth children. We had three hypotheses: first, we expected to find gradients in child development based on family wealth in a low-income country. Second, based on existing parenting literature, we expected to find that there are different types of parenting, such as cognitively-stimulating and emotionally-responsive, with unique implications for child development, and that parenting quality is associated with both family wealth and child development. Third, we hypothesized that parenting and nutrition status would jointly influence child development, and that parenting could have different impacts on child development depending on nutrition status. In sum, we hypothesized that family assets would exert a strong and consistent impact on child development, but that parenting quality and nutrition status could account for a portion of the variance attributable to family wealth, demonstrated by a reduction in the impact of family assets on child development when parenting and nutrition status are 6 included in the model. Our analysis is informed by two rounds of data, collected from a large sample of children and their caregivers in rural Cambodia. 2) Methods We investigated our hypotheses using data from a large panel of children in rural Cambodia. We rely on a unique panel data set of almost 7,000 children (age 2-4 in 2016) and their families, collected in rural Cambodia in 2016 and 2017. The child survey comprises a comprehensive battery of cognitive tests (executive function, language and early numeracy) as well as anthropometric measures (height and weight). The caregiver survey includes basic socio- economic characteristics, questions about the socio-emotional development of the child (Strengths and Difficulties Questionnaire), as well as 25 questions about parental practices, measuring three key dimensions of parenting (cognitive parenting, emotional parenting, negative parenting). Thus, the caregiver survey allows for distinguishing the role of socio-economic characteristics and the role of specific parental practices in explaining children’s development. Participants During the baseline data collection exercise in 2016, a survey firm1 used an adapted version of the EPI walk to select 6,934 eligible households in 305 villages.2 Each sampled household included children aged between 24 and 59 months.3 Approximately 11 months later, 1 Angkor Research©, under the close supervision of the Principal Investigators. 2 EPI refers to the Expanded Programme on Immunization of the World Health Organization (see e.g. Henderson & Sundaresan (1982)). 3 The study was approved by the institutional review board of the Paris School of Economics and the National Ethics Committee for Health Research by the Cambodian Ministry of Health. The data were collected as part of a randomized controlled trial investigating the effects of preschools on school readiness of children in rural Cambodia. Our main outcomes for this study, collected in 2017, stem from a first follow-up on the trial intervention. Between baseline and the first follow-up, community preschools were constructed in some of the sample villages. An awareness campaign was conducted to stimulate enrollment into the new preschools. Further, a home-based program, in which trained core-mothers address a range of topics around parenting, was implemented. These 7 the same survey firm revisited 93% of the sampled households.4 The follow-up sample, which forms the main sample for the analysis in this study, includes 6,452 households with 6,914 tested children at the age 3.91), of 2 ( 184), 3 ( 2271), 4 ( 2459), 5 ( 1965) and 6( 35). The average number of followed up households per village is 21.2 and the average number of eligible children per household is 1.07. Using the follow-up sample instead of baseline allows to control for baseline test scores in regressions and to observe outcomes during an age period which is critical for the development of language and early numeracy skills. Upon arrival at households’ home, the enumerator first identified the main caregiver and the eligible children,5 and administered the household and caregiver surveys as well as completed the anthropometric measures and development tests. The child’s testing session lasted an average of 31 minutes ( 50 30). Measures We developed and adapted the battery of child tests in close cooperation with the researchers of the Measuring Early Learning Quality and Outcomes project (MELQO; UNESCO, 2017) and the local survey firm. The survey firm recruited interviewers based on their familiarity with data collection and their experience with young children. To ensure that children and parents correctly understood the questionnaire and that the instruments were reliable, the research team pretested every instrument at least three times before collecting data on the interventions were randomized between villages and still at an early stage during the time of the first follow-up. We therefore ignore the interventions in this analysis. All our regression estimates are robust to the inclusion of treatment group variables to account for the randomized treatment. 4 The attrition rate of 7% can almost entirely be explained by seasonal and permanent relocation of households. The study does not follow up on households that move beyond the boundaries of sample villages. 5 The caregiver is defined as the direct relative that takes care of the child most of the time. In our sample, the primary caregiver is most often the mother (52% of the children), followed by the grandmother (35%), father (6%), grandfather (2%), aunt or uncle (3%), or an older sibling (1%). 8 sampled villages. The survey firm translated the questionnaires into Khmer and an independent third party back-translated them into English. The research team participated in the interviewer training conducted by the local firm’s fieldwork manager. Field staff was organized into six groups, each comprising of four interviewers, one supervisor and one field editor. All supervisors had several years of data collection experience in Cambodia and were responsible for household sampling and quality control procedures. Editors supported supervisors doing spot-checks, interviewer observations, and conducted independent re-interviews of at least 20% of interviews in each village.6 Household survey. The household survey includes information about family structure and socio-economic background. We constructed a dwelling quality index and household assets index.7 Both indexes were then standardized and summed up to create a wealth index. We use quintiles and quartiles of this index to analyze the cognitive gap in rural Cambodia. To construct a binary measurement of household wealth, we define low living standard as belonging to the lowest 2 quintiles of the wealth index. We use household size 4.48, 1.68 and the binary measure of wealth 0.40 as measures of socio-economic background at the household level. Caregiver survey. The educational background of the caregiver is measured by years of schooling 3.48, 3.15 . For better interpretability, we use a binary indicator for completion of primary school (6 years) .266 in our regressions. Caregiver non-verbal 6 Initial re-interviews were limited to a set of 20 questions to reduce the burden to respondents. In case of inconsistencies between the initial interview and the re-interview, households were called or visited again to obtain clarity. 7 The direction of each categorical variable for dwelling quality (quality of the floor, quality of the roof, electricity, toilet…) was first assessed using a Multiple Correspondence analysis. Then, each variable was coded using the sign of the coordinate of the first MCA dimension. We then took the average of the standardized version of each coded variable. The Cronbach alpha is 49%. For the household asset, we simply standardized each asset and took the standardized version of the sum. The Cronbach alpha is 65%. 9 reasoning is measured using a set of 12 Raven’s Progressive Matrices (Carpenter et al. 1990; Raven 2000). This test was also administered at baseline. When administered to the same caregiver again 11 months later, non-verbal reasoning scores showed a correlation of 0.46. To account for the multidimensionality of parent-child interaction (Bornstein, 2012), we developed a questionnaire that covers multiple aspects of parenting and that fits well the Cambodian context. Caregivers answered 25 items related to regular activities, parent-child bonding, punitive behavior and different child-rearing methods. Possible responses ranged from 1 (never or almost never) to 5 (every day or almost every day). The factor structure of these parenting questions is reported in the Results section. Child survey and assessment. We control for child age, measured as the number of days between birth and follow-up, divided by 365.25. Stunted growth is defined as having a height- for-age below two standard deviations of the growth standards by the World Health Organization (2006) .380, .706 . We use cutoff points other than -2 SD in the robustness analysis below. The child assessment includes different tests for the domains of language development, early numeracy and executive function. Most of the child tests are based on the Measuring Early Learning Quality and Outcomes (MELQO) toolkit (see UNESCO, 2017 for a description of the measures development process). We implemented multiple pretests to obtain reliable measures. Some tests from MELQO were removed during the adaptation process. For example, we removed a backward-digit span test and a mental transformation tasks because very few children were able to understand them. 10 Likewise, as children in our sample were unable to write Khmer letters, we removed tests that included writing. Sound (phoneme) discrimination tests were also too difficult for our sample of children. Some of the MELQO tests included in our battery were adapted to the local context or modified to fit the competences of the sampled children. For example, as some children rarely visit markets or stores, a question asking for things that can be bought and eaten at a market was changed to things that can be eaten. Likewise, letters and numbers used in some of the tests were selected to fit the local preschool curriculum. A number identification test was replaced with a quantity comparison test where children had to compare quantities of printed objects instead of printed numbers. A simple addition test was merged with the same quantity comparison test by asking children questions about the total number of objects obtained after adding or subtracting a given number of objects. Since the overall purpose was to generate a well-functioning test for Cambodia (as opposed to maintaining consistency for other purposes, such as international comparisons), adaptations were prioritized to ensure adequate fit to local context. To obtain a rich set of tests that capture different constructs of the main domains, the MELQO tests were complemented with other assessments in order to derive a test scale that covers these different constructs. Specifically, we added the Dimensional Change Card Sort test (Zelazo, 2006), a receptive vocabulary test similar to the PPVT (Dunn & Dunn, 2007), a test for familiarity with print objects and a sustained attention test. The measure of socio-emotional problems is based on the caregiver-reported SDQ instead of a child assessment. Note that the SDQ results are not reversed: a higher score means more behavioral difficulties. Before constructing the composite scores of child test domains, individual tests were first scored and standardized, thus ensuring that all tests contributed equal variance to their composite 11 score. Then, tests of one category were summed into a composite score and standardized again for better interpretability. After these steps, we obtained the following composite scores: Executive function (α .74, .60 .8 The final executive function score includes four tests, of which the first two listed here are from the MELQO toolkit. The construct inhibitory control is assessed with the head-knee task (McClelland et al., 2014). Working memory (short-term auditory memory) is assessed with a forward digit span test. The Dimensional Change Card Sort test is used as a measure of cognitive flexibility. We use a self- developed cancellation task to measure sustained attention. In this test, children see a printed matrix with different symbols and are asked to cross-out all symbols that match the given one (e.g. cross out all flowers). When completed, a larger matrix is given and a new symbol has to be crossed-out. The test continues until a child has completed 4 matrices, crossed out more wrong than correct images in a matrix, until the child loses attention, or states that it is done. Language development ( .71, .65 . Language and literacy skills are assessed with a letter name knowledge test in which children have to identify common letters of Khmer script. Expressive language skills are assessed with a self-developed test, similar to the PPVT (Dunn & Dunn, 2007) where children are asked to match a word to one out of four pictures. Receptive language is assessed with a listening comprehension test. Knowledge of reading concepts is assessed by showing a children’s storybook and asking questions such as where and in which direction one should start reading the story. A drawing test, where children copy shapes, like circles or squares, is used to assess fine-motor skills. 8 Cronbach’s for composite test scores are based on the covariance between the standardized test scores of the individual tests. The rank correlation coefficient for composite test scores gives the correlation between baseline and follow-up composite scores with 11 months in between on average. 12 Early numeracy ( .73, .56 . The tests for early numeracy includes a test for measurement concepts, e.g. if the child understands concepts such as tallest/shortest. In a test for verbal counting, children had to count up to 30. Numbers and operations are also administered with a quantitative comparison test, a number identification test analogous to the letter name knowledge test and a shape recognition test where children have to identify basic geometric shapes. Socio-emotional problems .69, .26 . To test for behavioral problems, we use the total difficulties score of the Strengths and Difficulties Questionnaire (SDQ). The score is composed of 5 questions for each of the four categories: emotional symptoms, conduct problems, hyperactivity/inattention, and peer relationship problems. This test for socio-emotional problems is the only child development measure that is based on the caregiver survey instead of the direct child assessment. 3) Results The analyses plan was designed to test the three main research questions: 1) are there gradients in child development based on family wealth for children three to five years of age; 2) does parenting show unique components and if so, how do they differentially affect child development, and how is parenting associated with family wealth and child development; 3) how do parenting and nutrition status jointly contribute to child development, and do they account for a portion of the variance in child development due to wealth gradients? We are also interested in looking across these three questions to estimate the proportion of child development accounted for by variations in family wealth, parenting quality, and nutrition status. 13 Are there gradients in child development based on family wealth for children aged three to five years? We start our analysis by looking at the wealth gradient for our four child development measures in Table 1 (“No control” columns). The data show large differences between children from households with low and high levels of wealth. The top wealth quintile performs between 0.44 and 0.70 SD higher than the lower quintile and between 0.42 and 0.60 SD separate the top and bottom quartile in cognitive competences (language, early numeracy and executive function). Since one of our language tests is fairly similar to the TVIP used by Schady et al. (2015), we use the language score to compare the size of these gaps with others documented in Latin America. In rural Cambodia, the cognitive gap lies at the lower end of the wealth gradients of the five Latin America countries studied by Schady et al. (2015) which range from 0.57 SD to 1.21 SD. The gap for socio-emotional problems (-0.22 between top and bottom quintiles) is smaller than for cognitive competences. A larger socio-emotional problem gap might arise later, potentially as a consequence of earlier cognitive difficulties. Importantly, the socio-emotional problems (SDQ) are reported by parents while other tests are directly administered to children, which means that measurement error may contribute to dampening the size of the gap. The last three rows (age 3, age 4 and age 5) of Table 1 and column ‘No controls’ clearly indicate that the wealth gaps significantly widen with time. Between 3 and 5 years old, the gap increases, especially for early numeracy and language. From a small 0.11 SD at age 3 in early numeracy (0.31 in language), the gap reaches a large magnitude of 0.87 SD in early numeracy (1.12 in language) by age 5. This result is consistent with the fact that the 3 to 5 years of age period is critical for both competences. The gap increases to a lesser extent in executive function from 0.37 to 0.54 SD. The gap for socio-emotional problems increases from -0.18 to -0.26 SD. 14 To better visualize the differences in child development by socio-economic status over age, we used separate nonparametric regressions for the lowest and highest quintiles of the SES index (Figure 1). Gaps in executive function, language and early numeracy skill are large and already evident at the age of 3, the earliest enrollment age for public preschools in Cambodia. The gradient in executive function skills stabilizes at the age of 4, with a size of approximately 0.5 standard deviations. This is in stark contrast to gradients in language and early numeracy skills. The initial gap in language at the age of three is comparable to executive function. However, the gap continues to grow until primary school-age to a magnitude of 0.9 SD (early numeracy) and 1.1 (language) standard deviations. This gap is of substantial magnitude and equivalent to a delay of 14-18 months (a child progressing of .75 SD per year in our sample). Figure 1. Nonparametric regressions of test scores on age using a bandwidth of 0.5 years. 15 Does parenting show unique components, and associations with family wealth and child development? We conducted explanatory factor analysis in order to determine if the parenting scores could be grouped by distinct dimensions of parenting. Sampling adequacy was supported by a Kaiser-Meyer-Olkin statistic of 0.868 (Kaiser, 1974), confirming a compact pattern of correlations between the 25 parenting items. Bartlett’s test for sphericity rejected the null hypothesis of zero covariance, p<.001. Parallel analysis suggested that there are 11 factors, equaling the number of factors with eigenvalue greater than zero. We followed a similar procedure as Jones et al. (2017) to reduce the number of factors retained. We rejected the large number of factors as only 3 of them had more than 1 item with a factor loading greater than 0.3 and all but 3 factors had poor reliability .6 . We gradually reduced the number of factors until all factors had reliability of at least .6 and at least two items with loadings greater than 0.3. With this procedure we obtained three parenting factors, which are summarized in Table 2. The three identified parenting factors capture three dimensions of parenting. First, “cognitive parenting” measures how parents actively interact with the child in ways that are likely to develop cognitive competencies: games, reading books and playing with toys or objects. The second dimension, “socio-emotional parenting” captures the emotional support and responsiveness of the child and the parent. Third, “negative parenting”, in which higher values indicate more negative parenting, refers to the degree to which parents use harsh or punitive statements to discipline the child. It is important to note that the constructs of emotional support and responsiveness and the use of harsh or punitive statements are not two extremes of the same dimension; instead, in our sample, they appear to be unique constructs. To assess the extent to which the three parenting scores capture meaningful differences in parenting behavior, we correlate them with other family characteristics. Table 3 indicates that 16 parenting quality is strongly correlated with socio-economic characteristics of the household and of the caregiver: wealthier and more educated families tend to adopt parenting behavior that is more cognitively stimulating and positive towards children. Cognitive parenting is positively correlated with caregiver education 0.26, .001 and caregiver non-verbal reasoning 0.18, .001 , but negatively correlated with a low household living standard ‐0.10, .001 and the caregiver’s age ‐0.15, .001 . A similar pattern was observed for socio-emotional parenting, although the correlations are smaller. Interestingly, while cognitive parenting and socio-emotional parenting are positively correlated 0.38, .001 , negative parenting is negatively correlated with cognitive parenting ‐0.05, .001 but positively correlated with socio-emotional parenting 0.14, .001 , i.e. parenting behaviors which often rely on punitive behavior are not inconsistent with other behaviors linked to positive socio- emotional interactions. Cognitive parenting is significantly positively correlated with scores on the child tests, negative parenting is significantly negatively correlated with these scores, and the direction of the effect of socio-emotional parenting is ambiguous in this correlation analysis. While socio- emotional parenting is negatively correlated with socio-emotional problems ‐0.04, .001 , it is also negatively correlated with the executive function score ‐0.03, .040 . Since parenting is significantly associated with several other household characteristics, isolating the independent role of parenting is important. We show below that after controlling for household and caregiver characteristics and even baseline test scores, the effect of parenting remains significant and that socio-emotional parenting does not negatively affect executive function skills. 17 What are the joint roles of parenting and nutrition in explaining variation in child development? The size of the raw wealth gap is substantially reduced after conditioning on parenting (the three separate parenting composite scores) or anthropometrics (height for age z- scores). Both sets of control variables reduce the gap in executive function by the same amount, about 8% (Table 1). The gap in language skills is reduced by about 10% when controlling for parenting and 6% when controlling for height for age. Similarly, the gap in early numeracy skills is reduced by 14% when controlling for parenting, and by 7% when controlling for nutrition. The biggest difference is in the case of socio-emotional problems where our parenting measures can account for 44% of the gap and height for age cannot explain any of it (0%).9 To investigate further the variation accounted for by parenting and other factors, as well as to test the hypothesis that the wealth gap is partially mediated by parenting and nutrition, in Table 4, we use four specifications in a hierarchical regression model. In a first step (column (1)), we regress the dependent variable on child age and child gender to determine how much of the total variation is explained by exogenous characteristics. The R-squareds in these regressions vary from 2.1% for socio-emotional problems, to 38% for early numeracy and approximately 43% for executive function and language skills. Note that the three child-administered test scores increased by approximately 0.72 to 0.78 SD by year of age. This means that the cognitive gap that separates a child from the bottom quartile and the top quartile (between 0.42 and 0.60 SD, see Table 1) corresponds to between 60 and 80 percent of one year of average development. 9 As with all these results, the estimates should not necessarily be interpreted as causal impacts. In this case in particular, for example, it is possible that socio-emotional problems could be the cause rather than result of negative parenting practices (i.e. reverse causation). 18 Girls perform significantly better than boys across all domains (9-12% of a SD). The gender difference is particularly pronounced for executive function. Common caregiver and household control variables added in column 2 explain additional variation for the domains of executive function Δ .021 , language Δ .039 , early numeracy Δ .023 , and socio-emotional problems Δ .016 . Joint significance of the added set of variables was confirmed with an -test ( .001 for all 4 outcomes). Low living standard, household size, caregiver education and caregiver non-verbal reasoning affect children’s cognition in the expected directions, with all coefficients highly significant and of relatively large magnitude. To test if the three parenting variables explain additional variation in child outcomes, they are added to the regression model in column 3 of Table 4. Parenting explains additional variance for all four child outcomes. The change in the R-squared is by far the highest for socio- emotional problems Δ .122 and of relatively smaller size for executive function Δ .003 , language Δ .010 , and early numeracy Δ .009 . Joint significance of the three parenting variables was confirmed with an -test ( .001 for all 4 outcomes). For the cognitive domains, the parameters of all household and caregiver controls are reduced once parenting is added to the model. The largest change is for caregiver education and household living standard for the socio-emotional problems outcome. Of our three parenting measures, cognitive parenting had the strongest effect on the three cognitive scores. A one SD increase in our measure of cognitive parenting is associated with a 0.092-0.098 SD increase in language and early numeracy scores, an increase that is comparable to half to two-thirds of the negative effect of stunting, low living standard, or low caregiver 19 education. Based on the magnitudes of these associations, to compensate for the cumulative detrimental effect of low caregiver education, low living standard and stunting, the cognitive parenting behavior score would need to move from the bottom 19 percentile to the top 95 percentile (3 SD), i.e. from a parent with poor cognitive parenting behavior to one of the best parents in our sample. Given this large magnitude, it is perhaps unreasonable to expect that cognitive stimulation alone might be able to make up for the effects of negative shocks or a negative environment. Socio-emotional parenting affects executive function skills positively but does not affect the other cognitive competences. Higher scores on socio-emotional parenting are also significantly associated with less socio-emotional problems. While smaller in size than for cognitive parenting, the effect of negative parenting is significantly negative for all cognitive competences. The size of the relationship is especially large for the socio-emotional problem score (0.346 SD). To test if nutrition also explains additional variance we performed a similar exercise by adding it to the model in column 4 instead of the parenting variables. Stunting cannot explain variance in socio-emotional problems Δ .000 , similar to Table 1. It explains a larger share of variance in executive function than our measures of parenting do Δ .008 , a similar share for language Δ .008 and a slightly smaller share for early numeracy Δ .006 . The coefficients of household and caregiver variables change only marginally when stunting is added to the model. Thus, while the relationships with outcomes are significant—equivalent to 20 about 3 months of cognitive development, for example—stunting seems to be a weaker mediator of caregiver education than parenting behaviors are. 10 We explore additional associations and mediating factors in Table 5. First, we include baseline test scores collected 11 months prior to the main outcomes of this study (which would also have the effect of controlling for unobserved factors already captured in those scores). While the magnitudes fall, the relationships with parenting typically remain significant, suggesting that the effect of parenting is cumulative (and not “one-shot”).11 We also hypothesized that parenting may have differential effects on child development, depending on children’s nutrition status. To test the hypothesis that the association of cognitive parenting and child outcomes is weaker for stunted children, we looked in column 3 of Table 5 at the effect of parenting, separately for stunted and non-stunted children and with baseline controls. We focused on cognitive parenting because of the notable role of nutrition status in cognitive development; we thus hypothesized that cognitive parenting and subsequent child development would be especially sensitive to nutrition status. While no significant interactions between parenting and child nutrition status were observed for executive function and socio- emotional problems, effects of cognitive parenting differed by stunting for language and early numeracy outcomes. The magnitude of the cognitive parenting effect is about 35% weaker for effects on language and 54% weaker for effects on early numeracy, showing that the positive effect of cognitive parenting is weaker for stunted children. We repeated the same exercise 10 The fact that wealth gradients remain large even after controlling for anthropometrics and parenting should be interpreted with care. As our measure of parenting is only self-reported and based on recall, it is possible that it does not capture the full dynamic of parenting throughout the children’s life time. Our measure of parenting could for instance be imperfectly correlated with past parenting (during pregnancy or breastfeeding). Our measure of parenting may also imperfectly measure the quality of parenting. Longitudinal and external measures of parenting would be necessary to enhance our measures. 11 Note that in Table 5 we control for both parenting and stunting. The effects of each are relatively unchanged. 21 without baseline controls in column 2 of Table 6 and observed similar interaction effect sizes. However, the interaction term for the language outcome is only significant in Table 5 where baseline test scores were added to the model.12 While the cognitive parenting x stunting interaction is statistically significant for our language and early numeracy outcomes, this finding can also be explained by an unobserved heterogeneous effect by parenting quality, since parenting quality is correlated with stunting. As a proxy for parenting quality, we therefore chose caregiver education, which is also negatively correlated with stunting. We then added a cognitive parenting x caregiver education to the models in the last columns of Table 5 and Table 6 to see if the coefficient of the cognitive parenting x stunting interaction changes. The inclusion of caregiver education, and its interaction with nutrition status, did not change the significance of the previous interaction. We thus conclude that cognitive parenting has a stronger impact on some domains of child development for non-stunted vs. stunted children. 4) Discussion The objective of this study was to quantify the role of parenting and nutrition in explaining the cognitive gap that we observed between children from low and high wealth backgrounds. This gap exists in Cambodia and is large: the difference between the highest and 12 We also tested if our results are robust to different cutoff points for the stunting indicator. As Perumal et al. (2018) point out, there is no biological basis for the commonly used cutoff at -2 SD in the height-for-age score (HAZ) and children slightly below and above the cutoff might be equally affected by undesirable health outcomes due to growth impairments. Using a binary indicator equal to one if HAZ<-1.5 M=0.592 does yield smaller and insignificant estimates for the cognitive parenting-stunted interaction term of column 3 in Table 5 for language -0.0233, .183 and early numeracy development -0.0165, .396 . Yet, using a lower cutoff score of HAZ<-2.5 M=0.201 yields slightly larger and significant interaction effects on language -0.0311, .057 and early numeracy development -0.0433, .014 than found in Table 5. The same pattern is observed without baseline control variables. This demonstrates that the weaker association of cognitive parenting and test scores is predominantly driven by children in the lower tail of the HAZ distribution. 22 lowest quintiles amounts to about 7-9 months of child development. For pre-academic competencies (language and early numeracy) the gap widens significantly between age 3 to 5: confirming that this period is crucial for development of early numeracy and language skills and is a time when gaps emerge and may become entrenched. We document four main findings. First, building on previous research on socio-economic gradients in developing countries, we confirm that a substantial cognitive gap exists between children from different backgrounds in rural Cambodia (confirming earlier findings in Naudeau et al., 2011). In our sample, 0.60 SD separate the children in the bottom quartile of a wealth distribution from the children in the top quartile (language test). The gap appears larger for cognitive competences (executive function, language, early numeracy) than for non- cognitive/socio-emotional outcomes (Strengths and Difficulties Questionnaire). While significant and large, these gaps in rural Cambodia are at the low end of those documented in rural Latin America by Schady et al. (2015) (between 0.57 to 1.12 SD on a similar language test). Second, the gap in child development due to family wealth widens with age (e.g. Paxson & Schady, 2007; Fernald et al., 2012). Between three and five years of age, the gap especially widens in language (from 0.31 SD at age 3 to 1.12 SD at age 5) and early numeracy, but less for executive function (from 0.37 SD to 0.54 SD). It also increases strongly for socio-emotional problems. The fact that the executive function gap is already present at age 3 and only increases by 46% between age 3 and 5 suggests that before age 3, children are already experiencing substantial growth in executive function that is affected by family socio-economic status. Third, we confirm that cognitive-stimulating and emotionally-responsive parenting is one of the important explanatory factors of cognitive differences due to family wealth. Our measures of parenting behaviors are strongly correlated with children’s performance, accounting for about 23 8-14% of the gap in cognitive competences observed among 3 and 5 year olds. While the explanatory power of parenting is robust, the magnitude of the parental effect is relatively small in comparison with other explanatory factors: We show that a large 1 SD improvement in the main parental competences--an improvement beyond the scope of most parenting programs-- would only compensate about a half of the effect of stunting. Fourth, we provide evidence that the positive effects of parenting may be stronger for children who are not stunted. This is an important finding for two reasons. First, it demonstrates the importance of designing and testing integrated program models of child development that anticipate the complex interactions that will arise between child characteristics and environmental stimuli. It is unlikely that “one size fits all” (Shonkoff, Radner & Foote, 2017). Second, this interaction – albeit a small one, and not consistent across all domains of parenting or child development – adds to the growing body of research documenting the variable effects of environments on child development, by adding undernutrition as a condition that may moderate the effects of later environmental stimuli on child development. Encouraging stimulating and responsive parenting and improving children’s educational environments are key to reducing cognitive inequality and improving children’s performance during that critical time when the cognitive gap widens. Yet, to fully compensate for initial inequality, improvements in the measured dimensions of parenting are not a panacea: other factors (nutrition, poverty level) are important. Poverty and undernutrition must be targeted, along with parenting, and program designers should expect that the effects of various interventions may vary based on children’s nutrition status, especially if undernutrition is already evident. As with all elements of early development, prevention and early intervention are more effective than amelioration – and in this case, that may pertain to addressing undernutrition 24 before conception and during pregnancy. At a minimum, the findings from this study confirm that a combination of interventions may be needed to effectively support young children’s development in low-income countries. The wealth gradient in child development is well-documented in many countries. We show that parenting behaviors explain a share of this—but other factors matter as well. More work is needed to examine the broader range of environmental conditions that contribute to this gap, and perhaps more critically, how they interact with one another. Such analyses would help guide the design of effective interventions which would exploit complementarities between parenting, nutrition and other factors for the promotion of child development in challenging contexts. 25 References Alkire, S., & Santos, M. E. (2010). Acute Multidimensional Poverty: A New Index for Developing Countries (No. HDRP-2010-11). Human Development Report Office (HDRO), United Nations Development Programme (UNDP). Belsky, J., Bakermans-Kranenburg, M. J., & Van IJzendoorn, M. H. (2007). For better and for worse: Differential susceptibility to environmental influences. Current directions in psychological science, 16(6), 300-304. Bentley, M. E., Wasser, H. M., & Creed‐Kanashiro, H. M. (2011). Responsive Feeding and Child Undernutrition in Low-and Middle-Income Countries, 2. The Journal of Nutrition, 141(3), 502-507. Black, M. M., Walker, S. P., Fernald, L. C., Andersen, C. T., DiGirolamo, A. M., Lu, C., ... & Devercelli, A. E. (2017). Early childhood development coming of age: science through the life course. The Lancet, 389(10064), 77-90.. Blair, C., & Raver, C. C. (2012). Child development in the context of adversity: experiential canalization of brain and behavior. American Psychologist, 67(4), 309. Bornstein, M. H. (Ed.). (2005). Handbook of parenting: volume 4 social conditions and applied parenting. Psychology Press. Bornstein, M. H., & Putnick, D. L. (2012). Cognitive and socioemotional caregiving in developing countries. Child development, 83(1), 46-61. Bradley, R. H., & Corwyn, R. F. (2002). Socioeconomic status and child development. Annual review of psychology, 53(1), 371-399. Britto, P. R., Lye, S. J., Proulx, K., Yousafzai, A. K., Matthews, S. G., Vaivada, T., ... & 26 MacMillan, H. (2017). Nurturing care: promoting early childhood development. The Lancet, 389(10064), 91-102. Carpenter, P. A., Just, M. A., & Shell, P. (1990). What one intelligence test measures: a theoretical account of the processing in the Raven Progressive Matrices Test. Psychological review, 97(3), 404. Chang, L., Schwartz, D., Dodge, K. A., & McBride-Chang, C. (2003). Harsh parenting in relation to child emotion regulation and aggression. Journal of family psychology, 17(4), 598. Cuevas, K., Deater‐Deckard, K., Kim‐Spoon, J., Watson, A. J., Morasch, K. C., & Bell, M. A. (2014). What's mom got to do with it? Contributions of maternal executive function and caregiving to the development of executive function across early childhood. Developmental Science, 17(2), 224-238. Dunn, L. M., Padilla, E. R., Lugo, D. E., & Dunn, L. M. (1986). Test de vocabulário en imágenes Peabody: adaptación hispanoamericana. American Guidance Service. Dunn, L. M., & Dunn, D. M. (2007). PPVT-4: Peabody picture vocabulary test. Pearson Assessments. Ellis, B. J., Boyce, W. T., Belsky, J., Bakermans-Kranenburg, M. J., & Van IJzendoorn, M. H. (2011). Differential susceptibility to the environment: An evolutionary– neurodevelopmental theory. Development and psychopathology, 23(1), 7-28. Galasso, E., Wagstaff, A., with Naudeau, S. and Shekar, M. (2017) The economic costs of stunting and how to reduce them. World Bank Group Policy Research Note. The World Bank. http://pubdocs.worldbank.org/en/536661487971403516/PRN05-March2017- Economic-Costs-of-Stunting.pdf 27 Grantham-McGregor, S. M., Powell, C. A., Walker, S. P., & Himes, J. H. (1991). Nutritional supplementation, psychosocial stimulation, and mental development of stunted children: the Jamaican Study. The Lancet, 338(8758), 1-5. Fernald, L. C., Kariger, P., Hidrobo, M., & Gertler, P. J. (2012). Socioeconomic gradients in child development in very young children: Evidence from India, Indonesia, Peru, and Senegal. Proceedings of the National Academy of Sciences, 109(Supplement 2), 17273- 17280. Heckman, J. J. (2006). Skill formation and the economics of investing in disadvantaged children. Science, 312(5782), 1900-1902. Heckman, J. J., Stixrud, J., & Urzua, S. (2006). The effects of cognitive and noncognitive abilities on labor market outcomes and social behavior. Journal of Labor economics, 24(3), 411-482. Henderson, R. H., & Sundaresan, T. (1982). Cluster sampling to assess immunization coverage: a review of experience with a simplified sampling method. Bulletin of the World Health Organization, 60(2), 253. Hughes, C., & Ensor, R. (2006). Behavioural problems in 2‐year‐olds: links with individual differences in theory of mind, exective function and harsh parenting. Journal of Child Psychology and Psychiatry, 47(5), 488-497. Jones, P. C., Pendergast, L. L., Schaefer, B. A., Rasheed, M., Svensen, E., Scharf, R., ... & Seidman, J. C. (2017). Measuring home environments across cultures: Invariance of the HOME scale across eight international sites from the MAL-ED study. Journal of school psychology, 64, 109-127. Kaiser, H. F. (1974). An index of factorial simplicity. Psychometrika, 39(1), 31-36. 28 McClelland, M. M., Cameron, C. E., Duncan, R., Bowles, R. P., Acock, A. C., Miao, A., & Pratt, M. E. (2014). Predictors of early growth in academic achievement: The head-toes-knees- shoulders task. Frontiers in psychology, 5, 599. McLoyd, V. C. (1998). Socioeconomic disadvantage and child development. American psychologist, 53(2), 185. Naudeau, S., Martinez, S., Premand, P., and Filmer, D. (2011). Cognitive development among young children in low-income countries. In Alderman, H., ed. No Small Matter: The impact of poverty, shocks, and human capital investments in early childhood development. The World Bank. Washington, DC. Nguyen, P. H., Saha, K. K., Ali, D., Menon, P., Manohar, S., Mai, L. T., ... & Ruel, M. T. (2014). Maternal mental health is associated with child undernutrition and illness in Bangladesh, Vietnam and Ethiopia. Public health nutrition, 17(6), 1318-1327. Perumal, N., Bassani, D. G., & Roth, D. E. (2018). Use and Misuse of Stunting as a Measure of Child Health. The Journal of nutrition, 148(3), 311-315. Paxson, C., and Schady, N. (2007). Cognitive development among young children in Ecuador: The roles of wealth, health, and parenting. Journal of Human Resources, 42(1):49-84. Phillips, D. A., & Shonkoff, J. P. (Eds.). (2000). From neurons to neighborhoods: The science of early childhood development. National Academies Press. Prado, E. L., & Dewey, K. G. (2014). Nutrition and brain development in early life. Nutrition reviews, 72(4), 267-284. Raven, J. (2000). The Raven's progressive matrices: change and stability over culture and time. Cognitive psychology, 41(1), 1-48. 29 Schady, N., Behrman, J., Araujo, M. C., Azuero, R., Bernal, R., Bravo, D., ... & Vakis, R. (2015). Wealth gradients in early cildhood cognitive development in five Latin American countries. Journal of Human Resources, 50(2), 446-463. Shonkoff, J. P., Radner, J. M., & Foote, N. (2017). Expanding the evidence base to drive more productive early childhood investment. The Lancet, 389(10064), 14-16. Surkan, P. J., Kennedy, C. E., Hurley, K. M., & Black, M. M. (2011). Maternal depression and early childhood growth in developing countries: systematic review and meta-analysis. Bulletin of the World Health Organization, 89, 607-615. UNICEF. (2017). Overview: MELQO: Measuring Early Learning Quality and Outcomes. UNESCO Publishing. Veena, S. R., Gale, C. R., Krishnaveni, G. V., Kehoe, S. H., Srinivasan, K., & Fall, C. H. (2016). Association between maternal nutritional status in pregnancy and offspring cognitive function during childhood and adolescence; a systematic review. BMC pregnancy and childbirth, 16(1), 220. World Bank. (2018). World Development Report 2018: Learning to Realize Education’s Promise. Washington, DC: World Bank. World Health Organization. (2006). WHO child growth standards: length/height for age, weight- for-age, weight-for-length, weigh-for-height and body mass index-for-age, methods and development. World Health Organization. Zelazo, P. D. (2006). The Dimensional Change Card Sort (DCCS): A method of assessing executive function in children. Nature protocols, 1(1), 297. 30 Table 1: Results of linear regressions predicting different domains of child development from quantile/quintile dummies and parenting variables or child anthropometrics. Executive function Language Early numeracy Socio‐emotional problems No Height  No Height  No Height  No Height  Parenting Parenting Parenting Parenting  control for age  control for age  control for age  control for age Quartile gap ‐ Full sample 0.435*** 0.398*** 0.396*** 0.601*** 0.537*** 0.562*** 0.418*** 0.359*** 0.386*** ‐0.212*** ‐0.121*** ‐0.214*** Quintile gap ‐ Full sample 0.476*** 0.436*** 0.436*** 0.692*** 0.624*** 0.653*** 0.449*** 0.385*** 0.416*** ‐0.224*** ‐0.125*** ‐0.227*** … age 3 0.372*** 0.350*** 0.340*** 0.311*** 0.283*** 0.288*** 0.113** 0.083* 0.101** ‐0.183*** ‐0.107* ‐0.194*** … age 4 0.523*** 0.471*** 0.490*** 0.725*** 0.645*** 0.693*** 0.439*** 0.374*** 0.419*** ‐0.234*** ‐0.148** ‐0.234*** … age 5 0.544*** 0.505*** 0.481*** 1.120*** 1.024*** 1.043*** 0.874*** 0.774*** 0.788*** ‐0.261*** ‐0.120 ‐0.256*** Notes. Table only shows parameter estimates for the highest quantile (first row) or highest quintile (row 2-5). Column ‘No control’ shows results using only quantile/quintile dummies as explanatory variables. Columns ‘Parenting’ add the three parenting composite scores to explanatory variables. Columns ‘Height for age’ add height for age z-scores to control variables. Dependent variables are standardized residuals of a linear regression of original test scores on the first three polynomials of age. Two-tailed p-values for t-tests using robust standard errors clustered on the village level.*p < .1 **p < .05 ***p < .01 31 Table 2: Pattern matrix for exploratory factor analysis (iterated principal factor) with promax oblique rotation of the parenting measures Items (abbreviated) Mean SD Cognitive parenting Socio‐emotional parenting Negative parenting Read book or magazine and look at pictures 1.78 1.25 0.57 ‐0.02 ‐0.03 Play games with words 1.70 1.23 0.50 0.00 0.01 Tell stories 1.57 1.09 0.50 ‐0.01 0.01 Sing songs 2.16 1.51 0.46 0.12 ‐0.01 Play games with numbers 2.64 1.57 0.51 0.10 ‐0.03 Play active games, e.g. with running or ball 1.42 1.00 0.42 0.01 ‐0.02 Draw or paint 1.49 1.01 0.56 ‐0.07 ‐0.01 Spend time with child so he/she can feel love and care 4.18 1.20 ‐0.04 0.47 ‐0.02 Show love by hugging, kissing, caring 4.44 1.07 0.02 0.48 ‐0.08 Teach to become self‐suficient in daily routines 4.50 0.98 ‐0.01 0.46 ‐0.01 Talk to child when he/she did something wrong 4.02 1.16 ‐0.01 0.45 0.17 Comfort child when he/she feels sad 3.89 1.22 ‐0.03 0.49 0.07 Encourage, compliment when child does something 3.77 1.20 0.11 0.51 ‐0.08 Tell child that he/she makes you happy 2.66 1.50 0.21 0.30 ‐0.02 Talk about reasons when forbidding something 4.18 1.11 ‐0.01 0.40 0.14 Feel that it is hard taking care of child 3.06 1.65 0.03 0.00 0.44 Physical punishment (kicking, slapping, beating) 2.04 1.20 0.04 ‐0.04 0.56 When angry, throw anger at child by shouting 2.84 1.38 ‐0.05 0.04 0.63 Feel annoyed with what child has done 2.68 1.37 0.01 ‐0.03 0.72 Feel annoyed or angry when child cries 2.59 1.40 ‐0.01 0.05 0.59 Shout or speak loudly to child 3.00 1.41 ‐0.07 0.04 0.65 Call child dumb, lazy or similar 1.59 1.11 0.05 ‐0.12 0.37 Talk, e..g name objects, describe daily routines 3.29 1.53 0.23 0.24 0.03 Ignore when child cries 1.47 0.97 0.08 ‐0.12 0.24 Go for a walk with child 3.36 1.51 0.01 0.29 0.06 Cronbach's alpha of composite score 0.72 0.69 0.77 Notes. Factor loadings >.3 in bold are the items of the three parenting composite scores. Possible responses ranged from 1 (never or almost never) to 5 (every day or almost every day). 32 Table 3: Pairwise correlations between follow-up child test scores, child, household and caregiver control variables. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 1 Executive function 1.00 2 Language development 0.75*** 1.00 3 Early numeracy 0.66*** 0.77*** 1.00 4 Socio‐emotional problems ‐0.15*** ‐0.16*** ‐0.16*** 1.00 5 Age (years) 0.65*** 0.66*** 0.61*** ‐0.14*** 1.00 6 Female 0.06*** 0.04*** 0.04*** ‐0.05*** ‐0.01 1.00 7 Stunted ‐0.13*** ‐0.13*** ‐0.11*** 0.02 ‐0.03*** ‐0.03** 1.00 8 Household size ‐0.01 ‐0.05*** ‐0.02 ‐0.02 0.01 ‐0.01 0.02* 1.00 9 Low living standard ‐0.11*** ‐0.14*** ‐0.12*** 0.06*** ‐0.01 0.02 0.08*** ‐0.12*** 1.00 10 Caregiver age 0.01 0.03*** 0.03** 0.04*** 0.01 ‐0.01 0.02 0.05*** ‐0.15*** 1.00 11 Caregiver primary school 0.07*** 0.09*** 0.07*** ‐0.11*** ‐0.03** 0.00 ‐0.06*** ‐0.05*** ‐0.15*** ‐0.32*** 1.00 12 Caregiver non‐verbal reasoning 0.08*** 0.10*** 0.07*** ‐0.06*** ‐0.02 ‐0.00 ‐0.06*** ‐0.04*** ‐0.09*** ‐0.37*** 0.38*** 1.00 13 Cognitive parenting 0.05*** 0.11*** 0.10*** ‐0.12*** ‐0.03*** 0.01 ‐0.03** ‐0.07*** ‐0.10*** ‐0.15*** 0.26*** 0.18*** 1.00 14 Socio‐emotional parenting ‐0.03** ‐0.02 ‐0.01 ‐0.04*** ‐0.11*** ‐0.01 ‐0.03** ‐0.06*** ‐0.06*** ‐0.03** 0.11*** 0.08*** 0.38*** 1.00 15 Negative parenting ‐0.07*** ‐0.07*** ‐0.08*** 0.35*** ‐0.04*** ‐0.10*** 0.00 0.02 0.05*** ‐0.04*** ‐0.05*** ‐0.01 ‐0.05*** 0.14*** Notes. Table shows pairwise correlations. Two-tailed p-values: *p < .1 **p < .05 ***p < .01 33 Table 4: Results of hierarchical regression analysis predicting different domains of child development from parental stimulation, nutrition and other household, caregiver and child variables. Notes. Changes in the R² (R² and F‐tests for joint significance of additional variables F refer to a comparison of (2) vs. (1) in column (2), (3) vs. (2) in column (3) and (4) vs. (2) in column (4). Two-tailed p-values for t-tests and F-tests using robust standard errors clustered on the village level. *p < .1 **p < .05 ***p < .01 34 Table 5: Results of hierarchical regression analysis predicting different domains of child development using baseline controls, a cognitive parenting x stunted and a cognitive parenting x caregiver education interaction term. Notes. Changes in the R² (R² and F‐tests for joint significance of additional variables F refer to a comparison of (2) vs. (1) in column (2), (3) vs. (2) in column (3) and (4) vs. (3) in column (4). Two-tailed p-values for t-tests and F-tests using robust standard errors clustered on the village level. *p < .1 **p < .05 ***p < .01 35 Table 6: Results of hierarchical regression analysis predicting different domains of child development using a cognitive parenting x stunted and cognitive parenting x caregiver education interaction term. Notes. Changes in the R² (R² and F‐tests for joint significance of additional variables F refer to a comparison of (2) vs. (1) in column (2) and (3) vs. (2) in column (3). Two-tailed p-values for t-tests and F-tests using robust standard errors clustered on the village level. *p < .1 **p < .05 ***p < .01 36