 Research
 Open Access
 Published:
Classification of caesarean section and normal vaginal deliveries using foetal heart rate signals and advanced machine learning algorithms
BioMedical Engineering OnLine volume 16, Article number: 89 (2017)
Abstract
Background
Visual inspection of cardiotocography traces by obstetricians and midwives is the gold standard for monitoring the wellbeing of the foetus during antenatal care. However, inter and intraobserver variability is high with only a 30% positive predictive value for the classification of pathological outcomes. This has a significant negative impact on the perinatal foetus and often results in cardiopulmonary arrest, brain and vital organ damage, cerebral palsy, hearing, visual and cognitive defects and in severe cases, death. This paper shows that using machine learning and foetal heart rate signals provides direct information about the foetal state and helps to filter the subjective opinions of medical practitioners when used as a decision support tool. The primary aim is to provide a proofofconcept that demonstrates how machine learning can be used to objectively determine when medical intervention, such as caesarean section, is required and help avoid preventable perinatal deaths.
Methods
This is evidenced using an open dataset that comprises 506 controls (normal virginal deliveries) and 46 cases (caesarean due to pH ≤ 7.20—acidosis, n = 18; pH > 7.20 and pH < 7.25—foetal deterioration, n = 4; or clinical decision without evidence of pathological outcome measures, n = 24). Several machinelearning algorithms are trained, and validated, using binary classifier performance measures.
Results
The findings show that deep learning classification achieves sensitivity = 94%, specificity = 91%, Area under the curve = 99%, Fscore = 100%, and mean square error = 1%.
Conclusions
The results demonstrate that machine learning significantly improves the efficiency for the detection of caesarean section and normal vaginal deliveries using foetal heart rate signals compared with obstetrician and midwife predictions and systems reported in previous studies.
Introduction
Worldwide, over 130 million babies are born each year. 3.6 million will die due to perinatal complication and 1 million of these will be intrapartum still births [1]. In the USA, the number of deliveries in 2012 was 3952,841; one in every 164 of these resulted in stillbirth.^{Footnote 1} In the UK, in the same year, there were 671,255 with one in every 200 being stillbirth^{Footnote 2} and 300 that died in the first four weeks of life [2].
Cardiotocography (CTG) is the most common method used to monitor the foetus during the early stages of delivery [3] and clinical decisions are made using the visual inspection of CTG traces. However, the main weakness with this approach is poor human interpretation which leads to high inter and intraobserver variability [4]. While significant pathological outcomes like hypoxia are uncommon, false alarms are not, which can lead to serious abnormalities, such as cardiopulmonary arrest, brain and vital organ damage, cerebral palsy, hearing, visual and cognitive defects and in severe cases, death, being overlooked [5]. Conversely, over interpretation of CTG is common and the direct cause of unnecessary caesarean sections (CS). In such cases, between 40 and 60% of babies are born without any evidence to support pathological outcomes, such as hypoxia and metabolic acidosis [6].
This paper aims to address this problem by incorporating a proofofconcept system alongside existing gold standard methods in antenatal care. Using foetal heart rate signals and machine learning an objective measure of foetal state is used to detect the onset of pathological cases. This will provide obstetricians and midwives with an additional level of foetal state interpretation and help decide if and when surgical intervention is required. The results show that the approach has superior predictive capacity when compared with the 30% positive predictive value produced by obstetricians and midwives when classifying normal vaginal and caesarean section deliveries [7].
The remainder of this paper is organized as follows: “Background” provides a brief description of CTG as a screening tool for foetal hypoxia and its causes, and the findings of related work. “Methods” describes the dataset adopted in this paper and various steps involved. “Results” presents the classification outcomes for both control and case records, and the findings are discussed in “Discussion”. The paper is concluded in “Conclusions and future work”.
Background
CTG was initially developed as a screening tool to predict foetal hypoxia [7]. However, there is no evidence to suggest that there has been any improvement in perinatal deaths since the introduction of CTG into clinical practice 45 years ago. It is generally agreed that 50% of birthrelated brain injuries are preventable, with incorrect CTG interpretation leading the list of causes [8]. Equally, over interpretation of CTG is common and the direct cause of unnecessary caesarean sections, which costs the NHS £1700 for each caesarean performed compared with £750 for a normal vaginal delivery. It is therefore generally agreed that predicting adverse pathological outcomes and diagnosing pathological outcomes earlier clearly have important consequences, for both health and the economy. One interesting approach is machine learning.
Warrick et al. [9] developed a system for the classification of normal and hypoxic foetuses by modelling the FHR and uterine contraction (UC) signal pairs as an input–output system to estimate their dynamic relation in terms of an impulse response function [10]. The authors report that their system can detect almost half of the pathological cases 1 h and 40 min prior to delivery with a 7.5% false positive rate. Kessler et al. [11] on the other hand, using 6010 high risk deliveries, combined CTG with ST waveform to apply timely intervention for caesarean or vaginal delivery, which they report, reduced foetal morbidity and mortality [12].
In comparative studies, Huang et al. [13] compared three different classifiers; a decision tree (DT), an artificial neural network (ANN), and discriminant analysis (DA). Using the ANN classifier, it was possible to obtain a 97.78% overall accuracy. This was followed by the DT and DA with 86.36 and 82.1% accuracy respectively. The sensitivity and specificity values were not provided making accuracy alone an insufficient performance measure for binary classifiers. This is particularly true in evaluations where datasets are skewed in favour of one class with significant differences between prior probabilities.
In a similar study, Ocak et al. [14] evaluated an SVM and genetic algorithm (GA) classifier and reported 99.3 and 100% accuracies for normal and pathological cases respectively. Similar results were reported in [15, 16]. Again, Sensitivity and Specificity values where not provided in these studies. Meanwhile Menai et al. [17] carried out a study to classify foetal state using a naive bayes (NB) classifier with four different feature selection (FS) techniques: Mutual information, correlationbased, ReliefF, and information gain. The study found that the NB classifier in conjunction with features produced using the ReliefF technique produce the best results when classifying foetal state with 93.97, 91.58, and 95.79% for accuracy, sensitivity and specificity, respectively. While the results are high, the dataset is multivariate and highly imbalanced. Alternative model evaluation metrics for multiclass data, such as micro and macroaveraging, and micro and macroFmeasure, would provide a more informed account of model performance. Furthermore, an appropriate account of how the class skew problem was addressed is missing.
The adaptive boosting (AdaBoost) classifier was adopted in a study by Karabulut et al. [18] who report an accuracy of 95.01%—again no sensivity or specificity values were provided. While Spilka et al. who are the current forerunners of pioneering work in machine learning and CTG classification [6], used a random forest (RF) classifier in conjunction with latent class analysis (LCA) [19] and reported sensitivity and specificity values of 72 and 78% respectively using the CTGUHB dataset [3]. Producing slightly better results in [20] using the same dataset, Spilka et al. attempted to detect hypoxia using a C4.5 decision tree, naive bayes, and SVM. The SVM produced the best results using a tenfold cross validation method achieving 73.4% for sensitivity and 76.3% of specificity.
Methods
This section describes the dataset adopted in this study and discusses the steps taken to preprocess the data and extract the features from raw FHR signals. The section is then concluded with a discussion on the feature selection technique and dimensionality reduction.
CTG data collection
Chudacek et al. [3] conducted a comprehensive study that captured intrapartum recordings between April 2010 and August 2012. The recordings were collected from the University Hospital in Brno (UHB), in the Czech Republic by obstetricians with the support of the Czech Technical University (CTU) in Prague. These records are publically available from the CTUUHB database, in Physionet [3].
The CTUUHB database contains 552 CTG recordings for singleton pregnancies with a gestational age less than 36 weeks that were selected from 9164 recordings. The STAN S21/S31 and Avalon FM 40/50 foetal monitors were used to acquire the CTG records. The dataset contains no prior known development factors (i.e. they are ordinary clean obstetrics cases); the duration of stage two labour is less than or equal to 30 min; foetal heart rate signal quality is greater than 50% in each 30 min’ window; and the pH umbilical arterial blood sample is available. In the dataset, 46 caesarean section deliveries are included and the rest are ordinary clean vaginal deliveries. Figure 1 shows a scatter plot of the dataset with eclipses defining the separation between both case and control groups. Note that in this study 46 cases are classified as a caesarean delivery due to pH ≤7.20—acidosis, n = 18; pH >7.20 and pH ≤7.25—foetal deterioration, n = 4; or clinical decision without evidence of pathological outcome measures, n = 24 as defined in [14]. Table 1 provides details of the outcome measures used in the CTUUHB database.
The ID column in Table 1 maps to the file number in the CTUUHB database, while the second column provides details about the mothers age. The pH column describes the umbilical artery pH value for each case and BDef provides information on the base deficit in extracellular fluid. The PCO_{2} describes the partial pressure of carbon dioxide. BE gives values for the base excess and finally Apgar scores are a subjective evaluation of the outcome of the delivery. For a more indepth discussion of the dataset and these parameters please refer to [3].
Each recording begins no more than 90 min before delivery. Each CTG record contains the FHR time series (measured in beats per minute) and uterine contraction (UC) signal—each sampled at 4 Hz. The FHR was obtained from an ultrasound transducer attached to the abdominal wall. In this study only the FHR signal is only considered in this study since it provides direct information about the foetal state.
Preprocessing
Each of the 552 FHR signal recordings were filtered using a 6th order lowpass butterworth filter with fc = 4 Hz and a cutoff frequency of 0.034 Hz. To correct the phase distortion introduced by a onepass filter, a twopass filter (forwards and reverse) was used to filter each of the signals. Noise, and missing values were removed using cubic Hermite spline interpolation [21].
FHR features
This section describes the statistical, higherorder statistical and higherorder spectral features extracted from the FHR signals.
Morphological features
The initial set of features considered are those defined by the international federation of gynecology and obstetrics^{Footnote 3} (FIGO) and the National Institute for Health and Care Excellence^{Footnote 4} (NICE). Consider a raw FHR time series signal X with length N, where X = {x_{n}, n = 1, 2…, N}, in which the virtual baseline mean (VBM), \( \overline{x} \) is defined as:
Such that \( \overline{\text{x}} \) can be used to remove accelerations and decelerations (\(if\,x_{n} > 10 + \overline{\text{x}} \;{\text{then}}{:}\,x_{n} = {\bar{\text{x}}} + 10;\;if\,x_{n} <  10 + \overline{\text{x}} \;{\text{then}}{:}\, x_{n} = {\bar{\text{x}}} +  10\) from the FHR signal so that the real baseline FHR (RBL) can be derived [22]:
where H and L are the upper and lower limits of the time series signal respectively, X is the signal and N is the length of the signal.
Using the RBL, FIGO accelerations and decelerations can be extracted. These are features commonly used by obstetricians to monitor the interplay between the sympathetic and parasympathetic systems. Accelerations and decelerations within X represent the transient increases and decreases (±15 bpm) that last for 15 s or more [23]. In the case of accelerations, this typically indicates adequate blood delivery and is reassuring for the obstetrician. Calculating accelerations in the signal is defined by:
where X is the signal, \( x_{i} \) is the ith element of X, RBL is the real baseline defined in (2), and D is the duration of time in which x _{ i } remains above RBL +15.
In contrast decelerations represent temporary decreases (−15 bpm) in FHR below the RBL that last for 15 s or more, which can indicate the presence of possible pathological outcomes such as, umbilical cord compression, hypoxia or acidosis [18]. The decelerations in the signal are calculated as:
where \( x_{i} \) is the ith element of signal X, RBL is the real baseline, and D is the time duration in which x _{ i } remains below RBL15.
Short and longterm variability (STV and LTV respectively) are further indicators used by obstetricians. The presence of both suggests an intact neuromodulation of the FHR and normal cardiac function and is one of the most reassuring measures in neonatal care [24]. When STV or LTV decreases or is absent, it can be a significant indicator for the presence of hypoxia or acidosis. Therefore, they are both considered to be important predictors. STV is calculated according to the following equation:
where M is the number of minutes contained in the X signal and \( R_{t} \) is defined as:
where H is the number of subintervals in 60 s (in this case H = 60/K), K is the sample frequency (4 Hz) multiplied by 2.5 s and \( \overline{S}_{j} \) is the average value of 2.5 s for a subinterval j = {1, 2,…,H}.
In contrast, LTV is defined as the difference between the minimum and maximum value in a 60s block and is averaged to the duration of the signal if it is more than 1 min long. LTV is defined as:
where N is the length of the X signal, b is 240 samples (60s blocks for a 4 Hz sample frequency).
Collectively, RBL, accelerations, decelerations, STV and LTV define the five main FIGO/NICE features used by obstetricians and midwives and are subsequently consider as predictors for separating caesarean section and normal vaginal deliveries in this study.
Time series features
FIGO feature sets are often extended in automated CTG analysis to include patterns in the signal that are not easily identifiable through visual inspection. Two useful timeseries features that have been heavily utilized in medical signal processing are root mean squares (RMS) and sample entropy (SampEn). RMS is a useful feature for estimating short term variability between accelerations and deceleration [25] and is commonly described for a signal X with length N as:
This feature is particularly good at estimating sympathetic/parasympathetic dominance where the later, in a similar way to decelerations, can indicate the presence of possible pathological incidences, such as hypoxia and acidosis.
Whereas, sample entropy, quantifies the nonlinear dynamics of the FHR and the loss of complexity in the FHR signal. Previous studies have reported that it is a worthwhile feature for determining if the foetus is deprived of oxygen [26]. Sample entropy is the negative natural logarithm of the conditional probability that a dataset of length N, having repeated itself for m samples within a tolerance of r, will also repeat itself for m + 1 samples. Based on the calculation in [27] the time series X that contains N points, \( x_{i} ,\;x_{2} \ldots ,\;x_{N} \) subsequences can be defined by length m, and given by: \( y_{i} = (x_{i} , \;x_{i + 1} , \ldots .,x_{i + m  1} ) \) where i = 1, 2, …, N – m + 1. This allows the following quantity to be defined: \( B_{i}^{m} \left( r \right) \) as \( (N  m  1)^{  1} \) times the number of vectors \( V_{j}^{m} \) within r of \( V_{i}^{m} \), where j ranges from 1 to N−m, and \( j \ne i \), to exclude selfmatches, followed by:
Similarly, \( A_{i}^{m} \left( r \right) \) is defined as \( (N  m  1)^{  1} \) times the number of vectors \( V_{j}^{m + 1} \) within r of \( V_{i}^{m + 1} \), where j ranges from 1 to N−m, and \( j \ne i \), and set:
The parameter SampEn(m, r) is then defined as:
Which can be estimated by the statistic:
where N is the length of the X signal, m is the length of sequences to be compared, and r is the tolerance for accepted matches.
Frequency domain features
To overcome signal quality variations in the FHR signal, due to electrode placement and the physical characteristics of subjects [28], frequency domain features have been studied using power spectral density (PSD) computed using fast fourier transform (FFT):
where X(f) contains the information for the signal and x(t) is obtained from X(f) using the inverse of the fourier transformation:
The most notable feature is the peak frequency (FPeak) within the PSD, which has been used extensively in heart rate variability studies [29]. It is regarded as a useful measure of variability and normal sympathetic and parasympathetic function. It describes the dominant frequency in the PSD that has the maximum spectral power. In this study, peak frequency is derived using Welch’s method [30]:
where \( {\text{s}}_{i} \left( i \right) \) is the power of the spectrum at bin i. As shown later in the paper, this feature has good discriminative capacity as a confounding coefficient.
Nonlinear features
Poincare plots are a geometrical representation of a time series that is also used extensively to measure heart rate variability [20]. This paper shows that it has excellent discriminatory capacity in CTG analysis. Unlike HRV where it is commonly used, in FHR the difference between two beats is given as NN rather than the RR interval. A line of identity is used as a 45° imaginary diagonal line on the plot and the points falling on the line have the property \( NN_{n} = NN_{n + 1} \) [31]. Three coefficients of the Poincare plot, SD1 (the standard deviation of points perpendicular to the axis of lineofidentity), SD2 (the standard deviation of points along the axis of lineofidentity) and SDRatio are used as features to describe the cloud of points in the plot. Fundamentally, SD1 and SD2 are directly related to the standard deviation of NN interval (SDNN) and the standard deviation of the successive difference of the NN interval (SDSD) that is given by:
where \( Y_{NN} \left( 0 \right) \) and \( Y_{NN} \left( 1 \right) \) describe the autocorrelation function for lag0 and lag1 of the NN interval, respectively. The mean of NN intervals is \( \overline{NN} \). Equation 16 shows that SD1 and SD2 measures are derived from the correlation and mean of the NN intervals time series with lag0 and lag1.
The SD1 feature is an index of instantaneous recording of the beattobeat shortterm variability (the parasympathetic action) and SD2 describes the longterm variability (the sympathetic action). SD1 and SD2 are combined to form the ratio of SD1/SD2 that shows the relation between short and longterm variations of NN intervals:
It is also possible to detect the existence of chaos in the FHR signal since the foetal heartbeat fluctuates on different time scales and has the property of being selfsimilar. In this study, the boxcounting dimension is used to estimate the dynamics of the FHR [32]. It is a quantitative measure of the morphological properties of a signal and its capacity that is determined by covering the signal with N boxes of side length r. The minimal number of optimally sized boxes required to cover the complete signal describes the boxcounting dimension coefficient such that:
where D is the box counting fractal dimension of the object, r is the side length of the box, and N(r) is the smallest number of boxes of side r to cover the time series signal.
Longterm timecorrelations or selfaffinity measures of the FHR signal have also proven in previous studies to be useful for separating normal and pathological cases [33]. In this study, detrend fluctuation analysis (DFA) is performed where the returned exponent value indicates the presence or absence of fractal properties, i.e. selfsimilarity. The DFA probes the signal at different time scales and provides a fractal scaling exponent x. First the times series is integrated as follows:
where y(k) is the cumulative sum of the ith sample and X_{avg} is the mean value of the entire signal. Windows are derived from y(k) of equal length n and linear approximations \( y_{n} \) are found using least squares fit (this represents a trend in a given window). The average fluctuation F(n) of the signal around the trend is given by:
The calculations are repeated for all values of n. In this instance the primary focus is the relation between F(n) and the size of the window n. In general F(n) increases with the size of window n.
Feature selection
Feature selection is performed using the recursive feature eliminator algorithm (RFE) [34]. In this study a feature set was derived from the raw FHR signals based on the feature definitions described and a model fit generated using the RFE algorithm (refer to Algorithm 1) [34].
Each feature within this set is ranked using its importance to the model where S is a sequence of ordered numbers, which are candidate values for the number of features to retain \( (S_{1} > S_{2, \ldots } ) \). This process is repeated and the \( S_{i} \) top ranked features are retained. The model is refit and the performance is reassessed. The top \( S_{i} \) features are used to fit the final model.
Synthetic minority oversampling technique
In a two class balanced dataset the prior probabilities will be equal for each. This is not the case for the CTUUHB dataset given there are 506 controls (majority class) and 46 cases (minority class). Classifiers are more sensitive to detecting the majority class and less sensitive to the minority class and this leads to biased classification [35]. Therefore, given a random sample taken from the dataset, the probability of a classifier classifying a foetus observation as a control will be much higher (91.6%–506/552) than the probability of it classifying a foetus observation as a case (8.3%–46/552). This imposes a higher cost for misclassifying the minority (predicting that a foetus is normal and the outcome being pathological) than the majority class, (predicting a foetus is pathological and the outcome being normal).
In order to address this problem, it is necessary to resample the dataset [36]. Various resampling techniques are available, and these include under sampling and over sampling. Under sampling reduces the number of records from the majority class to make it equal to the minority class—in this instance it would mean removing 460 records leaving us with a very small dataset. In contrast, data in the minority class can be increased using oversampling. In this study, the synthetic minority oversampling technique (SMOTE) as defined in Algorithm 2 is used rather than reducing the dataset further [37].
Several studies have shown that the SMOTE technique effectively solves the class skew problem [35, 38,39,40,41,42]. Using SMOTE the minority class (cases) is oversampled using each minority class record, in order to generate new synthetic records along line segments joining the k minority class nearest neighbours. This forces the decision region of the minority class to become more general and ensures that the classifier creates larger and less specific decision regions, rather than smaller specific ones. In [37] the authors indicated that this approach is an accepted technique for solving problems related to unbalanced datasets.
Machine learning classifiers
Deep learning classifier
Deep learning neural network architectures have recently proven to be very powerful classifiers [43]. To the best of our knowledge, this algorithm has not been used in CTG studies, and this paper is thus the first to consider its use in automated CTG analysis. A multilayer feedforward neural network architecture is used based on theoretical proofs in [44]. The supervised training phase is based on uniform adaptive optimized initialization that is determined by the size of the network. A Tansigmoid nonlinear activation function f is utilized and defined as:
where \( x_{i} \) and \( w_{i} \) represent the firing neuron’s input values and their weights, respectively; and α denotes the weighted combination.
The multinomial distribution is adopted with the crossentropy loss function, which is typically used for classification in deep learning. For each training example j, the objective is to minimize a loss function:
where W is the collection \( \{ W_{i} \}_{1:N  1} \), \( W_{i} \) denotes the weight matrix connecting layers i and i + 1 for a network of N layers. Similarly B is the collection \( \{ b_{i} \}_{1:N  1} \), where \( b_{i} \) denotes the column vector of biases for layer i + 1. In the case of cross entropy, the loss function can be calculated by:
where \( t^{(j)} \) and \( O^{(j)} \) are the predicted and actual outputs, respectively, training example j, y represents the output units, and O the output layer.
The process used in this study to minimize the loss function defined in (22) is stochastic gradient descent (SGD) (refer to Algorithm 3) [45].
To address the problem of overfitting the dropout regularization technique proposed in [45] is used. This ensures that during forward propagation, when a given training example is used, the activation of each neuron in the network is suppressed within probability P. This coefficient is typically <0.2 for input neurons and <=0.5 for hidden neurons. Dropout allows an exponentially large number of models to be averaged as an ensemble, which helps prevent overfitting and improve generalization.
Momentum and learning rate annealing are used to modify backpropagation to allow prior iterations to influence the current version. In particular a velocity vector, v, is defined to modify the updates:
where \( \theta \) describes the parameters W and B, \( \mu \) the momentum coefficient, and \( \alpha \) the learning rate. Using the momentum parameter helps to avoid local minima and any associated instability [46]. Learning rate annealing is used to gradually reduce the learning rate \( \alpha_{t} \) to “freeze” into a local minima in the optimized landscape and is based on the principles described in [47].
Fishers linear discriminant analysis classifier
Before the more advanced random forest classification model is considered this section discusses the fishers linear discriminant analysis (FLDA) classifier as a baseline classification model. FLDA finds a linear combination of features that determines the direction along which the two classes are best separated. In this study the criterion proposed by Fisher is used which is the ratio of betweenclass to withinclass variances. The data is projected onto a line, and the classification is performed in this onedimensional space. The projection maximizes the distance between the means of the two classes while minimizing the variance within each class:
where \( \alpha \) is the bias, W is calculated using Fishers LDA, and X is the training data without class labels such that \( f\left( y \right) \ge 0 \) for normal records and <0 for pathological records. W is derived from X such that the within class scatter matrix \( S_{W} \) is minimized by:
where C is the number of classes, \( X_{i} \) is the set of all points that belong to class i, \( \mu_{i} \) is the mean of class i, and \( X_{k} \) is the kth point of \( X_{i} \). The between class scatter matrix \( S_{B} \) is maximized by:
where C is the number of classes, \( N_{i} \) is the total number of points that belong to class i, \( \mu_{i} \) is the mean of class i, and \( \mu \) is the overall mean, i.e. the mean of the data when all classes are considered together.
Random forest classifier
Random forest (RF) classifiers have featured widely in biomedical research [11, 48,49,50]. They are based on an ensemble of many randomized decisiontrees that are used to vote on the classification outcome. Many studies have shown that they give classification accuracies comparable with the best current classifiers on many datasets. They are able to handle data with a large number of features. Those features that are important for classification are determined through the calculation of an importance score for each feature. Each decisiontree is randomized using a bootstrap statistical resampling technique, with random feature selection [51].
Given an M feature set, trees are constructed using m features randomly selected from the feature set at each node of the tree. The best split is calculated using these m features, which continues until the tree is fully grown without pruning. The procedure is repeated for all trees in the forest using different bootstrap samples of the data. Classifying new samples can then be achieved using a majority vote. The approach combines bagging with decision tree classifiers to achieve this (refer to Algorithm 4).
Performance measures
kfold cross validation is used as a prediction metric with fivefolds and 1 and 30 repetitions, respectively. The average performance obtained from 30 simulations is utilized. This number is considered, by statisticians, to be an adequate number of iterations to obtain an acceptable average. Let \( C_{k} \) denote the indices of the observations in part k, and \( n_{k} \) the number of observations in k: if n is a multiple of K, then \( n_{k} = n/K \). Compute:
where
and \( \widehat{x}_{i} \) is the fit for observation i, obtained from the data with part k removed.
Sensitivity (true positives) and specificity (true negatives) measure the predictive capabilities of classifiers in binary classification tests. Sensitivities refer to the true positive rate or recall rate (pathological cases). Specificities measure the proportion of true negatives (normal cases). Sensitivities are considered higher priority than specificities, in this study. It is important to predict a pathological case rather than missclassify a normal case. To evaluate the performance of classifiers fitted to imbalanced datasets the FMeasure is a useful metric that combines precision and recall into a single value with equal weighting on both measures [52].
The area under the curve (AUC) is an accepted performance metric that provides a value equal to the probability that a classifier will rank a randomly chosen positive instance higher than a randomly chosen negative one (this obviously assumes that positive ranges higher than negative). This has been chosen, as it is a suitable evaluation method for binary classification. Consider a classifier that gives estimates according to \( p(C_{i} x) \), it is possible to obtain values \( \{ a_{1} , \ldots , a_{n1} ; a_{i} = p\left( {C_{1} x} \right), x_{i} \in C_{1} \} \) and \( \{ b_{1} , \ldots , b_{n2} ; b_{i} = p\left( {C_{2} x} \right), x_{i} \in C_{2} \} \) and use them to measure how well separated the distributions of \( \widehat{p}(x) \) for class \( C_{1} \) and \( C_{2} \) patterns are [53].
Using the estimates, \( \{ a_{1} , \ldots , a_{n1} ,b_{1} , \ldots , b_{n1} \} \) they can be ranked in increasing order. The class \( C_{1} \) test points can be summed to see that the number of pairs of points, one from class \( C_{1} \) and one from \( C_{2} \) with \( \widehat{p}(x) \) smaller for class \( C_{2} \) than the \( \widehat{p}(x) \) value for class \( C_{1} \), is:
where \( r_{i} \) is the ranked estimate, \( S_{o} \) is the sum of the ranks of the class \( C_{1} \) test patterns. Since there are \( n_{1} n_{2} \) pairs, the estimate of the probability that a randomly chosen class \( C_{2} \) pattern has a lower estimated probability of belonging to class \( C_{1} \) than a randomly chosen class \( C_{1} \) is:
This is equivalent to the area under the ROC which provides an estimate obtained using the rankings alone and not threshold values to calculate it [51].
Results
This section presents the classification results for control and case records using the CTUUHB dataset. The features extracted from the FHR signals are used to model each of the classifiers. The performance is measured using sensitivity, specificity, AUC, FMeasure and MSE values.
Using all features from original data
In the first evaluation, all the features in the feature set are used to train the FLDA, RF and DL classifiers.
Classifier performance
The results in Table 2 show that the sensitivities for all classifiers are very low, while corresponding specificities are high. This is expected given that there are a limited number of case records from which the classifiers can learn a suitable fit. The FMeasure is a good metric when using imbalanced datasets and provides a better indication of classifier performance than sensitivity, specificity and AUC. As can be seen the FMeasure for the FLDA and RF are low with slightly better results produced by the DL (only slight better than chance).
It is clear that the models are capable of classifying control records but not case records. This is because there are 506 controls and only 46 cases from which the classifiers can learn, which is significantly lower. The AUC values are relatively low for the FLDA with slightly higher values for the RF and higher values again for the DL. Yet, the Sensitivities, which are considered more important in this study, are all low. Table 3 shows the error estimate for fivefold crossvalidation using both 1 and 30 repetitions.
The errors are consistent with the expected MSE baserate of 8.3% (46 pathological/552 FHR records) with the exception of the DL which produced an MSE = 3%.
Model selection
The receiver operator characteristic (ROC) curve is a useful graphic that shows the cutoff values for the false and true positive rates. It is particularly useful in binary classification to illustrate classifier performance. In the current evaluation, Fig. 2 shows that the FLDA performed poorly. The RF and DL classifiers produced slightly better results, which reflect the sensitivity, specificity and AUC values in Table 2.
The primary reason for the low sensitivities (despite the AUC for the RF and DL being relatively high) is that there are insufficient case records to model the class. This is in contrast to the classification of control records that are skewed in its favour. This causes significant problems in machine learning. As such, resampling the classes in the absence of real pathological cases is a conventional way of addressing this problem [54].
Using all features from synthetic minority oversampling technique data
The 46 case records are resampled using the SMOTE algorithm. The SMOTE algorithm generates a new dataset containing extra cases derived from the minority class while reducing the majority class samples accordingly. Figure 3 shows the separation of classes following oversampling. Compared with Fig. 1 it is clear that both case and control data are now evenly distributed between the two groups. There is significant overlap between case and controls and no two feature combinations were able to increase this decision boundary any further than that presented in Fig. 3.
Classifier performance
Using the new SMOTEd feature set (300 cases and 300 controls—empirically this distribution produced the best Sensitivity, Specificity, AUC, FMeasure and MSE results), Table 4 indicates that the Sensitivities for all models improved (90% in most cases). This is however at the expense of lower specificities (10% decreases). The results are encouraging given that accurately predicting cases is more important than predicting controls. The FMeasure acts as a support metric in this evaluation and produces encouraging results in the RF and DL classifiers.
Table 5 shows a marked improvement in error rates in all classifiers except the FLDA, which has increased by 12%. In the case of the DL classifier, the results indicate a 1.7% error rate, which is significantly less than the expected MSE baserate of 50% (300 cases/600 FHR records).
Model selection
The ROC curve in Fig. 4 illustrates that all the models have significantly improved with the exception of the FLDA where the overall performance remained more or less the same.
The results show that adopting the SMOTE oversampling technique improves classifier performance. While oversampling data is not ideal, it is an accepted technique within many clinical studies when skewed datasets need to be normally distributed [54,55,56].
The remaining evaluations build on these results with a particular focus on dimensionality reduction.
Using RFE selected features from SMOTE data
Recursive feature extraction (RFE)
Using the RFE algorithm, each feature is assessed to determine their discriminatory capacity. Figure 5 shows the crossvalidation results using various feature combinations.
The results indicate that the optimal number of features is eight as can been seen in Table 6.
The eight ranked features are DFA, RMS, FPeak, Acc, SD2, SDRatio, SAMPEN, and Dec. The following evaluation determines whether this reduced feature set can improve on or maintain the previous set of results.
Classifier performance
Looking at Table 7 it can be seen that most of the classifiers perform slightly worse using the eight features in terms of sensitivity. This is with exception to the RF classifier, which can maintain similar results using the reduced feature set.
The MSE values, reported in Table 8, for all but the RF classifier (whose error more or less stayed the same) were slightly worse than in the previous evaluation.
Model selection
In this final evaluation, the ROC curve in Fig. 6 illustrates that there are no real improvements on the previous evaluation for the FLDA and DL, but that the RF performs very well with a reduced set of features.
Discussion
Obstetricians and midwives visually inspect CTG traces to monitor the wellbeing of the foetus during antenatal care. However, inter and intraobserver variability and low positive prediction is accountable for the 3.6 million babies that die each year. This paper, presented a proofofconcept using machine learning and FHR signals as an ambulatory decision support to antenatal care. The results indicate that it is possible to provide high predictive capacity when separating normal vaginal deliveries and caesarean section deliveries and in many cases produce much better results than those reported in previous studies (see Table 9).
Using the original unbalanced dataset the best classifier (DL classifier) achieves SE = 0%, SP = 99%, AUC = 87%, and FMeasure = 52%. While the specificity values are high, all sensitivity values are below 3%. The low Sensitivity is attributed to the disproportionate number of normal records compared with pathological records and the fact that unbalanced datasets in general cause bias in favour of the majority class. The minimum error rate MSE = 3% was achieved by the DL using 30 repetitions. This relatively small MSE appeared to be a good error rate. However, the classifiers were simply classifying by minimizing the probability of error, in the absence of sufficient evidence to help them to classify otherwise.
The SMOTE algorithm using all 13 features significantly improved the Sensitivity values for all classifiers. While oversampling is not ideal, it is a way to solve the class skew problem that is widely used in medical data analysis [40, 61,62,63,64]. The best classification algorithm is again the DL classifier, which achieves SE = 94%, SP = 91%, AUC = 100%, FMeasure = 100% and MSE = 2%. The reason for this is that the algorithm has the ability to extract complex nonlinear patterns generally observed in physiological data like FHR signals. Through the extraction of these patterns, the DL algorithm uses relatively simpler linear models for data analysis tasks, such as classification. The DL generalizes, and finds the global minima, which allows it to generate learning patterns and relationships beyond immediate neighbours in data. It is able to provide much more complex representations of data by extracting representations directly from unsupervised data without domain knowledge or inference.
Using the RFE algorithm as a feature selection technique the algorithm eliminated five features from the original 13 that were considered to have very low discriminatory capacity. The remaining eight features were used to fit the models and the results show that the RF achieved the best overall results with SE = 91%, SP = 91%, AUC = 98%, FMeasure = 91% and MSE = 11%. The primary reason for these good results is that the RF algorithm is based on an ensemble of many randomized decisiontrees that are used to vote on the classification outcome. They are able to handle data with a very large number of features (although the feature set in this study is not particularly large) and those features that are important for classification can be determined through the calculation of an importance score for each feature. The score metric based on voting is similar to the approach adopted in knearest neighbour classification and the voting mechanism to classify new data points based on the majority surrounding data points of a particular class. The DL classifier performed worse on the reduced dataset but still produces better results than several studies discussed in this paper [6, 20, 57, 60].
Conclusions and future work
The primary aim in this paper was to evaluate a proofofconcept approach to separating caesarean section and normal vaginal deliveries using FHR signals and machine learning. The results show that using a deep learning classifier it is possible to achieve 94% for Sensitivity, 91% for Specificity, 99% for AUC, 100% for Fscore, and 1% for Mean Square Error. This shows significant improvements over the 30% positive predictive value achieved by obstetricians and midwives and warrants further investigation as a potential decision support tool for use alongside the current CTG gold standard.
Nonetheless, despite the encouraging results reported, the study needs further evaluation using truly independent data to fully assess its value. In future work this will be made possible by soliciting support for clinical trials and utilising other open datasets that have adopted a similar study design. Other important work will include regression analysis, using a larger number of classes to predict the expected pathological event, in terms of the number of hours or days to delivery, not just whether the outcome is likely to be a caesarean section or a normal vaginal delivery. We also need to integrate and use the clinical data provided with this study in future analysis tasks.
It will also be important to evaluate different parameter adjustment settings, particularly in the case of the DL algorithm to determine if the results can be further improved. Automatic feature detection will also be explored using the DL to extract features from the raw FHR signals.
It is less than ideal to use oversampled data. Therefore, another direction for future work will explore opportunities to obtain data through funded clinical trials. This will also help provide a much more indepth account of the value of machine learning and its perceived benefits on predicting caesarean section and normal vaginal deliveries.
While only the FHR signal is considered in this paper, since it provides direct information about the foetus’s state, it would be useful to combine this signal with the UC signal, which has been studied in previous work [64].
Overall, the proposed methodology is robust, contributes to the biomedical data analytics field, and provides new insights into the use of deep learning algorithms when analysing FHR traces that warrants further investigation.
Abbreviations
 CTG:

cardiotocography
 ST:

ST segment connects the QRS complex and the T wave
 UHB:

University Hospital in Brno
 CTU:

Czech Technical University
 CTUUHB:

Czech Technical University–University Hospital in Brno
 STAN S21/S31:

product name for CTG analysis
 Avalon FM 40/50:

foetal monitor
 FHR:

foetal heart rate
 UC:

uterine contraction
 FIGO:

international federation of gynecology and obstetrics
 NICE:

National Institute for Health and Care Excellence
 VBM:

virtual baseline
 RBL:

real baseline
 STV:

short term variability
 LTV:

long term variability
 RMS:

root mean squares
 SampEn:

sample entropy
 FFT:

fast fourier transform
 PSD:

power spectral density
 FPeak:

peak frequency
 HRV:

heart rate variability
 SD:

standard deviation
 DFA:

detrend fluctuation analysis
 SVM:

support vector machine
 DL:

deep learning
 RFE:

recursive feature eliminator
 SMOTE:

synthetic minority oversampling technique
 SGD:

stochastic gradient descent
 FLDA:

fishers linear discriminant analysis
 RF:

random forest
 AUC:

area under the curve
 ROC:

receiver operator curve
 MSE:

mean square errors
References
 1.
Warren JB, Lambert WE, Fu R, Anderson JM, Edelman AB. Global neonatal and perinatal mortality: a review and case study for the Loreto Province of Peru. Res Rep Neonatol. 2012;2:103–13.
 2.
Brown R, Wijekoon JHB, Fernando A, Johnstone ED, Heazell AEP. Continuous objective recording of fetal heart rate and fetal movements could reliably identify fetal compromise, which could reduce stillbirth rates by facilitating timely management. Med Hypotheses. 2014;83(3):410–7.
 3.
Chudacek B, Spilka J, Bursa M, Janku P, Hruban L, Huptych M, Lhotska L. Open access intrapartum CTG database. BMC Pregnancy Childbirth. 2014;14(16):1–12.
 4.
Ugwumadu A. Are we (mis)guided by current guidelines on intrapartum fetal heart rate monitoring? case for a more physiological approach to interpretation. Int J Obstet Gynaecol. 2014;121(9):1063–70.
 5.
Sola A, Golombek SG, Bueno MTM, LemusVarela L, Auluaga C, Dominquez F, Al E. Safe oxygen saturation targeting and monitoring in preterm infants: can we avoid hypoxia and hyperoxia? Acta Paediatr. 2014;103(10):1009–18.
 6.
Spilka J, Georgoulas G, Karvelis P, Chudacek V. Discriminating normal from ‘abnormal’ pregnancy cases using an automated FHR evaluation method. Artif Intell Methods Appl. 2014;8445:521–31.
 7.
Pinas A, Chadraharan E. Continuous cardiotocography during labour: analysis, classification and management. Best Pract Res Clin Obstet Gynaecol. 2016;30:33–47.
 8.
Warrick PA, Hamilton EF, Precup D, Kearney RE. Classification of normal and hypoxic fetuses from systems modeling of intrapartum cardiotocography. IEEE Trans Biomed Eng. 2010;57(4):771–9.
 9.
Menai ME, Mohder FJ, Almutairi F. Influence of feature selection on naïve bayes classifier for recognizing patterns in cardiotocograms. J Med Bioeng. 2013;2(1):66–70.
 10.
Koop G, Pesaran MH, Potter SM. Impulse response analysis in nonlinear multivariate models. J Econom. 1996;74(1):119–47.
 11.
Hand J, Till RJ. A Simple generalisation of the area under the ROC curve for multiple class classification problems. Mach Learn. 2001;45(2):171–86.
 12.
Kessler J, Moster D, Albrechfsen S. Delay in intervention increases neonatal morbidity in births monitored with cardiotocography and STwaveform analysis. Acta Obs Gynecol Scand. 2014;93(2):175–81.
 13.
Spilka J, Chudacek V, Koucky M, Lhotska L, Huptych M, Janku P, Georgoulas G, Stylios C. Using nonlinear features for fetal heart rate classification. Biomed Signal Process Control. 2012;7(4):350–7.
 14.
OmoAghoja L. Maternal and fetal acidbase chemistry: a major determinant of perinatal outcome. Ann Med Health Sci Res. 2014;4(1):8–17.
 15.
Kreyszig E. Advanced engineering mathematics. Hoboken: Wiley; 2005. p. 816.
 16.
Goldberger AL, Amaral LAN, Glass L, Hausdorff JM, Ivanov PC, Mark RG, Mietus JE, Moody GB, Peng CK, Stanley HE. PhysioBank, PhysioToolkit, and PhysioNet: components of a new research resource for complex physiologic signals. Circulation. 2000;101(23):e215–20. doi:10.1161/01.CIR.101.23.e215.
 17.
Mantel R, van Geijn HP, Caron FJ, Swartjies JM, van Woerden EE, Jongsma HW. Computer analysis of antepartum fetal heart rate: 2. Detection of accelerations and decelerations. Int J Biomed Comput. 1990;25(4):273–86.
 18.
Pinto P, Bernardes J, CostaSantos C, AmorimCosta C, Silva M, AyresdeCampos D. Development and evaluation of an algorithm for computer analysis of maternal heart rate during labor. Comput Biol Med. 2014;49(1):30–5.
 19.
Camm J. Heart rate variability: standards of measurement, physiological interpretation and clinical use. task force of the European society of cardiology and the North American society of pacing and electrophysiology. Circulation. 1996;93(5):1043–65.
 20.
Schiermeier S, Van Leeuwen P, Lange S, Geue D, Daumer M, Reinhard J, Gronemeyer D, Hatzmann W. Fetal heart rate variation in magnetocardiography and cardiotocography–a direct comparison of the two methods. Z Geburtshilfe Neonatol. 2007;211(5):179–84.
 21.
Signorini MG, Fanelli A, Magenes G. Monitoring fetal heart rate during pregnancy: contributions from advanced signal processing and wearable technology. Comput Math Methods Med. 2014;2014(707581):1–10.
 22.
Radomski D, Grzanka A, Graczyk S, Przelaskowski A. Assessment of uterine contractile activity during a pregnancy based on a nonlinear analysis of the uterine electromyographic signal. In: Pietka E, Kawa J, editors. Information technologies in biomedicine. Heidelberg: Springer; 2008. p. 325–31.
 23.
Maner WL, Garfield RE, Maul H, Olson G, Saade G. Predicting term and preterm delivery with transabdominal uterine electromyography. Obstet Gynecol. 2003;101(6):1254–60.
 24.
Williams DP, Jarczok MN, Ellis RJ, Hillecke TK. Twoweek test–retest reliability of the Polar^{®} RS800CXTM to record heart rate variability. Clin Physiol Funct Imaging. 2016. doi:10.1111/cpf.12321.
 25.
Welch PD. The use of fast fourier transform for the estimation of power spectra: a method based on time aver. aging over short, modified periodograms. IEEE Trans Audio Electoacoustics. 1967;15(2):70–3.
 26.
Karmakar CK, Khandoker AH, Gubbi J, Palaniswami M. Complex correlation measure: a novel descriptor for poincare plot. Biomed Eng Online. 2009. doi:10.1186/1475925X817.
 27.
Sarkar N, Chaudhuri BB. An efficient differential boxcounting approach to compute fractal dimension of image. IEEE Trans Syst Man Cybern. 1994;24(1):115–20.
 28.
Abry P, Roux SG, Chudacek V, Borgnat P, Goncalves P, Doret M. Hurst exponent and IntraPartum fetal heart rate: impact of decelerations. In: 26th IEEE international symposium on computerbased medical systems; 2013. p. 131–136.
 29.
Granitto PM, Bohorquez AB. Feature selection on wide multiclass problems using OVARFE. Intel Artif. 2009;44(2009):27–34.
 30.
Blagus R, Lusa L. SMOTE for highdimensional classimbalanced data. BMC Bioinform. 2013;14(106):1–16.
 31.
Khoshgoftaar TM, van Hulse J, Napolitano A. Comparing boosting and bagging techniques with noisy and imbalanced data. IEEE Trans Syst Man Cybern Syst. 2011;41(3):552–68.
 32.
Chawla NV, Bowyer KW, Hall LO, Kegelmeyer WP. SMOTE: synthetic minority oversampling technique. J Artif Intell Res. 2002;16:321–57.
 33.
Sun T, Zhang R, Wang J, Li X, Guo X. Computeraided diagnosis for earlystage lung cancer based on longitudinal and balanced data. PLoS ONE. 2013;8(5):63559.
 34.
Taft LM, Evans RS, Shyu CR, Eggar MJ, Chawla N. Countering imbalanced datasets to improve adverse drug event predictive models in labor and delivery. J Biomed Inform. 2009;42(2):356–64.
 35.
Lin W, Chen JJ. Classimbalanced classifiers for highdimensional data. Brief Bioinform. 2013;14(1):13–26.
 36.
Nahar J, Imam T, Tickle KS, Ali ABMS, Chen YP. Computational intelligence for microarray data and biomedical image analysis for the early diagnosis of breast cancer. Expert Syst Appl. 2012;39(16):12371–7.
 37.
Wang Y, Simon M, Bonde P, Harris BU, Teuteberg JJ. Prognosis of right ventricular failure in patients with left ventricular assist device based on decision tree with SMOTE. Trans Inf Technol Biomed. 2012;16(3):383–90.
 38.
Silver D, Huang A, Maddison CJ, Guez A, Sifre L, Al E. Mastering the game of go with deep neural networks and tree search. Nature. 2016;529(7587):484–9.
 39.
Candel A, Lanford J, LeDell E, Parmar V, Arora A. Deep learning with H20. 2015.
 40.
Srivastava N, Hinton G, Krizhevsky A, Sutskever I, Salakhutdinov R. Dropout: a simple way to prevent neural networks from overfitting. J Mach Learn Res. 2014;15:1929–58.
 41.
Sutskever I, Vinyals O, Le QV. Sequency to Sequency Learning with Neural Networks. In: 27th annual conference on neural information processing systems (NIPS); 2014. p. 1–9.
 42.
Zeiler MD. ADADELTA: an adaptive laerning rate method. arXiv.org; 2012.
 43.
Tomas P, Krohova J, Dohnalek P, Gajdos P. Classification of cardiotocography records by random forest. In: 36th IEEE international conference on telecommunications and signal processing; 2013. p. 620–623.
 44.
Peterek T, Gajdoš P, Dohnálek P, Krohová J. Human fetus health classification on cardiotocographic data using random forests. In: Pan JS, Snasel V, Corchado ES, Abraham A, Wang SL, editors. Intelligent data analysis and its applications, Volume II: proceeding of the first EuroChina conference on intelligent data analysis and applications, June 13–15, 2014, Shenzhen, China. Switzerland: Springer; 2014. p. 189–98.
 45.
Tetschke F, Schneider U, Schleussner E, Witte OW, Hoyer D. Assessment of fetal maturation age by heart rate variability measures using random forest methodology. Comput Biol Med. 2016;70(1):157–62.
 46.
Vressler R, Kreisberg RB, Bernard B, Niederhuber JE, Vockley JG, Shmulevich I, Knijnenburg TA. CloudForest: a scalable and efficient random forest implementation for biological data. PLoS ONE. 2015;10(12):e0144820.
 47.
Webb AR, Copsey KD. Statistical pattern recognition. Chichester: Wiley; 2011. p. 1–642.
 48.
Ocak H, Ertunc HM. Prediction of fetal state from the cardiotocogram recordings using adaptive neurofuzzy inference systems. Neural Comput Appl. 2013;23(6):1583–9.
 49.
Lopez V, Fernandez A, Garcia S, Palade V, Herrera F. An insight into classification with imbalanced data: empirical results and current trends on using data intrinsic characteristics. Inf Sci (Ny). 2013;250(20):113–41.
 50.
Tong L, Change Y, Lin S. Determining the optimal resampling strategy for a classification model with imbalanced data using design of experiments and response surface methodologies. Expert Syst Appl. 2011;38(4):4222–7.
 51.
Spilka J, Georgoulas G, Karvelis P, Vangelis P, Oikonomou P, Chudacek V, Stylios C, Lhotska L, Janku P. Automatic evaluation of FHR recordings from CTUUHB CTG database. Inf Technol Bio Med Inform. 2013;8060:47–61.
 52.
Blagus R, Lusa L. SMOTE for highdimensional classimbalanced data. BMC Bioinform. 2013;14(106):1–16.
 53.
Spilka J. Complex approach to fetal heart rate analysis: a hierarchical classification model. PhD thesis, Czech Technical University in Prague; 2013.
 54.
Georgieva A, Payne SJ, Moulden M, Redman CWG. Artificial neural networks applied to fetal monitoring in labour. Neural Comput Appl. 2013;22(1):85–93.
 55.
Czabanski R, Jezewski J, Matonia A, Jezewski M. Computerized analysis of fetal heart rate signals as the predictor of neonatal acidemia. Expert Syst Appl. 2012;39(15):11846–60.
 56.
Chudacek V, Spilka J, Huptych M, Georgoulas G, Janku P, Koucky M, Stylios C, Lhotska L. Comparison of linear and nonlinear features for intrapartum cardiotocography evaluation–clinical usability vs. contribution to classification. In: Biosginal; 2010. p. 369–372.
 57.
Taft LM, Evans RS, Shyu CR, Egger MJ, Chawla N, Mitchell JA, Thornton SN, Bray B, Varner M. Countering imbalanced datasets to improve adverse drug event predictive models in labor and delivery. J Biomed Inform. 2009;42(2):356–64.
 58.
Sun T, Zhang R, Wang J, Li X, Guo X. Computeraided diagnosis for earlystage lung cancer based on longitudinal and balanced data. PLoS ONE. 2013;8(5):e63559.
 59.
Lin W, Chen JJ. Classimbalanced classifiers for highdimensional data. Brief Bioinform. 2013;14(1):13–26.
 60.
Sun T, Zhang R, Wang J, Li X, Guo X. Computeraided diagnosis for earlystage lung cancer based on longitudal and balanced data. PLoS ONE. 2013;8(5):e63559.
 61.
Nahar J, Imam T, Tickle KS, Ali AS, Chen YP. Computational intelligence for microarray data and biomedical image analysis for the early diagnosis of breast cancer. Expert Syst Appl. 2012;39(16):12371–7.
 62.
Blagus R, Lusa L. SMOTE for highdimensional classimbalanced data. BMC Bioinform. 2013;14(1):106.
 63.
Wang Y, Simon M, Bonde P, Harris BU, Teuteberg JJ, Kormos RL, Antaki JF. Prognosis of right ventricular failure in patients with left ventricular assist device based on decision tree with SMOTE. Trans Inf Technol Biomed. 2012;16(3):383–90.
 64.
Fergus P, Cheung P, Hussain A, AlJumeily D, Dobbins C, Iram S. Prediction of preterm deliveries from EHG signals using machine learning. PLoS ONE. 2013;8(10):e77154.
Authors’ contributions
PF did all of the data processing and analysis, AH conducted the background research, DA managed and conducted editorial reviews, DH and NB evaluated the mathematical formulas and overall methodology. PF wrote the discussion and conclusion sections. All authors read and approved the final manuscript.
Acknowledgements
None.
Competing interests
The authors declare that they have no competing interests.
Availability of data and materials
The datasets generated during the current study are available in Physionet, https://physionet.org/physiobank/database/ctuuhbctgdb/.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Author information
Affiliations
Corresponding author
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Fergus, P., Hussain, A., AlJumeily, D. et al. Classification of caesarean section and normal vaginal deliveries using foetal heart rate signals and advanced machine learning algorithms. BioMed Eng OnLine 16, 89 (2017). https://doi.org/10.1186/s129380170378z
Received:
Accepted:
Published:
Keywords
 Classification
 Feature extraction and selection
 Deep learning
 Intrapartum cardiotocography
 Machine learning
 Random forest