Next Article in Journal
3D Printing Polymeric Materials for Robots with Embedded Systems
Next Article in Special Issue
Visual Robotic Perception System with Incremental Learning for Child–Robot Interaction Scenarios
Previous Article in Journal
Lifetime of Catalyst under Voltage Cycling in Polymer Electrolyte Fuel Cell Due to Platinum Oxidation and Dissolution
Previous Article in Special Issue
Image-Label Recovery on Fashion Data Using Image Similarity from Triple Siamese Network
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multiclass Confusion Matrix Reduction Method and Its Application on Net Promoter Score Classification Problem

by
Ioannis Markoulidakis
,
Ioannis Rallis
,
Ioannis Georgoulas
,
George Kopsiaftis
*,
Anastasios Doulamis
and
Nikolaos Doulamis
School of Rural and Surveying Engineering, National Technical University of Athens, 9th, Heroon Polytechniou Str., 15773 Athens, Greece
*
Author to whom correspondence should be addressed.
Technologies 2021, 9(4), 81; https://doi.org/10.3390/technologies9040081
Submission received: 29 August 2021 / Revised: 13 October 2021 / Accepted: 21 October 2021 / Published: 2 November 2021
(This article belongs to the Collection Selected Papers from the PETRA Conference Series)

Abstract

:
The current paper presents a novel method for reducing a multiclass confusion matrix into a 2 × 2 version enabling the exploitation of the relevant performance metrics and methods such as the receiver operating characteristic and area under the curve for the assessment of different classification algorithms. The reduction method is based on class grouping and leads to a special type of matrix called the reduced confusion matrix. The developed method is then exploited for the assessment of state of the art machine learning algorithms applied on the net promoter score classification problem in the field of customer experience analytics indicating the value of the proposed method in real world classification problems.

1. Introduction

1.1. Data Classification Problem Analysis

The problem of classification in the field of machine learning is an instance of pattern recognition [1]. The objective of the classification problem is to identify the class to which an observation belongs based on the set of associated features (explanatory variables). Typically, a set of training dataset of observations with known class membership is available, and in this context, classification is a special case of supervised learning. A classification problem can be further categorized into (a) binary classification problem in case that the class label can take only two values (e.g., A or B), (b) multiclass classification problem in case the class label can take more than two different values (e.g., A, B, C), and (d) multi-label classification problem in case that each observation is associated with multiple classes.
A wide variety of algorithms has been proposed for the solution of classification problems [1]. Example algorithms for binary classification include: linear regression, logistic regression, naïve Bayes, support vector machines, decision tree, random forest, k-nearest neighbors, convolutional neural networks, artificial neural networks, etc. It appears that depending on the specific problem type a different algorithm may be more appropriate. For that reason the performance of different algorithms should be evaluated before selecting the one that best fits the problem characteristics.
At the core of the performance evaluation of different classification algorithms we find the so called confusion matrix. The confusion matrix is defined as the matrix providing the mix of predicted vs. actual class instances. It allows for the definition of a wide range of performance metrics (e.g., accuracy, precision, recall, Mathews correlation coefficient, etc.) [2] or techniques such as the receiver operating characteristic (ROC) and area under the curve (AUC) [3]. In binary classification all performance metrics and ROC analysis are applicable. However, in multiclass classification only a limited set of performance metrics is available (typically accuracy, recall, precision, and F1-score). As mentioned in [3] ROC analysis is also applicable in multiclass classification problems, however, the complexity of analysis increases with the number of classes. Assuming a problem with N classes, a method for applying ROC analysis in this case is to derive a ROC chart for each class (considering the selected class as positive and the union of the rest classes as negative prediction) leading to N ROC charts. An analysis for a 3-dimensional classification problem is provided in [4].
The current paper provides a method which allows the application of the performance analysis available for a binary classification problem to a multiclass problem. The proposed method allows for the reduction in the dimension of a multiclass confusion matrix eventually reaching an equivalent form of a binary classification confusion matrix. The method is based on the ability to group the available class labels into sub-groups leading to a reduced confusion matrix which, as shown, has a specific structure depending on the way that classes are grouped. The paper then provides the formulation of the performance metrics for the reduced version of the confusion matrix and the way that ROC analysis can be performed. As a final step, the proposed method is applied to a specific multiclass classification problem associated with a customer experience metric (Net promoter score) indicating the method ability to support the algorithm performance evaluation in real world problems.

1.2. Customer Experience and Associated Metrics

Customer experience (CX) has been introduced as a concept a few decades ago [5,6,7,8,9] and since then it has evolved as a pillar of corporate strategy primarily aiming at customer loyalty and customer value. In this context, CX strategies are applied across industries and market sectors ranging from hotels, banks, retail stores, telecoms, and utility companies to health care, business to business, and e-commerce. Over time a set of best practices have been developed and shared across companies on how to build and implement CX strategies. Such practices include the introduction of CX programs with dedicated organizational roles, the inclusion of CX targets in the targets of all company units, the adoption of specific action plans with projects aiming at CX improvement, differentiation, and innovation.
The CX measurement framework is a key element of any CX program. Such a framework refers to a set of methods and metrics for measuring or characterizing the customer perceptions and the satisfaction level from the customer relation with a company. A number of CX metrics have been adopted for collecting customer feedback in an effective manner. Such metrics include net promoter score (NPS) [9], the customer satisfaction score (CSAT) [10] or the customer effort score (CES) [11]. These metrics are characterized by the question the customer is requested to respond to. In particular, the NPS question has the following form: “How likely is that you would recommend this company to your friends or colleagues?” The customer response in this case is an integer in a range from 0 to 10. Depending on the range of customer response customers are categorized as shown in Table 1a.
CSAT question has the following form: “How would you rate your overall satisfaction with this company?”. Customers may be requested to provide an integer score ranging from 1 to 5, 1 to 7, or 1 to 10 or select in multiple choice response ranging from “very dissatisfied” to “very satisfied”, as shown in Table 1b. CES question is in the form: “How easy was for you to complete this task or action?” with customers providing their response in the form of an integer ranging, e.g., from 1 to 10.
The measurement of CX metrics is performed based on two main methods:
  • Market surveys: The survey is performed typically in a random sample of the market population. This method has the advantage of measuring CX for all market competitors. Moreover, apart from the main CX metric the survey also measures the customer satisfaction for a number CX attributes, such as the product experience, the value perception, the touchpoint experience (call center, website, mobile app, shops, etc.), or key customer journeys (e.g., billing, product purchase);
  • Customer feedback during or right after a transaction: this method is used so as to measure the customer satisfaction at different stages of a transaction (or more generally a customer journey) or to measure the CX of a specific touchpoint (e.g., shop, call center, website, mobile app, etc.) Such measurements capture only the feedback of own customers, however, they can reveal significant customer insights.
Both methods are necessary for a successful CX program as they provide essential information on the direction that a company should follow in order to improve its competitive position in the market. The current paper focuses on the analysis applied on the market survey data. In particular, this analysis includes as a first step the monitoring of the CX metric and the CX attribute satisfaction trends vs. time and in comparison with the competing companies. The next step is the so called key drivers’ analysis which provides the link between CX attributes and the CX metric as described in the following section.

1.3. The CX Metric Classification Problem

The CX metric drivers’ analysis problem has the objective to identify the impact of the improvement of a CX attribute (e.g., website experience) to the overall CX Metric (e.g., NPS). This analysis allows companies to issue targeted actions of CX improvement with maximized impact on the overall CX metric. This problem is a multiclass classification problem as the CX metrics and the CX attribute satisfaction scores have integer scores within a specific range so they can be considered to be “labels”. The classification problem has the following formulation:
C X ^ M ( i ) = f ( C X a t t ( i , 1 ) , C X a t t ( i , 2 ) , , C X a t t ( i , K ) )
where, C X ^ M ( i ) is the CX metric score of survey responder i, ( i = 1 , 2 , , n ) , C X a t t ( i , j ) , (j = 1, 2,…, K) is the satisfaction score of the same survey responder i, for CX attribute j and f is the function that allows for the association between these metrics.
The problem of CX metric classification can be addressed based on machine learning algorithms as shown in [12]. The confusion matrix of this problem is a multiclass confusion matrix with a dimension equal to the number of different class labels (i.e., the CX metric score labels). For example for NPS classification an 11 × 11 confusion matrix will occur (class labels: 0, 1, 2,…,10) and for CSAT classification a 5 × 5 , 7 × 7 , or 10 × 10 confusion matrix will occur depending on the score range applied in the relevant questionnaire (1–5, 1–7, 1–10, accordingly). The following section provides an overview of the state of the art regarding the way that the performance of the various algorithms applied for solving a multiclass classification problem.

2. Classification Algorithm Performance Analysis

As mentioned in the introduction, there is a set of performance metrics which can be applied on the confusion matrix of a classification problem so as to assess an algorithm or compare the performance of different algorithms. The confusion matrix for a binary and a multiclass classification problem are presented in Figure 1a,b accordingly. Each column of the matrix represents the instances of a predicted class, while each row represents the instances of an actual class. An element of the confusion matrix at row i and column j provides the number of instances for which the predicted class is j and the actual class is i. In this context, the confusion matrix shows the ways in which the classification model is “confused” when it makes predictions. It can provide the insight not only into the errors being made by a classifier, but more importantly the type of the resulting errors.

2.1. Algorithm Performance for Binary Classification Problems

In the case of a binary classification problem, the confusion matrix has a dimension of 2 × 2 (Figure 1a) where one of the labels is considered to be “Positive” and the other one “Negative”. The matrix elements are characterized based on the predicted label (positive, negative) and the result of the comparison of the predicted with the actual class label (true, false) [13]: true positives (TP), true negatives (TN), false positives (FP), and false negatives (FN). A wide number of performance metrics have been defined for the binary classification confusion matrix as shown in Table 2 [14,15].
In principle, not a single metric can judge the performance of an algorithm and for that reason multiple metrics should be taken into account. Moreover, when data instances are imbalanced (e.g., the number of negative instances is substantially lower vs. the positive instances) some metrics may provide a misleading picture.
An advanced method for assessing the performance of a binary classification algorithm is the so called receiver operating characteristic (ROC) [3]. This is a methodology which provides a chart of true positive rate or recall (TPR) vs. false positive rate or fall out rate (FPR). The resulting chart corresponds to the metrics occurring when a specific threshold θ ranging from 0 to 1 is applied as a criterion for selecting the positive or the negative Prediction based on the probabilities assigned to the positive and negative states as an outcome of the algorithm. The ROC chart allows to compare the performance of an algorithm with the random selection algorithm (an algorithm that selects the positive or the negative class randomly), as well as with the performance of another algorithm. In the ROC space (i.e., pairs of TPR, FPR) the random selection algorithm is characterized by the straight line connecting the points (0,0) and (1,1). Any algorithm performing better than the random selection algorithm should have its ROC curve points lying above this straight line and vice versa.
The area under the curve (AUC) is a metric defined in a ROC chart as the area under the ROC curve of an algorithm. Hence, the higher the AUC value the better the algorithm performance. The AUC (combined with the observation of the ROC curve) allows for the comparison of different algorithms by selecting the algorithms with the highest AUC metric as the best performing ones.

2.2. Multiclass Confusion Matrix and Metrics

In the case of multiclass classification, the metrics defined for binary classification, do not apply in full extend. The multiclass confusion matrix (see Figure 1b) has a dimension N × N where N is the number of different class labels C 0 , C 1 , , C N (e.g., 11 for NPS). Therefore the characterization of TP, TN, FP, FN instances is not applicable in this case. Instead, it is feasible to perform an analysis focusing on a specific class based on the characterization provided in the example of Figure 1b. A set of metrics for each class can be defined based on this approach. Then, based on the proper combination of these metrics it is feasible to provide metrics for the entire confusion matrix.
Table 3 provides an overview of the metrics defined for a multiclass confusion matrix and, in particular, accuracy, recall, precision, and F 1 -score. As it can be seen there are two main approaches in defining performance metrics called “micro” and “macro”.

3. Multiclass Confusion Matrix Reduction Methods

Let us consider a multiclass classification problem with set A containing the N different class labels C i ( i = 1 , 2 , , N ) of the dependent parameter of the problem: A = { C 1 , C 2 , , C N } . The confusion matrix for that problem is an N × N matrix with the form of Figure 1b. The concept of the reduction process of such a multiclass confusion matrix originates from the fact that often from the problem analysis viewpoint, classes could be grouped into sets of classes. For example, in the problem of NPS or CSAT classification, customers with a certain score range can be assigned to specific categories (e.g., NPS between 0 and 6 are categorized as “detractors”, as shown in Table 1a).
Since the class labels can be re-ordered without affecting the resulting performance metrics of a multiclass confusion matrix, without loss of generality we may consider a group of classes as a sub-set of the set of classes defined as follows:
G = { C 1 , C 2 , , C m } , G A
An example of such a grouping of classes is presented in Figure 2. As shown in this figure there is a clear definition of false positives (FP) and false negatives (FN) instances. False positives correspond to instances for which the predicted class belongs to Grouped class G and the actual class does not belong to group G. Vice versa, false negatives correspond to instances for which the predicted class does not belong to group G while the actual class belongs to group G. Regarding true positives (TP) there are some options analyzed in the following section. The total FP and FN instances of the group G are provided as follows:
F P ( G ) = C i G C j G C i , j , F N ( G ) = C i G C j G C i , j

3.1. Class Grouping Options

Regarding the definition of TP for group G the following three options were identified and presented in Figure 3:
  • Relaxed Grouping of Classes [ G ] R : As shown in Figure 3a, in this case any prediction of a class C i G with actual class C j G is considered to be a true positive instance. In the example of NPS classification, assuming that we are interested in the group of detractors (score from 0 to 6) then a prediction of score 1 with an actual score 3 is considered to be true positive since both the predicted and the actual group is “detractor”;
  • Strict Grouping of Classes [ G ] S : As shown in Figure 3b, in this case only the predictions C i G which are identical to the actual class are considered to be true positive instances, i.e., only the instances C i , i with C i G . For example, assuming the grouped class of detractors in NPS problem, if the predicted class is 3 then a TP instance occurs only if the actual class is 3;
  • Hybrid-RS Grouping of Classes [ G ] R S : As shown in Figure 3c, in this case apart from the instances where the predicted class is identical to the actual class, there is an additional set of combinations of predicted and actual classes which are considered to be a TP instance. For example, assuming the group of detractors in NPS problem, assume that we are interested in an algorithm that predicts the scores that are equal of better than the actual scores.
The selection of the proper class grouping option is subject to the nature and the requirements of the specific classification problem under study.

3.2. The Grouped Class Formal Definition

Based on the above class grouping options we define the “Grouped Class” as a pair of two sets as follows:
[ G ] H = ( G , G T P ( G , H ) )
where, G is the subset of A containing the classes that are grouped as defined in Equation (2). H indicates the applied grouping option with the following values: (a) R (for Relaxed), (b) S (for Strict) and (c) R S (for Hybrid-RS). G T P is defined as the set of True Positive combinations for this specific grouping as follows:
G T P ( G , H ) = { ( C i , C j ) : C i , C j G and G is a TP instance based on H option }
In this case, the number of TP instances for any grouping scheme is the following:
T P ( [ G ] H ) = ( C i , C j ) G T P ( G , H ) C i , j

3.3. The Intragroup Mismatch Instances

As also observed from Figure 3, in the case of strict and hybrid-RS grouping there is a set of instances which cannot be characterized as TP. These instances are called intragroup mismatch (IM) instances because both the predicted and the actual class belong to the same group G (thus, intragroup) however, the predicted class is different than the actual one and do not correspond to a TP instance (i.e., a mismatch). According to the definition of G T P set they can be calculated as follows ( H = R , S , R S ) :
I M ( [ G ] H ) = C i , C j G , ( C i , C j ) G T P ( G , H ) C i , j
From Equations (6) and (7) we achieve:
T P ( [ G ] H ) + I M ( [ G ] H ) = C i , C j G C i , j
As a consequence of these definitions, the total number of predicted instances that belongs to G and the total number of actual instances that belong to G will be ( H : R , S , R S ) :
n p r e d i c t e d ( [ G ] H ) = T P ( [ G ] H ) + F P ( [ G ] H ) + I M ( [ G ] H )
n a c t u a l ( [ G ] H ) = T P ( [ G ] H ) + F N ( [ G ] H ) + I M ( [ G ] H )
The above equations indicate clearly that IM of group G cannot be considered to be either FP or FN. The way that IM instances are treated is covered in the following section.

4. The Concept of the Reduced Confusion Matrix

Based on the definitions of the previous section, it is now possible to define a grouping of classes of the original multiclass confusion matrix as a number of M grouped classes. To achieve this, the original set of classes A is split into a number of M disjoint subsets (groups) with the following properties:
G j A , j = 1 , 2 , , M , for i j , G i G j = and i = 1 M G i = A
Then, for each subset of classes Gj a grouping option Hj is being selected ( R , S , R S ) leading to a set of grouped classes GA defined as follows:
G A = { [ G 1 ] H 1 , [ G 2 ] H 2 . , [ G M ] H M , H j = R , S o r R S }
Based on the above definition of the set of grouped classes we can now define the reduced confusion matrix shown in Figure 4. The representation of the reduced version of the confusion matrix should also take into account the new category of intragroup mismatch (IM) instances. To achieve this, we consider a specific format for the reduced confusion matrix called “ M × M + I M ” (M indicates the new reduced dimension of the confusion matrix and I M indicates the presence of intragroup mismatch instances), as presented in the example of Figure 4. In this new proposed format the I M category is introduced as an extra row and column with the column being the transposed version of the I M row. This matrix format leads to Equations (9) and (10) when summing up the instances in column or the row of group Gj so as to calculate the number of actual or the number of predicted instances accordingly.

4.1. Consecutive Confusion Matrix Reduction Steps

Based on the rationale developed in the previous section, it is now possible to define a multi-step confusion matrix reduction process where in the general case either classes or the original matrix or groups of classes can further be grouped according to the problem specific requirements. Based on Equations (11) and (12) the set of grouped classes of step k can be defined as follows:
G A ( k ) = [ G k , 1 ] H k , 1 , [ G k , 2 ] H k , 2 , , [ G k , M k ] H k , M k
G k , i G A ( k 1 ) , i = 1 , 2 , , M k
For i j , G k , i G k , j = and i = 1 M k G k , i = G A ( k 1 )
G A ( 0 ) = A , G 0 , i = C i , H 0 , i = R , i = 1 , 2 , , M
where, G A ( k ) is set of grouped classes at step k ( k 0 ) , M k is the number of grouped classes ( M k 2 ), G k , i , i = 1 , 2 , , M k is the i th grouped class of step k, H k , i is the applied grouping option (R, S or RS) for the i th grouped class of step k ( G k , i ) . Equation (16) provides the initial conditions of the process, i.e., the link to the original set of problem classes A.
In each step of this process, the confusion matrix of the previous step is considered as the basis. According to this approach, the number of FP, FN, TP, and IM instances for a specific grouped class [ G k , M k ] H k , M k can be derived from the generalization of the Equations (3), (7) and (8) as follows:
F P [ G k , M k ] H k , M k = G k 1 , i G k , m G k 1 , j G k , m G k 1 , i , j
F N [ G k , M k ] H k , M k = G k 1 , i G k , m G k 1 , j G k , m G k 1 , i , j
T P [ G k , M k ] H k , M k = ( G k 1 , i , G k 1 , j ) G T P ( G k , m , H k , m ) G k 1 , i , j
I M [ G k , M k ] H k , M k = G k 1 , i G k , m I M ( G k 1 , i ) + G k 1 , i , G k 1 , j G k , m , G k 1 , i , j G T P ( G k , m , H k , m ) G k 1 , i , j
where, G k 1 , i , j is the number of instances where the actual class belongs to grouped class G k 1 , i and the predicted class belongs to grouped class G k 1 , j .
It is easy to show that when the applied grouping option in a consecutive reduction process is either relaxed (R) or strict (S) then the following property applies:
j G j H H = c i [ j G j ] { C i } H , for H either R or S
In the case that different grouping methods are applied in each step then the above property does not hold.

4.2. Performance Metrics for a Reduced Confusion Matrix

The performance metrics for an “ M × M + I M ” reduced confusion matrix ( M > 2 ) are defined according to the formulas of Table 4. Based on the calculation of recall and precision metrics for a specific group of classes, as shown in Table 4, it is feasible to estimate the reduced confusion matrix metrics based on macro average and micro average according to the formulas of Table 3.
Moreover, it is very easy to show that for the same set of classes G, the following properties apply regarding the applied grouping options:
T P ( [ G ] S ) T P ( [ G ] H ) , H = R , S , R S
A c c ( [ G ] S ) A c c ( [ G ] H ) , H = R , S , R S
It is also easy to observe that when strict grouping is applied the accuracy of the reduced confusion matrix equals to the accuracy of original confusion matrix.
In the case that the confusion matrix reduction process leads to a “ 2 × 2 + I M ” confusion matrix, then the confusion matrix has the form of Figure 5.
In this case the typical binary classification confusion matrix metrics can be also calculated based on the following observation regarding the number of positive and negative instances:
Actual Positives :   P a c t u a l = T P + F N + I M P
Predicted Positives :   P p r e d i c t e d = T P + F P + I M P
Actual Negatives :   N a c t u a l = T N + F P + I M N
Actual Negatives :   N p r e d i c t e d = T N + F N + I M P
The performance metrics presented in Table 2 for a 2 × 2 confusion matrix have been properly adapted for the “ 2 × 2 + I M ” confusion matrix and presented in Table 5. Note, that the formulas of Table 2 still apply for the metrics (using, however, the revised version of the required input parameters as provided in Table 5): F 1 -score, Fowlkes–Mallows index, balanced accuracy, prevalence threshold, informedness, markedness, and threat score (critical success index).
The Mathews correlation coefficient based on the proof provided in Appendix A is provided by the following equation:
M C C = I M P · I M N + ( T P + I M P ) T N F N · F P ( T P + F N + I M P ) ( T P + F P + I M P ) ( F N + T N + I M N ) ( F P + T N + I M N )
Moreover, we define the following new metrics related to intragroup mismatch rates:
Positive IM Rate : P I M R = I M P T P + F N + I M P
Negative IM Rate : N I M R = I M N T N + F P + I M N
Positive Predictive IM Rate : P P I M R = I M P T P + F P + I M P
Negative Predictive IM Rate : N P I M R = I M N T N + F N + I M N
Based on the above definitions (Equation (29) to Equation (32)) we achieve:
T P R + P I M R + F N R = 1
T N R + N I M R + F P R = 1
P P V + P P I M R + F D R = 1
N P V + N P I M R + F O R = 1

4.3. Receiver Operating Characteristic for a Reduced 2 × 2 + I M Confusion Matrix

Since ROC chart and the related metric of AUC is a commonly used method to assess the performance of different classification algorithms, it would be advantageous to define a method to derive ROC chart and AUC metric for a reduced 2 × 2 + I M confusion matrix.
A two-step approach is followed in this case for the identification of the predicted class:
  • Step 1: As in the ordinary binary classification, the prediction of positive vs. negative grouped class is based on a threshold θ ;
  • Step 2: The prediction of a specific class from the set of grouped positive or negative classes is based on maximum likelihood.
Based on this approach the ROC chart can be derived based on precision (TPR) and fall-out rate (FPR) metrics based on the equations provided in Table 5. In particular, as the threshold approaches zero ( θ 0 ) all predictions will be negatives. So we obtain:
lim θ 0 T P = 0 , lim θ 0 F P = 0 , lim θ 0 I M P = 0
lim θ 0 T N = T N M L , lim θ 0 F N = P a c t u a l , lim θ 0 I M N = N a c t u a l T N M L
lim θ 0 T P R = 0 , lim θ 0 F P R = 0
where, N a c t u a l is the number of actually negative instances, P a c t u a l is the number of actual positive instances and T N M L is the number of true negative instances based on maximum likelihood selection of the predicted class as θ approaches 0.
When the threshold approaches one ( θ 1 ) all predictions will be positives. So we obtain:
lim θ 1 T P = T P M L , lim θ 1 F P = N a c t u a l , lim θ 1 I M P = P a c t u a l T P M L
lim θ 1 T N = 0 , lim θ 1 F N = 0 , lim θ 1 I M N = 0
lim θ 1 T P R = T P M L / P a c t u a l , lim θ 1 F P R = 1
where P a c t u a l is the number of actually positive instances and T P M L is the number of true positive instances based on maximum likelihood selection of the predicted class when θ approaches 1.
Based on the above analysis, it appears that the ROC space for a reduced confusion matrix is the following: T P R in [ 0 , T P M L / P a c t u a l ] and F P R in [ 0 , 1 ] . This range is the same for all the applied algorithms. However, the ratio T P M L / P a c t u a l may depend on the applied algorithm. Moreover, it is very easy to show that the random selection process corresponds to the straight line between the points: ( 0 , 0 ) and ( T P M L / P a c t u a l , 1 ) .
The area under the curve (AUC) can be estimated based on the following definition [14]:
A U C = 0 1 T P R ( F P R 1 ( θ ) ) d θ

5. Confusion Matrix Reduction for NPS Classification

In this section, we apply the proposed method to the NPS classification problem for a grouping scenario based on specific business requirements.

5.1. NPS Classification Dataset

The analysis has been performed using real NPS survey data from the Greek mobile communication market. The dataset corresponds to a set of approximately 2500 samples per operator collected in a period of 6 consecutive months based on market surveys. For confidentiality purposes the set of data used for the analysis presented in this paper correspond to random selection of approximately 2500 responses from customers of the three different operators.
Each NPS survey was based on a questionnaire which apart from the NPS question it also included questions on the customer level of satisfaction from a set of 9 customer experience attributes: tariff plan, network voice service, network data service, shops, call center, website, mobile application, billing, and roaming service.

5.2. Machine Learning Algorithms for NPS Classification

The NPS classification problem can then be formulated as shown in Equation (1). In our NPS survey data analysis, the NPS classification problem refers to the consideration of the customer scores on the 9 customer experience attributes as a basis to predict the NPS score of each individual responder. Such an analysis allows for the identification of the level of contribution or the importance of each one of the 9 attributes in the formulation of the NPS score.
In this paper, the following algorithms have been tested for the NPS classification problem [12]: logistics regression (LR), k-nearest neighbor (k-NN), naïve Bayes (NB), support vector machines (SVM), decision trees (DT), random forest (RF), convolutional neural networks (CNN), and artificial neural networks (ANN).
The available NPS survey data of approximately 2500 samples were split into training and validation data with an 80–20% ratio. Further details on the applied algorithms are provided in Appendix B.

5.3. Confusion Matrix for the NPS Classification Problem

The original set of classes in this case corresponds to the 11 different classes of NPS score (0–10):
A ( N P S ) = 0 , 1 , 2 , 3 , 4 , 5 , 6 , 7 , 8 , 9 , 10
The NPS classification problem has an 11 × 11 confusion matrix, as shown in Figure 6a resulting from the application of logistic regression algorithm to the NPS survey data. Based on the definition of the confusion matrix reduction process the following grouping steps can be applied:
Step 1: Relaxed grouping based on the definition of customer categories: detractor, passive and promoter (as shown in Figure 6a):
A r e d u c e d ( N P S ) = [ 0 , 1 , 2 , 3 , 4 , 5 , 6 ] R , [ 7 , 8 ] R , [ 9 , 10 ] R = C d e t r a c t o r s , C p a s s i v e s , C p r o m o t e r s
The above reduced confusion matrices has the form of 3 × 3 + I M , as shown in Figure 6. Note that since the applied grouping is relaxed (R) all I M values are zero.
Step 2: In order to compare the performance of different classification algorithms we are going to consider the following grouping step:
A r e d u c e d ( N P S ) = { [ C d e t r a c t o r s ] S , [ C p a s s i v e s , C p o r m o t e r s ] S }
In the reduced set of classes, a customer may belong to one of the two available grouped classes depending on the NPS score: (a) detractor (considered to be the negative case) and (b) either a promoter or a passive (considered to be the positive case).
The result of this step is a reduced confusion in the form of 2 × 2 + I M (see Figure 7) which allows the extensive performance analysis of different classification algorithms.

5.4. Performance Results

The analysis presented above has been followed for all the machine learning algorithms listed in Section 5.2. Just like in the example of logistic regression presented in Figure 6 and Figure 7, every machine learning algorithm led to the production of an original 11 × 11 confusion matrix which was then reduced according to the steps presented in the previous section into a 3 × 3 + I M and then into a 2 × 2 + I M reduced confusion matrix.
The following tables provide the performance metrics derived for each tested machine learning algorithm. Specifically, Table 6 provides the performance metrics for the original 11 × 11 confusion matrix, Table 7 provides the metrics for the reduced 3 × 3 + I M confusion matrix, and Table 8 provides the results for the reduced 2 × 2 + I M confusion matrix.
The first issue considered is the impact of the confusion matrix reduction steps on the performance metrics. Figure 8 compares the metrics of 11 × 11, 3 × 3 + I M and 2 × 2 + I M confusion matrices of the NPS classification problem based on logistic regression. This figure indicates a substantial improvement of the metrics when relaxed grouping is applied (from 11 × 11 to 3 × 3 + I M reduction). This is expected as in relaxed grouping the cases considered to be true positive increases compared to the original confusion matrix. This can be observed in the example of Figure 6a,b. In the second grouping step (from 3 × 3 + I M to 2 × 2 + I M reduction) the application of strict grouping leads to the same accuracy while differences of smaller scale are observed in the rest metrics (precision, recall, F 1 -score). This is also explained by the fact that in strict grouping the number of true positives remains the same so accuracy metric is not affected (this can be also derived from Equations (22) and (23)).
The second issue considered is the ability to compare the performance of machine learning algorithms for the same classification problem (i.e., NPS classification). As it can be seen from Table 6 where the metrics for the 11 × 11 confusion matrix are provided, decision trees appear to outperform the rest of the algorithms for all available metrics (accuracy, precision, recall, and F 1 -score). However, the actual level of these metrics is rather low (e.g., accuracy 0.39). This analysis compares the algorithms in their ability to correctly predict the NPS score (i.e., from 0 to 10). Considering the results of the 3 × 3 + I M reduced confusion matrix in Table 7 we observe that CNN algorithm outperforms the rest algorithms in accuracy, recall, and F 1 -score metrics while considering precision metric ANN algorithm outperforms CNN. These results indicate that CNN is probably the most appropriate algorithm in predicting the customer NPS class, i.e., detractor, passive, promoter with an accuracy level of 0.63. Regarding the results of the 2 × 2 + I M reduced confusion matrix presented in Table 8, we observe that CNN algorithm outperforms the rest algorithms for the majority of the available metrics. The table indicates the ability to deep dive in the comparison of algorithms with the availability of metrics such as the Mathews correlation coefficient. Moreover, metrics such as PIMP and PPIMP indicate the performance of the algorithm in terms of the resulting portion of cases under Intra-group mismatch category.
To complete the comparison of the performance of the applied algorithms we may apply ROC analysis for the 2 × 2 + I M matrix of each algorithm. Figure 9 provides the ROC charts of the 2 × 2 + I M matrix and Figure 10 provides the relevant AUC values. The ROC analysis leads to the conclusion that the best performing algorithms are CNN, ANN, logistic regression, and SVM which clearly outperform the rest of the algorithms. The figures demonstrate the value of the proposed method in assessing the performance of the applied machine learning algorithms.

6. Conclusions

The paper deals with the performance analysis of multi-class classification algorithms based on the analysis of the associated multi-class confusion matrix. Triggered from the fact that the available performance metrics for a multi-class confusion matrix are limited compared to a binary classification problem, the current paper proposes a method that reduces the dimensions of a multi-class confusion matrix based on class grouping process. The key method outcomes include: (a) the options of class grouping (relaxed, strict, and hybrid RS grouping) which are selected based on the classification problem requirements (b) the definition of the “reduced confusion matrix” with an extra row and a column for the intragroup mismatch (IM) instances (c) the performance metrics formulas for a reduced confusion matrix and (d) the application of ROC, AUC analysis for a reduced 2 × 2 + I M confusion matrix.
The paper has applied the proposed method in the case of NPS classification problem demonstrating the capability of the method to handle an 11 × 11 confusion matrix. Based on the NPS classification problem specific characteristics, a two-step grouping approach was followed allowing for the performance analysis of a wide range of machine learning algorithms. The results of ROC analysis indicate that CNN, ANN, logistic regression, and SVM outperform random forest, naïve Bayes, decision trees, and k-NN algorithms for this specific problem. Based on these results and taking into account that NPS scoring patterns may change dynamically depending on the market conditions the next research steps include the investigation of techniques that can further improve the performance of neural networks such as the ones proposed in [16] (retraining of neural networks) and [17] (deep learning framework applied to dynamically changing conditions).
The proposed method can be applied to any multiclass classification problem since the applied classification algorithms always lead to the generation of multiclass confusion matrix. As shown with the example of NPS classification, to exploit the proposed method, the appropriate grouping of classes should be applied depending on the problem specific characteristics and the scope of the analysis. A topic for future research will be the investigation of potential generic grouping approaches that would allow the analysis of multiclass classification problems for which the grouping of classes may not be an obvious task (e.g., optical character recognition, animal species classification, etc.).

Author Contributions

Conceptualization, I.M., A.D. and N.D.; methodology, I.M., I.R., G.K., A.D. and N.D.; software, I.M., I.G. and I.R.; validation, I.M. and G.K.; formal analysis, I.M., I.R. and G.K.; investigation, I.M., I.G., I.R. and G.K.; resources, I.M., G.K. and N.D.; writing—original draft version, I.M., I.R. and G.K.; writing—revies and editing, I.M., G.K., A.D. and N.D.; visualizations, I.M. and I.G.; supervision, A.D. and N.D.; project administration, G.K., A.D. and N.D., funding acquisition, I.M., A.D. and N.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research has been co-financed by the European Union and Greek national funds through the Operational Program Competitiveness, Entrepreneurship and Innovation, under the call RESEARCH–CREATE–INNOVATE (project code: T1EDK-05063).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Mathews Correlation Coefficient for Reduced Multiclass Confusion Matrix

Based on [2] the Mathews correlation coefficient can be derived as a “binarization” of the Pearson coefficient of determination based on the following equation:
R 2 = n ( S n S ¯ ) ( P n P ¯ ) n ( S n S ¯ ) 2 n ( P n P ¯ ) 2
where, S n is the predicted binary class with value 1 for positive and 0 for negative, P n is the actual binary class with value 1 for Positive and 0 for negative and S ¯ and P ¯ are the mean values of S n and P n values accordingly.
Then, for a reduced matrix 2 × 2 + I M the following equations hold:
Total Sample : n t o t a l = T P + F P + I M P + T N + F N + I M N
Predicted Positives P P = n S n = n t o t a l · S ¯ = T P + F N + I M P
Actual Positives : P = n P n = n t o t a l · P ¯ = T P + F P + I M P
Moreover, it is easy to observe that:
S n 2 = S n , P n 2 = P n and n S n P n = T P + I M P
Then, the numerator of Equation (A1) becomes:
n ( S n S ¯ ) ( P n P ¯ ) = T P + I M P n t o t a l S ¯ P ¯
The components of the denominator of Equation (A1) can be derived as follows:
n ( S n S ¯ ) 2 = n S n 2 2 · S ¯ n S n + n t o t a l · S ¯ 2 = n t o t a l S ¯ ( 1 S ¯ )
n ( P n P ¯ ) 2 = n t o t a l · P ¯ ( 1 P ¯ )
So MMC becomes:
M C C = T P + I M P n t o t a l · S ¯ · P ¯ n t o t a l S ¯ · P ¯ ( 1 S ¯ ) ( 1 P ¯ )
Based on Equations (A2)–(A4) we obtain:
M C C = I M P · I M N + ( T P + I M P ) T N F N · F P ( T P + F N + I M N ) ( T P + F P + I M P ) ( F N + T N + I M N ) ( F P + T N + I M N )
From the above equation we observe that as expected for I M P = I M N = 0 we achieve the M C C for a regular 2 × 2 confusion matrix.

Appendix B. The Applied Machine Learning Algorithms

Appendix B.1. Decision Trees

Decision tree learning is a predictive modeling approach used in statistics, data mining, and machine learning [18]. It uses a decision tree to go from observations about an item to conclusions about the item’s target value (represented in the leaves) [19]. In these tree structures, “leaves” represent class labels (i.e., the labels y Y = f ( C 1 , C 2 , , C M ) and “branches” represent conjunctions of features that lead to those class labels. Decision trees where the target variable can take continuous values, typically real numbers, are called regression trees. Table A1 contains the decision trees parameters that we used in the current paper.
Table A1. Parameters of the decision-trees classifiers.
Table A1. Parameters of the decision-trees classifiers.
ParameterValues
Function measuring quality of splitEntropy
Maximum depth of tree3
Weights associated with classes1

Appendix B.2. k-Nearest Neighbors

The k-nearest neighbors (k-NN) algorithm is a non-parametric method used for a wide range of classification problems [20,21,22]. A majority vote of its neighbors classifies an object, with the object being assigned to the class most common among its k-nearest neighbors; it is, therefore, a type of instance-based learning, where the function is only approximated locally and all computation is deferred until classification. Often a fuzzy variation of the k-NN algorithm is used [23]. Table A2 presents the main parameters of the method applied in this paper.
Table A2. Parameters of the k-NN classifiers.
Table A2. Parameters of the k-NN classifiers.
ParameterValues
Number of neighbors5
Distance metricMinkowski
Weight functionuniform

Appendix B.3. Support Vector Machines

Support vector machines (SVM) are supervised learning models with associated learning algorithms [24,25,26] so SVM requires a training dataset. An SVM model is a representation of the examples as points in space, mapped so that the examples of the separate categories are divided by a clear margin that is as wide as possible. The parameters of the SVM method are presented in Table A3.
Table A3. Parameters of the support vector machines classifiers.
Table A3. Parameters of the support vector machines classifiers.
ParameterValues
Kernel typeLinear
Degree of polynomial kernel function3
Weights associated with classes1

Appendix B.4. Random Forest

The random forest (RF) classifier consists of a combination of tree classifiers where each of them is generated using a random vector sampled independently from the input vector, and each tree casts a unit vote for the most popular class to classify an input vector [27]. Table A4 contains the parameters of the RF classifier used in the current work.
Table A4. Parameters of the random forest classifiers.
Table A4. Parameters of the random forest classifiers.
ParameterValues
Number of trees100
Measurements of quality of splitGini index

Appendix B.5. Artificial Neural Networks

Artificial neural networks (ANNs) are highly non-linear classifiers with many applications to extensive domains [16]. Their structures try to resemble how human brains work with neurons and synapses. In particular, an ANN consists of one input layer that receives the input signals as data, one or more hidden layers of neurons that process these data under a non-linear way and one output layer that yields the final classification outcome. Each neuron sums the input signals under a weighted (parametrize) way and then passes them to a non-linear activation functions [28]. The activation functions, actually non-linearly transforms the input signals to the outputs and they consist of function bases, if they are increasing, bounded, and almost satisfy the continuity property (through the Kolmogorov theorem [29]). Especially, their significance in modeling time-varying signals is of particular interest [30]. In our implementation, an ANN with one input, one hidden, and one output layer is employed. Figure 5 indicates the specific parameters of our ANN. The ReLu function is adopted as activation for the input and hidden layers, while the softmax for the output layer. This is due to the fact that these activation functions seem to work better than other approaches as is proven through the deep learning paradigm [31]. The ANN parameters are presented in Table A5.
Table A5. Parameters of the artificial neural networks.
Table A5. Parameters of the artificial neural networks.
ParameterValues
Number of hidden neurons6
Activation function applied for the input and hidden layerReIU
Activation function applied for the output layerSoftmax
Optimizer network functionAdam
Calculated lossSparse categorical cross-entropy
Epochs used100
Batch size10

Appendix B.6. Convolutional Neural Networks

Convolutional neural networks (CNNs) exploit machine learning paradigms on deep structures. It firstly extracts a set of appropriate features from the raw data, by applying convolutions on the input signals propagating them into deep layers while at the last layer a classification is carried out to assign the input data into classes but on the use of the deep features identified by the convolutional layers. CNNs utilize trainable filters and pooling operations on their input resulting in a hierarchy of increasingly complex features [31,32,33]. Convolutional layers consist of a rectangular grid of neurons (filters), each of which takes inputs from rectangular sections of the previous layer. Each convolution layer is followed by a pooling layer in which subsamples block-wise the output of the precedent convolutional layer and produce a scalar output for each block. Table A6 provides the parameters of the CNN method.
Table A6. Parameters of the convolutional neural networks.
Table A6. Parameters of the convolutional neural networks.
ParameterValues
ModelSequential (array of Keras Layers)
Kernel size3
Pool size4
Activation function appliedReIU
Calculated losscategorical cross-entropy
Epochs used100
Batch size128

Appendix B.7. Naïve Bayes

Naïve Bayes classifiers are a family of probabilistic classifiers based on applying Bayes’ theorem with strong independence assumptions between the features. These classifiers are highly scalable, requiring a number of parameters linear in the number of variables (features/predictors) in a learning problem [34]. Maximum-likelihood training can be completed by evaluating a closed-form expression, which takes linear time, rather than by expensive iterative approximation as used for many other types of classifiers.

Appendix B.8. Logistic Regression

Logistic regression is a statistical model that in its basic form uses a logistic function to model a binary dependent variable [35]. Mathematically, a binary logistic model has a dependent variable with two possible values, such as pass/fail, which is represented by an indicator variable, where the two values are labeled 0 and 1. The defining characteristic of the logistic model is that increasing one of the independent variables multiplicatively scales the odds of the given outcome at a constant rate, with each independent variable having its own parameter. Table A7 provides the main logistic regression parameters.
Table A7. Parameters of the logistic regression.
Table A7. Parameters of the logistic regression.
ParameterValues
Maximum number of iterations300
Algorithm used in optimizationL-BFGS
Weights associated with classes1

References

  1. Alpaydin, E. Introduction to Machine Learning; The MIT Press: Cambridge, MA, USA, 2020. [Google Scholar]
  2. Matthews, B.W. Comparison of the Predicted and Observed Secondary Structure of T4 Phage Lysozyme. Biochim. Biophys. Acta-(Bba)-Protein Struct. 1975, 405, 442–451. [Google Scholar] [CrossRef]
  3. Fawcett, T. An Introduction to ROC Analysis. Pattern Recognit. Lett. 2006, 27, 861–874. [Google Scholar] [CrossRef]
  4. Mossman, D. Three-Way ROCs. Med Decis. Mak. 1999, 19, 78–89. [Google Scholar] [CrossRef]
  5. Bolton, R.N.; Drew, J.H. A Multistage Model of Customers’ Assessments of Service Quality and Value. J. Consum. Res. 1991, 17, 375–384. [Google Scholar] [CrossRef]
  6. Aksoy, L.; Buoye, A.; Aksoy, P.; Larivière, B.; Keiningham, T.L. A Cross-National Investigation of the Satisfaction and Loyalty Linkage for Mobile Telecommunications Services across Eight Countries. J. Interact. Mark. 2013, 27, 74–82. [Google Scholar] [CrossRef]
  7. Rageh Ismail, A.; Melewar, T.C.; Lim, L.; Woodside, A. Customer Experiences with Brands: Literature Review and Research Directions. Mark. Rev. 2011, 11, 205–225. [Google Scholar] [CrossRef] [Green Version]
  8. Gentile, C.; Spiller, N.; Noci, G. How to Sustain the Customer Experience:: An Overview of Experience Components That Co-Create Value With the Customer. Eur. Manag. J. 2007, 25, 395–410. [Google Scholar]
  9. Reichheld, F.F.; Covey, S.R. The Ultimate Question: Driving Good Profits and True Growth; Harvard Business School Press: Boston, MA, USA, 2006; Volume 211. [Google Scholar]
  10. Fornell, C.; Johnson, M.D.; Anderson, E.W.; Cha, J.; Bryant, B.E. The American Customer Satisfaction Index: Nature, Purpose, and Findings. J. Mark. 1996, 60, 7–18. [Google Scholar] [CrossRef] [Green Version]
  11. de Haan, E.; Verhoef, P.C.; Wiesel, T. The Predictive Ability of Different Customer Feedback Metrics for Retention. Int. J. Res. Mark. 2015, 32, 195–206. [Google Scholar] [CrossRef]
  12. Markoulidakis, I.; Rallis, I.; Georgoulas, I.; Kopsiaftis, G.; Doulamis, A.; Doulamis, N. A Machine Learning Based Classification Method for Customer Experience Survey Analysis. Technologies 2020, 8, 76. [Google Scholar] [CrossRef]
  13. Stehman, S.V. Selecting and Interpreting Measures of Thematic Classification Accuracy. Remote. Sens. Environ. 1997, 62, 77–89. [Google Scholar] [CrossRef]
  14. Powers, D. Evaluation: From Precision, Recall and F-Measure to ROC, Informedness, Markedness & Correlation. arXiv 2020, arXiv:2010.16061. [Google Scholar]
  15. Tharwat, A. Classification Assessment Methods. Appl. Comput. Inform. 2020. [Google Scholar] [CrossRef]
  16. Doulamis, A.; Doulamis, N.; Kollias, S. On-Line Retrainable Neural Networks: Improving the Performance of Neural Networks in Image Analysis Problems. IEEE Trans. Neural Netw. 2020, 11, 137–155. [Google Scholar] [CrossRef]
  17. Doulamis, N.; Voulodimos, A. FAST-MDL: Fast Adaptive Supervised Training of Multi-Layered Deep Learning Models for Consistent Object Tracking and Classification. In Proceedings of the 2016 IEEE International Conference on Imaging Systems and Techniques (IST), Chania, Greece, 4–6 October 2016; pp. 318–323. [Google Scholar] [CrossRef]
  18. Quinlan, J.R. Induction of Decision Trees. Mach. Learn. 1986, 1, 81–106. [Google Scholar] [CrossRef] [Green Version]
  19. Rokach, L.; Maimon, O.Z. Data Mining with Decision Trees: Theory and Applications; World Scientific: Singapore, 2007. [Google Scholar]
  20. Bhatia, N. Survey of nearest neighbor techniques. arXiv 2010, arXiv:1007.0085. [Google Scholar]
  21. Murphy, K.P. Machine Learning: A Probabilistic Perspective; MIT Press: Cambridge, MA, USA, 2012. [Google Scholar]
  22. Protopapadakis, E.; Voulodimos, A.; Doulamis, A.; Camarinopoulos, S.; Doulamis, N.; Miaoulis, G. Dance Pose Identification from Motion Capture Data: A Comparison of Classifiers. Technologies 2018, 6, 31. [Google Scholar] [CrossRef] [Green Version]
  23. Keller, J.M.; Gray, M.R.; Givens, J.A. A Fuzzy K-Nearest Neighbor Algorithm. IEEE Trans. Syst. Man Cybern. 1985, 580–585. [Google Scholar] [CrossRef]
  24. Abe, S. Support Vector Machines for Pattern Classification; Advances in Pattern Recognition; Springer: London, UK, 2010. [Google Scholar] [CrossRef]
  25. Smola, A.J.; Schölkopf, B. A Tutorial on Support Vector Regression. Stat. Comput. 2004, 14, 199–222. [Google Scholar] [CrossRef] [Green Version]
  26. Kopsiaftis, G.; Protopapadakis, E.; Voulodimos, A.; Doulamis, N.; Mantoglou, A. Gaussian Process Regression Tuned by Bayesian Optimization for Seawater Intrusion Prediction. Comput. Intell. Neurosci. 2019, 2019, e2859429. [Google Scholar] [CrossRef]
  27. Pal, M. Random Forest Classifier for Remote Sensing Classification. Int. J. Remote Sens. 2005, 26, 217–222. [Google Scholar] [CrossRef]
  28. Haykin, S.; Network, N. A comprehensive foundation. Neural Netw. 2004, 2, 41. [Google Scholar]
  29. Hecht-Nielsen, R.; Drive, O.; Diego, S. Kolmogorov’s Mapping Neural Network Existence Theorem. In Proceedings of the International Conference on Neural Networks; IEEE Press: New York, NY, USA, 1987; Volume 3, pp. 11–14. [Google Scholar]
  30. Doulamis, N.; Doulamis, A.; Varvarigou, T. Adaptable Neural Networks for Modeling Recursive Non-Linear Systems. In Proceedings of the 2002 14th International Conference on Digital Signal Processing Proceedings. DSP 2002 (Cat. No.02TH8628), Santorini, Greece, 1–3 July 2002; Volume 2, pp. 1191–1194. [Google Scholar] [CrossRef]
  31. Voulodimos, A.; Doulamis, N.; Doulamis, A.; Protopapadakis, E. Deep Learning for Computer Vision: A Brief Review. Comput. Intell. Neurosci. 2018, 2018, 1–13. [Google Scholar] [CrossRef]
  32. Protopapadakis, E.; Voulodimos, A.; Doulamis, A. On the Impact of Labeled Sample Selection in Semisupervised Learning for Complex Visual Recognition Tasks. Complexity 2018, 2018, e6531203. [Google Scholar] [CrossRef]
  33. Doulamis, A.; Doulamis, N.; Protopapadakis, E.; Voulodimos, A. Combined Convolutional Neural Networks and Fuzzy Spectral Clustering for Real Time Crack Detection in Tunnels. In Proceedings of the 2018 25th IEEE International Conference on Image Processing (ICIP), Athens, Greece, 7–10 October 2018; pp. 4153–4157. [Google Scholar] [CrossRef]
  34. Haouari, B.; Ben Amor, N.; Elouedi, Z.; Mellouli, K. Naïve Possibilistic Network Classifiers. Fuzzy Sets Syst. 2009, 160, 3224–3238. [Google Scholar] [CrossRef]
  35. Hosmer, D.W.; Lemeshow, S. Applied Logistic Regression; John Wiley & Sons: New York, NY, USA, 2000. [Google Scholar]
Figure 1. Confusion matrix examples. (a) Binary classification problem confusion matrix. (b) Multiclass classification problem confusion matrix.
Figure 1. Confusion matrix examples. (a) Binary classification problem confusion matrix. (b) Multiclass classification problem confusion matrix.
Technologies 09 00081 g001
Figure 2. Confusion matrix class grouping example.
Figure 2. Confusion matrix class grouping example.
Technologies 09 00081 g002
Figure 3. True positive instances definition examples. (a) TP for relaxed grouping. (b) TP for strict grouping. (c) Example TP for hybrid RS-grouping.
Figure 3. True positive instances definition examples. (a) TP for relaxed grouping. (b) TP for strict grouping. (c) Example TP for hybrid RS-grouping.
Technologies 09 00081 g003
Figure 4. The reduced confusion matrix “ M × M + I M ” format ( H = R , S , R S ).
Figure 4. The reduced confusion matrix “ M × M + I M ” format ( H = R , S , R S ).
Technologies 09 00081 g004
Figure 5. The concept of reduced confusion matrix in “ 2 × 2 + I M ” form.
Figure 5. The concept of reduced confusion matrix in “ 2 × 2 + I M ” form.
Technologies 09 00081 g005
Figure 6. The reduction in the 11 × 11 NPS classification confusion matrix into a 3 × 3 + I M reduced confusion matrix. (a) The 11 × 11 NPS confusion matrix. (b) The reduced 3 × 3 + I M confusion matrix.
Figure 6. The reduction in the 11 × 11 NPS classification confusion matrix into a 3 × 3 + I M reduced confusion matrix. (a) The 11 × 11 NPS confusion matrix. (b) The reduced 3 × 3 + I M confusion matrix.
Technologies 09 00081 g006
Figure 7. The reduction in the 3 × 3 + I M confusion matrix into a 2 × 2 + I M .
Figure 7. The reduction in the 3 × 3 + I M confusion matrix into a 2 × 2 + I M .
Technologies 09 00081 g007
Figure 8. Logistic regression performance for the original and reduced confusion matrices.
Figure 8. Logistic regression performance for the original and reduced confusion matrices.
Technologies 09 00081 g008
Figure 9. The ROC chart for NPS classification problem.
Figure 9. The ROC chart for NPS classification problem.
Technologies 09 00081 g009
Figure 10. The AUC vs. the applied algorithm.
Figure 10. The AUC vs. the applied algorithm.
Technologies 09 00081 g010
Table 1. The two most widely used customer categorizations.
Table 1. The two most widely used customer categorizations.
(a) The NPS customer catergorization
NPS ResponseNPS Label
9–10Promoter
7–8Passive
0–6Detractor
(b) The CSAT customer catergorization
CSAT ResponseCSAT Label
5Very Satisfied
4Satisfied
3Neutral
2Dissatisfied
1Very Dissatisfied
Table 2. The performance metrics of the binary classification confusion matrix.
Table 2. The performance metrics of the binary classification confusion matrix.
MetricFormula
Accuracy A c c = T P + T N T P + T N + F P + F N
True Positive Rate (Recall) T P R = T P T P + F N
True Negative Rate (Specificity) T N R = T N T N + F P
Positive Predictive Value (Precision) P P V = T P T P + F P
Negative Predictive Value N P V = T N T N + F N
F 1 -Score F 1 = 2 · T P R · P P V T P R + P P V
False Negative Rate (Miss Rate) F N R = F N T P + F N
False Positive Rate (Fall Out Rate) F P R = F P T N + F P
False Discovery Rate F D R = F P T P + F P
False Omission Rate F O R = F N T N + F N
Fowlkes-Mallows index F M = P P V · T P R
Balanced Accuracy B A = T P R + T N R 2
Mathews Correlation coefficient M C C = T P · T N F P · F N ( T P + F N ) ( T P + F P ) ( F N + T N ) ( F P + T N )
Prevalence Threshold P T = T P R ( 1 T N R ) + T N R 1 T P R + T N R 1
Informedness B M = T P R F P R
Markedness M K = P P V F O R
Threat Score (Critical Success Index) T S = T P T P + F N + F P
Table 3. Performance metrics for a multiclass confusion matrix.
Table 3. Performance metrics for a multiclass confusion matrix.
MetricFormula
Accuracy A c c ( A r e d u c e d ) = i = 1 N T P ( C i ) i = 1 N j = 1 N C i , j
Recall of Class C i T P R ( C i ) = T P ( C i ) T P ( C i ) + F N ( C i )
Precision of Class C i P P V ( C i ) = T P ( C i ) T P ( C i ) + F P ( C i )
F 1 -Score of Class C i F 1 ( C 1 ) = 2 · T P R ( C i ) · P P V ( C i ) T P R ( C i ) · P P V ( C i )
Recall (macro average) T P R ( m a c r o ) = 1 N i = 1 N T P R ( C i )
Precision (macro average) P P V ( m a c r o ) = 1 N i = 1 N P P V ( C i )
F 1 -Score (macro average) F 1 ( m a c r o ) = 2 · T P R ( m a c r o ) · P P V ( m a c r o ) T P R ( m a c r o ) + P P V ( m a c r o )
Recall (micro average) T P R ( m i c r o ) = i = 1 N T P ( C i ) i = 1 N [ T P ( C i ) + F P ( C i ) ]
Precision (micro average) P P V ( m i c r o ) = i = 1 N R P ( C i ) i = 1 N [ T P ( C i ) + F P ( C i ) ]
F 1 -Score (micro average) F 1 ( m i c r o ) = 2 · T P R ( m i c r o ) · P P V ( m i c r o ) T P R ( m i c r o ) + P P V ( m i c r o )
Table 4. Performance metrics of reduced confusion matrix.
Table 4. Performance metrics of reduced confusion matrix.
MetricFormula
Accuracy of Reduced Confusion Matrix A c c ( A r e d u c e d ) = j T P ( [ G j ] H j ) j T P ( [ G j ] H j ) + i j , i j G i , j + j I M ( [ G j ] H j )
Recall of Group G j T P R ( [ G j ] H j ) = T P ( [ G j ] H j ) T P ( [ G j ] H j ) + F N ( [ G j ] H j ) + I M ( [ G j ] H j )
Precision of Group G j P P V ( [ G j ] H j ) = T P ( [ G j ] H j ) T P ( [ G j ] H j ) + F P ( [ G j ] H j ) + I M ( [ G j ] H j )
Table 5. Performance metrics of reduced confusion matrix.
Table 5. Performance metrics of reduced confusion matrix.
MetricFormula
Accuracy A c c = T P + T N T P + T N + F P + F N + I M P + I M N
True Positive Rate (Recall) T P ( T P + F N + I M P )
True Negative Rate (Specificity) T N T N + F P + I M N
Positive Predictive Value (Precision) P P V = T P T P + F P + I M P
Negative Predictive Value N P V = T N T N + F N + I M N
False Negative Rate (Miss Rate) F N R = F N T P + F N + I M P
False Positive Rate (Fall Out Rate) F P R = F P T N + F P + I M N
False Discovery Rate F D R = F P T P + F P + I M P
False Omission Rate F O R = F N T N + F N + I M N
Table 6. The performance metrics of the 11 × 11 confusion matrix for each applied algorithm.
Table 6. The performance metrics of the 11 × 11 confusion matrix for each applied algorithm.
Logistic Regr.SVMk-NNDecision TreesRandom ForestNaïve BayesCNNANN
Accuracy0.370.380.340.390.330.310.380.37
Precision0.160.170.170.230.190.170.210.14
Recall0.150.170.160.220.170.150.190.17
F1-score0.130.140.160.210.180.150.190.15
Table 7. The performance metrics of the “ 3 × 3 + I M ” reduced confusion matrix for each applied algorithm.
Table 7. The performance metrics of the “ 3 × 3 + I M ” reduced confusion matrix for each applied algorithm.
Logistic Regr.SVMk-NNDecision TreesRandom ForestNaïve BayesCNNANN
Accuracy0.580.560.600.510.540.560.630.57
Precision0.680.650.620.500.550.540.660.69
Recall0.510.480.560.500.500.580.600.52
F1-score0.500.530.580.500.510.570.620.51
Table 8. The performance metrics of the “ 2 × 2 + I M ” reduced confusion matrix for each one of the applied algorithms.
Table 8. The performance metrics of the “ 2 × 2 + I M ” reduced confusion matrix for each one of the applied algorithms.
Logistic Regr.SVMk-NNDecision TreesRandom ForestNaïve BayesCNNANN
Accuracy0.580.560.600.510.540.560.630.57
Precision0.570.550.600.520.540.620.620.56
Recall0.640.610.630.520.580.540.650.62
F1-Score0.600.580.610.520.560.580.630.59
Specificity0.090.070.160.150.110.270.190.10
Miss Rate0.010.010.050.100.050.190.040.02
Negative Predictive Value0.110.090.180.150.130.210.210.13
Fall Out Rate0.140.100.100.100.120.060.090.13
False Discovery Rate0.120.100.100.100.120.070.090.12
False Omission Rate0.030.010.110.180.110.320.080.04
Fowlkes-Mallows index0.600.580.610.520.560.580.630.59
Mathews Correlation Coefficient0.330.370.350.370.360.340.380.36
PIMR0.350.390.320.380.370.270.310.36
PPIMR0.310.350.300.380.350.310.300.32
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Markoulidakis, I.; Rallis, I.; Georgoulas, I.; Kopsiaftis, G.; Doulamis, A.; Doulamis, N. Multiclass Confusion Matrix Reduction Method and Its Application on Net Promoter Score Classification Problem. Technologies 2021, 9, 81. https://doi.org/10.3390/technologies9040081

AMA Style

Markoulidakis I, Rallis I, Georgoulas I, Kopsiaftis G, Doulamis A, Doulamis N. Multiclass Confusion Matrix Reduction Method and Its Application on Net Promoter Score Classification Problem. Technologies. 2021; 9(4):81. https://doi.org/10.3390/technologies9040081

Chicago/Turabian Style

Markoulidakis, Ioannis, Ioannis Rallis, Ioannis Georgoulas, George Kopsiaftis, Anastasios Doulamis, and Nikolaos Doulamis. 2021. "Multiclass Confusion Matrix Reduction Method and Its Application on Net Promoter Score Classification Problem" Technologies 9, no. 4: 81. https://doi.org/10.3390/technologies9040081

APA Style

Markoulidakis, I., Rallis, I., Georgoulas, I., Kopsiaftis, G., Doulamis, A., & Doulamis, N. (2021). Multiclass Confusion Matrix Reduction Method and Its Application on Net Promoter Score Classification Problem. Technologies, 9(4), 81. https://doi.org/10.3390/technologies9040081

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop