MINI REVIEW article
Multivariate Analysis and Machine Learning in Cerebral Palsy Research
- Department of Neurology, Washington University in St. Louis, St. Louis, MO, United States
Cerebral palsy (CP), a common pediatric movement disorder, causes the most severe physical disability in children. Early diagnosis in high-risk infants is critical for early intervention and possible early recovery. In recent years, multivariate analytic and machine learning (ML) approaches have been increasingly used in CP research. This paper aims to identify such multivariate studies and provide an overview of this relatively young field. Studies reviewed in this paper have demonstrated that multivariate analytic methods are useful in identification of risk factors, detection of CP, movement assessment for CP prediction, and outcome assessment, and ML approaches have made it possible to automatically identify movement impairments in high-risk infants. In addition, outcome predictors for surgical treatments have been identified by multivariate outcome studies. To make the multivariate and ML approaches useful in clinical settings, further research with large samples is needed to verify and improve these multivariate methods in risk factor identification, CP detection, movement assessment, and outcome evaluation or prediction. As multivariate analysis, ML and data processing technologies advance in the era of Big Data of this century, it is expected that multivariate analysis and ML will play a bigger role in improving the diagnosis and treatment of CP to reduce mortality and morbidity rates, and enhance patient care for children with CP.
Cerebral palsy (CP) is the most common movement disorder in children (1) and causes the most severe physical disability in neurodevelopmental disorders (2). Among children with CP, around 33% of them can not walk, 25% can not talk, 25% have epilepsy, 50% have an intellectual disability, and most of them are in pain (3). Spastic CP is the most common subtype of this disorder, often shown as muscle stiffness that causes movement difficulties in a hand, arm, foot, or leg on one or both sides of the body, affecting the majority (>85%) of children with CP (4). Other subtypes of CP include dyskinetic (athetoid or dystonic) CP, ataxic CP (e.g., with tremors), and mixed CP.
Currently, there is no cure for CP, but medications (such as baclofen and botulinum toxin), supportive treatments (such as physical therapy), and surgical procedures [such as orthopedic surgery and selective dorsal rhizotomy (SDR)] can help patients alleviate symptoms and improve motor skills (5). The signs and symptoms of CP usually appear in the early months of life, but the average age for diagnosis of CP is around 2 years (2). Therefore, early identification and intervention is crucial for patients with CP because infants have higher potential for recovery from neural lesions than adults (1, 2). Neuroimaging, motor assessment (such as general movement assessment), and neurological examinations can help identify high-risk infants, monitor neurodevelopment, and detect or predict CP. Neuroimaging such as magnetic resonance imaging (MRI) and cranial ultrasound are useful to detect structural changes [intraventricular hemorrhage, periventricular leukomalacia (PVL), etc.] in the newborn brain, monitor lesion progression, and assess treatment effects, although compared with MRI, cranial ultrasound is less sensitive to lesions in the gray matter or malformations. Severe CP (caused by severe brain lesions such as PVL) can be identified by MRI or cranial ultrasound as soon as the lesions become recognizable on imaging after birth. However, 12–14% of children with CP have negative MRI scans due to subtle lesions in the brain (6). Thus, an integrated approach (imaging, motor assessment and neurological examinations) is needed to predict mild or moderate CP.
To predict CP in infants, Prechtl has described a general movement assessment method as a clinical assessment approach to identify CP motor impairments in infants by evaluating their spontaneous general movements (7). In particular, two atypical motor development features [(1) the presence of cramped-synchronized general movements at a preterm or term age and (2) the absence of small smooth movements or fidgety movements at 3–5 months] have been defined (8–10), which can identify CP in high-risk infants reliably (11). However, only well-trained physicians can perform such assessment, and general movement assessment based on visual observation by physicians is often influenced by subjective impressions and observer fatigue. Therefore, there is growing interest in developing multivariate and machine learning (ML)-based movement assessment tools for a more objective and quantitative motor assessment to detect movement impairments in high-risk infants (12, 13).
Multivariate analysis is a statistical analytical approach that simultaneously evaluates multiple variables, which compared with univariate analysis, may have more advantages (e.g., free from restrictions of various assumptions in univariate analysis) in identifying the associations between multiple data variables (e.g., variables associated with CP outcomes), grouping data into different groups or subgroups (e.g., different CP subtypes), and developing new diagnostic tests (e.g., differentiate CP subtypes with key feature variables). Multivariate analysis includes statistical methods such as principal components analysis (PCA), canonical variate analysis, independent components analysis, and multivariate regression. ML (or statistical learning) is a group of multivariate analytic methods that first identify the most significant data features or patterns that can best separate the data into different classes in the training dataset, and then apply these data features or patterns to the test dataset for data classification or prediction. ML has been increasingly applied to the biomedical field (14, 15), and examples of ML methods include linear discriminate analysis (16), support vector machine (SVM) (17), artificial neural networks (ANN) (18), random forest (19), and cluster analysis (20).
With growing interests, multivariate analysis has been increasingly employed in CP research in recent years, and research with multivariate analyses in CP is in infancy (14). To provide an overview of this relatively young field, PubMed search was performed with keywords “multivariate analysis cerebral palsy pediatric,” “machine learning cerebral palsy,” or “multivariate analysis cerebral palsy imaging.” The search yielded 126 articles. Articles were excluded if their subjects were not pediatric or the statistical methods used were not multivariate or the article was published before year 1990. This paper assessed the studies that used multivariate analysis in CP research and found that multivariate studies in CP are mainly in four categories: (1) risk factor identification; (2) detection of CP and identification of CP abnormalities; (3) movement assessment for CP prediction; and (4) outcome evaluation.
Multivariate Analysis in Risk Factor Identification
Early work on CP risk factor identification started from birth certificates. In a large population-based cohort study, data from birth certificates for 192 children with CP in four counties in California were compared with 155,636 healthy children in the same regions and the study found that low birth weight and (early or late) gestational age at birth were associated with high prevalence of CP, but early prenatal care and delivery at a hospital (for low birth weight neonates) were not associated prevalence of CP (21). Using multivariate analysis on clinical data of 113 CP infants (identified from 1,105 infants), Pinto-Martin et al. found that in low birth weight infants, cranial ultrasound imaging abnormalities such as parenchymal echodensities/lucencies (or ventricular enlargement) and germinal matrix/intraventricular hemorrhage were strong risk factors for disabling CP, but factors such as birth weight, gestational age, and Apgar score were not associated with it (22). In addition, a multicenter, large sample study (across eight European study centers, n = 585) revealed that there was a high rate of infection in mothers of CP children during their pregnancy and major CP abnormalities on structural MRI included white-matter damage due to immaturity (e.g., PVL) (42.5%), lesions in the basal ganglia (12.8%), cortical or subcortical lesions (9.4%), and malformation (9.1%) (23). A number of studies that identified CP risk factors have performed both univariate and multivariate analyses (24–26) where the risk factors identified by multivariate analyses were a subset of those identified by univariate analyses (24, 25), and the results of multivariate analyses were more rigorous and valid.
Studies that applied multivariate analysis to CP risk factor identification are summarized in Table 1 (22, 24–31). Multivariate logistic regression has often been used to identify risk factors for CP (22, 26, 29–31). Some of the risk factors identified by multivariate analytic studies include: premature birth, low birth weight, severe birth asphyxia, preterm rupture of membrane, abnormal cranial ultrasound, or structural MRI imaging findings (e.g., parenchymal echodensities/lucencies, or ventricular enlargement), intraventricular hemorrhage, PVL, neonatal sepsis, hypoxia-ischemic encephalopathy, hypoglycemia, neonatal jaundice, etc. (22–31). These risk factors are useful to understand the causes of CP, identify high-risk infants, and aid in the diagnosis of CP.
Table 1. Summary of studies with multivariate analyses in identification of risk factors and detection of CP.
Further, the risk factors for CP revealed by the multivariate studies are useful to prevent CP. Several CP risk factors such as brain injury and infection can be managed and avoided by preventing their causative mechanisms, and preventive efforts such as rubella vaccination, anti-D vaccination, and preventing methylmercury contamination are effective in preventing CP (2). In addition, meta-analysis has indicated that CP may be reduced by 30% in premature infants (<32 weeks gestation) by providing mothers of imminent-labor with magnesium sulfate for neuroprotection of their babies (43, 44). Further, early interventions such as hypothermia have prevented CP in 12.5% of infants with neonatal encephalopathy following an acute intrapartum hypoxic event (45). Since, currently, there is no cure for CP, CP prevention is critical to reduce the prevalence of CP and save children from CP and CP-caused life-long disabilities. Multivariate analysis may help identify significant risk factors and early interventions in order to prevent CP.
Taken together, multivariate analysis is important in identification of risk factors for CP and the risk factors identified such as premature birth and abnormal (cranial ultrasound or structural MRI) imaging findings are useful not only for CP cause identification and diagnosis but also for CP prevention. Further research is needed to identify more manageable and avoidable risk factors and early interventions (such as neuroprotective drugs or therapies) to prevent CP and reduce CP morbidity rate.
Multivariate Analysis in Detection of CP and Identification of CP Abnormalities
Since the current average age for diagnosis of CP is around 2 years and infants have higher potential for neural recovery (1, 2), early detection of CP is critical to make early intervention possible. To identify CP in high-risk infants, neuroimaging such as cranial ultrasound and MRI is important for lesion detection and deciding the timing of the lesion. Studies applied multivariate analytic methods to detection of CP and identification of CP abnormalities are summarized in the latter part of Table 1 (32–42). A multicenter study of very-low birth weight infants (n = 381, survival rate = 87%, 36, or 9.4% with CP) indicated that cranial ultrasound findings such as grade 3–4 intraventricular hemorrhage and PVL were useful in predicting CP; in particular, PVL and ventriculomegaly were related to high detection rate (≥30%) for CP (24). Further, to differentiate CP subtypes, Griffiths et al. examined the T2 MRI images of children with spastic or dyskinetic CP (n = 20 in each group), and found that patients with spastic CP had more severe injury to white matter near the paracentral lobule, while patients with dyskinetic CP had more injury to the subthalamic nucleus (STN) (33). Multivariate logistic regression further identified the associated factors (i.e., lesions in distinctly different anatomical locations) for differentiation of spastic and dyskinetic CP (33).
When brain injuries in patients with CP are subtle, advanced imaging such as diffusion tensor imaging (DTI) and diffusion-weighted imaging is useful to detect CP injuries with subtle abnormalities. DTI metrics such as mean diffusion (MD) and fractional anisotropy (FA) are often used to identify injuries in white matter tracts. The value of DTI in identifying degenerative changes in patients with spastic CP due to periventricular white matter injury has been demonstrated by an early study, which reported that children with spastic hemiparetic CP (caused by periventricular white matter injury) had reduced DTI fiber count on the ipsilateral (the same side as the lesion) side of the corticospinal tract (CST), corticobulbar tract (CBT), and superior thalamic radiation, and had MD and FA changes reflected neurodegeneration of the motor and sensory pathways (n = 5) (46). Further DTI studies found that white matter damage in the posterior thalamic radiation pathways was more severe than that in the CSTs in children with CP (n = 28) (47), and DTI abnormalities in several white matter tracts such as posterior thalamic radiation pathways or superior regions of the thalamocortical and corticomotor tracts correlated with motor function measured by, e.g., Gross Motor Function Classification System (GMFCS) level (n = 28–34) (34, 41, 47). A review paper summarized the results of 22 DTI studies in CP and reported common findings of decrease FA (or increased MD) in the corticomotor and sensorimotor pathways, which correlated with clinical measures (48). Some research findings suggested that the CP injury in the somatosensory circuits might be more severe than that in the motor circuits, which may contribute more to motor impairment in CP (49). Further research is needed to unfold the mechanisms underlie sensorimotor impairment in CP and to improve detection of CP through neuroimaging.
Apart from neuroimaging, the abnormalities of CP have been identified via multivariate analysis using patients’ clinical data (perinatal history, CP type, CP frequency, motor function, degree of disability, etc.) and data from other sources such as electromyography (EMG) and bone mineral density (BMD) (32, 35–40, 42). For example, multivariate analyses have indicated that 85% of CP patients have oropharyngeal dysphagia (36), 13% have sleep disturbance (37), 70% have abnormal optic nerve head (40), and 50% have low BMD (42). Factors associated with these CP abnormalities have also been identified (35, 37, 39).
In addition, since neonatal encephalopathy can cause CP, detection of neonatal encephalopathy helps detect potential CP. Structural brain connectivity networks of infants with neonatal encephalopathy have been examined using diffusion tractography extracted from DTI images, and ML methods such as SVM have been applied to structural connectivity features to detect neonatal encephalopathy (50). Moreover, since epilepsy and seizure disorders are common in children with CP, electroencephalography (EEG) is used to detect co-occurring seizures in high-risk infants or children. ML approaches such as SVM and ANN have been applied to EEG features to identify ictal and interictal spikes and achieved high detection rate for seizures (51). Further, multivariate analysis has found that children with CP after perinatal or neonatal stroke are more likely to have severe disability, cognitive impairment or epilepsy than CP children after delayed stroke (32).
Taken together, multivariate analytic studies in CP detection have identified imaging markers such as intraventricular hemorrhage and PVL on cranial ultrasound (24), injury to white matter near the paracentral lobule or to the STN on T2 MRI images (33), and injury in the CSTs and the posterior thalamic radiation pathways on DTI images (34, 41, 46, 47), which are useful in detecting CP and differentiating CP subtypes. Multivariate analyses have also identified CP abnormalities and their associated factors from non-imaging data (32, 35–40, 42). Further research is needed to identify biomarkers at the early stages of the disease to improve the diagnosis of CP, reduce diagnosis delay, and allow early identification and intervention for CP.
ML in Movement Assessment for CP Prediction
Studies reported applications of ML in movement assessment for CP prediction are summarized in Table 2 (52–61). Early identification of motor impairments in high-risk infants enables early detection of CP. The two atypical movement features (related to the cramped-synchronized general movements and the absence of fidgety movements) in general movement assessment are strong predictors for CP diagnosis (8–10). Based on these key motor impairment features, their movement characteristics and associated movement variables have been identified to detect movement impairments in high-risk infants (53–55, 57–59). ML approaches have made it possible to analyze recorded movement data and identify motor impairments automatically.
Table 2. Summary of studies with multivariate analytic and machine learning approaches in movement assessment and outcome evaluation in CP.
As a pioneer study, Meinecke et al. analyzed the 3D movement data of infants (n = 22, seven with CP) from video recordings, extracted an optimal combination of movement features with cluster analysis, and identified CP motor impairments with quadratic discriminant analysis (overall detection rate: 73%) (52). Further, the characteristics of fidgety movements and associated movement measurements have been identified to distinguish infants with movement impairments from those without (53, 54, 58). For example, extracting motion features related to fidgety movements (such as motion distance and relative frequency) from video recordings with an optical flow-based method, Stahl et al. examined the motion patterns of 82 infants (15 with CP), applied SVM classifier to detection of CP movement impairments, and achieved a good classification accuracy (93.7 ± 2.1% with features of relative frequency; sensitivity: 85.3 ± 2.8%; specificity: 95.5 ± 2.5%) (58).
Apart from video recordings, other movement recording systems such as accelerometers and electromagnetic movement tracking system have been employed for movement assessment and CP prediction. Heinze et al. examined the general movements of a group of newborns and infants (n = 23, 4 with CP) with accelerometers, selected optimal (combined) movement features with genetic algorithm, classified CP motor impairments with a decision tree-based classifier, and obtained overall detection rates of 88–92% (55). To distinguish the gait patterns between patients with CP (n = 4), patients with multiple sclerosis (n = 4), and healthy controls (n = 12), Alaqtash et al. extracted gait features from 3D ground reaction force data, compared the gait patterns of the three groups, and applied nearest-neighbor classifier and ANN to gait feature classification, which led to overall classification accuracies of 85% (ANN, with a combination of gait features) and 95% (after optimizing the gait features to an optimal set of six gait features) (56). Moreover, using electromagnetic movement tracking recordings, Karch et al. studied the general movements from 63 infants (10 with CP), extracted movement features such as joint centers, and computed stereotype scores with dynamic time warping, yielding a high CP classification accuracy with stereotype score of upper lime movement (sensitivity: 90%; specificity: 96%) (57). For a review on movement recognition techniques in general movement assessment for CP prediction in high-risk infants, see Ref. (12).
In addition, multivariate and ML approaches have been used in the assessment of physical therapy, and the effect of orthotic devices such as ankle foot orthosis on CP patients (55, 56). To evaluate the quality of exercises in CP physical therapy, Parmar and Morris (n = 5) applied four classifiers (SVM, neural networks, AdaBoosted decision tree, and dynamic time warping) to movement feature (joint and angle data in time or frequency domain) classification to identify correct or wrong exercise, and found that among the four classifiers, AdaBoosted decision tree performed the best with high classification accuracies (94.68% for joint data; 90.3% for angle data) (56).
However, these machine-learning-based movement assessment studies are at the early stage of research. For example, the subjects in the study of Alaqtash et al. (56) were healthy subjects (n = 5), the movement data of wrong exercises (with errors) were simulated data, thus, the results were preliminary. Further research is needed to apply machine-learning methods to real movement data of patients with CP. In addition, the sample sizes of patients with CP in these machine-learning-based movement assessment studies are small (n = 4–15), studies with large samples are needed to further verify and improve these machine-learning methods. Further, the performance of these machine-learning studies or classification systems may be improved by optimizing the data processing chains (feature extraction, feature selection, classification, and verification). For details on optimizing data processing chains to improve classification performance, see Ref. (79).
Recently, Marschik et al. proposed an integrated system called a fingerprint model that monitors the movement and speech-language development of new-born babies and infants at risk and automatically detects neurodevelopmental disorders such as CP by multidimensional data analysis and machine-learning approaches (13). Although it is challenging (in technical details), the fingerprint model enables neurological assessment of at-risk infants in an objective and quantitative manner and facilitates early detection of CP and other neurodevelopmental disorders, which may be the future direction of pediatric clinical practice.
Multivariate Analysis in CP Outcome Evaluation
Although there is no cure for CP, currently, treatment effects and outcomes in CP patients have been studied extensively. Multivariate approaches have been applied to outcome assessment (including survival analysis) in CP, and Table 2 (the latter part) provides a summary of these studies (62–78). The majority of the outcome studies employed a two-step approach: first, univariate analysis is used to identify variables that are associated with outcome; second, multivariate analysis is used to further examine the variables indicated by the univariate analysis and identify outcome predictors. A large sample of children with CP (n = 4,007) in UK were studied, and multivariate survival analysis indicated that the death rate was ~8, 85–94% of the children survived to age 20 years old, and the best predictors of CP survival were the number and severity of impairments (62). The multivariate outcome studies in CP fall into three categories: (1) outcome evaluation of medication and supportive treatments; (2) surgical outcome evaluation; and (3) quality of life (QOL) evaluation.
Outcome Evaluation of Medication and Supportive Treatments
The effect of commonly used medication oral baclofen on children with CP has been assessed with a multivariate model for the population pharmacokinetics analysis (n = 61), and it has been found that baclofen dosage based on body weight was appropriate to treat patients (≥2 years old), and determinants of apparent clearance in these children included body weight, a possible genetic factor, and age (70). Plasticity (shown as increased FA in the CSTs on DTI and improved motor function measures) induced by combined therapy (botulinum followed by physiotherapy) in children with spastic quadriplegia (n = 8) has been reported (80), while a later DTI study indicated that the addition of botulinum to physiotherapy did not influence the outcome at 6 months in children with spastic diplegic CP (n = 18) (81). DTI has also been used to evaluate the motor function outcomes of hemiplegic CP patients after rehabilitation treatment and DTI measurements such as the fiber number and FA of bilateral CSTs were correlated with functional level of hemiplegia scale (82). The quality of exercises in CP physical therapy has been evaluated with several classifiers, and AdaBoosted decision tree obtained good detection rate of exercise errors (61). In addition, multivariate outcome evaluation of therapies such as transcranial direct current stimulation (tDCS) has been performed (75). Grecco et al. investigated the functional outcome of tDCS in children with CP (n = 56), and multivariate logistic regression analyses identified that the presence of motor evoked potential was a predictor for walk test and gait speed, subcortical injury was a predictor for gross motor function, and both of them were predictors of motor function gain arise from tDCS combined with gait training in these patients (75). However, there are few outcome studies of medications and supportive treatments in CP using multivariate analyses, and multivariate analysis may play a bigger role in such outcome evaluation to reveal the true therapeutic effects of these treatments and their outcomes in CP patients.
Surgical Outcome Evaluation
A number of multivariate studies have investigated the outcomes of surgical procedures in CP (63, 64, 67, 71, 74–77). For example, the outcomes of SDR surgery have been studied in children with CP and outcome-associated factors have been identified (63, 64). Kim et al. examined factors associated with poor outcome of SDR surgery in pediatric patients with CP (n = 174) using multivariate logistic regression, and found that the poor-outcome rate was 6.3%, and the type of CP disability (diplegia, quadriplegia, etc.) was the predictor of poor outcome after SDR surgery (63). Further, Golan et al. evaluated the risk of spinal deformity in children with CP after SDR surgery (n = 98) and multivariate regression analysis identified several risk factors for spinal deformity: CP severity, ambulatory function, age at surgery, and gender (64). In addition, Long et al. investigated intraoperative opioid dosing and associated outcomes with multivariate regression in children with CP (n = 71) who underwent orthopedic surgery and reported that less intraoperative opioid was administered to CP children than non-CP children, and intraoperative opioid dosing was the outcome predictor of postoperative oxygen desaturation and ICU admission (67). Shore et al. examined the surgical outcome of children with CP (n = 320) who underwent proximal femoral varus derotation osteotomy (VDRO) for treatment of hip displacement, and found that success rate was 92% for GMFCS levels I and II vs. 76% for GMFCS level V, and multivariate analysis indicated that soft-tissue release during the VDRO procedure was the predictor of surgical success (73). The surgical outcome predictors identified by these studies are useful for outcome prediction for individual CP patients.
Complications after CP surgical treatments have also been investigated. Using multivariate regression analysis, Sponseller et al. studied deep wound infection after spinal fusion surgery at seven institutions in pediatric patients with CP (n = 204) and found that 6.4% of patients developed deep wound infection following surgery, and presence of a gastrostomy/gastrojejunostomy tube was the factor associated with infection (69). In addition, Kato et al. investigated cervical spine in patients with athetoid CP who underwent posterior decompression surgery (n = 31) and multivariate analysis showed that pedicle sclerosis associated with a higher risk of breach of cervical pedicle screws (71). Further, Minhas et al. evaluated the effect of body mass index class on complications after orthopedic surgery in children with CP (n = 1,746) and multivariate logistic regression analysis revealed that underweight status was the risk factor for complications in osteotomies and spine surgery (76). The risk factors identified by these studies are helpful to avoid the surgical complications and improve surgical treatments in CP.
In children with CP who underwent surgery, intraoperative neuromonitoring (IONM) often fail (failure rate 61%) (74). Mo et al. studied IONM in children with CP who underwent surgical scoliosis correction (n = 206) and multivariate logistic regression analysis revealed that PVL, hydrocephalus, and encephalomalacia were the predictors of poor IONM signals, while moderate or marked hydrocephalus and encephalomalacia were the predictors of no signals (74). Further, outcome prediction of CP surgical procedures has been explored in a recent study. Galarraga et al. examined children with CP who underwent (hip, ankle, foot, etc.) surgery (n = 115), and multi-regression analysis revealed that preoperative and surgical data could predict postoperative kinematics, and mean prediction errors (varying from 4° to 10°) were smaller compared with the variability of gait parameters (77). These results are encouraging because they indicated that the postsurgical kinematics of patients with CP could be predicted (relatively accurately with small mean prediction errors) using presurgical and surgical data, which allows an estimate of postsurgical outcome ahead of time.
Quality of life in physical ability, intellectual ability, self-care, and other aspects of life is an important outcome in CP. Multivariate analysis has been frequently used to assess QOL in patients with CP (65, 66, 68, 72, 78), and factors associated with physical QOL and self-care have been identified. For example, a multivariate analysis on QOL data of infants with CP (n = 92) identified GMFCS and intellectual capacity as the associated factors of self-care activity development, and GMFCS as the associated factors of mobility activities development (72). Further, a recent multivariate analysis showed that physical activity was positively associated with physical and total QOL in patients with CP (n = 128), and walking performance was positively associated with physical QOL (78). The factors identified by these studies may improve the QOL of patients with CP.
Taken together, since there is no cure for CP yet, and the death rate of CP is high (~8%), there is much to do to improve the outcomes of CP, and multivariate analytic approaches may play a bigger role in meeting such clinical demands. Surgical outcome predictors and risk factors for complications in CP surgical treatments have been identified by a number of multivariate outcome studies (63, 64, 67, 69, 71, 73, 76), which are useful not only for outcome evaluation and prediction but also for avoiding complications and improving surgical treatments in CP. However, there are few outcome studies for medications and supportive treatments (such as physical therapy) in CP using multivariate analysis. Thus, further research is needed to evaluate the outcomes of medications and supportive treatments, and multivariate analysis may play a bigger role in such outcome evaluation to reveal the true therapeutic effects of these treatments and their outcomes in CP patients, and help improve the outcomes of these treatments for patients with CP.
Multivariate analysis has been applied to several areas in CP research such as identification of risk factors for CP, detection of CP and identification of CP abnormalities, movement assessment for CP prediction, and outcome assessment. The studies reviewed in this paper have demonstrated that multivariate analytic and ML approaches have made it possible to analyze movement recordings and identify CP movement impairments automatically. In addition, outcome predictors for surgical treatments have been identified by multivariate outcome studies. To make the multivariate analytic and ML approaches useful in clinical settings, further research with large samples is needed to verify and improve these methods in CP detection, movement assessment, and outcome evaluation/prediction. As multivariate analysis, ML and data processing technologies advance in the era of Big Data, it is expected that multivariate analysis and ML will play a bigger role in improving the diagnosis and treatment of CP to reduce mortality and morbidity rates, and enhance patient care for children with CP.
JZ reviewed the multivariate analytic studies in cerebral palsy and wrote up the manuscript.
Conflict of Interest Statement
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This work was inspired by the clinical team led by Dr. Manish Shah at the University of Texas Health Science Center (in Houston) and affiliated hospitals who are dedicated to patient care for children with cerebral palsy. Proofread was kindly provided by Drs. Gary and Carla Brandenburger.
5. NIH (National Institute of Neurological Disorders and Stroke). Cerebral Palsy: Hope through Research. (2017). Available from: https://www.ninds.nih.gov/Disorders/Patient-Caregiver-Education/Hope-Through-Research/Cerebral-Palsy-Hope-Through-Research#3104_2
6. Krageloh-Mann I, Horber V. The role of magnetic resonance imaging in elucidating the pathogenesis of cerebral palsy: a systematic review. Dev Med Child Neurol (2007) 49:144–51. doi:10.1111/j.1469-8749.2007.00144.x
8. Einspieler C, Prechtl HFR. Prechtl’s assessment of general movements: a diagnostic tool for the functional assessment of the young nervous system. Ment Retard Dev Disabil Res Rev (2005) 11:61–7. doi:10.1002/mrdd.20051
9. Einspieler C, Marschik PB, Bos AF, Ferrari F, Cioni G, Prechtl HFR. Early markers for cerebral palsy: insights from the assessment of general movements. Future Neurol (2012) 7(6):709–17. doi:10.2217/fnl.12.60
12. Marcroft C, Khan A, Embleton ND, Trenell M, Plötz T. Movement recognition technology as a method of assessing spontaneous general movements in high risk infants. Front Neurol (2015) 5:284. doi:10.3389/fneur.2014.00284
13. Marschik PB, Pokorny FB, Peharz R, Zhang D, O’Muircheartaigh J, Roeyers H, et al. A novel way to measure and predict development: a heuristic approach to facilitate the early detection of neurodevelopmental disorders. Curr Neurol Neurosci Rep (2017) 17(5):43. doi:10.1007/s11910-017-0748-8
22. Pinto-Martin JA, Riolo S, Cnaan A, Holzman C, Susser MW, Paneth N. Cranial ultrasound prediction of disabling and nondisabling cerebral palsy at age two in a low birth weight population. Pediatrics (1995) 95(2):249–54.
24. Allan WC, Vohr B, Makuch RW, Katz KH, Ment LR. Antecedents of cerebral palsy in a multicenter trial of indomethacin for intraventricular hemorrhage. Arch Pediatr Adolesc Med (1997) 151(6):580–5. doi:10.1001/archpedi.1997.02170430046010
25. Kim JN, Namgung R, Chang W, Oh CH, Shin JC, Park ES, et al. Prospective evaluation of perinatal risk factors for cerebral palsy and delayed development in high risk infants. Yonsei Med J (1999) 40(4):363–70. doi:10.3349/ymj.19184.108.40.2063
28. Zhong Y, Wu J, Wu K, Wen R, Hou G, Peng D, et al. Evaluation of risk factors associated with cerebral palsy in children of Leshan Prefecture, Sichuan: a case control study on 308 cases. Zhonghua Yu Fang Yi Xue Za Zhi (2002) 36(5):323–6.
29. Miamoto CB, Pereira LJ, Paiva SM, Pordeus IA, Ramos-Jorge ML, Marques LS. Prevalence and risk indicators of temporomandibular disorder signs and symptoms in a pediatric population with spastic cerebral palsy. J Clin Pediatr Dent (2011) 35(3):259–63. doi:10.17796/jcpd.35.3.738t75v74l1m1p22
31. Yu T, Rong L, Wang Q, You Y, Fu JX, Kang LM, et al. Influence of neonatal diseases and treatments on the development of cerebral palsy in preterm infant. Sichuan Da Xue Xue Bao Yi Xue Ban (2013) 44(2):270–3.
32. Golomb MR, Saha C, Garg BP, Azzouz F, Williams LS. Association of cerebral palsy with other disabilities in children with perinatal arterial ischemic stroke. Pediatr Neurol (2007) 37(4):245–9. doi:10.1016/j.pediatrneurol.2007.06.003
33. Griffiths PD, Radon MR, Crossman AR, Zurakowski D, Connolly DJ. Anatomic localization of dyskinesia in children with “profound” perinatal hypoxic-ischemic injury. AJNR Am J Neuroradiol (2010) 31(3):436–41. doi:10.3174/ajnr.A1854
34. Yoshida S, Hayakawa K, Yamamoto A, Okano S, Kanda T, Yamori Y, et al. Quantitative diffusion tensor tractography of the motor and sensory tract in children with cerebral palsy. Dev Med Child Neurol (2010) 52(10):935–40. doi:10.1111/j.1469-8749.2010.03669.x
35. Coppola G, Fortunato D, Mainolfi C, Porcaro F, Roccaro D, Signoriello G, et al. Bone mineral density in a population of children and adolescents with cerebral palsy and mental retardation with or without epilepsy. Epilepsia (2012) 53(12):2172–7. doi:10.1111/j.1528-1167.2012.03639.x
36. Benfer KA, Weir KA, Bell KL, Ware RS, Davies PS, Boyd RN. Oropharyngeal dysphagia and gross motor skills in children with cerebral palsy. Pediatrics (2013) 131(5):e1553–62. doi:10.1542/peds.2012-3093
37. Romeo DM, Brogna C, Musto E, Baranello G, Pagliano E, Casalino T, et al. Sleep disturbances in preschool age children with cerebral palsy: a questionnaire study. Sleep Med (2014) 15(9):1089–93. doi:10.1016/j.sleep.2014.05.008
38. Adler C, Berweck S, Lidzba K, Becher T, Staudt M. Mirror movements in unilateral spastic cerebral palsy: specific negative impact on bimanual activities of daily living. Eur J Paediatr Neurol (2015) 19(5):504–9. doi:10.1016/j.ejpn.2015.03.007
39. Tao W, Zhang X, Chen X, Wu D, Zhou P. Multi-scale complexity analysis of muscle coactivation during gait in children with cerebral palsy. Front Hum Neurosci (2015) 9:367. doi:10.3389/fnhum.2015.00367
40. Ghate D, Vedanarayanan V, Kamour A, Corbett JJ, Kedar S. Optic nerve morphology as marker for disease severity in cerebral palsy of perinatal origin. J Neurol Sci (2016) 368:25–31. doi:10.1016/j.jns.2016.06.029
41. Reid LB, Cunnington R, Boyd RN, Rose SE. Surface-based fMRI-driven diffusion tractography in the presence of significant brain pathology: a study linking structure and function in cerebral palsy. PLoS One (2016) 11(8):e0159540. doi:10.1371/journal.pone.0159540
42. Tosun A, Erisen Karaca S, Unuvar T, Yurekli Y, Yenisey C, Omurlu IK. Bone mineral density and vitamin D status in children with epilepsy, cerebral palsy, and cerebral palsy with epilepsy. Childs Nerv Syst (2017) 33(1):153–8. doi:10.1007/s00381-016-3258-0
44. Crowther CA, Middleton PF, Voysey M, Askie L, Duley L, Pryde PG, et al. Assessing the neuroprotective benefits for babies of antenatal magnesium sulphate: an individual participant data meta-analysis. PLoS Med (2017) 14(10):e1002398. doi:10.1371/journal.pmed.1002398
46. Thomas B, Eyssen M, Peeters R, Molenaers G, Van Hecke P, De Cock P, et al. Quantitative diffusion tensor imaging in cerebral palsy due to periventricular white matter injury. Brain (2005) 128(Pt 11):2562–77. doi:10.1093/brain/awh600
47. Hoon AH Jr, Stashinko EE, Nagae LM, Lin DD, Keller J, Bastian A, et al. Sensory and motor deficits in children with cerebral palsy born preterm correlate with diffusion tensor imaging abnormalities in thalamocortical pathways. Dev Med Child Neurol (2009) 51(9):697–704. doi:10.1111/j.1469-8749.2009.03306.x
48. Scheck SM, Boyd RN, Rose SE. New insights into the pathology of white matter tracts in cerebral palsy from diffusion magnetic resonance imaging: a systematic review. Dev Med Child Neurol (2012) 54(8):684–96. doi:10.1111/j.1469-8749.2012.04332.x
50. Ziv E, Tymofiyeva O, Ferriero DM, Barkovich AJ, Hess CP, Xu D. A machine learning approach to automated structural network analysis: application to neonatal encephalopathy. PLoS One (2013) 8(11):e78824. doi:10.1371/journal.pone.0078824
51. Cabrerizo M, Ayala M, Goryawala M, Jayakar P, Adjouadi M. A new parametric feature descriptor for the classification of epileptic and control EEG records in pediatric population. Int J Neural Syst (2012) 22:1250001. doi:10.1142/S0129065712500013
52. Meinecke L, Breitbach-Faller N, Bartz C, Damen R, Rau G, Disselhorst-Klug C. Movement analysis in the early detection of newborns at risk for developing spasticity due to infantile cerebral palsy. Hum Mov Sci (2006) 25(2):125–44. doi:10.1016/j.humov.2005.09.012
54. Adde L, Helbostad JL, Jensenius AR, Taraldsen G, Grunewaldt KH, Støen R. Early prediction of cerebral palsy by computer-based video analysis of general movements: a feasibility study. Dev Med Child Neurol (2010) 52(8):773–8. doi:10.1111/j.1469-8749.2010.03629.x
55. Heinze F, Hesels K, Breitbach-Faller N, Schmitz-Rode T, Disselhorst-Klug C. Movement analysis by accelerometry of newborns and infants for the early detection of movement disorders due to infantile cerebral palsy. Med Biol Eng Comput (2010) 48(8):765–72. doi:10.1007/s11517-010-0624-z
56. Alaqtash M, Sarkodie-Gyan T, Yu H, Fuentes O, Brower R, Abdelgawad A. Automatic classification of pathological gait patterns using ground reaction forces and machine learning algorithms. Conf Proc IEEE Eng Med Biol Soc (2011):453–7. doi:10.1109/IEMBS.2011.6090063
57. Karch D, Kang K-S, Wochner K, Philippi H, Hadders-Algra M, Pietz J, et al. Kinematic assessment of stereotypy in spontaneous movements in infants. Gait Posture (2012) 36(2):307–11. doi:10.1016/j.gaitpost.2012.03.017
58. Stahl A, Schellewald C, Stavdahl O, Aamo OM, Adde L, Kirkerod H. An optical flow-based method to predict infantile cerebral palsy. IEEE Trans Neural Syst Rehabil Eng (2012) 20(4):605–14. doi:10.1109/TNSRE.2012.2195030
59. Kanemaru N, Watanabe H, Kihara H, Nakano H, Nakamura T, Nakano J, et al. Jerky spontaneous movements at term age in preterm infants who later developed cerebral palsy. Early Hum Dev (2014) 90(8):387–92. doi:10.1016/j.earlhumdev.2014.05.004
60. Wahid F, Begg R, Sangeux M, Halgamuge S, Ackland DC. The effects of an ankle foot orthosis on cerebral palsy gait: a multiple regression analysis. Conf Proc IEEE Eng Med Biol Soc (2015) 5509–12. doi:10.1109/EMBC.2015.7319639
63. Kim HS, Steinbok P, Wickenheiser D. Predictors of poor outcome after selective dorsal rhizotomy in treatment of spastic cerebral palsy. Childs Nerv Syst (2006) 22(1):60–6. doi:10.1007/s00381-005-1160-2
64. Golan JD, Hall JA, O’Gorman G, Poulin C, Benaroch TE, Cantin MA, et al. Spinal deformities following selective dorsal rhizotomy. J Neurosurg (2007) 106(6 Suppl):441–9. doi:10.3171/ped.2007.106.6.441
66. White-Koning M, Arnaud C, Dickinson HO, Thyen U, Beckung E, Fauconnier J, et al. Determinants of child-parent agreement in quality-of-life reports: a European study of children with cerebral palsy. Pediatrics (2007) 120(4):e804–14. doi:10.1542/peds.2006-3272
68. Smits DW, Ketelaar M, Gorter JW, van Schie P, Dallmeijer A, Jongmans M, et al. Development of daily activities in school-age children with cerebral palsy. Res Dev Disabil (2011) 32(1):222–34. doi:10.1016/j.ridd.2010.09.025
69. Sponseller PD, Jain A, Shah SA, Samdani A, Yaszay B, Newton PO, et al. Deep wound infections after spinal fusion in children with cerebral palsy: a prospective cohort study. Spine (Phila Pa 1976) (2013) 38(23):2023–7. doi:10.1097/BRS.0b013e3182a83e59
70. He Y, Brunstrom-Hernandez JE, Thio LL, Lackey S, Gaebler-Spira D, Kuroda MM, et al. Population pharmacokinetics of oral baclofen in pediatric patients with cerebral palsy. J Pediatr (2014) 164(5):1181–8. doi:10.1016/j.jpeds.2014.01.029
71. Kato S, Shoda N, Chikuda H, Seichi A, Takeshita K. Morphological characteristics of cervical spine in patients with athetoid cerebral palsy and the accuracy of pedicle screw placement. Spine (Phila Pa 1976) (2014) 39(8):E508–13. doi:10.1097/BRS.0000000000000234
72. Kruijsen-Terpstra AJ, Ketelaar M, Verschuren O, Smits DW, Jongmans MJ, Gorter JW. Determinants of developmental gain in daily activities in young children with cerebral palsy. Phys Occup Ther Pediatr (2015) 35(3):265–79. doi:10.3109/01942638.2014.957429
73. Shore BJ, Zurakowski D, Dufreny C, Powell D, Matheney TH, Snyder BD. Proximal femoral varus derotation osteotomy in children with cerebral palsy: the effect of age, gross motor function classification system level, and surgeon volume on surgical success. J Bone Joint Surg Am (2015) 97(24):2024–31. doi:10.2106/JBJS.O.00505
74. Mo AZ, Asemota AO, Venkatesan A, Ritzl EK, Njoku DB, Sponseller PD. Why no signals? Cerebral anatomy predicts success of intraoperative neuromonitoring during correction of scoliosis secondary to cerebral Palsy. J Pediatr Orthop (2015) 37(8):e451–8. doi:10.1097/BPO.0000000000000707
75. Grecco LA, Oliveira CS, Galli M, Cosmo C, Duarte Nde A, Zanon N, et al. Spared primary motor cortex and the presence of MEP in cerebral palsy dictate the responsiveness to tDCS during gait training. Front Hum Neurosci (2016) 10:361. doi:10.3389/fnhum.2016.00361
76. Minhas SV, Chow I, Otsuka NY. The effect of body mass index on postoperative morbidity after orthopaedic surgery in children with cerebral palsy. J Pediatr Orthop (2016) 36(5):505–10. doi:10.1097/BPO.0000000000000475
78. Mann K, Tsao E, Bjornson KF. Physical activity and walking performance: influence on quality of life in ambulatory children with cerebral palsy (CP). J Pediatr Rehabil Med (2016) 9(4):279–86. doi:10.3233/PRM-160395
80. Chaturvedi SK, Rai Y, Chourasia A, Goel P, Paliwal VK, Garg RK, et al. Comparative assessment of therapeutic response to physiotherapy with or without botulinum toxin injection using diffusion tensor tractography and clinical scores in term diplegic cerebral palsy children. Brain Dev (2013) 35(7):647–53. doi:10.1016/j.braindev.2012.10.012
81. Trivedi R, Gupta RK, Shah V, Tripathi M, Rathore RK, Kumar M, et al. Treatment-induced plasticity in cerebral palsy: a diffusion tensor imaging study. Pediatr Neurol (2008) 39(5):341–9. doi:10.1016/j.pediatrneurol.2008.07.012
Keywords: multivariate analysis, machine learning, cerebral palsy, early diagnosis, outcome assessment
Citation: Zhang J (2017) Multivariate Analysis and Machine Learning in Cerebral Palsy Research. Front. Neurol. 8:715. doi: 10.3389/fneur.2017.00715
Received: 07 September 2017; Accepted: 11 December 2017;
Published: 21 December 2017
Edited by:Andrea Guzzetta, University of Pisa, Italy
Reviewed by:D. Mishra, Maulana Azad Medical College, India
Finn Lennartsson, Lund University, Sweden
Copyright: © 2017 Zhang. 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) or licensor 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: Jing Zhang, email@example.com