%
%
%Solution to Assignment 7, part 2, exercise 3
%
%F. Gabbiani 11/25/08
%


rate1 = 30; %in spike/sec
rate2 = 50; %in spike/sec
interval = 200; %in msec
max_count = 20; %in spikes


[cnt1, pd1, pd2] = poisscnt2(rate1,rate2,interval,max_count);

h_fig = figure(1);

%sets decorations and plots the isi distribution
plot(cnt1,pd1,'k--');
hold on; 
plot(cnt1,pd2,'k');
titt = sprintf('Spike count distributions for an interval of %g [msec]',...
               interval);
title(titt);
xlabel('spike count');
ylabel('probability');


[pd, pf, err] = roc(cnt1,pd1,pd2);

h_fig = figure(2);
set(h_fig,'Position',[296 420 560 800]);


%sets decorations and plots the isi distribution
subplot(2,1,1);
plot(pf,pd,'k',(0:0.1:1),(0:0.1:1),'k--');
titt = sprintf('ROC curve');
title(titt);
xlabel('probability of false alarm');
ylabel('probability of correct detection');
subplot(2,1,2);
plot(pf,err,'k',(0:0.1:1),ones(11,1)*0.5,'k--');
title('Error rate');
xlabel('probability of false alarm');
ylabel('probability of error');