{'fraction_non_alphanumeric': 0.026078651218803206, 'fraction_numerical': 0.014202068020431046, 'mean_word_length': 5.058364779874214, 'pattern_counts': {'":': 0, '<': 5, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 2, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Risk factors for recurrent injuries in victims of suspected non-accidental trauma: a retrospective cohort study

Background Many children who are victims of non-accidental trauma (NAT) may be repeatedly evaluated for injuries related to maltreatment. The purpose of this study was to identify risk factors for repeated injuries in children with suspected NAT. Methods We conducted a retrospective cohort study using claims data from a pediatric Medicaid accountable care organization. Children with birth claims and at least one non-birth related claim indicating a diagnosis of NAT or skeletal survey in 2007–2011 were included. Recurrent events were defined as independent episodes of care involving an urgent/emergent care setting that included a diagnosis code specific for child abuse, a CPT code for a skeletal survey, or a diagnosis code for an injury suspicious for abuse. Cox proportional hazards models were used to examine risk factors for recurrent events. Results Of the 1,361 children with suspected NAT, a recurrent NAT event occurred in 26% within 1 year and 40% within 2 years of their initial event. Independent risk factors for a recurrent NAT event included a rural residence, age < 30 months old, having only 1 or 2 initially detected injuries, and having a dislocation, open wound, or superficial injury at the previous event (p ≤ 0.01 for all). Conclusions Over 25% of children who experienced a suspected NAT event had a recurrent episode within one year. These children were younger and more likely to present with “minor” injuries at their previous event.


Background
Non-accidental trauma (NAT) is a leading cause of injury and death throughout early childhood [1,2]. Repeated evaluations in the medical setting for traumatic injuries should raise concerns that these injuries may be caused by either negligent behavior on the part of the caretaker or by recurrent intentional mechanisms.
Previously reported predictors of recurrent NAT include prior child protective services involvement, history of domestic violence, chronicity of maltreatment, child's age, parental history of maltreatment as a child, and parental substance abuse, criminal record, and mental health issues, or after specific injuries [5,6,[9][10][11][12][13]. These previous studies are limited in that they either do not assess risk factors related specifically to trauma, such as sentinel traumatic events, or they do not address recurrence of maltreatment. The purpose of this study was to identify patterns of injuries and factors associated with suspected episodes of recurrent NAT in a cohort of young children enrolled in a Medicaid managed care program who had at least one highly suspicious encounter for NAT.

Data source
Partners for Kids (PFK) is Nationwide Children's Hospital's pediatric accountable care organization. PFK contracts with the Medicaid Managed Care Organizations in Central and Southeastern Ohio to manage the care of almost 300,000 children across 37 counties, from urban Columbus to rural Appalachia. At the time of this study, over 2,000 physicians were submitting claims to PFK. The PFK claims database includes information on all billable medical care, procedures, and encounters for its enrollees, allowing for tracking of patients over time, across institutions, and across both inpatient and outpatient encounters. Access to this claims database is available to researchers at our institution, though is not freely available to individuals outside of our institution, and was granted by the PFK accountable care organization.

Study population
This study used enrollment data and facility and professional claims data from January 2007 to December 2011 for children born during this time period. We identified all children with a birth record claim who also had at least one claim indicating a diagnosis of abuse (physical, emotional, or neglect) or a skeletal survey at a non-birth related episode of care. (Figure 1) Suspected NAT events were defined as episodes of care in which a claim contained either (a) an International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM) discharge diagnosis code specific for child abuse, (b) a Current Procedural Terminology (CPT) coded skeletal survey, or (c) ICD-9 coded injuries suspicious for abuse; these events could be the event that brought the child into the study cohort, or they could occur either before or after that event. We excluded events that had an ICD-9 E-code for a trauma mechanism that could explain the injury or an ICD-9 code for an underlying medical illness that could explain the injury or need for skeletal survey. Episodes of care coded as follow-up care were excluded. Episodes of care with only a diagnosis of minor cutaneous injury from a specific mechanism and no other codes indicative of suspected NAT were also excluded. In order to include all claims for care related to a single incident of suspected NAT, an episode of care encompassed all claims for service provided concurrently or within two days of the care documented in the claim. In order to minimize the risk of defining claims for follow-up care as new events, only episodes of care that included encounters in the emergency department, urgent care, or inpatient setting were considered for inclusion as recurrent events. Figure 1 outlines cohort development and includes all ICD-9 and CPT codes used to define the cohort and events.

Independent variables
Variables determined at the time of each event included age, sex, days since last event, the presence of symptoms or diseases of the respiratory system, digestive system, nervous system and sense organs, skin and subcutaneous tissue, endocrine, nutritional, metabolic, or immunity disorders, vaccination during the episode of care, location, type, and mechanism of injuries, number of injuries, injury severity (evaluated as the probability of death based on the trauma mortality prediction model, TMPM-ICD9) [14,15], and death during the episode of care. The type and mechanism of injury were defined using ICD-9 diagnosis codes and E-codes respectively. The location of each injury was categorized into one of six body regions based on the Abbreviated Injury Scale [16]. The number of injuries was defined as the number of unique injury diagnosis codes listed during the episode of care of the event. Because family socioeconomic status (SES) indicators were not available, zip code level SES variables (median household family income, percentage of the population over age 25 with a Bachelor's degree or higher) and urban vs. rural residence were determined from the 2000 U.S. Census based on each child's zip code at their first event [17]. Enrollment duration in months and enrollment continuity were determined for each child. Other independent variables were determined according to their presence on any claim submitted prior to the suspected NAT event including musculoskeletal disease, congenital anomalies, and prematurity.

Statistical analysis
Characteristics at each suspected NAT event and in children with and without recurrent events were summarized using descriptive statistics (medians and inter-quartile ranges (IQR) or frequencies and percentages). Kaplan-Meier curves were used to display the proportion of children with recurrent events over time after the initial event. To determine risk factors for recurrent events, we used an extension of the Cox proportional hazards model for recurrent event data, the Prentice, Williams and Peterson gap time (PWP-GT) model [18]. Events beyond the fifth event were not included due to their insufficient number for analysis. Predictor variables in these models were the independent variables as measured at the previous suspected NAT event, with the exception of the zip code based variables, which were determined at the first event only. The reported hazard ratios (HR) estimate the relative hazard rates of having an event in those with and without the characteristic being examined.
For the examination of multivariable associations between the predictors and the time to the next suspected NAT event, Cox proportional hazards PWP-GT recurrent event models were used. All variables with bivariate associations significant at p < 0.20 were included with subsequent variable elimination until all remaining variables had p < 0.10. The final multivariable model revealed the overall associations of factors measured at any particular event with the risk of a subsequent suspected NAT event, after adjustment for other measured risk factors. We included all children regardless of their duration of follow-up in our analyses in order to minimize selection bias; in all of the survival analyses performed, children were included in the pool of patients at risk for subsequent events from the time of their initial event until the end of their last month of enrollment in PFK during the study period. Subsequently, several sensitivity analyses were performed to evaluate how the inclusion of patients with short follow-up, discontinuous follow-up, or without birth records in the database affected the results. The sensitivity analyses involved repeating the multivariable modeling excluding those children who did not maintain continuous enrollment in PFK, then excluding those children with less than 60 days of follow-up after their first event, and finally performing the analyses in those children who did not have birth records in the PFK database. All statistical analyses were performed using SAS (Statistical Analysis Software v9.3, Cary, NC). The conduct of this study was approved by Nationwide Children's Hospital Institutional Research Board with a waiver of informed consent. This research study has adhered to the STROBE guidelines for observational studies as outlined at http:// www.strobe-statement.org. Additional file 1.

Identification of cohort
Of the 140,828 children born and enrolled in PFK from 2007-2011, 2,362 had a claim with either a diagnosis of child maltreatment or a skeletal survey. Sixty-one percent of these children (n = 1,434) had birth records in the PFK database. After removing events with diagnosis codes for a medical illness or trauma mechanism that could potentially explain the injuries, the cohort was further refined to 1,361 children who were included in the main analysis ( Figure 1).

Frequency of recidivism
Three hundred and seventy-three (27.4%) patients in our cohort had more than one episode of care for a suspected NAT event during the study period ( Table 1). The incidence of suspected NAT events in the total cohort was 49 events per 1000 person-months. Two hundred and sixty-one children had 2 events, 74 children had 3 events, 22 children had 4 events, 13 children had 5 events, and 3 children had 6 events during the study period. Of all of these events, 51% had documentation of a skeletal survey, 35% had an abuse diagnosis, and 65% had an injury. Thirty percent of all events had only injuries, with no evidence of a skeletal survey or abuse diagnosis; at these events, the most common injuries were open wounds (32%) and contusions (27%). These potentially accidental injuries equated to an injury rate of 177/1000 person-years. Based on Kaplan-Meier analysis, 26% of the children had ≥1 recurrent event within 1 year of their initial event and 40% had ≥1 recurrence within 2 years of their initial event. The time between events decreased significantly with each subsequent event ( Figure 2, p < 0.0001). It is important to note that the duration of follow-up after the initial event widely varied (median (IQR) 383 days (145, 773)). However, the finding of significantly decreased time between events with increasing event number held in the subsample of 476 children with at least 600 days of follow-up after their first event (p = 0.005).

Demographics, comorbidities and injury characteristics of patient population
Characteristics and injuries of the children with a single event were compared to children with recurrent events (Table 1). Among those children with multiple events during the time they remained in PFK, the median time between the first and second events was 191 days (IQR 69, 389). The median probability of death, according to the trauma mortality prediction model (TMPM-ICD9), was higher at the first event (3.1%) than at subsequent events (<2% for all).

Risk of recidivism
Results of univariate comparisons of demographics, comorbidities and injuries between children with a single event and children with recurrent events are shown in Table 1. Factors independently associated with the risk of suspected recurrent NAT based on multivariable modeling are shown in Table 2. Living in a rural area and being less than 30 months of age at any event were associated with a higher risk of having a subsequent event ( Table 2). Having a dislocation, open wound, or superficial injury (p ≤ 0.01 for all) was associated with an increased risk of having a subsequent event. In addition, children who had 1-2 injuries at any event were more likely to have a subsequent event, whereas children with 3 or more injuries were not at increased risk for another event compared to children with no injury diagnoses. The most common body locations of dislocations, open wounds, and superficial injuries at the initial event among children with suspected recurrent NAT events were as follows: 14

Sensitivity analyses
These analyses were repeated in the subgroup of children who maintained continuous enrollment in PFK from their birth until the end of 2011 (n = 891, 65.5%). The results were similar, with the addition of injuries due to a fall now becoming a significant independent predictor of an increased risk of suspected recurrent NAT events in multivariable modeling (HR 1.42, 95% CI 1.07-1.88, p = 0.02). In order to account for potential selection bias due to informative censoring caused by the removal of children from their home after abuse, the analyses were repeated in only those children followed for at least 60 days after their first event. All results were similar. Finally, all analyses were also repeated in those children who did not have birth records in the PFK database. Compared to children with birth records in PFK, children without birth records tended to be older and to live in a zip code with a slightly higher median family income at their first event documented in PFK ( Table 3). Because of these differences and the likelihood that suspected NAT events were missed in children who entered into PFK after birth, these children were not included in the main analyses. There were several differences in the final multivariable models between those with and without birth records in PFK. Specifically, in addition to the risk factors previously identified, children who had a musculoskeletal disease (HR 1.42, 95% CI 1.03-1.92, p = 0.03) or a congenital anomaly (HR 1.44, 95% CI 1.04-2.00, p = 0.03) were more likely to experience a subsequent event.

Discussion
Many children who are victims of NAT may not experience it as a single isolated event, but rather as part of a pattern of recurrent violence that represents the normative structure of their social environment. This study used administrative claims data from a pediatric Medicaid accountable care organization to identify children with repeated medical encounters for injuries that are suspicious for NAT. We have identified several demographic and injury characteristics that are associated with an increased risk for suspected recurrent NAT events. These include living in a rural area, younger age at an event, fewer injuries  at an initial event, and specific injury categories including dislocations, superficial injuries, and open wounds. In addition, suspected recurrent NAT events were often observed months after the initial event and the time to a next event decreased with each subsequent event.
Missing child abuse at initial presentation can lead to significant subsequent morbidity [6,19]. With regards to NAT related to TBI, 30% of children hospitalized with abusive head injuries had a sentinel injury [12]. Data from our group using the Ohio State Trauma Registry suggests that victims of recurrent NAT who are hospitalized for their injuries have higher mortality rates compared to victims of single episodes of NAT (25% vs. 10%) [8]. By gaining a better understanding of the types and timing of injuries that portend risk to a child for recurrent NAT, we may be able to develop targeted screening tools and appropriate interventions that can be used to prevent recurrent NAT and its associated morbidity and mortality.
Previously identified risk factors for recurrent NAT include prior child protective services involvement, chronicity of maltreatment, child's age, and parental history including domestic violence, substance abuse, criminal   record, mental health issues, and being maltreated as a child [9][10][11]. In addition, several case series have described recurrence of maltreatment following specific injuries [5,6,12,13]. On a population level, Friedlaender et al., used Medicaid claims data to demonstrate that victims of maltreatment changed ambulatory care providers with greater frequency in the year before their first episode than those children who were not abused [20]. The current study is the first to utilize system-level administrative data to identify patterns of injuries and factors associated with suspected episodes of recurrent NAT. This populationbased approach allows us to examine all medical encounters for a patient, including episodes of care that occur outside the patient's usual hospital or health care system. Using this approach, we identified several trauma-related risk factors for suspected recurrent NAT. In this study, more than a quarter of children had a recurrent event within just one year of their first event. Risk factors a recurrent event included having fewer injuries (≤2 injuries) or having a dislocation, open wound, or and superficial injury at the previous event. These data potentially identify a bias in either the diagnosis of abuse and/or the variable response of child protective services to children based on the number and severity of physically evident injuries. Children with fewer or less severe injuries may not be reported to child protective services or are not removed from the unsafe environment leading to subsequent events. Identification of these more minor injuries as potential targets for additional screening or referral to child abuse specialists warrants further prospective study.
This study also found that the median length of time between the first and second suspected NAT events was 191 days (IQR 69, 389). This is important to note because the average length of child protective services involvement with a family may be significantly shorter. In addition, the risk of having a subsequent event increased with each event; 26% of children who experienced a first event proceeded to experience a second event within a year, whereas 60% children who experienced a 4th event proceeded to experience a fifth event within a  year. Understanding both the prolonged length of time between a first and second event, as well as the increasing risks with recurrent events may inform secondary prevention strategies for both medical and child welfare staff. Several limitations inherent in using system-level administrative claims data are relevant to this study. First, approximately 35% of patients had at least one break in Medicaid enrollment. In our analysis, we included these children as if they had remained continuously in the cohort throughout the study period. With this approach, there is the potential that children suffered a recurrent event during the time of non-enrollment, and therefore our data would be an underestimate of the number of recurrent events. However, the appropriateness of this assumption is increased by the finding of similar results in the subgroup of children with continuous enrollment in PFK from birth until the end of the study period. Second, some children who were removed from their home by child welfare after their initial event were lost to follow-up in this study. Whether or not a child remains in PFK after out-of-home placement varies by county in Ohio. Thus, we are unsure of the exact number of children lost to follow-up for this reason. However, when analyses were repeated in children who remained in PFK for at least 2 months after their first event, the results were very similar. Third, we are limited in the sensitivity and specificity of the ICD-9 coding practices used to identify key variables. In particular, ICD-9 coding performed after discharge is likely to underestimate the actual prevalence of abusive injuries in part because physicians may be reticent to assign intentional causality without confirmation from a multi-disciplinary team of social workers and law enforcement agents whose consensus is not often available until after discharge. In addition, ICD-9 codes provide limited ability to distinguish between different types of abuse. In this study, we aimed to focus on suspected physical abuse, but some of the codes chosen to define abuse could have certainly represented instances of emotional or sexual abuse, or child neglect. Furthermore, we were fairly liberal in our definition of potentially abusive injuries. Although some of the injury-only events could have involved accidental injuries, it is important to note that the rate of injuryonly episodes was remarkably high in this population (177 events per 1000 person-years), a rate more than 40 times the rate of 3.17 events per 1000 person-years that was reported in a general population of 0-3 year olds [21]. This extraordinarily high injury rate is concerning, regardless of whether the injuries were purposefully inflicted or represent neglect. An additional limitation of this study is that, administrative datasets provide limited data on covariates of interest. For example, this study would have benefited from additional data on race, parental characteristics, and family-level SES characteristics. By integrating US Census data, however, it was possible to evaluate zip-code level SES characteristics. Although the above limitations were unavoidable in the use of this administrative database, it is likely that they mainly resulted in under-identification of suspected NAT events and therefore minimized, rather than exaggerated our findings.

Conclusion
Factors associated with an increased risk for suspected recurrent NAT events in this study include rural residence, younger age, fewer initially detected injuries, and specific injury types at a previous event. Recurrent events often occur months after the initial event. These findings potentially identify a bias in either the diagnosis of NAT or the response of child protective services to children who present with less severe or less numerous injuries.

Additional file
Additional file 1: STROBE Statement-Checklist of items that should be included in reports of cohort studies.

Competing interests
The authors declare that they have no competing interests.
Authors' contributions KJD conceptualized and designed the study, interpreted the data, drafted parts of the initial manuscript, reviewed and revised the manuscript, and approved the final manuscript as submitted. JT interpreted the data, drafted parts of the initial manuscript, reviewed and revised the manuscript, and approved the final manuscript as submitted. JIG interpreted the data, critically reviewed and revised the manuscript, and approved the final manuscript as submitted. JNC acquired and analyzed the data, interpreted the data, drafted parts of the initial manuscript, reviewed and revised the manuscript, and approved the final manuscript as submitted. PCM conceptualized and designed the study, interpreted the data, drafted parts of the initial manuscript, reviewed and revised the manuscript, and approved the final manuscript as submitted. All authors read and approved the final manuscript.