RESEARCH ARTICLE Year : 2017  Volume : 8  Issue : 1  Page : 24 A reduced set of features for chronic kidney disease prediction Rajesh Misir^{1}, Malay Mitra^{2}, Ranjit Kumar Samanta^{2}, ^{1} Department of Computer Science, Vidyasagar University, Medinipur, India ^{2} Department of Computer Science and Application, Expert Systems Laboratory, University of North Bengal, Darjeeling, West Bengal, India Correspondence Address: Chronic kidney disease (CKD) is one of the lifethreatening diseases. Early detection and proper management are solicited for augmenting survivability. As per the UCI data set, there are 24 attributes for predicting CKD or nonCKD. At least there are 16 attributes need pathological investigations involving more resources, money, time, and uncertainties. The objective of this work is to explore whether we can predict CKD or nonCKD with reasonable accuracy using less number of features. An intelligent system development approach has been used in this study. We attempted one important feature selection technique to discover reduced features that explain the data set much better. Two intelligent binary classification techniques have been adopted for the validity of the reduced feature set. Performances were evaluated in terms of four important classification evaluation parameters. As suggested from our results, we may more concentrate on those reduced features for identifying CKD and thereby reduces uncertainty, saves time, and reduces costs.
Introduction Disorders affecting kidney structure and function in heterogeneous form are generally termed chronic kidney disease (CKD).[1] National Kidney Foundation of America in 2002 created guidelines for a more clear definition and classification system for CKD.[2] CKD is classified into Stages I–V according to the estimated glomerular filtration rate (GFR) shown in [Table 1].[2] GFR is estimated having mathematical equations using serum creatinine, age, sex, body size, ethnic origin, etc.[1],[3] If the normal functionality of a kidney is degraded to an extent, wastes can build up to high levels in your blood making you feel sick. Unfortunately, this degradation is noted at the later stages of CKD making the matter complicated in much of the cases. Sometimes, it is too late when we consult a physician. Hence, the early detection is solicited for longterm survivability. It is argued that once in a year or once in 2 years, CKDrelated investigations may be done. General awareness of the people has to be increased along with providing less number of pathological tests with less cost and less time.{Table 1} [Table 2] represents the CKD data set from UCI that contains 24 attributes plus one attribute for class (binary).[4] It contains 400 samples to two different classes (“CKD”  250 cases; “NOTCKD”  150 cases). Out of 24 attributes, 11 are numeric and 13 are nominal. The data set contains a number of missing values. After excluding tuples with missing values, 158 samples were used in this work.{Table 2} This paper aims to determine the important features for CKD predictions. Features were selected using one correlationbased algorithm with eight different searching techniques. Two different intelligent classification algorithms based on artificial neural networks were applied to test the effectiveness of the feature set. Classification results were compared in terms of four modelvalidating parameters. Methods Correlationbased feature subset selection Correlationbased feature subset selection (CFS) is based on the following hypothesis: A feature of a subset is considered good which are highly correlated with the class but may be uncorrelated with other features of the class.[5] The feature evaluation mathematical formula provides an operational definition of the above hypothesis as follows: [INLINE:1] Where rfc is the correlation between the summed features and the class variable; k is the number of features rfc and above line is the average of the correlation between the features and class variable; and rff and above line is the average intercorrelation between features.[6] Summarily, a feature will be admitted to the subset if its correlation with the class is higher than the highest correlation between it and anyone of the already selected features. To accommodate all categories of features in equation 1, continuous features are transformed to categorical features using a discretization method.[7] The degree of associations between nominal features was estimated by a specific method.[8] For searching the feature subset in a reasonable amount of time, heuristic search strategies are often used.[9] Further, different types of feature selection methods use the correlationbased approach in different applications.[10],[11] Eight search methods were applied in the study, namely, Best First, Exhaustive Search, Genetic Search, Greedy Stepwise, Linear Forward Selection, Random Search, Scatter Search, and Subset Size Forward Selection.[12] Incremental back propagation learning networks classifier The normal back propagation network is not an incremental by its nature.[13] The network learns by the back propagation rule of Rumelhart et al. under the constraint that the change to each weight for each instance is bounded.[14] It is accomplished by introducing scaling factors which scale down all weight adjustments so that all of them are within bounds. The learning rule is now [INLINE:2] Where Wij is the weight from unit i to unit j, η (0<η<1) is a trialindependent learning rate, δj is the error gradient at unit j, Oi is the activation level at unit i, and the parameter k denotes the kth iteration. The incremental back propagation learning networks proceed as follows: Given a single misclassified instance: Begin Repeatedly apply the bounded weight adaptation learning rule (2) on the instance until stopping criteria are met. If the instance can be correctly learned, then restore the old weights and apply the bounded weight adaptation learning rule once; Else restore the old weights and apply the structural adaptation learning rules. End The stopping criteria are the instance can be correctly learned or the output error fluctuates in a small range.[13] Levenberg–Marquardt classifier The Levenberg–Marquardt (LM) algorithm is basically an iterative method that locates the minimum of a multivariate function that is expressed as the sum of squares of nonlinear realvalued functions.[15],[16] LM can be thought of as a combination of steepest descent and the Gauss–Newton (GN) method. LM algorithm is more robust than GN algorithm which essentially means that it finds a solution even if it starts far off the final minimum. During the iterations, the new configuration of weights in step k + 1 is calculated as follows [INLINE:3] Where J  the Jacobian matrix, λ  adjustable parameter, and μ  error vector. The parameter λ is modified based on the development of error function E. If the step causes a reduction of E, we accept it. Otherwise, λ is changed; reset the original value and recalculate w (k + 1). Applications This study consists of three stages: feature extraction and reduction by CFS using eight different search algorithms finding the most relevant and reduced set of features and then classification by two important classification algorithms. The schematic view of our proposed system is shown in [Figure 1].{Figure 1} Data Preprocessing It is one of the important steps for the development of any model. We completely randomize the data sets with missing records. Missing tuples were excluded leaving 158 data sets for use in this work without missing values. We apply CFS with the aforementioned eight search techniques using a free software named WEKA.[12] Out of the eight search algorithms, six algorithms, a majority, suggested eight common reduced attribute set as shown in [Table 3].{Table 3} Next, we apply the two discussed intelligent techniques on the data set to testify the effectiveness of the above reduced feature set as shown in [Table 3]. Network Architecture The artificial neural network model selection includes a choice of network architecture and feature selection. Theoretically, a network with one hidden layer and logistic function as the activation function at the hidden and output nodes is capable of approximating any function arbitrarily closely provided that the number of hidden nodes is large enough.[17] Hence, we used one input layer, one hidden layer, and one output layer. The number of hidden nodes is evaluated from the following formula as proposed by Huang et al.[18] [INLINE:4] In the present study, m = number of input nodes = 8, n = number of output = 2, and hence s = 6 after round off. Hence, in our study, we use six neurons at the hidden layer for all combinations. Modeling Results The classification algorithms were implemented in Alyuda NeuroIntelligence.[19] We used Intel Core 2 Duo Processor E7400 CPU (2.8 GHz Dual Core, 1066 MHz FSB, 3 MB L2 cache) with 2048 MB DDR2 RAM for implementation. [Table 4] shows the network structure, epochs, number of retrains, and number patterns used in training, validation, and testing phases. As overtraining control measure, we retain the copy of the network with the lowest validation error.{Table 4} Performance evaluation To predict the performance of the system, we computed correct classification accuracy (CCR), specificity, sensitivity, and receiver operating characteristic area under the curve (AUC) as these are very important parameters to predict the performance of the system without knowing the distribution of data. We computed true positive (TP), true negative (TN), false positive (FP), and false negative (FN) to further compute other performance parameters as discussed below: [INLINE:5] Further, we had drawn curves using sensitivity and specificity and computed the AUC. The high value of AUC indicates highperformance matching. The compiled results from 100 simulations are shown in [Table 5]. We observe that both the classification methods show excellent performance justifying the validity of the reduced set of features. However, out of two methods, CFS + LM shows slightly better performance in terms of CCR and specificity.{Table 5} Conclusions Our findings suggest a reduced set of eight features: specific gravity, albumin, serum creatinine, hemoglobin, packed cell volume, white blood cell count, red blood cell count, and hypertension as more significant for investigating CKD. To justify the validity of the reduced feature set, we deployed two different binary classifiers. The classifiers show excellent performance. Hence, we suggest that these reduced feature set might be worthwhile to scrutiny when the final decision is made by the doctors. Reduced parameters reduce laboratory costs and time. At the same time, reduced features reduce uncertainty in decisionmaking. We suggest that the techniques used here could be applied to other diseases. Financial support and sponsorship Nil. Conflicts of interest There are no conflicts of interest. References


