Skip to main content

ORIGINAL RESEARCH article

Front. Public Health, 12 March 2024
Sec. Public Mental Health
This article is part of the Research Topic Artificial Intelligence and Mental Health Care View all 13 articles

Factors and pathways of non-suicidal self-injury in children: insights from computational causal analysis

  • 1Department of Health Policy and Management, West China School of Public Health and West China Fourth Hospital, Sichuan University, Chengdu, China
  • 2College of Computer Science, Sichuan University, Chengdu, Sichuan, China
  • 3Teaching and Research Section of General Practice, The General Practice Medical Center, West China Hospital of Sichuan University, Chengdu, China

Background: Non-suicidal self-injury (NSSI) has become a significant public health issue, especially prevalent among adolescents. The complexity and multifactorial nature of NSSI necessitate a comprehensive understanding of its underlying causal factors. This research leverages the causal discovery methodology to explore these causal associations in children.

Methods: An observational dataset was scrutinized using the causal discovery method, particularly employing the constraint-based approach. By integrating machine learning and causal inference techniques, the study aimed to determine direct causal relationships associated with NSSI. The robustness of the causal relationships was evaluated using three methods to construct and validate it: the PC (Peter and Clark) method, Fast Causal Inference (FCI) method, and the GAE (Graphical Autoencoder) method.

Results: Analysis identified nine nodes with direct causal relationships to NSSI, including life satisfaction, depression, family dysfunction, sugary beverage consumption, PYD (positive youth development), internet addiction, COVID-19 related PTSD, academic anxiety, and sleep duration. Four principal causal pathways were identified, highlighting the roles of lockdown-induced lifestyle changes, screen time, positive adolescent development, and family dynamics in influencing NSSI risk.

Conclusions: An in-depth analysis of the factors leading to Non-Suicidal Self-Injury (NSSI), highlighting the intricate connections among individual, family, and pandemic-related influences. The results, derived from computational causal analysis, underscore the critical need for targeted interventions that tackle these diverse causative factors.

1 Introduction

Non-suicidal self-injury (NSSI) has become a pressing public health issue, with rising prevalence both in developed and developing countries (1). Rates of NSSI fluctuate between 11.5 and 33.8%, contingent on sample type and study design, but there is an undeniable uptrend worldwide, even in developing nations. Adolescence is the peak risk period for NSSI (2), with about 23% of adolescents, 13.4% of young adults, and 5.5% of adults being affected (3). Alarmingly, up to 24.7% of Chinese adolescents report experiencing NSSI (4), which necessitates further attention (5).

The repercussions of NSSI in children and adolescents are severe and long-lasting (6, 7). It is closely linked to heightened suicidal ideation and attempts (8). Even when accounting for suicidal thoughts, NSSI remains a potent predictor of suicidal actions (911). Specifically, 39.6% of those who have self-harmed report suicidal behaviors, and of that cohort, 66.3% have a history of NSSI (12). Besides, NSSI correlates with several psychological challenges, like depression, anxiety, and post-traumatic stress disorder(PTSD), and negatively impacts familial and interpersonal bonds (13, 14).

The risk factors for non-suicidal self-injury in children and adolescents are diverse and multifaceted. While much emphasis has been placed on the psychological, family, and school levels—including psychological disorders and symptoms, adverse family experiences, and victimization (4, 15, 16)—the underlying mechanisms driving NSSI behaviors are complex. Many studies to date have been constrained by methodological limitations, failing to grasp the problem from a comprehensive and systematic viewpoint. Moreover, the emergence of the COVID-19 pandemic and its associated lockdown measures has further complicated the scenario (17, 18). Factors such as changes in sleep time (19), physical activity, screen time (20), and increased stress due to isolation and academic pressures have heightened the mental and emotional issues among children and adolescents, possibly leading to a surge in NSSI behaviors (2123).

Protective factors against Non-Suicidal Self-Injury (NSSI), including life satisfaction and Positive Youth Development (PYD), have been highlighted as promising avenues for intervention. Research highlights that higher levels of life satisfaction can serve as a significant buffer against self-injurious behaviors, while the presence of PYD qualities has been shown to not only reduce the risk of NSSI but also lessen the impact of depressive symptoms on such behaviors (2426). This underscores the PYD perspective's shift from focusing on youth deficits or psychopathology to emphasizing their strengths, skills, and assets, which can be nurtured and improved (27, 28). These factors are theorized to buffer individuals from the deleterious effects of risk factors like depression and family dysfunction. Enhancing these protective mechanisms may mitigate the risk of NSSI, providing a critical strategy for supporting at-risk adolescents.

Based on the literature review above and the four-function model of self-injury (29, 30), we developed a comprehensive hypothetical framework (Figure 1) for the emergence of NSSI in adolescents. we have formulated a series of causal pathway hypotheses to be explored within our comprehensive hypothetical framework.

• Changes in sleep times, physical activity, consumption of sugar-sweetened beverages, and screen time during COVID-19 lockdowns are hypothesized to augment the probability of developing familial and psychological dysfunctions, subsequently enhancing the risk of NSSI.

• COVID-19-related PTSD is presumed to escalate the risk of psychological issues, further contributing to NSSI.

• Family dysfunction may indirectly lead to NSSI by exacerbating psychological distress.

• PYD qualities could potentially modulate the relationship between psychological distress and NSSI.

Figure 1
www.frontiersin.org

Figure 1. NSSI comprehensive hypothetical framework.

Building on the comprehensive review and the hypothetical framework outlined above, this research adopts causal discovery methods within a robust theoretical framework to elucidate the complex web of factors influencing NSSI among adolescents. Traditional research methodologies, which are predominantly observational and correlational, have contributed valuable insights yet frequently encounter limitations in establishing causality (31). Moreover, there is a tendency within existing research to adopt a narrow focus on isolated variables, neglecting the broader constellation of contributing factors. Recognizing the limitations of conventional research methodologies in capturing the multifaceted nature of NSSI, our study utilizes a causal discovery approach applied to a unique observational dataset, aiming to identify causal factors associated with NSSI in children.

Causal discovery, a methodological paradigm dedicated to unearthing cause-and-effect relationships among variables, emerges as a formidable tool in scientific inquiry and data analysis. It enables researchers to transcend the boundaries of correlation, probing the underlying “why” behind observed phenomena (3234). In the domain of NSSI, discerning the causal dynamics behind this complex behavior is crucial for crafting effective prevention and intervention measures. Through this approach, our study aims to contribute a nuanced understanding of NSSI, fostering the development of more precise and efficacious strategies to combat this pressing issue among youth.

2 Materials and methods

2.1 Study design and participants

Data for this research was sourced from the Chengdu Positive Child Development (CPCD) survey (35). Launched in December 2019 in Chengdu, this school-based longitudinal study targeted students aged 6–16 years, drawn from five primary and secondary schools. Information was collected using a mix of paper and electronic questionnaires completed by children and caregivers, as well as through direct interviews. The study aimed to explore connections between students' sociodemographic factors, health, lifestyle behaviors, and academic performance, with supplementary data obtained from caregiver health evaluations and school records. The first follow-up of the survey, conducted between June 16, 2020, and July 8, 2020, aligned with the resumption of classroom activities following the COVID-19 lockdowns.

A meticulous vetting process was utilized to refine the dataset, ultimately including 5,807 students to ensure high-quality data for robust causal inference. This process involved identifying and eliminating statistical outliers, removing variables with a missing rate exceeding 90%, and excluding samples lacking caregiver questionnaires or COVID-19 related responses. The “Sample Selection and Refinement Process for NSSI Study” is illustrated in Figure 2.

Figure 2
www.frontiersin.org

Figure 2. Sample selection and refinement process for NSSI study.

2.2 Ethics statement

The confidentiality of all data was meticulously maintained by the research team, with no disclosure of any personal information within the study's findings. This research adhered to the principles set forth in the Helsinki Declaration, secured approval from the Ethics Committee of Sichuan University (Approval No. K2020025), and obtained consent from the relevant school authorities, parents, and students.

2.3 Measures

2.3.1 Non-suicidal self-injury

The assessment of Non-suicidal self-injury (NSSI) in our study was conducted using the 9-item Intentional Self-Injury Scale, developed and validated by Gratz (36). This scale evaluates eight specific NSSI behaviors, including cutting, burning, biting, stabbing, hitting, pinching, and ingesting non-food items, along with a self-assessment for NSSI-related hospitalizations. Each behavior is rated on a 4-point scale, ranging from 1 (never) to 4 (three or more times). Participants reporting a score of 1 or higher on any item were identified as exhibiting NSSI behaviors. This measurement approach has been previously validated in studies involving Chinese adolescents (37).

2.3.2 Explanatory variables

In our research, we carefully selected 21 variables collected during the first follow-up period of the CPCD Survey. This selection was strategically guided by the existing literature (4, 810, 15, 38), and the aim of our study was to explore causal relationships with non-suicidal self-injury (NSSI). We focus on various sociodemographic, behavioral, and psychological factors that have been previously associated with NSSI behaviors. This approach aimed to create a manageable and comprehensive analytical framework that aligned our study with established findings and theory in NSSI research.

Sociodemographic data, including gender, age, grade, and BMI, were collected, prioritizing “Grade” for its relevance to educational context and peer interactions. Mental health status and psychological characteristics were evaluated through scales measuring depression, academic anxiety, positive child development, life satisfaction, family dysfunction, and internet addiction, aiming to create a comprehensive analytical framework that aligns with established NSSI research findings and theories.

Moreover, the study evaluated COVID-19-related behavioral factors like PTSD, changes in sleep patterns, physical activity, screen time, and sugary beverage consumption to understand the pandemic's impact on students' routines and its potential link to NSSI behaviors. These evaluations used precise binary coding for behavioral changes, offering detailed insights into lifestyle adjustments during the pandemic. Additional considerations included internet addiction and caregivers' mental health status, with a comprehensive overview of all variables, definitions, and measurement methods provided in the Supplementary material.

2.4 Data pre-processing

We conducted essential data preprocessing steps to ensure that we had a well-prepared dataset for subsequent causal discovery analysis. To harmonize the numerical and categorical features within the dataset, we employed standardized procedures. For the categorical variables, which included ordinal data with a discernible degree of ordering, we applied the LabelEncoder technique to transform them into numerical representations while preserving the ordinal relationships. Simultaneously, we standardized the numerical features using the StandardScaler method, which centered the variables around a mean of zero and scaled them to unit variance. This preprocessing not only alleviated potential issues arising from varying scales but also facilitated the compatibility of our data with the causal discovery algorithm, ensuring that it operated effectively in uncovering causal relationships within our dataset.

2.5 Causal discovery model application

It is crucial to recognize that causality and correlation are distinct concepts with profound implications for scientific research and data analysis. While correlation measures the strength and direction of statistical association between variables, it does not elucidate the direction of influence or establish a cause-and-effect relationship. Causality, on the other hand, delves into the intricate web of cause-and-effect connections, emphasizing that changes in one variable lead to changes in another. In this study, we aimed to identify the complex causality hidden in NSSI based on social-demographic characteristics, psychological characteristics, COVID-related behavioral factors, behavioral factors, mental health status and caregiver's mental health status.

The primary methods for learning the causal graph structure among variables include constraint-based, score-based, and gradient-based approaches. Constraint-based methods, such as PC (Peter and Clark) and FCI (39), identify causal relationships by recognizing patterns of conditional independence. Score-based methods, including GES (40), assign numerical scores to different causal structures and search for the structure with the highest score. Gradient-based methods, including GAE (41), have emerged as notable alternatives, exhibiting superior accuracy and computational efficiency. In comparison to constraint-based and score-based methods, gradient-based methods directly learn the causal graph through end-to-end optimization, eliminating the need for explicit detection of conditional independence or the use of scoring functions in intermediate steps.

Consequently, having recognized the strengths of each approach, we integrate a combination of these methods to enhance the accuracy and robustness of our causal Directed Acyclic Graph (DAG) learning. We employed three approaches for constructing and validating the DAG: the PC method (34), FCI (42), and the GAE. This is aimed at fortifying the reliability and credibility of our DAG, ensuring its resilience and accuracy across diverse analytical frameworks. The framework of our prediction work is presented in Figure 3.

Figure 3
www.frontiersin.org

Figure 3. The workflow for computational casual discovery.

PC algorithm represents a fundamental approach to causality analysis that is grounded in conditional independence testing. Its fundamental logic hinges upon the principle that genuine causal relationships manifest as conditional independence relationships in observational data. FCI algorithm is an extension of PC, designed to handle latent variables more effectively. As both PC and FCI are Constraint-based methods, they first use conditional independence tests to learn the skeleton of the underlying casual graph, and then orient the edges based on a series of orientation rules. It's important to note that the order in which variables are considered can influence the outcomes of these two algorithms.

To mitigate the impact of order-dependency, we conducted multiple runs of both PC and FCI with different random orders of variables. This ensures that the algorithm are exposed to various variable orders, mitigating the influence of a specific sequence on the final results. In our study, for FCI and the three variations of PC (original, stable, and parallel), we performed 1,000 runs for each, employing random variable orders. We identified the edges that appear most frequently across the multiple runs. The final structure, represented by edges that consistently appear, can be considered as the more robust DAG. The general steps for PC and FCI are described as follows:

1. Prior knowledge on variable selection. It is crucial to note that the PC and FCI algorithm assume causal sufficiency, and therefore, all potential confounding factors should be present in the data for a comprehensive causal analysis. We conscientiously leveraged prior knowledge in the face of a large number of total features to identify the variables relevant to NSSI.

2. Prior knowledge on variable relationship. We added prior knowledge to guide our causal modeling by informing the anticipated forbidden causal pathways. The rationale behind the prohibition of each path is detailed below.

• Prohibition of paths leading to “Grade” and “Gender” variables. Within our causal modeling framework, we have enforced a stringent constraint by prohibiting any causal pathways that lead to the “Grade” and “Gender” variables. This restriction is grounded in the fundamental premise that “Grade” and “Gender” are considered intrinsic attributes of individuals, impervious to external causal influences.

• Prohibition of paths leading to “COVID-Related Behavioral Factors”. We also prohibited any causal pathways leading to the “COVID-Related Behavioral Factors”. This constraint arose from the distinctive temporal characteristics and underlying nature of these factors. Unlike other variables in our analysis, which were sourced from a single survey conducted during the COVID-19 pandemic, the “COVID-Related Behavioral Factors” represented the dynamic changes observed between this COVID-19 period survey and another before-COVID-19 period survey. The timeline of these behavioral factors is inherently distinct, encompassing the evolving responses of individuals to the unique circumstances brought about by the pandemic. Thus, we enforced the prohibition of all paths leading to “COVID-Related Behavioral Factors”.

3. Causal discovery process. Our application of the PC and FCI algorithms followed these essential steps:

• Skeleton discovery: Begin with a fully connected undirected graph that includes all variables. Utilized conditional independence tests to uncover the initial skeleton of the causal graph, representing potential pairwise relationships between variables. We performed independence test using Fisher-z's test.

• Edge Orientation: PC orients edges based on conditional independence tests. FCI orients edges within Markov equivalence classes, considering latent variables and capturing more complex causal relationships.

• Refinement: PC refines the graph iteratively by applying conditional independence tests and edge removals. FCI refines the graph by considering additional conditional independence tests and making decisions about edge orientations within Markov equivalence classes.

• Multiple runs and final structure: For FCI and three variations of PC, we performed 1,000 runs with random variable orders. Identified edges that consistently appeared across multiple runs to establish a final structure.

GAE algorithm represents Graph Autoencoder for causal structure learning. It is an alternative generalization of NOTEARS (43) to handle nonlinear causal relations. After decoder, GAE can generate an adjacency matrix that captures the relationships between nodes. The adjacency matrix typically includes positive values for positive edges and negative values for negative edges. Positive value might represent a positive relationship or interaction, while negative value might represent a negative or inhibitory relationship. We used the learned latent representations from the GAE to identify factors associated with protective effects and those contributing to an increased risk, shown in Figure 3.

3 Results

3.1 Sample characteristics

Among the 5,307 participants included in the final analysis, 1,394 individuals (26.27%) reported engaging in non-suicidal self-injury (NSSI) in the past year. Among the 1,394 participants who reported engaging in NSSI, the NSSI detection rate was 24.82% for males (n = 667) and 27.75% for females (n = 727). The NSSI detection rate among primary school students was 23.51% (n = 886), while among middle school students, it was 33% (n = 508). The characteristics of the sample are presented in Table 1.

Table 1
www.frontiersin.org

Table 1. Sample characteristics.

3.2 Direct causes and effects of NSSI (local causal network)

We utilize a directed acyclic graph (DAG) to exhibit the complexity of variables influencing adolescent and child non-suicidal self-injury (NSSI) behaviors during the COVID-19 pandemic. Within this illustrative graph, the structure of the local causal network model, fundamental to our study, is evident. It includes an array of variables, each having its significance in the broader scope. The DAG aptly captures the “upstream” and “downstream” variables (31). Here, “upstream” variables act as forerunners or primary triggers that might affect succeeding variables. Conversely, the “downstream” variables are those likely influenced or modified due to the preceding ones.

The DAG offers a comprehensive visualization, revealing both direct neighbors and secondary neighbors within two “steps” of NSSI behaviors. This graphical representation paints an in-depth portrait of the relationships between various factors during the COVID-19 pandemic's acute phase for adolescents. Our application of the GAE algorithm not only corroborates these findings but also enriches them by providing a nuanced understanding of the factors contributing to NSSI. This advanced method, represented in Figure 3, facilitated the quantification of each variable's impact by calculating their polarity and magnitude. The GAE algorithm's unique ability to model nonlinear relationships allowed us to further delineate factors with protective effects (denoted in green) from those associated with increased risk (indicated in red). The results, as depicted in Figure 4, confirm and strengthen our causal network model, providing a nuanced understanding of the variables that directly or indirectly contribute to NSSI behaviors.

Figure 4
www.frontiersin.org

Figure 4. The DAG for local variable in NSSI.

The nodes with direct causal relationship with NSSI are life satisfaction, depression, family dysfunction, sugary beverage consumption, PYD, Internet addiction, COVID-19 related PTSD, academic anxiety and sleep duration. The analysis identified four key causal pathways influencing NSSI: (1) Lockdown-induced reductions in physical activity escalate academic anxiety, potentially leading to PTSD, heightening internet addiction and NSSI risks. (2) Enhanced screen time and sugary beverage consumption are linked to elevated depression risks and increased NSSI likelihood. (3) Positive adolescent development acts as a buffer, mitigating the adverse effects of family dysfunction and internet addiction on NSSI. (4) Family dysfunction negatively impacts life satisfaction and fosters depression, directly contributing to NSSI. The detailed DAG illustrating these relationships is provided in Supplementary Figure S1.

The graph further enlightens us about individual attributes that could influence this dynamic, underscoring the phenomena's multifaceted nature. The focus on NSSI behavior in the DAG, and its potential links with these various factors, is particularly noteworthy. Moreover, the graph showcases the intricate ties between these variables, elucidating the ripple effects and possible feedback mechanisms within the system. Ultimately, the DAG serves as a comprehensive visual tool, elucidating the intricate connections and potential causal pathways pertaining to NSSI behavior in adolescents amid the COVID-19 pandemic.

4 Discussion

Our research harnessed machine learning and causal discovery techniques to understand the intricate web of causal links surrounding NSSI in children and adolescents. Aiming to shed light on influential factors related to individual characteristics, family dynamics, and the repercussions of the COVID-19 pandemic, we strived to provide actionable insights that could spearhead innovative prevention and intervention approaches. Given the complex nature of NSSI, dictated by numerous variables, and the limitations of current research methodologies, we sought a fresh approach to discern causation, ultimately guiding more effective strategies. The PC and FCI algorithms exhibit sensitivity to the order of variables, implying that the obtained results may vary based on the sequence in which variables are considered during execution. To address this, we conducted multiple runs of these algorithms using different orders of features. Additionally, to enhance the robustness of the causal structure, we incorporated the GAE method into our approach. Subsequently, the most frequently identified edges were retained for further in-depth analysis.

Our study has delineated a network of nine pivotal variables that exhibit direct causal associations with NSSI behavior in children, as visualized in Figure 4 (for the full names and assignments corresponding to the abbreviations of the variables shown in Supplementary Table S1). Central among these findings is the pervasive influence of the COVID-19 pandemic, which manifest as a marked decrease in sugary beverage intake, diminished sleep quality, and the emergence of COVID-19-triggered PTSD. Individual psychological drivers encompassed themes of academic stress and manifestations of depression. In tandem, key psychological attributes were identified: the holistic life satisfaction measure, the nuanced interplays of family dynamics, the embodiment of Positive Youth Development (PYD) (26), and the grip of internet addiction. The coherence of these findings with existing literature underscores the robustness of our methodological approach (4, 44, 45).

Significantly, the COVID-19 pandemic has altered the sleep patterns of children and adolescents, a crucial aspect of mental health. Prolonged lockdowns have disrupted regular sleep schedules, heightening stress and anxiety levels. These sleep disruptions, combined with the stress of new learning modes and social isolation, have intensified academic anxieties, affecting family dynamics and increasing household tensions. Moreover, the pandemic has triggered PTSD in young people, a concerning development given its long-term mental health implications. This emergence of PTSD, fueled by constant pandemic-related news and personal experiences, adds a critical dimension to our understanding of the pandemic's impact on youth mental health.

Digging deeper into these variables, the interconnectedness within this network is profound. Academic stress finds its roots in the diminished sleep and curtailed physical activity patterns during the pandemic, as well as the challenges posed by family dynamics. The haunting presence of COVID-19-related PTSD is sculpted by the restrictive physical activity regimes, academic stressors, and family dysfunction. Depression's emergence is further amplified by excessive screen time, both for academic and recreational purposes, interwoven with academic anxieties, the nurturing aspects of Positive Youth Development, and overarching life satisfaction. Completing this intricate web, internet addiction bears the imprints of PTSD linked to the pandemic, the shadows of depression, family struggles, and the protective or exacerbating elements of Positive Youth Development.

Emerging from our findings is a comprehensive understanding that provides invaluable insights for safeguarding the mental wellbeing of children and adolescents in the wake of unforeseen public health crises, such as the COVID-19 pandemic. This understanding emphasizes the potential of our data in guiding practical interventions and preventive measures. The notable reduction in sugary beverage consumption and sleep duration, as well as the emergence of pandemic-induced PTSD, highlight the profound physiological and psychological shifts induced by prolonged lockdowns and associated societal changes (46, 47). Our findings advocate for prevention strategies that are not only trauma-informed but also adaptive to the evolving public health landscape.

The established causal relationship between sugary beverage consumption and NSSI echoes previous studies associating unhealthy diets with increased depressive symptoms in adolescents (48, 49). Inflammatory diets can intensify mental health problems, possibly through obesity and inflammation (50, 51). Thus, the importance of balanced diets is underscored. Emphasizing a diet rich in fruits, vegetables, and anti-inflammatory foods, combined with strategies like restoring regular sleep patterns and trauma-informed interventions, may offer a comprehensive approach to improving the mental wellbeing of children and adolescents.

Furthermore, understanding the intricate network of causal factors—including academic anxieties (52), family dynamics (53), and individual psychological attributes—means that interventions can be more targeted and precise. This precision is vital in the clinical setting, where tailored interventions can lead to more effective outcomes. Rather than applying generic measures, strategies can be devised to specifically counteract or augment identified causal agents. This approach enhances the practical utility of our findings, offering a roadmap for clinicians and policymakers in developing targeted interventions. The identified drivers, such as academic stress and depression, could be targeted through school-based programs emphasizing coping mechanisms, emotional regulation, and peer support. Family-centered interventions might focus on strengthening familial bonds and improving communication, reducing the chance of family dysfunction exacerbating mental health issues (54). The correlation between internet addiction and NSSI emphasizes the importance of digital literacy programs that equip adolescents with skills to navigate the online world safely.

A key distinction of our study lies in the use of computational causal discovery. Traditional methodologies often restrict themselves to observational correlations, which, although informative, don't offer a genuine window into the underlying causative structures. Causal discovery goes beyond merely identifying these associations, allowing us to pinpoint the drivers of adverse outcomes such as NSSI. The superiority of this approach lies in its potential to tailor prevention and intervention strategies based on the actual causes, rather than mere symptoms or correlated factors. This means that initiatives informed by our findings can be significantly more effective, as they strike at the heart of the issue, directly addressing and mitigating the root causes. As we move forward in our collective endeavor to nurture the mental health of our younger generation, leveraging advanced methodologies like computational causal discovery will be paramount in ensuring our strategies are not only well-informed but also impactful.

5 Strengths and limitations

A significant strength of our study is the use of multiple causal discovery algorithms, enhancing the robustness and interpretability of the results. This approach marks a departure from traditional approaches that predominantly rely on correlations. This novel approach facilitates a nuanced understanding of the intricate “cause-and-effect” dynamics underlying NSSI behaviors in children and adolescents, particularly in the unique context of the COVID-19 pandemic. By pinpointing fundamental causative elements, our study lays the foundation for more targeted and efficacious interventions that address root causes, providing both immediate and sustained psychological health benefits.

Nevertheless, the study has limitations, notably the influence of variable ordering on model outcomes. We've addressed this by running numerous iterations, lending consistency to our findings despite potential variable ordering effects.

While the adoption of causal algorithms marks an advancement in our analysis, they cannot fully negate the impact of unobserved confounders or bidirectional relationships. The cross-sectional data limits temporal causality claims, necessitating further validation with longitudinal studies.

Additionally, while AI methods enhance efficiency and aid in the application of causal inference to health data, they are not infallible. The necessity for human judgment remains a key component in the interpretative process, not least because different AI algorithms might yield varying interpretations or conclusions (55). This involves using different AI algorithms for mutual validation and conducting repeated experiments to test their robustness. This dual approach of leveraging both AI and human judgment facilitates a more nuanced and robust analysis than could be achieved by either one alone.

In terms of academic recommendations and future directions, there is a compelling need for further exploration of the identified mechanisms underpinning NSSI behaviors. Using longitudinal datasets with time series information, in conjunction with computational causal discovery, can offer more robust and definitive insights into causality. forging interdisciplinary collaborations that meld psychological, societal, and technological insights could provide a more holistic understanding and usher in innovative intervention strategies tailored to the multifaceted challenges of the modern era.

6 Conclusions

Drawing on unique computational causal discovery and machine learning methods, this study illuminated the intricate causal network of factors influencing NSSI in children during the COVID-19 pandemic. Our findings underscore nine critical variables intricately interwoven, reflecting the profound effects of the pandemic, academic stress, family dynamics, and individual psychological attributes. The study's insights offer a fresh perspective for devising impactful interventions, emphasizing the significance of addressing root causes, particularly in the wake of unprecedented global challenges.

Data availability statement

The data analyzed in this study is subject to the following licenses/restrictions: the dataset can only be used for non-commercial, academic research purposes. Requests to access these datasets should be directed to Rui Hu, aHUyODU3OTExODk2JiN4MDAwNDA7MTYzLmNvbQ==.

Ethics statement

The studies involving humans were approved by Ethics Committee of Sichuan University (Approval No. K2020025). The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants' legal guardians/next of kin.

Author contributions

XG: Conceptualization, Data curation, Investigation, Writing – original draft. LW: Formal analysis, Methodology, Writing – original draft. ZF: Writing – review & editing. ZL: Data curation, Formal analysis, Visualization, Writing – original draft. LL: Supervision, Writing – review & editing. LJ: Supervision, Writing – review & editing. LZ: Funding acquisition, Supervision, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Natural Science Foundation of China (NSFC) (No. 82273748).

Acknowledgments

We would like to express our gratitude to ChatGPT, a language model based on the GPT-4 architecture developed by OpenAI, for assisting with the translation of portions of this manuscript.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2024.1305746/full#supplementary-material

References

1. Mannekote Thippaiah S, Shankarapura Nanjappa M, Gude JG, Voyiaziakis E, Patwa S, Birur B, et al. Non-suicidal self-injury in developing countries: a review. Int J Soc Psychiatry. (2021) 67:472–82. doi: 10.1177/0020764020943627

Crossref Full Text | Google Scholar

2. Lim KS, Wong CH, McIntyre RS, Wang J, Zhang Z, Tran BX, et al. Global lifetime and 12-month prevalence of suicidal behavior, deliberate self-harm and non-suicidal self-injury in children and adolescents between 1989 and 2018: a meta-analysis. Int J Environ Res Public Health. (2019) 16:4581. doi: 10.3390/ijerph16224581

PubMed Abstract | Crossref Full Text | Google Scholar

3. Rodriguez-Blanco L, Carballo JJ, Baca-Garcia E. Use of ecological momentary assessment (EMA) in non-suicidal self-injury (NSSI): a systematic review. Psychiatry Res. (2018) 263:212–9. doi: 10.1016/j.psychres.2018.02.051

Crossref Full Text | Google Scholar

4. Qu D, Wen X, Liu B, Zhang X, He Y, Chen D, et al. Non-suicidal self-injury in Chinese population: a scoping review of prevalence, method, risk factors and preventive interventions. Lancet Reg Health West Pac. (2023) 37:100794. doi: 10.1016/j.lanwpc.2023.100794

PubMed Abstract | Crossref Full Text | Google Scholar

5. Wang C, Zhang P, Zhang N. Adolescent mental health in China requires more attention. Lancet Public Health. (2020) 5:e637. doi: 10.1016/S2468-2667(20)30094-3

PubMed Abstract | Crossref Full Text | Google Scholar

6. Borschmann R, Becker D, Coffey C, Spry E, Moreno-Betancur M, Moran P, et al. 20-year outcomes in adolescents who self-harm: a population-based cohort study. Lancet Child Adolesc Health. (2017) 1:195–202. doi: 10.1016/S2352-4642(17)30007-X

PubMed Abstract | Crossref Full Text | Google Scholar

7. Mars B, Heron J, Crane C, Hawton K, Lewis G, Macleod J, et al. Clinical and social outcomes of adolescent self harm: population based birth cohort study. BMJ. (2014) 349:g5954. doi: 10.1136/bmj.g5954

PubMed Abstract | Crossref Full Text | Google Scholar

8. Poudel A, Lamichhane A, Magar KR, Khanal GP. Non suicidal self injury and suicidal behavior among adolescents: co-occurrence and associated risk factors. BMC Psychiatry. (2022) 22:96. doi: 10.1186/s12888-022-03763-z

Crossref Full Text | Google Scholar

9. Mars B, Heron J, Klonsky ED, Moran P, O'Connor RC, Tilling K, et al. Predictors of future suicide attempt among adolescents with suicidal thoughts or non-suicidal self-harm: a population-based birth cohort study. Lancet Psychiatry. (2019) 6:327–37. doi: 10.1016/S2215-0366(19)30030-6

PubMed Abstract | Crossref Full Text | Google Scholar

10. Hawton K, Bale L, Brand F, Townsend E, Ness J, Waters K, et al. Mortality in children and adolescents following presentation to hospital after non-fatal self-harm in the multicentre study of self-harm: a prospective observational cohort study. Lancet Child Adolesc Health. (2020) 4:111–20. doi: 10.1016/S2352-4642(19)30373-6

PubMed Abstract | Crossref Full Text | Google Scholar

11. Groschwitz RC, Kaess M, Fischer G, Ameis N, Schulze UM, Brunner R, et al. The association of non-suicidal self-injury and suicidal behavior according to DSM-5 in adolescent psychiatric inpatients. Psychiatry Res. (2015) 228:454–61. doi: 10.1016/j.psychres.2015.06.019

PubMed Abstract | Crossref Full Text | Google Scholar

12. Voss C, Hoyer J, Venz J, Pieper L, Beesdo-Baum K. Non-suicidal self-injury and its co-occurrence with suicidal behavior: an epidemiological-study among adolescents and young adults. Acta Psychiatr Scand. (2020) 142:496–508. doi: 10.1111/acps.13237

PubMed Abstract | Crossref Full Text | Google Scholar

13. Bentley KH, Cassiello-Robbins CF, Vittorio L, Sauer-Zavala S, Barlow DH. The association between nonsuicidal self-injury and the emotional disorders: a meta-analytic review. Clin Psychol Rev. (2015) 37:72–88. doi: 10.1016/j.cpr.2015.02.006

Crossref Full Text | Google Scholar

14. Taylor PJ, Jomar K, Dhingra K, Forrester R, Shahmalak U, Dickson JM, et al. meta-analysis of the prevalence of different functions of non-suicidal self-injury. J Affect Disord. (2018) 227:759–69. doi: 10.1016/j.jad.2017.11.073

Crossref Full Text | Google Scholar

15. Liang K, Zhao L, Lei Y, Zou K, Ji S, Wang R, et al. Nonsuicidal self-injury behaviour in a city of China and its association with family environment, media use and psychopathology. Compr Psychiatry. (2022) 115:152311. doi: 10.1016/j.comppsych.2022.152311

PubMed Abstract | Crossref Full Text | Google Scholar

16. Wang Y, Luo B, Hong B, Yang M, Zhao L, Jia P. The relationship between family functioning and non-suicidal self-injury in adolescents: a structural equation modeling analysis. J Affect Disord. (2022) 309:193–200. doi: 10.1016/j.jad.2022.04.124

PubMed Abstract | Crossref Full Text | Google Scholar

17. Murata S, Rezeppa T, Thoma B, Marengo L, Krancevich K, Chiyka E, et al. The psychiatric sequelae of the COVID-19 pandemic in adolescents, adults, and health care workers. Depress Anxiety. (2021) 38:233–46. doi: 10.1002/da.23120

PubMed Abstract | Crossref Full Text | Google Scholar

18. De Luca L, Giletta M, Nocentini A, Menesini E. Non-suicidal self-injury in adolescence: the role of pre-existing vulnerabilities and COVID-19-related stress. J Youth Adolesc. (2022) 51:2383–95. doi: 10.1007/s10964-022-01669-3

PubMed Abstract | Crossref Full Text | Google Scholar

19. Wright KP, Linton SK, Withrow D, Casiraghi L, Lanza SM, de la Iglesia H, et al. Sleep in university students prior to and during COVID-19 stay-at-home orders. Curr Biol. (2020) 30:R797–8. doi: 10.1016/j.cub.2020.06.022

PubMed Abstract | Crossref Full Text | Google Scholar

20. Brito LMS, da Silva Boguszewski MC, de Souza MTR, Martins F, Mota J, Leite N. Indoor physical activities, eating and sleeping habits among school adolescents during COVID-19 pandemic. Rev Bras Atividade Física Saúde. (2020) 25:1–6. doi: 10.12820/rbafs.25e0117

Crossref Full Text | Google Scholar

21. Schwartz-Mette RA, Duell N, Lawrence HR, Balkind EG. COVID-19 distress impacts adolescents' depressive symptoms, NSSI, and suicide risk in the rural, northeast US. J Clin Child Adolesc Psychol. (2022) 52:2042697. doi: 10.1080/15374416.2022.2042697

PubMed Abstract | Crossref Full Text | Google Scholar

22. John A, Eyles E, Webb RT, Okolie C, Schmidt L, Arensman E, et al. The impact of the COVID-19 pandemic on self-harm and suicidal behaviour: update of living systematic review. F1000Research. (2020) 9:1097. doi: 10.12688/f1000research.25522.2

Crossref Full Text | Google Scholar

23. Xiao J, Wang R, Hu Y, He T, Ruan Z, Chen Q, et al. Impacts of the psychological stress response on nonsuicidal self-injury behavior in students during the COVID-19 epidemic in China: the mediating role of sleep disorders. BMC Psychol. (2022) 10:87. doi: 10.1186/s40359-022-00789-6

PubMed Abstract | Crossref Full Text | Google Scholar

24. Kress VE, Newgent RA, Whitlock J, Mease L. Spirituality/religiosity, life satisfaction, and life meaning as protective factors for nonsuicidal self-injury in college students. J Coll Counsel. (2015) 18:160–74. doi: 10.1002/jocc.12012

Crossref Full Text | Google Scholar

25. Chen Y, Luo J, Jiang L, Shi W, Jia P, Zhang J, et al. Association between positive youth development and non-suicidal self-injury: a longitudinal survey of children and adolescents in southwest China, 2019-21. J Affect Disord. (2024) 350:755–60. doi: 10.1016/j.jad.2024.01.072

PubMed Abstract | Crossref Full Text | Google Scholar

26. Zhu X, Shek D. The predictive effect of depression on self-injury: positive youth development as a moderator. Appl Res Qual Life. (2023) 18:2877–94. doi: 10.1007/s11482-023-10211-x

Crossref Full Text | Google Scholar

27. Shek DT, Dou D, Zhu X, Chai W. Positive youth development: current perspectives. Adolesc Health Med Ther. (2019) 10:131–41. doi: 10.2147/AHMT.S179946

PubMed Abstract | Crossref Full Text | Google Scholar

28. Tolan P, Ross K, Arkin N, Godine N, Clark E. Toward an integrated approach to positive development: Implications for intervention. Appl Dev Sci. (2016) 20:214–36. doi: 10.1080/10888691.2016.1146080

Crossref Full Text | Google Scholar

29. Nock MK. Why do people hurt themselves? New insights into the nature and functions of self-injury. Curr Direct Psychol Sci. (2009) 18:78–83. doi: 10.1111/j.1467-8721.2009.01613.x

PubMed Abstract | Crossref Full Text | Google Scholar

30. Hepp J, Carpenter RW, Störkel LM, Schmitz SE, Schmahl C, Niedtfeld I. A systematic review of daily life studies on non-suicidal self-injury based on the four-function model. Clin Psychol Rev. (2020) 82:101888. doi: 10.1016/j.cpr.2020.101888

PubMed Abstract | Crossref Full Text | Google Scholar

31. Saxe GN, Ma S, Morales LJ, Galatzer-Levy IR, Aliferis C, Marmar CR. Computational causal discovery for post-traumatic stress in police officers. Transl Psychiatry. (2020) 10:233. doi: 10.1038/s41398-020-00910-6

PubMed Abstract | Crossref Full Text | Google Scholar

32. Zanga A, Ozkirimli E, Stella F. A survey on causal discovery: theory and practice. Int J Approx Reason. (2022) 151:101–29. doi: 10.1016/j.ijar.2022.09.004

Crossref Full Text | Google Scholar

33. Glymour C, Zhang K, Spirtes P. Review of causal discovery methods based on graphical models. Front Genet. (2019) 10:524. doi: 10.3389/fgene.2019.00524

PubMed Abstract | Crossref Full Text | Google Scholar

34. Zhang K, Zhu S, Kalander M, Ng I, Ye J, Chen Z, et al. gcastle: a python toolbox for causal discovery. arXiv. (2021). doi: 10.48550/arXiv.2111.15155

Crossref Full Text | Google Scholar

35. Zhao L, Shek DT, Zou K, Lei Y, Jia P. Cohort profile: Chengdu positive child development (CPCD) survey. Int J Epidemiol. (2022) 51:e95–e107. doi: 10.1093/ije/dyab237

PubMed Abstract | Crossref Full Text | Google Scholar

36. Gratz KL. Measurement of deliberate self-harm: preliminary data on the deliberate self-harm inventory. J Psychopathol Behav Assess. (2001) 23:253–63. doi: 10.1037/t04163-000

Crossref Full Text | Google Scholar

37. Wu J, Liu S, Luo J, Li X, You J. The effects of childhood abuse, depression, and self-compassion on adolescent nonsuicidal self-injury: a moderated mediation model. Child Abuse Neglect. (2023) 136:105993. doi: 10.1016/j.chiabu.2022.105993

PubMed Abstract | Crossref Full Text | Google Scholar

38. Hu R, Peng Ll, Du Y, Feng Yw, Xie Ls, Shi W, et al. Reciprocal effect between depressive symptoms and adolescent non-suicidal self-injury before and after COVID-19: a longitudinal study. Preprint. (2023). doi: 10.21203/rs.3.rs-2529545/v1

Crossref Full Text | Google Scholar

39. Spirtes P, Glymour CN, Scheines R. Causation, Prediction, and Search. Cambridge, MA: MIT Press (2000).

Google Scholar

40. Chickering DM. Optimal structure identification with greedy search. J Mach Learn Res. (2002) 3:507–54. doi: 10.1162/153244303321897717

Crossref Full Text | Google Scholar

41. Ng I, Zhu S, Chen Z, Fang Z. A graph autoencoder approach to causal structure learning. arXiv [preprint]. (2019). doi: 10.48550/arXiv.1911.07420

Crossref Full Text | Google Scholar

42. Spirtes PL, Meek C, Richardson TS. Causal inference in the presence of latent variables and selection bias. arXiv [preprint]. (2013). doi: 10.48550/arXiv.1302.4983

Crossref Full Text | Google Scholar

43. Zheng X, Aragam B, Ravikumar PK, Xing EP. Dags with no tears: Continuous optimization for structure learning. Adv Neur Inf Process Syst. (2018) 31:9492–503. doi: 10.48550/arXiv.1803.01422

Crossref Full Text | Google Scholar

44. Xiong A, Liao S, Luo B, Luo S, Tong Y, Li Z. Associations between problematic internet use, life satisfaction, and deliberate self-harm among Chinese adolescents: a multi-centered longitudinal study. Addict Behav. (2023) 147:107808. doi: 10.1016/j.addbeh.2023.107808

PubMed Abstract | Crossref Full Text | Google Scholar

45. Zhang Y, Jin Z, Li S, Xu H, Wan Y, Tao F. Relationship between chronotype and mental behavioural health among adolescents: a cross-sectional study based on the social ecological system. BMC Psychiatry. (2023) 23:1–12. doi: 10.1186/s12888-023-04879-6

PubMed Abstract | Crossref Full Text | Google Scholar

46. Du N, Ouyang Y, Xiao Y, Li Y. Psychosocial factors associated with increased adolescent non-suicidal self-injury during the COVID-19 pandemic. Front Psychiatry. (2021) 12:743526. doi: 10.3389/fpsyt.2021.743526

PubMed Abstract | Crossref Full Text | Google Scholar

47. Fegert JM, Vitiello B, Plener PL, Clemens V. Challenges and burden of the Coronavirus 2019 (COVID-19) pandemic for child and adolescent mental health: a narrative review to highlight clinical and research needs in the acute phase and the long return to normality. Child Adolesc Psychiatry Ment Health. (2020) 14:1–11. doi: 10.1186/s13034-020-00329-3

PubMed Abstract | Crossref Full Text | Google Scholar

48. Orlando L, Savel KA, Madigan S, Colasanto M, Korczak DJ. Dietary patterns and internalizing symptoms in children and adolescents: a meta-analysis. Aust New Zeal J Psychiatry. (2022) 56:617–41. doi: 10.1177/00048674211031486

PubMed Abstract | Crossref Full Text | Google Scholar

49. Chopra C, Mandalika S, Kinger N. Does diet play a role in the prevention and management of depression among adolescents? A narrative review. Nutr Health. (2021) 27:243–63. doi: 10.1177/0260106020980532

Crossref Full Text | Google Scholar

50. Sureda A, Bibiloni MdM, Julibert A, Bouzas C, Argelich E, Llompart I, et al. Adherence to the mediterranean diet and inflammatory markers. Nutrients. (2018) 10:62. doi: 10.3390/nu10010062

Crossref Full Text | Google Scholar

51. Oddy WH, Allen KL, Trapp GS, Ambrosini GL, Black LJ, Huang RC, et al. Dietary patterns, body mass index and inflammation: pathways to depression and mental health problems in adolescents. Brain Behav Immun. (2018) 69:428–39. doi: 10.1016/j.bbi.2018.01.002

PubMed Abstract | Crossref Full Text | Google Scholar

52. Chen H, Guo H, Chen H, Cao X, Liu J, Chen X, et al. Influence of academic stress and school bullying on self-harm behaviors among Chinese middle school students: the mediation effect of depression and anxiety. Front Public Health. (2023) 10:1049051. doi: 10.3389/fpubh.2022.1049051

PubMed Abstract | Crossref Full Text | Google Scholar

53. Cassels M, van Harmelen AL, Neufeld S, Goodyer I, Jones PB, Wilkinson P. Poor family functioning mediates the link between childhood adversity and adolescent nonsuicidal self-injury. J Child Psychol Psychiatry. (2018) 59:881–7. doi: 10.1111/jcpp.12866

PubMed Abstract | Crossref Full Text | Google Scholar

54. Kelada L, Hasking P, Melvin G. Adolescent NSSI and recovery: the role of family functioning and emotion regulation. Youth Soc. (2018) 50:1056–77. doi: 10.1177/0044118X16653153

Crossref Full Text | Google Scholar

55. Hunter DJ, Holmes C. Where medical statistics meets artificial intelligence. N Engl J Med. (2023) 389:1211–9. doi: 10.1056/NEJMra2212850

Crossref Full Text | Google Scholar

Keywords: NSSI, causal discovery, mental health, artificial intelligence, risk factors, COVID-19

Citation: Guo X, Wang L, Li Z, Feng Z, Lu L, Jiang L and Zhao L (2024) Factors and pathways of non-suicidal self-injury in children: insights from computational causal analysis. Front. Public Health 12:1305746. doi: 10.3389/fpubh.2024.1305746

Received: 02 October 2023; Accepted: 23 February 2024;
Published: 12 March 2024.

Edited by:

Jorge Piano Simoes, University of Twente, Netherlands

Reviewed by:

Paola Longo, University of Turin, Italy
Dan Qiu, Central South University, China

Copyright © 2024 Guo, Wang, Li, Feng, Lu, Jiang and Zhao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Lihua Jiang, bGhqaWFuZyYjeDAwMDQwO3NjdS5lZHUuY24=; Li Lu, bHVsaSYjeDAwMDQwO3NjdS5lZHUuY24=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.