% selection of departments arount Paris
R=COVID19Modelingv2('ParisSuburb');
nR = length(R);
%% plot
clf
col = [ 0 0 0; .5 .5 .5; tooclear(cbrewer('qual','Set2',nR-2)) ];
datefmt = 'dd/mm/yyyy';
E = datenum('11/05/2020',datefmt);
formatfig(gcf,'figname',datestr(today,'ParisSuburb_yyyy-mm-dd'),'position',[956 312 730 950],'paperposition',[1.5000    4.1750   18.0000   21.3500])
hs = subplots(1,[.4 1],0,0); hax=axes('position',[.4936 .2846 .3813 .20]);
hp = zeros(nR,1);
for i=1:nR
    subplot(hs(1)), hold on
    ymax = R(i).data.ypred(end);
    plot(R(i).data.date,R(i).data.ypredCI(:,2)/ymax,'-','color',col(i,:),'linewidth',3)
    plot(R(i).data.sampledDate,min(1,R(i).data.infected/ymax),'o','markerfacecolor',col(i,:),'markeredgecolor','w','markersize',12)
    subplot(hax), hold on
    hp(i) = plot(R(i).data.date,R(i).data.pdfpred/max(R(i).data.pdfpred),'-','color',col(i,:),'linewidth',3);
    ind = (R(i).data.date>=today) & (R(i).data.date<=E);
    plot(R(i).data.date(ind),R(i).data.pdfpred(ind)/max(R(i).data.pdfpred),'-','color',col(i,:),'linewidth',6)

    subplot(hs(2)), hold on
    hp(i) = plot(R(i).data.date,R(i).data.pdfpred/max(R(i).data.pdfpred),'-','color',col(i,:),'linewidth',3);
    ind = (R(i).data.date>=today) & (R(i).data.date<=E);
    plot(R(i).data.date(ind),R(i).data.pdfpred(ind)/max(R(i).data.pdfpred),'-','color',col(i,:),'linewidth',6)
end
formatax(hax,'fontsize',10,'ylim',[-.002 .17],'xtick',E+(-14:7:21))
datetick(hs(1),'x','dd/mm')
datetick(hs(2),'x','dd/mm')
datetick(hax,'x','dd-mm-yy')
formatax(hs,'fontsize',14,'ylim',[-.02 1.1],'xlim',[datenum('20/03/2020',datefmt) datenum('20/06/2020',datefmt)])
set(hax,'xlim',[datenum('20/04/2020',datefmt) datenum('01/06/2020',datefmt)])
subplot(hs(2))
legend(hp,regexprep([R.entity],'^DPT ',''),'location','best','fontsize',12,'box','off')