This articles discusses about various model validation techniques of a classification or logistic regression model. The below validation techniques do not restrict to logistic regression only. It can be used for other classification techniques such as decision tree, random forest, gradient boosting and other machine learning techniques. These validation techniques are considered as benchmarks for comparing predictive models in marketing analytics and credit risk modeling domain. Model validation is a crucial step of a predictive modeling project.

Primarily there are three methods of validation. They are listed below -

- Split Sample Validation
- Cross Validation
- Bootstrapping Validation

*The detailed explanation of these methods are listed below -*

**1. Split Sample Validation**

- Randomly split data into two samples: 70% = training sample, 30% = validation sample.
- Score (predicted probability) the validation sample using the response model under consideration.
- Rank the scored file, in descending order by estimated probability
- Split the ranked file into 10 sections (deciles)
- Number of observations in each decile
- Number of actual events in each decile
- Number of cumulative actual events in each decile
- Percentage of cumulative actual events in each decile. It is called Gain Score.
- Divide the gain score by % of data used in each portion of 10 bins. For example, in second decile, divide gain score by 20.
- Calculate
**KS statistics**(It measure of the degree of separation between the positive and negative distributions. In other words, it checks the maximum difference between distribution of cumulative events and cumulative non-events)

**Model Validation Metrics**

**1. KS Statistics**

KS Test measures to check whether model is able to separate events and non-events. In probability of default (bank defaulters) model, it checks whether the credit risk model is able to distinguish between good and bad customers. The calculation of KS test is explained below -

KS = Maximum difference between Cumulative % Event and Cumulative % Non-Event

KS Test |

**Important Note -**In this case, KS is maximum at third decile and KS score is 59.1. Ideally, it should be in first three deciles and score lies between 40 and 70. And there should not be more than 10 points (in absolute) difference between training and validation KS score. Score above 70 is susceptible and might be overfitting so rigorous validation is required.

**Calculating KS Test with SAS**

**2. Rank Ordering**

To see rank ordering, calculate the percentage of events (defaults) in each decile group and check the event rate should be monotonically decreasing. It means the model predicts the highest number of events in the first decile and then goes progressively down. You can check the rank ordering in the table shown above. The rank ordering is maintained in this example.

**3. Area under curve**

It explains the trade-off between true positive rate (Sensitivity) and false positive rate (1-Specificity).

**It is calculated by summing Concordance value and (0.5 times of Tied Percent).**

ROC Curve |

AUC should be more than 0.7 in both the training and validation samples. Should not be a significant difference between AUC score of both these samples. If it is more than 0.8, it is considered as an excellent score.

**4. Hosmer Lemeshow Test**

**5. Lift Chart**

It measures how much better one can expect to do with the predictive model comparing without a model.

**Understand Gain and Lift Chart**

**Model Validation Rules : Summary**

- Same significant variables should come in both the training and validation sample.
- The behavior of variables should be same in both the samples (same sign of coefficients)
**Beta coefficients**should be close in training and validation samples**KS statistics**should be in top 3 deciles.- KS statistics should be between 40 and 70. Should not be significantly different from Training KS (more than 10 points in absolute)
**Rank Ordering**- There should not be any break in rank ordering.**Lift Curve**- The larger the cumulative lift value the better the accuracy**Area under Curve (AUC) -**Should be more than 0.7.**Goodness of Fit Tests -**Model should fit the data well. Check Hosmer and Lemeshow Test and Deviance and Residual Test.

**SAS Code : Model Validation - Logistic Regression**

/* Split data into two datasets : 70%- training 30%- validation */

Proc Surveyselect data=finaldata out=split samprate=.7 outall;

Run;

Data training validation;

Set split;

if selected = 1 then output training;

else output validation;

Run;

/* Logistic Model*/

ods graphics on;

Proc Logistic Data = training;

Model Sbp_flag = age_flag sex bmi_flag/ lackfit ctable pprob =0.5 outroc=troc;

Output out= test p=ppred;

score data=validation out = Logit_File outroc=vroc;

Run;

Proc rank data= logit_file descending groups=10 out=predrank_Dev;

var P_1 ;

/****PREDRANK - variable name for ranks****/

ranks predrank;

run ;

PROC TABULATE Data= predrank_Dev;

CLASS predrank;

VAR Sbp_flag;

TABLE predrank=" " all, Sbp_flag*(N="Count" SUM="Number of Responses" COLPCTSUM="% of events") / box="Decile";

TITLE Creating a Lift Table for Model;

RUN;

Proc Freq data = predrank_Dev;

tables Sbp_flag;

run;

Proc npar1way data=logit_file edf;

class Sbp_flag;

var p_1;

run;

**2. Cross Validation**

**1. Jack-knife / Leave-one-out**: The model is fitted on all the cases except one observation and is then tested on the set-aside case. This procedure can be repeated as many times as the number of observations in the original sample (random without replacement sampling). It is implemented in PROC LOGISTIC with

**predprobs=crossvalidate.**

**Limitation :**

If the model is tested on a single observation, it is not possible to assess one of the most important dimensions of model’s performance, i.e. calibration (measure of how close the predicted probabilities are to the actual rate of events).

**2. K-fold cross-validation :**Splits the data into K subsets; each is held out in turn as the validation set. (Random without replacement technique)

**10-fold cross- validation :**

- Randomly divide your data into ten parts.
- Hold aside the first tenth of the data as a validation dataset; fit a logistic model using the remaining 9/10 (the training dataset).
- Using the resulting training model, calculate the predicted probability for each validation observation.
- Repeat this 9 more times (so that each tenth of the dataset becomes the validation dataset exactly once).
- Now, you have a predicted probability for each observation from a model that was not based on that observation.
- An AUC score is calculated for each of the 10 runs, and then calculate average AUC.

**SAS Code : Jack-knife / Leave-one-out Validation**

proc logistic data=fil.liver descending;

model complications = age_at_op comorb lobeormore_code bilat_resec_code numsegs_resec;

output out=preds predprobs=crossvalidate;

run;

ods select none;

ods output WilcoxonScores=WilcoxonScore;

proc npar1way wilcoxon data= preds;

where complications^=.;

class complications;

var XP_1;

run;

ods select all;

data AUC;

set WilcoxonScore end=eof;

retain v1 v2 1;

if _n_=1 then v1=abs(ExpectedSum - SumOfScores);

v2=N*v2;

if eof then do;

d=v1/v2;

Gini=d * 2;

AUC= d + 0.5;

put AUC= GINI=;

keep AUC Gini;

output;

end;

run;

proc logistic data=preds descending;

model complications = age_at_op comorb lobeormore_code bilat_resec_code numsegs_resec;

roc pred=xp_1;

roccontrast;

run;

%macro xval(dsn=,outcome=,covars=,k=10,sel=stepwise,outdsn=_xval_,outdsn2=comparison);

data _modif;

set &dsn;

unif=&k*ranuni(20052905);

xv=ceil(unif);

run;

%do i=1 %to &k;

proc logistic data=_modif(where=(xv ne &i)) outmodel=_mod&i;

model &outcome (event="1") =&covars / selection=&sel;

ods output association=assoc&i;

run;

%if print^=0 %then %do;proc printto file='junk.txt';%end;

proc logistic inmodel=_mod&i;

score data=_modif(where=(xv=&i)) out=out&i;

run;

ods select none;

ods output KolSmir2Stats=KS&i;

proc npar1way data= out&i edf;

where &outcome^=.;

class &outcome;

var P_1;

run;

ods select all;

ods select none;

ods output WilcoxonScores=Wil&i;

proc npar1way wilcoxon data= out&i;

where &outcome^=.;

class &outcome;

var P_1;

run;

ods select all;

data AUC&i;

set Wil&i end=eof;

retain v1 v2 1;

if _n_=1 then v1=abs(ExpectedSum - SumOfScores);

v2=N*v2;

if eof then do;

d=v1/v2;

/*Gini=d * 2; */

Scoring_AUC= d + 0.5;

put Scoring_AUC=;

put "****Open work.results dataset to see results of training datasets....";

keep Scoring_AUC;

output;

end;

run;

%if print^=0 %then %do;proc printto;run;%end;

%end;

data &outdsn;

set %do j=1 %to &k;out&j %end;;

run;

data training (keep =label2 nvalue2 rename= (nvalue2=Training_AUC));

set %do j=1 %to &k;assoc&j %end;;

where label2= 'c';

if label2='c' then label2 ='AUC';

run;

data ks (keep =label2 nvalue2 rename= (nvalue2=Scoring_KS));

set %do j=1 %to &k;ks&j %end;;

where label2= 'D';

if label2='D' then label2 ='KS';

run;

data validation;

set %do j=1 %to &k;AUC&j %end;;

run;

data &outdsn2 (drop = label2);

merge training validation ks;

run;

ods select none;

ods output WilcoxonScores=WilcoxonScore;

proc npar1way wilcoxon data= &outdsn;

where &outcome^=.;

class &outcome;

var P_1;

run;

ods select all;

data AUC;

set WilcoxonScore end=eof;

retain v1 v2 1;

if _n_=1 then v1=abs(ExpectedSum - SumOfScores);

v2=N*v2;

if eof then do;

d=v1/v2;

Gini=d * 2;

AUC= d + 0.5;

put AUC= GINI=;

put "****Open work.results dataset to see results of training datasets....";

keep AUC Gini;

output;

end;

run;

%mend;

%xval(dsn=alldata,outcome=y,covars=entry,k=5,sel=stepwise, outdsn=kfold, outdsn2=comparison);

model complications = age_at_op comorb lobeormore_code bilat_resec_code numsegs_resec;

output out=preds predprobs=crossvalidate;

run;

ods select none;

ods output WilcoxonScores=WilcoxonScore;

proc npar1way wilcoxon data= preds;

where complications^=.;

class complications;

var XP_1;

run;

ods select all;

data AUC;

set WilcoxonScore end=eof;

retain v1 v2 1;

if _n_=1 then v1=abs(ExpectedSum - SumOfScores);

v2=N*v2;

if eof then do;

d=v1/v2;

Gini=d * 2;

AUC= d + 0.5;

put AUC= GINI=;

keep AUC Gini;

output;

end;

run;

proc logistic data=preds descending;

model complications = age_at_op comorb lobeormore_code bilat_resec_code numsegs_resec;

roc pred=xp_1;

roccontrast;

run;

**Rule :**Compare area under curve of both the samples.**SAS Macro : K-Fold Cross Validation***The following SAS program was written by Mithat Gonen. I modified it to calculate AUC of validation dataset and store modeling results of each fold in a dataset.*%macro xval(dsn=,outcome=,covars=,k=10,sel=stepwise,outdsn=_xval_,outdsn2=comparison);

data _modif;

set &dsn;

unif=&k*ranuni(20052905);

xv=ceil(unif);

run;

%do i=1 %to &k;

proc logistic data=_modif(where=(xv ne &i)) outmodel=_mod&i;

model &outcome (event="1") =&covars / selection=&sel;

ods output association=assoc&i;

run;

%if print^=0 %then %do;proc printto file='junk.txt';%end;

proc logistic inmodel=_mod&i;

score data=_modif(where=(xv=&i)) out=out&i;

run;

ods select none;

ods output KolSmir2Stats=KS&i;

proc npar1way data= out&i edf;

where &outcome^=.;

class &outcome;

var P_1;

run;

ods select all;

ods select none;

ods output WilcoxonScores=Wil&i;

proc npar1way wilcoxon data= out&i;

where &outcome^=.;

class &outcome;

var P_1;

run;

ods select all;

data AUC&i;

set Wil&i end=eof;

retain v1 v2 1;

if _n_=1 then v1=abs(ExpectedSum - SumOfScores);

v2=N*v2;

if eof then do;

d=v1/v2;

/*Gini=d * 2; */

Scoring_AUC= d + 0.5;

put Scoring_AUC=;

put "****Open work.results dataset to see results of training datasets....";

keep Scoring_AUC;

output;

end;

run;

%if print^=0 %then %do;proc printto;run;%end;

%end;

data &outdsn;

set %do j=1 %to &k;out&j %end;;

run;

data training (keep =label2 nvalue2 rename= (nvalue2=Training_AUC));

set %do j=1 %to &k;assoc&j %end;;

where label2= 'c';

if label2='c' then label2 ='AUC';

run;

data ks (keep =label2 nvalue2 rename= (nvalue2=Scoring_KS));

set %do j=1 %to &k;ks&j %end;;

where label2= 'D';

if label2='D' then label2 ='KS';

run;

data validation;

set %do j=1 %to &k;AUC&j %end;;

run;

data &outdsn2 (drop = label2);

merge training validation ks;

run;

ods select none;

ods output WilcoxonScores=WilcoxonScore;

proc npar1way wilcoxon data= &outdsn;

where &outcome^=.;

class &outcome;

var P_1;

run;

ods select all;

data AUC;

set WilcoxonScore end=eof;

retain v1 v2 1;

if _n_=1 then v1=abs(ExpectedSum - SumOfScores);

v2=N*v2;

if eof then do;

d=v1/v2;

Gini=d * 2;

AUC= d + 0.5;

put AUC= GINI=;

put "****Open work.results dataset to see results of training datasets....";

keep AUC Gini;

output;

end;

run;

%mend;

%xval(dsn=alldata,outcome=y,covars=entry,k=5,sel=stepwise, outdsn=kfold, outdsn2=comparison);

**Note:**Calculate AUC of each of the runs and then calculate average of all of the 10 runs.

**3. Bootstrapping Validation**

In other words, you randomly draw your observations with replacement, then calculates the logistic regression and stores the coefficients. This is repeated n times. So you'll end up with 10'000 different regression coefficients.

**Rule :**If a variable is truly representative of the model it will occur in the majority of the N fitted models. The c statistics is calculated for each iteration in order to examine the predicted probability of each model.The overall accuracy is the average of the N measures.

*The following SAS program was written by Mithat Gonen.*

%macro bval(dsn=,outcome=,covars=,B=10,sel=stepwise);

proc sql noprint;

select n(&outcome) into:_n from &dsn;

run;

proc surveyselect data=&dsn method=urs outhits rep=&B n=&_n out=bsamples noprint;

run;

%do i=1 %to &B;

proc logistic data=bsamples(where=(replicate=&i)) outmodel=_mod&i descending noprint;

model &outcome=&covars / selection=&sel;

run;

proc logistic inmodel=_mod&i;

score data=&dsn out=out1&i;

run;

score data=bsamples(where=(replicate=&i)) out=out2&i;

run;

proc printto;run;

%end;

%end;

data bval1;

set %do j=1 %to &B;out1&j(in=in&j) %end;;

%do j=1 %to &B; if in&j then bsamp=&j; %end;

run;

data bval2;

set %do j=1 %to &B;out2&j(in=in&j) %end;;

%do j=1 %to &B; if in&j then bsamp=&j; %end;

run;

proc printto file='junk.txt' new;

proc logistic data=bval1 descending;

by bsamp;

model &outcome=p_1;

ods output association=assoc1;

run;

proc logistic data=bval2;

by bsamp;

model &outcome=p_1;

ods output association=assoc2;

run;

set %do j=1 %to &B;out1&j(in=in&j) %end;;

%do j=1 %to &B; if in&j then bsamp=&j; %end;

run;

data bval2;

set %do j=1 %to &B;out2&j(in=in&j) %end;;

%do j=1 %to &B; if in&j then bsamp=&j; %end;

run;

proc printto file='junk.txt' new;

proc logistic data=bval1 descending;

by bsamp;

model &outcome=p_1;

ods output association=assoc1;

run;

proc logistic data=bval2;

by bsamp;

model &outcome=p_1;

ods output association=assoc2;

run;

proc logistic data=&dsn;

model &outcome=&covars / selection=&sel;

ods output association=assoc3;

run;

proc printto;

data assoc3;

set assoc3;

bsamp=1;

run;

data optim;

merge assoc1(where=(label2='c') keep=bsamp label2 nvalue2 rename=(nvalue2=auc1))

assoc2(where=(label2='c') keep=bsamp label2 nvalue2 rename=(nvalue2=auc2))

assoc3(where=(label2='c') keep=bsamp label2 nvalue2 rename=(nvalue2=auc3));

by bsamp;

run;

proc sql;

select mean(auc3) as OptimisticAUC, mean(auc2-auc1) as OptimisimCorrection,

mean(auc3)-mean(auc2-auc1) as CorrectedAUC from optim;

quit;

%mend;

%bval(dsn=liver,outcome=complications,covars=age_at_op comorb lobeormore_code bilat_resec_code numsegs_resec,B=10);

model &outcome=&covars / selection=&sel;

ods output association=assoc3;

run;

proc printto;

data assoc3;

set assoc3;

bsamp=1;

run;

data optim;

merge assoc1(where=(label2='c') keep=bsamp label2 nvalue2 rename=(nvalue2=auc1))

assoc2(where=(label2='c') keep=bsamp label2 nvalue2 rename=(nvalue2=auc2))

assoc3(where=(label2='c') keep=bsamp label2 nvalue2 rename=(nvalue2=auc3));

by bsamp;

run;

proc sql;

select mean(auc3) as OptimisticAUC, mean(auc2-auc1) as OptimisimCorrection,

mean(auc3)-mean(auc2-auc1) as CorrectedAUC from optim;

quit;

%mend;

%bval(dsn=liver,outcome=complications,covars=age_at_op comorb lobeormore_code bilat_resec_code numsegs_resec,B=10);

Thanks mate for all your hardwork.

ReplyDeleteKs statistics using npar1way is recorded in decimals(fragments) with p-value. How do u get 40-70?

After getting the ranked probabilities, how do u know wch customer belongs to a decile?

Thanks mate for all your hardwork.

ReplyDeleteKs statistics using npar1way is recorded in decimals(fragments) with p-value. How do u get 40-70?

After getting the ranked probabilities, how do u know wch customer belongs to a decile?

By 40-70, i meant 40% - 70% (0.4 - 0.7). Run PROC RANK with GROUPS = 10 to know the customer placement in a decile.

DeleteNice work, Deepanshu! The way you listed steps and SAS codes for model validation in logistic regression is really helpful. It would be more helpful if you have a one line statement regarding each SAS code stating what it is doing and where does it belong in the 10 steps split sample validation.

ReplyDeleteThanks,

VB

The bootstrap validation code seems to produce an error message. The statement below

ReplyDeleteproc logistic data=bval1 descending;

by bsamp;

model &outcome=p_1;

ods output association=assoc1;

run;

is where the problem is -- SAS complains that the "Variable "p_1 is not found". The same error occurs in PROC LOGISTIC for the data=BVAL2.

Please help!

Great work

ReplyDeleteReally helpful post Deepanshu, I wanted to know the thoery behind calculation of KS statistic, Gini Coefficient. Also how to calculate the GOF value?

ReplyDeleteIts good work

ReplyDeleteawesome job Deepanshu,

ReplyDeleteBut on the theoretical note, how to we score with bootstrapping logistic regression ??

Pardon,

ReplyDeletewhere can I find the dataset of the previous SAS codes?

Will you please do the same (or provide code) with the GLM model.

ReplyDeleteThanks in advance