Decision Tree (Gini)

LOADING DATA INTO R ENVIRONMENT

Reading the Data

# reading the data
library(data.table)
df <- fread("MCICreditCardDefault.csv")

Converting some data columns into factors

# convert 'Male' as a factor
df[, Male := as.factor(Male)]
# convert 'Education' as a factor
df[, Education := as.factor(Education)]
# convert 'MaritalStatus' as a factor
df[, MaritalStatus := as.factor(MaritalStatus)]
# convert 'Default' as a factor
df[, Default := as.factor(Default)]

Splitting data into training and test set

library(caret)
# data partition
set.seed(2341)
trainIndex <- createDataPartition(df$Default, p = 0.80, list = FALSE)
# 80% training data
train.df <- df[trainIndex, ]
# 20% testing data
test.df <- df[-trainIndex, ]
dim(train.df)

## [1] 23681     9

dim(test.df)

## [1] 5920    9

TRAINING THE DECISION TREE MODEL

Setting Control parameters

library(caret)
# control parameters
trctrl <- trainControl(method = "cv", n = 10, classProbs = TRUE)

Running the Training Model

## setting levels as "Yes" and "NO"
train.df$Default <- ifelse(train.df$Default == "1","Yes","No")

set.seed(3333)
# fitting decision tree classification model
DTModel <- train(Default ~ CreditLimit
                         + Male
                         + Education
                         + MaritalStatus
                         + Age
                         + BillOutstanding
                         + LastPayment, 
                         data = train.df, 
                         method = "rpart",
                         metric = "ROC",
                         parms  = list(split = "gini"), 
                         trControl = trctrl)
# model summary
DTModel

## CART 
## 
## 23681 samples
##     7 predictor
##     2 classes: 'No', 'Yes' 
## 
## No pre-processing
## Resampling: Cross-Validated (10 fold) 
## Summary of sample sizes: 21313, 21314, 21313, 21313, 21312, 21314, ... 
## Resampling results across tuning parameters:
## 
##   cp           Accuracy   Kappa     
##   0.001703255  0.7791055  0.07998385
##   0.002081756  0.7791899  0.08118771
##   0.003406510  0.7765722  0.02036714
## 
## Accuracy was used to select the optimal model using the largest value.
## The final value used for the model was cp = 0.002081756.

Visualizing Tree using package rpart.plot

# viasulaziation
library(rpart.plot)
prp(DTModel$finalModel, box.palette = "Reds", tweak = 1.2, varlen = 20)

Variable Importance in Decision Tree Model

# plotting variable importance
plot(varImp(DTModel))

TESTING THE DECISION TREE MODEL

Predicting Model on Test Data Set

# predicting the model on test data set
PredDTModel <- predict(DTModel, 
                    newdata = test.df,type = "prob")

Plotting the Predicted Probabilities

# plot of probabilities
plot(PredDTModel$Yes, 
     main = "Scatterplot of Probabilities of default (test data)", 
     xlab = "Customer ID", 
     ylab = "Predicted Probability of default")

Confusion Matrix at 50% Cut-Off Probability

# taking the cut-off probability 50%
Predicted <- ifelse(PredDTModel$Yes > 0.50, "Yes", "No")

# ordering the vectors
Actual <- ifelse(test.df$Default=="1","Yes","No")
# making confusion matrix
cm <-confusionMatrix(table(Predicted,Actual))
cm

## Confusion Matrix and Statistics
## 
##          Actual
## Predicted   No  Yes
##       No  4530 1230
##       Yes   69   91
##                                           
##                Accuracy : 0.7806          
##                  95% CI : (0.7698, 0.7911)
##     No Information Rate : 0.7769          
##     P-Value [Acc > NIR] : 0.2516          
##                                           
##                   Kappa : 0.0785          
##                                           
##  Mcnemar's Test P-Value : <2e-16          
##                                           
##             Sensitivity : 0.98500         
##             Specificity : 0.06889         
##          Pos Pred Value : 0.78646         
##          Neg Pred Value : 0.56875         
##              Prevalence : 0.77686         
##          Detection Rate : 0.76520         
##    Detection Prevalence : 0.97297         
##       Balanced Accuracy : 0.52694         
##                                           
##        'Positive' Class : No              
## 

ROC Curve using Package ROCR

# loading the package
library(ROCR)
DTPrediction <- predict(DTModel, test.df,type = "prob")
Prediction <- prediction(DTPrediction[2],test.df$Default)
performance <- performance(Prediction, "tpr","fpr")
# plotting ROC curve
plot(performance,main = "ROC Curve",col = 2,lwd = 2)
abline(a = 0,b = 1,lwd = 2,lty = 3,col = "black")

AUC (Area Under the Curve)

library(ROCR)
# area under curve
DTPrediction <- prediction(PredDTModel[2],test.df$Default)
aucDT <- performance(DTPrediction, measure = "auc")
aucDT <- aucDT@y.values[[1]]
aucDT

## [1] 0.5726524