Élaguer l'arbre de décision

Modélisation du risque de crédit en R

Lore Dirick

Manager of Data Science Curriculum at Flatiron School

Problèmes des grands arbres de décision

Trop complexe : moins lisible
Surapprentissage sur l'ensemble de test
Solution : utiliser printcp(), plotcp() pour élaguer

Printcp et tree_undersample

printcp(tree_undersample)

Classification tree:
rpart(formula = loan_status ~ ., data = undersampled_training_set, method = "class",
 control = rpart.control(cp = 0.001))
Variables actually used in tree construction:
age    annual_inc     emp_cat     grade    home_ownership   ir_cat     loan_amnt     
Root node error: 2190/6570 = 0.33333
n= 6570 
        CP    nsplit  rel error   xerror      xstd
1  0.0059361      0    1.00000   1.00000   0.017447
2  0.0044140      4    0.97443   0.99909   0.017443
3  0.0036530      7    0.96119   0.98174   0.017366
4  0.0031963      8    0.95753   0.98904   0.017399
               ...  
16 0.0010654     76    0.84247   1.02511   0.017554
17 0.0010000     79    0.83927   1.02511   0.017554

Plotcp et tree_undersample

Capture d'écran du 2020-06-22 à 17 h 57 min 10 s.png

Plotcp et tree_undersample

Capture d'écran du 2020-06-22 à 17 h 56 min 53 s.png

$CP = 0.003653$

Tracer l'arbre élagué

Capture d'écran du 2020-06-23 à 18 h 16 min 04 s.png

ptree_undersample=prune(tree_undersample,
                        cp = 0.003653)

plot(ptree_undersample,
     uniform=TRUE)

text(ptree_undersample)

Tracer l'arbre élagué

Capture d'écran du 2020-06-23 à 18 h 15 min 42 s.png

ptree_undersample=prune(tree_undersample,
                        cp = 0.003653)

plot(ptree_undersample,
     uniform=TRUE)

text(ptree_undersample,
     use.n=TRUE)

prp() dans le paquet rpart.plot

Capture d'écran du 2020-06-22 à 18 h 05 min 09 s.png

library(rpart.plot)
prp(ptree_undersample)

prp() dans le paquet part.plot

Capture d'écran du 2020-06-22 à 18 h 04 min 33 s.png

library(rpart.plot)
prp(ptree_undersample, extra = 1)

Passons à la pratique !

Modélisation du risque de crédit en R