function Evaluate_Probabilistic_Classifier(Y_tst,Y_hat,p_hat,w)

num_ex = length(Y_tst);

[p_hat_s,ind_s] = sort(p_hat);
Y_tst_s = Y_tst(ind_s);
Y_hat_s = Y_hat(ind_s);

for i = 1:num_ex-w
    
    predicted_prob(i) = sum(p_hat_s(i:i+w))/(w+1);
    
    t1 = Y_tst_s(i:i+w) - Y_hat_s(i:i+w);
    ind_corr = find(t1==0);
    actual_prob_est(i) = length(ind_corr)/(w+1);
    err_in_pr(i) = predicted_prob(i) - actual_prob_est(i);
end

figure;
plot(predicted_prob,actual_prob_est,'r*-');
hold on
xlabel('Predicted Probability of Correct Classification','FontSize',14);
ylabel('Actual Probability','FontSize',14);
t1 = [predicted_prob,actual_prob_est];
min_t1 = min(t1);
plot([min_t1 1], [min_t1 1], 'b--');
hold off;

t2 = 5;