The Brier score is an important measure of calibration i.e. the mean squared difference between the predicted probability and the actual outcome.
Lower the Brier score is for a set of predictions, the better the predictions are calibrated.
- If the predicted probability is 1 and it happens, then the Brier Score is 0, the best score achievable.
- If the predicted probability is 1 and it does not happen, then the Brier Score is 1, the worst score achievable.
- If the predicted probability is 0.8 and it happens, then the Brier Score is (0.8-1)^2 =0.04.
- If the predicted probability is 0.2 and it happens, then the Brier Score is (0.2-1)^2 =0.64.
- If the predicted probability is 0.5, then the Brier Score is (0.5-1)^2 =0.25, irregardless of whether it happens.
By specifying fitstat option in proc logistic, SAS returns Brier score and other fit statistics such as AUC, AIC, BIC etc.
proc logistic data=train;
model y(event="1") = entry/ outroc = rocstats;
score data=valid out=valpred fitstat;
run;
Hi Deepanshu,
ReplyDeleteI tried getting Brier score using fitstat option but seems I didnt get it.. Can you please place a snapshot here? Do i need to add any other option to get this score.
Which version of SAS you are using? The fitstat option in score statement is available in SAS versions starting from SAS 9.3
Deletehello,
ReplyDeletehow is it possible to perform easily a calibration plot (observed probabilities as a function of the expected one) for a logistic model ?
best regards