close all
clear all 
disp('Ddimer: ROC')
lw = 2;  
set(0, 'DefaultAxesFontSize', 16);
fs = 14;
msize = 10;
%Rudimentary ROC in D-dimer Example    
[a1 b1 c1 d1 e1 f1 g1 h1]=sesp(42,96,3,32);
[a2 b2 c2 d2 e2 f2 g2 h2]=sesp(31,33,14,95);

figure(1)
%Figure roc2.eps
plot([0, 1-b2, 1-b1,1],[0 a2 a1 1],'k-','linewidth',lw)
hold on
%  ylabel('\rm Sensitivity','Interpreter','LaTeX')
%  xlabel('\rm 1-Specificity','Interpreter','LaTeX')
   ylabel('Sensitivity')
   xlabel('1-Specificity')
 
fill([0, 1-b2, 1-b1,1,1,0],[0 a2 a1 1,0,0],'y')
plot([0, 1-b2, 1-b1,1], [0 a2 a1 1],'o',  ...
       'MarkerSize',msize,'MarkerEdgeColor','k',...
       'MarkerFaceColor',[.49 1 .63])  
plot([0, 1],[0,   1],'g--','LineWidth',2)             
     text('Interpreter','latex',...
	'String','$(0.2578, 0.6889)$',...
	'Position',[0.27 0.65],...
	'FontSize',17)

    text('Interpreter','latex',...
	'String','$(0.75, 0.9333)$',...
	'Position',[0.45 .95],...
	'FontSize',17) 

     text('Interpreter','latex',...
	'String','$\fbox{Area = 0.7297}$',...
	'Position',[0.55 .45],...
	'FontSize',17) 

   % print -depsc 'C:\BESTAT\ROC\ROCeps\roc2.eps'
    area = auc([0, 1-b2, 1-b1,1],[0 a2 a1 1])
    
    
    