Internalizing problems, which include symptoms of depression and anxiety, are among the most prevalent mental health difficulties (Epkins & Heckler, Reference Epkins and Heckler2011). Internalizing problems can emerge as young as early childhood, with the incidence increasing from childhood through adolescence (Brock & Kochanska, Reference Brock and Kochanska2015; Cyranowski, Frank, Young, & Shear, Reference Cyranowski, Frank, Young and Shear2000; Gilliom & Shaw, Reference Gilliom and Shaw2004; Hollenstein & Lougheed, Reference Hollenstein and Lougheed2013). Internalizing symptoms tend to present similarly in children and adolescents, and the same criteria are used to diagnose internalizing disorders in both developmental age ranges (McLellan & Hudson, Reference McLellan, Hudson, Goldstein and DeVries2017; Powell, Ocean, & Stanick, Reference Powell, Ocean, Stanick, Goldstein and De Vries2017). Several family factors are associated with increased risk for youths’ internalizing problems, including family emotional climate, parental socialization of emotion, and relationship quality (Epkins & Heckler, Reference Epkins and Heckler2011; Lougheed, Reference Lougheed2020; Morris, Cui, Criss, & Simmons, Reference Morris, Cui, Criss, Simmons, Cole and Hollenstein2018). Conflictual relationships between parents and children specifically may increase the risk for children's internalizing problems because of their negative implications for children's emotional development (Smith, Nelson, & Adelson, Reference Smith, Nelson and Adelson2019). It is important to understand how associations between these aspects of parental relationships and children's internalizing problems unfold over time, for example, whether relationship characteristics and internalizing problems are mutually influential, or whether one tends to precede the other. The goal of this study was to disentangle the longitudinal associations between one aspect of parent–child relationships (i.e., parent–child conflict) and children's internalizing problems during mid-childhood.
Multiple theoretical approaches such as developmental psychopathology (Cicchetti & Cohen, Reference Cicchetti and Cohen2006), dynamic systems perspectives (Lougheed, Reference Lougheed2020), and interpersonal theories (Epkins & Heckler, Reference Epkins and Heckler2011; Joiner, Coyne, & Blalock, Reference Joiner, Coyne, Blalock, Joiner and Coyne1999) point to the importance of parent–child relationships in the emergence and maintenance of internalizing problems. Common among these perspectives is the idea that the interactions between parents and children shape children's ability to manage emotions, which in turn (and in combination with other factors such as genetic propensity and event exposures) is associated with the likelihood that children develop internalizing problems. The developmental psychopathology perspective (Cicchetti & Cohen, Reference Cicchetti and Cohen2006) is especially helpful when thinking about these complex developmental processes. This perspective emphasizes that gaining a complete understanding of psychopathology involves (a) research on both typical and atypical development, (b) an emphasis on individual differences, given that similar developmental outcomes may be reached via different trajectories and that similar trajectories can lead to different outcomes, and (c) disentangling complex associations between multiple levels of analysis (e.g., individual, social; Cicchetti & Cohen, Reference Cicchetti and Cohen2006; Cicchetti & Rogosch, Reference Cicchetti and Rogosch1996). We draw on the developmental psychopathology approach to frame our examination of longitudinal, bidirectional associations between parent–child conflict and typically developing children's internalizing problems, accounting for individual differences.
Associations between parent–child conflict and children's internalizing problems
In middle childhood (approximately age 6 years until the onset of puberty, although specific definitions of this developmental period vary across research groups), the maturing capacity for self-regulation enables children to be less reliant on primary caregivers than they were in early childhood (Bosmans & Kerns, Reference Bosmans and Kerns2015). One implication of this maturation is that children may start to experience the early stages of autonomy development (Stuart Parrigon & Kerns, Reference Stuart Parrigon and & Kerns2016). In day-to-day life, parent–child conflicts are an important context for children's autonomy development, are a typical part of family life, and can be a feature of healthy relationships when family members assert their perspectives while working towards solutions (Laursen & Hafen, Reference Laursen and Hafen2010). Parent–child conflict tends to peak in early adolescence (i.e., ages 10 to 13 years of age; Smetana, Campionne-Barr, & Metzger, Reference Smetana, Campionne-Barr and Metzger2006), and increases in conflict frequency and intensity may begin to manifest in middle childhood.
Conflicts in the context of supportive relationships that have a constructive outcome (e.g., greater mutual understanding, resolution of problems) may be associated with beneficial outcomes for individuals, such as improved relationship skills and ability to set boundaries (Laursen & Hafen, Reference Laursen and Hafen2010). However, some patterns of parent–child conflicts are associated with relationship difficulties and children's internalizing problems. For example, relationship quality may be relatively low in families where mother–child and father–child conflicts are frequent, intense, coercive, and unresolved (Laursen & Hafen, Reference Laursen and Hafen2010; Sheeber, Davis, Leve, Hops, & Tildesley, Reference Sheeber, Davis, Leve, Hops and Tildesley2007). Such families with conflictual relationships may tend to have poor communication and problem-solving skills, which have been associated with youths’ internalizing problems (Sheeber et al., Reference Sheeber, Davis, Leve, Hops and Tildesley2007). It is important to disentangle the direction of association between parent–child conflict and internalizing problems, and there are several hypotheses about the direction of effects. One hypothesis is that conflict contributes to children's internalizing problems. Conflictual relationships may contribute to the risk for children's internalizing problems, as a hostile family environment may trigger or exacerbate internalizing problems (Morris, Silk, Steinberg, Myers, & Robinson, Reference Morris, Silk, Steinberg, Myers and Robinson2007). A second hypothesis is that children's internalizing problems contribute to greater conflict in the parent–child relationship. Internalizing problems can include irritability, which can lead to greater likelihood of conflictual interpersonal interactions (Joiner et al., Reference Joiner, Coyne, Blalock, Joiner and Coyne1999). A third hypothesis is a combination of the first two, and is that the effects between children's internalizing problems and conflict are bidirectional (Lougheed, Reference Lougheed2020)—increased child irritability associated with internalizing symptoms may increase the risk for parent–child conflict and, in turn, parent–child conflict can reinforce children's internalizing problems.
A number of studies have attempted to address the question of directional associations between parent–child conflict and children's internalizing problems. Much of this research has focused on parent–adolescent relationships spanning youths' ages from 11 to 17 years. Findings from research on parent–adolescent conflict are mixed and likely vary by the research method used. Some studies have shown that mother–adolescent conflict (Allen et al., Reference Allen, Insabella, Porter, Smith, Land and Phillips2006) and parent–adolescent conflict (both mothers and fathers; Goodman, Samek, Wilson, Iacono, & McGue, Reference Goodman, Samek, Wilson, Iacono and McGue2019; Samek et al., Reference Samek, Wilson, McGue and Iacono2018) longitudinally predicts adolescents' internalizing problems. Other studies have shown that adolescent internalizing problems and conflictual relationships with parents generally (Brière, Archambault, & Janosz, Reference Brière, Archambault and Janosz2013) and with both mothers and fathers specifically (Lewis, Collishaw, Thapar, & Harold, Reference Lewis, Collishaw, Thapar and Harold2014) predict each other longitudinally, although the link from internalizing problems to mother–adolescent and father–adolescent relationship difficulties may be stronger than the link from relationship difficulties to internalizing problems (Branje, Hale, Frijns, & Meeus, Reference Branje, Hale, Frijns and Meeus2010; Epkins & Heckler, Reference Epkins and Heckler2011).
Family factors associated with increased risk for youths’ internalizing problems differ between childhood and adolescence (Collins & Laursen, Reference Collins and Laursen2004), and thus may show different associations with internalizing problems at different developmental stages. A few studies have examined the associations between parent–child conflict and internalizing problems in childhood. A longitudinal study of observed mother–child and father–child relationship quality and children's internalizing problems from age 4.5 to 10 years showed that negative aspects of the mother–child relationship (e.g., conflictual interactions) preceded increases in children's internalizing problems (Brock & Kochanska, Reference Brock and Kochanska2015). Yet another longitudinal study showed reciprocal associations between parent-reported warmth and harsh punishment (as reported by primarily biological mothers in the sample) and child-reported depressive symptoms from ages 7 to 12 years, although the parent–child relationship was a stronger predictor of depressive symptoms than were depressive symptoms of the parent–child relationship (Hipwell et al., Reference Hipwell, Keenan, Kasza, Loeber, Stouthamer-Loeber and Bean2008). Evidence from another longitudinal study spanning 4th to 8th grade showed that child internalizing problems were associated with an increase in conflictual parent–child relationships, which in turn contributed to the stabilization of internalizing problems in adolescence (Yong, Fleming, McCarty, & Catalano, Reference Yong, Fleming, McCarty and Catalano2014). Taken together, the body of research spanning childhood and adolescence shows mixed evidence regarding the direction of associations between parent–child conflict and internalizing problems, although the evidence to date suggests that relationship issues predict internalizing problems more strongly than the other way around (see Epkins & Heckler, Reference Epkins and Heckler2011 for a review; Yong et al., Reference Yong, Fleming, McCarty and Catalano2014). In addition, the majority of work on this topic has been conducted on adolescent samples, which may not generalize to childhood samples given the confluence of biological, psychological, and social developments that is unique to adolescent development (Hollenstein & Lougheed, Reference Hollenstein and Lougheed2013).
Disentangling within-family and between-family associations
Most of the studies on the longitudinal associations between parent–child relationship difficulties and children's internalizing problems to date have used a cross-lagged panel model (CLPM). This analytic approach is an intuitive application of path modeling that involves simultaneously predicting two dependent variables (in this case conflict and internalizing problems) from the same variables at the previous time point. The analysis yields an autoregressive effect for each variable (e.g., conflict at time t being predicted from conflict at time t − 1) and a directional (cross-lagged) effect for each variable from the other variable at the previous time point (e.g., conflict at time t being predicted from internalizing problems at time t − 1).
Recently, CLPMs have been criticized because cross-lagged effects are presumed and interpreted to index within-person (or in this case within-family) processes. However, these models can lead to grossly mis-specified cross-lagged paths because examination of within-family change processes requires accounting for the trait-like stability of a construct, and autoregressive paths on their own do not adequately capture trait-like stability in a construct (Hamaker, Kuiper, & Grasman, Reference Hamaker, Kuiper and Grasman2015). Failure to adequately control for stable between-family differences can inflate cross-lagged and autoregressive paths and can even change the sign of the cross-lagged paths. Adding a random intercept to CLPMs corrects this bias by explicitly modeling stable between-family differences in each construct (conflict and youth internalizing) and the association of the two stable between-family factors in addition to the autoregressive and cross-lagged paths. Thus, cross-lagged paths are interpreted as time-specific within-family associations in the random-intercept CLPM (RI-CLPM).
To date, none of the research studies examining longitudinal associations between parent–child conflict and internalizing problems with CLPMs have used random intercepts. Therefore, existing literature on this topic (Branje et al., Reference Branje, Hale, Frijns and Meeus2010; Brière et al., Reference Brière, Archambault and Janosz2013; Goodman et al., Reference Goodman, Samek, Wilson, Iacono and McGue2019; Samek et al., Reference Samek, Wilson, McGue and Iacono2018; Yong et al., Reference Yong, Fleming, McCarty and Catalano2014) may overestimate directional and/or transactional associations between conflict and internalizing problems. It is possible that these associations could be explained, at least in part, by between-family differences that persist over time rather than being within-family processes that unfold over time. It should be noted, however, that some studies to date have used other models to separate within- and between-family associations, such as multilevel models (e.g., Brock & Kochanska, Reference Brock and Kochanska2015), generalized estimating equations (e.g., Hipwell et al., Reference Hipwell, Keenan, Kasza, Loeber, Stouthamer-Loeber and Bean2008), and some behavioral genetic designs (e.g., Lewis et al., Reference Lewis, Collishaw, Thapar and Harold2014). Based on the findings from these studies, we may expect that cross-lagged associations from parent–child conflict to internalizing problems may be stronger than reciprocal associations, and that there will likely be a strong stable between-family correlation of the random intercepts (given that there is child-driven genetic covariation across these constructs; Lewis et al., Reference Lewis, Collishaw, Thapar and Harold2014).
The role of gender
Gender differences in internalizing problems are well established, with girls and women being more likely to experience them than boys and men (Cyranowski et al., Reference Cyranowski, Frank, Young and Shear2000). This gender difference typically first emerges in adolescence and remains stable into adulthood (Holsen, Kraft, & Vittersø, Reference Holsen, Kraft and Vittersø2000). Interpersonal perspectives suggest that gender differences in internalizing problems are related to gender differences in affiliative relationship patterns (e.g., girls' greater need for close emotional communication, intimacy, and responsiveness; Cyranowski et al., Reference Cyranowski, Frank, Young and Shear2000). Girls may be more likely to respond to relationship stressors in a way that exacerbates internalizing problems (Goodman et al., Reference Goodman, Samek, Wilson, Iacono and McGue2019; Rudolph, Flynn, & Abaied, Reference Rudolph, Flynn, Abaied, Abela and Hankin2008). For example, if girls respond to conflicts with cognitive processes that up-regulate rather than resolve negative emotions (i.e., rumination), the emotion difficulties associated with internalizing problems may be reinforced (Nolen-Hoeksema, Reference Nolen-Hoeksema2012). It may be these gendered affiliative tendencies, together with the unique biopsychosocial transitions of female adolescents and experience of negative life events, trigger the onset of internalizing problems in adolescence (Cyranowski et al., Reference Cyranowski, Frank, Young and Shear2000). However, evidence on gender differences in the longitudinal associations between relationship difficulties and internalizing problems is mixed. Some studies have found no significant differences between adolescent boys and girls on the longitudinal links between relationship difficulties (e.g., conflict) and internalizing problems (Goodman et al., Reference Goodman, Samek, Wilson, Iacono and McGue2019), whereas others suggest that girls may show a stronger association between maternal and paternal relationship difficulties and internalizing problems than boys (Finan, Ohannessian, & Gordon, Reference Finan, Ohannessian and Gordon2018; Lewis et al., Reference Lewis, Collishaw, Thapar and Harold2014).
Less is known about the role of gender in internalizing symptoms in childhood. One longitudinal study found no gender differences in the mediating effects of parent–child conflict from late childhood to early adolescence (Yong et al., Reference Yong, Fleming, McCarty and Catalano2014). There is a dearth of research exploring gender differences in the longitudinal associations among children's internalizing problems and difficulties in the parental relationship. Research in this area is needed because boys and girls show gender differences in affiliative tendencies in childhood (e.g., girls being more likely than boys to disclose emotions in relationships and to seek social support; Zeman & Shipman, Reference Zeman and Shipman1997), and internalizing problems may begin to increase over childhood (Brock & Kochanska, Reference Brock and Kochanska2015). Therefore, it is possible that gender plays an important role in the longitudinal associations between parent–child conflict and children's internalizing problems—but whether or not this is the case is currently unknown.
The current study
The overarching goal of the current study was to disentangle the longitudinal associations between mother–child conflict and children's internalizing problems in mid-childhood. Of the available longitudinal studies, the evidence is mixed regarding whether mother–child conflict precedes children's internalizing problems, whether children's internalizing problems precede conflictual relationships, or whether there are bidirectional effects. In addition, few studies to date have examined the extent to which mother–child conflict and children's internalizing problems reflect between-family differences versus within-family processes (for exceptions, see Brock & Kochanska, Reference Brock and Kochanska2015; Lewis et al., Reference Lewis, Collishaw, Thapar and Harold2014). It is also unknown whether longitudinal associations vary by child gender. We used data from a publicly available longitudinal study—the Study of Early Child Care and Youth Development (SECCYD)—to test the following three research questions focusing on relationships with mothers.
(1) To what extent do mother–child conflict and children's internalizing problems reflect stable between-family differences in middle childhood? In line with longitudinal studies examining both within- and between-family factors (Brock & Kochanska, Reference Brock and Kochanska2015; Lewis et al., Reference Lewis, Collishaw, Thapar and Harold2014), we expected between-family correlations such that families with higher conflictual mother–child relationships would tend to have children with higher internalizing problems.
(2) Are the longitudinal associations between mother–child conflict and children's internalizing problems directional—with one preceding the other—or are they bidirectional? Based on research to date (e.g., Epkins & Heckler, Reference Epkins and Heckler2011; Yong et al., Reference Yong, Fleming, McCarty and Catalano2014), we expected that mother–child conflict would more strongly predict children's internalizing problems over time than children's internalizing problems would predict mother–child conflict.
(3) Do the longitudinal associations between mother–child conflict and children's internalizing problems vary by child gender? We hypothesized that the longitudinal associations would be stronger for girls than for boys. We tested our research questions using RI-CLPMs (Hamaker et al., Reference Hamaker, Kuiper and Grasman2015).
Method
Participants and procedure
Data from the publicly available SECCYD were used in the present study. Detailed information on recruitment and the data achieve is available elsewhere (https://www.icpsr.umich.edu).
The SECCYD aimed to investigate various childcare factors related to health, academic, behavioral, emotional, and other developmental outcomes from infancy to middle adolescence. Participants were recruited in 1991 at 10 locations in the USA (Little Rock, AR; Irvine, CA; Lawrence, KS; Boston, MA; Hickory, NC; Philadelphia, PA; Pittsburgh, PA; Charlottesville, VA; Seattle, WA; Madison, WI) through hospital visits. During selected 24-hr intervals among the selected hospitals, all women giving birth (N = 8,986) were screened for eligibility for the research. Some families (N = 3,142) were excluded due to criteria such as not being proficient in English or planning to move within the next 3 years. An additional 1,353 families could not be reached or refused to participate in the study at a follow-up phone interviews at 2 weeks. When children were 1 month old, home visits were completed among eligible families (N = 1,364) that were randomly chosen from the pool of remaining participants and data collection continued in four phases. At Phase I of the SECCYD (birth to 3 years), 1,364 children participated; at Phase III (2nd through 6th grade), the retention rate was high (77.79%) with 1,061 children participating. At recruitment, the highest completed education level was high school for 26% of mothers, 22% of mothers were from ethnic/racial minority backgrounds (i.e., non-White), and 21% of mothers had incomes lower than 200% of the poverty line. For more details on sample characteristics, see the report by the NICHD Early Child Care Research Network (2005).
In the current study, approximately 51% of the children were male. Mothers’ mean number of years of education was 14.23 years at the time they were recruited into the study. Approximately 20% of the participants were from minority (Hispanic, African American and Other) ethnic groups in the current sample. Descriptive statistics of participant ethnicity and other demographic variables are provided in Table 1. Ethical approval for the current study was obtained from the Purdue University Institutional Review Board (#1908022540).
Descriptive statistics

Note. For gender and race/ethnicity, the N values are the number of responses on the categorical variables. The specific dummy-coded variables reflect the proportion of the sample that falls into that category (i.e., if multiplied by 100 they would be the percentages).
Measures
Mother–child conflict
The conflict with child subscale of the Child–Parent Relationship Scale was used to assess mother–child conflict. Mothers reported mother–child conflict during home visits at 3rd and 5th grades, and during lab visits at 4th and 6th grades. The Child–Parent Relationship Scale was adapted from the Student–Teacher Relationship Scale (Pianta, Reference Pianta1994, Reference Pianta2001), and captured children's attachment behaviors through parent-report. Mothers responded to the items on a 5-point Likert-type scale from 1 (definitely does not apply) to 5 (definitely applies). Example items from the conflict with child subscale include “My child and I always seem to be struggling with each other,” and “My child easily becomes angry at me.” Higher scores indicated higher conflict. The Cronbach alpha of the conflict with child subscale ranged between .83 and .85 through 3rd to 6th grade.
Internalizing problems
Mothers reported their children's internalizing problems during home visits at 3rd and 5th grades, and during lab visits at 4th and 6th grades. The internalizing problems subscale of the Child Behavior Checklist (CBCL) (Achenbach, Reference Achenbach1991) was used. The CBCL is a valid (Achenbach & Rescorla, Reference Achenbach and Rescorla2000) and widely used instrument to measure social competence and problem behaviors among children aged between 4 and 18 years. Children's behaviors (e.g., “Unhappy, sad, or depressed”) were scored on 3-point Likert-type scales from 0 (not true of the child) to 2 (very true of the child) (Achenbach, Reference Achenbach1991). Higher scores indicate greater internalizing problems.
Child, maternal, and family covariates
In addition to the primary measures of interest above, several aspects of child, mother, and family status were measured to include in the analyses. Gender (1 = male, 0 = female) and ethnicity (dummy codes for Black, Hispanic, and Other with White as the reference group) were included as child covariates. Maternal education, maternal vocabulary, and maternal depression were included as maternal covariates. Maternal education was assessed by years of completed schooling when the child was 1 month old. Maternal vocabulary was measured by the Peabody Picture Vocabulary Test (Revised) when the child was 3 years old (Dunn & Dunn, Reference Dunn and Dunn1981). Maternal depression was assessed by the Center for Epidemiological Studies Depression Scale (Radloff, Reference Radloff1977) by mother-report during the home visits at 3rd grade. Although the items were rated on 4-point Likert-type scale from 1 (less than once a week) to 4 (5–7 days a week), they were transformed to a scale of 0–3 to match the original scale. An example item was “I thought my life had been a failure.” The Cronbach alpha of the scale was .90. Family income was reported during the home visits at 3rd grade and included as a family-level covariate. It was measured by the income-to-needs ratio (i.e., where 1 represents the federal poverty level, values less than 1 represent families living below the poverty level, and values greater than 1 represent families living above the poverty level).
Analytic plan
We conducted data management and descriptive statistics using Stata 15 (StataCorp, 2017) and all structural equation models using Mplus 8 (Muthén & Muthén, Reference Muthén and Muthén2017). Next, we conducted a series of model-building steps to ensure an appropriate base model for hypothesis testing. First, we ran measurement models with random intercepts and autoregressive paths independently for mother–child conflict and child internalizing problems to determine whether autoregressive paths should be constrained to be equal across time or freely estimated. Next, we combined the selected models for mother–child conflict and child internalizing problems into a single RI-CLPM without any covariates to examine model fit. Finally, our primary model used to test the research questions included the single RI-CLPM specification with covariates predicting the random intercepts. Covariates included maternal vocabulary, education, and depression (from the 3rd grade wave), race/ethnicity, family income-to-needs ratio (from the 3rd grade wave), and gender.
Hypothesis 1 (that there would be stable between-family correlations of higher conflictual mother–child relationships with higher child internalizing problems) was tested by the included correlation between the two random intercepts.
Hypothesis 2 (that mother–child conflict would more strongly predict children's internalizing problems over time than children's internalizing problems would predict mother–child conflict) was tested first by observing whether the cross-lagged paths in each direction were significant. If paths were significant in one or both directions, differences in magnitude were tested by constraining them to equality in a nested model and observing whether there was a decrement in model fit in doing so.
Next, we conducted a multiple group analysis consisting of two models in order to explore gender differences in order to test hypothesis 3 (that the longitudinal associations would be stronger for girls than for boys). Model 1 was an unconstrained model: our primary model—excluding gender as a covariate—with random intercepts, factor loadings, and paths freely estimated for males and females. The only estimates constrained equal between males and females were the time-specific intercepts for mother–child conflict and child internalizing problems (this was required for the model to converge in Mplus). Model 2 was a nested model that constrained factor loadings and path estimates for males and females to be equal. A decrement in model fit (a significant chi-square test), comparing Model 2 to Model 1 would indicate the presence of gender differences; a non-significant chi-square would indicate no decrement in model fit.
Finally, as sensitivity analyses to provide a direct comparison with most of the literature, we ran two additional CLPMs without random intercepts. One model included one lag (the autoregressive path from the prior time point only, t − 1) and the other model included two lags (the autoregressive path from the prior time point, t − 1, and the autoregressive path from two time points before, t − 2). Differences in the model conclusions are presented in the Results section of this paper and interpreted in the Discussion section.
Missing data and attrition
Missing data was handled by using full information maximum likelihood (FIML) in all analyses. FIML uses all available information in the model and does not restrict analyses to cases with no missing data, leading to less biased estimates than list-wise deletion and performing similarly to multiple imputation methods (Acock, Reference Acock, Cooper, Camic, Long, Panter, Rindskopf and Sher2012). Relatively little data on our primary constructs were lost due to attrition, given the scope of the study as birth cohort design. Of the original 1,364 participants, 1,090 (80%) provided at least one observation on mother–child conflict or internalizing problems (911 provided complete data on both constructs from 3rd through 6th grade). Our primary model was run on the subsample of 1,090 that provided at least one observation on mother–child conflict or internalizing, and the results were nearly identical and came to the same conclusions (results available by request). Gender and race/ethnicity experienced no missing data. The number of observations for all variables included in the analyses are reported in Table 1.
Results
Descriptive statistics and correlations
In general, children's levels of internalizing problems and mother–child conflict were stable from 3rd to 6th grade (see Table 1). They were also moderately correlated across time points (r values ranging from .34 to .47 and p < .001, see Table 2). Of the covariates, maternal depression at 3rd grade was most consistently correlated with children's internalizing and mother–child conflict, followed by family income-to-needs ratios and maternal education.
Correlation matrix

Note. * p < .05, ** p < .01, *** p < .001.
Model building
All the model-fit statistics are presented in Table S1. For the measurement models, chi-squared difference tests suggested that the model with constrained autoregressive paths fitted significantly worse for mother–child conflict, but the model with constrained estimated autoregressive paths did not provide a worse fit to the data for child internalizing problems (see Table S1). Thus, the two measurement models (including freely estimated autoregressive paths for mother–child conflict, but constrained autoregressive paths for child internalizing problems) were combined into a RI-CLPM that fitted the data well. The primary model, additionally including covariates, also fitted the data very well (Table S1). In the multiple group analysis exploring gender differences, constraining estimates for males and females did not result in a decrement in model fit (Table S1), contrary to the third hypothesis, and thus results from the primary model—the most parsimonious model—are presented.
Longitudinal associations between internalizing and mother–child conflict
The results of the RI-CLPM are presented in Figure 1. As expected in the first hypothesis, the random intercepts between children's internalizing problems and mother–child conflict were moderately correlated (r = .49, p < .001), suggesting that relatively stable between-family differences account for part of the association between the two constructs. The constructs also had significantly correlated errors at each time point, although these were weaker in magnitude than the correlation between random intercepts. These correlated errors suggest that time-specific variations in mother–child conflict and child internalizing problems are related to each other (r = .12–.28, p < .01), such that positive deviations in children's internalizing problems are associated with positive deviations in mother–child conflict (and vice versa) at each specific wave (i.e., beyond relative stability in both constructs). Although autoregressive paths for both mother–child conflict and child internalizing problems were significant, these paths were relatively small in magnitude.
Random intercept cross-lagged panel model (CLPM) for internalizing and conflict. Only significant paths from covariates to random intercepts are shown. Cross-lagged paths included, but are not significant (shown as grey). * p < .05, ** p < .01, *** p < .001.

Our second hypothesis was that mother–child conflict would more strongly predict children's internalizing problems over time than children's internalizing problems would predict mother–child conflict. Contrary to expectations, no significant cross-lagged paths emerged. Specifically, prior time-point-specific variation in child internalizing problems did not predict subsequent time-point-specific variation in mother–child conflict, or vice versa.
Gender differences
Model 1, with paths freely estimated for males and females, is presented in the supplementary material (Figure S1) as a test of our third hypothesis that longitudinal associations would be stronger for females than males. Although some differences emerged in which associations were statistically significant, the majority of associations were similar in magnitude between males and females, which was contrary to our expectations. In addition, as noted above, in the group analyses, the model with constrained estimates between males and females estimated the data very well and did not result in a decrement in model fit relative to the unconstrained model, and was thus preferred to the model that freely estimated associations (see Table S1 in the supplementary material). Notably, in the model with freely estimated cross-lagged paths for each gender, none reached statistical significance, which is consistent with conclusions drawn from the primary model.Footnote 1
Covariates
A few covariates emerged as significant predictors of the random intercepts (see Table S2 in the supplementary material for all unstandardized coefficients, standard errors, and p-values from the model). Maternal depressive symptoms significantly predicted children's internalizing problems (β =0.37, p < .001) and mother–child conflict (β = 0.31, p < .001). Black children, compared to White children, had significantly fewer internalizing problems (β = −0.06, p = .04) and significantly less mother–child conflict (β = −0.09, p < .001). Families with higher income-to-needs ratios had children with significantly fewer internalizing problems (β = −0.11, p = .02), although higher maternal education was associated with significantly more internalizing problems (β = 0.08, p = .01).
Sensitivity analyses with CLPMs
Finally, we also estimated two CLPMs with no random intercepts: one using autoregressive paths (see Figure S2) with a single lag of one time point (t − 1) and the other using autoregressive paths with two lags (t − 1 and t − 2; see wure S3). The model with a single lagged path provided a significantly worse fit than the model with two lagged paths (see Table S1). In the model including two lags, all autoregressive paths from t − 1 and t − 2 were statistically significant. Statistically significant bidirectional paths also emerged. Internalizing problems at 3rd grade significantly predicted mother–child conflict at 4th grade (β = 0.07, p < .001); 4th grade mother–child conflict significantly predicted 5th grade internalizing (β = 0.06, p = .033); and 5th grade internalizing significantly predicted 6th grade mother–child conflict (β = 0.09, p = .001).
Discussion
Our primary aim in this study was to disentangle the longitudinal associations between mother–child conflict and children's internalizing problems in middle childhood. We examined the extent of between-family differences and within-family processes in these longitudinal associations by incorporating random intercepts into our CLPM analytic approach. We found that mother–child conflict and children's internalizing problems were relatively stable between 3rd and 6th grade, that the relatively stable factors of each construct were positively correlated, and that their error terms were positively correlated at each time point, which were in line with expectations. However, contrary to expectations, we found no cross-lagged associations between mother–child conflict and children's internalizing problems (i.e., no significant bidirectional paths). In addition, we did not observe gender differences in the longitudinal associations between mother–child conflict and children's internalizing symptoms, which was also contrary to our expectations.
Disentangling within- and between-family associations: considering time scale
The idea that developmental processes at multiple levels of analysis (e.g., individual and social levels) are transactional is central to the developmental psychopathology approach (Cicchetti & Cohen, Reference Cicchetti and Cohen2006). Some studies on parent–child relationships and children's internalizing problems provide evidence for such bidirectional associations (e.g., Hipwell et al., Reference Hipwell, Keenan, Kasza, Loeber, Stouthamer-Loeber and Bean2008; Yong et al., Reference Yong, Fleming, McCarty and Catalano2014); however, the majority of studies on this topic to date do not account for between-family differences in their tests of within-family processes. As a consequence, estimates of directional associations between mother–child conflict and children's internalizing symptoms may be inaccurate (Hamaker et al., Reference Hamaker, Kuiper and Grasman2015). The results of our study highlight the differences in conclusions drawn when random intercepts are incorporated into CLPMs versus when they are not. Had we not accounted for stable between-family differences with random intercepts, we would have concluded from the data that, in middle childhood, children's internalizing problems predict subsequent mother–child conflict, which in turn predicts greater children's internalizing problems. However, these directional associations were not found in the model incorporating random intercepts, which suggests that the bidirectional effects found in previous studies using CLPMs may be influenced by relatively stable between-family differences that were not accounted for.
It is important to note that both the RI-CLPM and the CLPM with two lags fitted the data well, but the CLPM with one lag did not. Well-fitting models are not proof of their accuracy (Tomarken & Waller, Reference Tomarken and Waller2003), but the key conceptual distinctions should be considered between the RI-CLPM and CLPM with two lags in comparing them. In the RI-CLPM, the assumption is that unobserved factors contribute to relative stability in both constructs across time (potentially due to stable environments or individual characteristics). In the CLPM with two lags, stability in the constructs is generated by both t − 1 and t − 2 autoregressive paths. The causal inferences for this model would imply that both mother–child conflict and internalizing are not only temporally related to subsequent mother–child conflict and internalizing (t − 1), but they also have delayed sleeper effects that operate independently of the intermediary behaviors (t − 2) and covariates. Thus, the RI-CLPM aligns better with current theory on developmental influences (i.e., both proximal construct-specific and broader stable influences), and the CLPM with two lags is largely inconsistent with most theoretical explanations.
There are several possible explanations of the null directional effects in the results of our RI-CLPM. One is that there is a true lack of bidirectional associations between mother–child conflict and children's internalizing problems in middle childhood. Another explanation is that the lack of significant associations are type II errors. Given these possible explanations, and in the context of our overall results, the lack of directional associations observed in the current study may also have to do with the time scale at which the data were collected. Our results showed that the error terms between mother–child conflict and children's internalizing problems covary over time—as one increases in a given year, so does the other (independent of relative stability), and vice versa. The extent of conflict in mother–child relationships can vary over time scales of months, weeks, days, and hours (Lougheed, Reference Lougheed, Kunnen, van der Gaag, de Ruiter-Wilcox and Jeronimus2019, Reference Lougheed2020). We speculate that, while the current study does not provide evidence that mother–child conflict and children's internalizing problems are bidirectionally associated over yearly development, potential causal bidirectional effects may manifest on shorter time scales (e.g., weeks, months). It will be important for future research to incorporate multiple time scales into studies such as with multiple burst designs (Ram & Diehl, Reference Ram, Diehl, Diehl, Hooker and Sliwinski2015). Multiple burst designs involve “bursts” of intensive longitudinal assessments at short time scales (e.g., behavioral observations of mother–child conflicts, daily experience sampling of conflict) that are repeated over longer time scales such as months or years (Lougheed, Reference Lougheed2020). This type of design would enable researchers to more thoroughly disentangle the direction of associations between conflict and internalizing problems that we observed in their correlated error terms. This approach would be more closely aligned to the theoretical tenets of developmental psychopathology that individual development is embedded in social processes at multiple levels (Cicchetti & Cohen, Reference Cicchetti and Cohen2006; Granic, Reference Granic2005).
The role of gender
Gender differences in socioemotional development can manifest in middle childhood (Chaplin & Aldao, Reference Chaplin and Aldao2013). It was thus contrary to our expectations that we did not observe significant gender differences in the associations between mother–child conflict and children's internalizing problems. Previous research has shown that pubertal development is linked to both parent relationships (Marceau, Ram, & Susman, Reference Marceau, Ram and Susman2015; Paikoff & Brooks-Gunn, Reference Paikoff and Brooks-Gunn1991) and adolescent internalizing problems (Cyranowski et al., Reference Cyranowski, Frank, Young and Shear2000). We speculate that any gender differences in the longitudinal associations between mother–child conflict and children's internalizing problems, if they exist, may manifest around the same time that gender differences in internalizing problems become apparent in adolescence (Cyranowski et al., Reference Cyranowski, Frank, Young and Shear2000).
Limitations and future directions
Our results should be interpreted in the context of a few limitations of the current study.
First, mother–child conflict and children's internalizing symptoms were measured once per year; future studies should use designs with a shorter time scale, as discussed above.
Second, it will be important that future research is conducted longitudinally from middle childhood through adolescence to gain a clearer understanding of if, and when, conflictual relationships and internalizing problems become linked. Longitudinal designs spanning the years before, during, and after puberty will have the greatest potential to determine if and when gender is associated with longitudinal links between mother–child conflicts and children's internalizing problems.
Third, all the measures used in the current study were mother-report questionnaires; shared method variance could thus explain some of the findings. In this study, the sole use of mother-report questionnaires may have inflated the correlations between the random intercepts. Our results are somewhat inconsistent with other studies on bidirectional associations between negative aspects of parent–child relationships and children's internalizing problems that have either obtained data from multiple reporters or examined observed relationship characteristics with multiple family members. For example, one study found bidirectional associations between mother-reported negative parenting and subsequent daughter-reported depressed mood (Hipwell et al., Reference Hipwell, Keenan, Kasza, Loeber, Stouthamer-Loeber and Bean2008). Another study showed that observed negativity (e.g., conflictual interactions) in mother–child but not father–child relationships was associated with subsequent increases in children's internalizing symptoms (Brock & Kochanska, Reference Brock and Kochanska2015). Future research on parent–child conflict and children's internalizing problems should incorporate measures from multiple reporters to reduce the potential influence of shared method variance.
There are several other important directions for future research. One is to examine children's externalizing symptoms in relation to parent–child conflict and children's internalizing symptoms. Internalizing and externalizing symptoms can be comorbid in middle childhood, and children's externalizing problems have been associated with parent–child conflicts (Granic & Lougheed, Reference Granic, Lougheed, Dishion and Snyder2016). Another future direction is to simultaneously examine change over time in parental symptoms (e.g., internalizing, externalizing). In our study, we found that maternal depression was positively associated with stable between-family differences in both child internalizing symptoms and parent–child conflict. However, parental symptoms may vary over time in ways that predict both children's symptoms and parent–child conflict (Lougheed, Reference Lougheed2020). It is also important for future research to examine how conflict and internalizing problems are associated in other family relationships (e.g., children's relationships with fathers, non-biological parents, and other caregivers) and in more diverse samples.
Conclusion
The results of this study show that mother–child conflict and children's internalizing problems covary longitudinally in middle childhood—as one increases or decreases on a year-to-year basis, so does the other. This suggests that both of these aspects of children's lives are related to each other, although not through cross-lagged associations at a yearly time scale when between-family differences have been accounted for. Future studies that incorporate methods that are more directly aligned with the developmental psychopathology approach (Cicchetti & Cohen, Reference Cicchetti and Cohen2006; Granic, Reference Granic2005) by incorporating multiple time scales (e.g., multiple burst designs) and that account for both within-family processes and between-family differences (i.e., random intercepts) will be powerful for disentangling the associations between family relationship processes and children's internalizing problems.
Supplementary Material
The supplementary material for this article can be found at https://doi.org/10.1017/S0954579420000863
Financial Support
A cooperative agreement (5 U10 HD027040) between the study investigators and the Eunice Kennedy Shriver National Institute of Child Health and Human Development supported the design and data collection of the Study of Early Childcare and Youth Development (SECCYD) from birth through age 15 years. The content of this paper is solely the responsibility of the authors and does not necessarily reflect the views of the NICHD or the NICHD Early Child Care Research Network.
Conflicts of interest
None.
Acknowledgements
We gratefully acknowledge the support of the Center for Families (Purdue University). Dr Marceau was supported by the National Institute on Drug Abuse (K01 DA039288).

