Participants
Data for the present study were collected in waves 1 to 3 of the Estudio Longitudinal para una Infancia Saludable (Longitudinal Study for a Healthy Childhood; [ELISA]), a prospective longitudinal study conducted in Galicia (NW Spain). Data collection started in 2017 (T1), encompassing preschool children who were born in 2011–2013, and with information provided by both parents and teachers. Only children with available data in some of the main study variables, namely psychopathic traits and conduct problems, were included in the present study (n = 2,470). Sixteen participants with an affirmed diagnosis of, or being assessed for, autism spectrum disorder were excluded, resulting in a final sample of 2,454 children (48.2% girls), aged 3 to 6 (Mage = 4.26; SD = 0.91). A total of 72 public (79.2%), charter (18.1%), and private (2.8%) schools participated in the study, which were located in predominantly working-class communities, with low diversity in terms of ethnicity (93.9% of children were Spanish). Information was collected through 2,250 parents’ reports (87.2% mothers), and 2,407 reports from preschool teachers. Regarding children’s family background, 23.7% of mothers and 39.8% of fathers completed compulsory education, 47.4% and 31.2% completed higher education, and 28.9% and 29% completed vocational training studies.
Two follow-ups were conducted within one-year intervals. The first follow-up (T2) was conducted one year later in a sample of 2,333 children (Mage = 5.35; SD = 0.92), with information provided by 1,993 parents (81.25% of the total sample) and 2,170 teachers (88.46%). The level of attrition between T1-T2 participants was 4.69% considering the total sample, 11.42% based on parent-reports and 9.85% based on teacher-reports. The second follow-up (T3) was conducted two years following the initial assessment in a sample of 2,272 children (Mage = 6.33; SD = 0.92), with information provided by 1,790 parents (72.98% of the total sample) and 2,024 teachers (82.51%). The level of attrition between T1-T3 participants was 7.38% considering the total sample, 20.44% based on parent-reports and 15.91% based on teacher-reports. Comparisons between children with complete follow-up data (i.e., participation in three waves; n = 2.218; 90.4%), children who missed one of the follow-up studies (n = 172; 7%) and children with no follow-up data (i.e., participation only in T1; n = 63; 2.6%) revealed no significant differences in terms of gender, χ² (2) = 4.92, p = .476; age F (2,450) = 0.006, p = .994, and baseline levels of CP reported by parents, F (2,213) = 0.763, p = .467. There were differences according to family’s SES, F (2235) = 13.03, p < .001, and the baseline levels of conduct problems reported by teachers, F (2409) = 4.42, p < .05, with lower levels of SES and higher levels of conduct problems for children who missed one of the follow-up studies.
Measures
Psychopathic traits. Both parents and teachers rated the 28 items of the CPTI (Colins et al., 2014) in all three waves of the study. Eight items intend to measure the interpersonal or Grandiose-deceitful (GD) psychopathy component (e.g., “Thinks that he or she is better than everyone on almost everything”), 10 items intend to measure the affective or Callous-unemotional (CU) psychopathy component (e.g., “Never seems to have bad conscience for things that he or she has done”), and 10 items intend to measure the behavioral or Impulsive-need of stimulation (INS) psychopathy component (e.g., “Provides himself or herself with different things very fast and eagerly”). The CPTI items were rated on the basis of how the child usually behaves rather than how he/she behaves at the moment, in a response scale ranging from 1 (does not apply at all) to 4 (applies very well). The optimal factor structure of CPTI was examined as part of the present study.
Conduct problems. Both parents and teachers rated The Conduct Problems Scale, composed of 10 items (e.g., “Has been very angry”, and “Has beaten, torn, shoved, kicked, or thrown something on others without a reason”) that is closely based on DSM-IV (APA, 1994) criteria of oppositional defiant disorder (ODD) and CD, and were relevant to preschool children as well as older children and adolescents (Colins et al., 2014). Items were scored using a 5-point response scale (1 = never to 5 = very often). Cronbach’s alpha (α) for the three waves ranged between .86-.88 for parent reports, and between .93-.94 for teacher reports. In line with prior work (López-Romero, Colins et al., 2022), children were classified as exhibiting stable conduct problems (CP) if they were 0.5 SD above the mean of the CP measure in T2 (4–6 years old) and T3 (5–7 years old).
Procedure
The ELISA study was approved by [BLINDED]. A total of 126 public, charter and private schools were initially contacted in order to ask for potential collaboration. The initial contacts were made by phone, and information letters were subsequently sent by email. Once the school accepted the conditions and agreed to participate, families were contacted and invited to enrol in the study via information letters and group meetings in the schools, where a member of the research lab explained the conditions of the study. An active consent form was filled out by the families (approximately 25–50% response rate per school), after which the preschool teachers could also complete the questionnaires. Preschool teachers, who handed out the information to the parents, collected the informed consents. One teacher could complete the questionnaires for as many children in his/her classroom as there were written parental consent forms. Only one parent (i.e., mother, father, or principal caregiver) was asked to complete the questionnaires. Data collections took place during the Spring to assure that teachers have spent at least six months with the child before rating the questionnaire items. In all waves of the study, participants were given one month to complete the questionnaires. After that period, reminders were sent to those who were late, firstly by the preschool teacher and then directly by the ELISA staff via email. Neither families nor teachers received any monetary compensation for their participation in the study. Nonetheless, as a reward for their participation, all the schools received a set of educational games for preschoolers in T1, whilst both families and schools participated in a draw of several sets of books and educational games, valued between 50€ and 100€, at the end of the third wave data collection (T3).
Statistical Analysis
Exploring the optimal factor structure of the CPTI. To maximize the possibility to find the most adequate factor structure of the CPTI, two exploratory statistical frameworks were used. All exploratory models were based on parent-reported data from T1. First, we estimated the partial correlation network of the CPTI items using Copula gaussian graphical model estimation implemented in the R library BGGM (missing data were handled using multiple imputation with chained equations and predictive mean matching). We then inspected which items were strongly associated (i.e., correlated), with strong item-item associations being considered an indicator of a broader dimension. Zero-order polychoric correlations were used to estimate associations among items and we pooled strongly correlated items and reconducted the correlations until no correlations above .60 emerged. A correlation of .60 was selected because it indicates a moderate to strong correlation according to most criteria.
Second, we used exploratory factor analysis (EFA) to explore possible factor structures. EFA was based on the polychoric correlation matrix and the Kaiser-Meyer-Olkin (KMO) test values were used to examine whether the items were suitable for EFA. KMO values indicate the proportion of variance in variables that might be explained by latent factors and values above .80 are considered to indicate that EFA is well suited. Bartlett’s test of sphericity was also used, where a significant test result (i.e., < .05) indicates that EFA is suitable. Horn’s parallel analysis was used to determine the number of factors to retain, and these factors were extracted using principal axis factoring and promax rotation.
Confirmatory tests of factor models. The proposed model(s) identified using the methods described above, using parent-reported data from T1, were tested with new data (parent ratings from T2 and T3; and teacher-ratings from T1, T2 and T3) using confirmatory factor analysis (CFA). Model/data fit was evaluated using the Comparative Fit Index (CFI), Root Mean Square Error of Approximation (RMSEA), Standardized Mean Square Residual (SRMR), and Tucker-Lewis fit Index (TLI). Adequate model fit is indicated by higher CFI/TLI (values > 0.90 are indicative of adequate fit and values above .95 of good fit), and lower RMSEA and SRMR (values < 0.06 and 0.08, respectively, indicate good fit) (Schermelleh-Engel et al., 2003). Model fit of all models was contrasted with the fit of the original 3-factor CPTI model (Colins et al., 2014). CFAs were run using the R library lavaan and because of the ordinal response scale, diagonally weighted least squares estimation and scaled fit indices were used and examined.
Internal Structure of the CPTI Factors and Associations with Conduct Problems.
When the best fitting factor model had been identified, we estimated the internal structure of the factors/dimensions by modeling them as a network. The R library BGGM and Copula gaussian graphical model estimation was used to identify edges among the dimensions. To control for false positive rate, we used 95% credible intervals (CIs) for the edges. All edges whose 95% CI did not include zero were considered statistically significant. The nodes and all significant edges were plotted as a network using the Fruchterman-Reingold algorithm implemented in the R-package qgraph. To examine whether any node was more strongly associated with other nodes in the network, we estimated the predictability (an R2-like measure) of each node. High predictability indicates that a node has many and strong edges with other nodes in the network. We compared the predictability of all nodes (i.e., factors/dimensions) and differences for which the 99% CI did not include zero were considered statistically significant. A 99% CI was used because of the large sample size and multiple comparisons.
To examine how the CPTI dimensions were associated with CP, we added a node to the network that indicated the degree of CP that the child exhibited at T1 (i.e., cross-sectional associations). To evaluate whether some CPTI dimension were more strongly related to CP than others, we compared all edges between the CPTI dimensions and CP. Differences for which the 99% CI did not include zero were considered statistically significant. Parent- and teacher-rated CPTI and CP data from T1 (ages 3–5) were used.
To examine which factors/dimensions were most important to predict later CP, we used regression models. CPTI dimensions were added as independent variables and later CP as the dependent variable. Two measures of later CP were used: (1) continuous parent- and teacher-rated CP scores at T3 and (2) stable CP defined as 0.5 SD above the mean of the CP measure at T2 and T3. We made inference based on the degree of explained variance of the full model and which independent variables were significantly associated with later CP. To make further inference, dominance analysis was used in which the unique contribution (in the form of explained variance) of each independent variable to later CP was estimated. For continuous CP scores, we used linear regression and for stable CP, we used logistic regression. For the logistic regression and the subsequent dominance analysis, Cox and Snell’s R2 were used to interpret explained variance. All predictive models were first conducted using only CPTI dimensions as independent variables and then by adding T1 CP as a covariate.
For comparative reasons, main analyses were replicated for the original 3-factor structure, with results presented as Supplemental material. All additional data and study materials are available upon request to the corresponding author.