function scoreStair2(t, scoreN, bestIndexList)
h1=line(t,log10(scoreN),'marker','.','color',[0.8 0.8 0.8],'markerfacecolor',[0.8 0.8 0.8], 'marker','o','markersize',3,'linestyle','none');
hold on
h2=stairs(t(bestIndexList),log10(scoreN(bestIndexList)),'k-','linewidth',2);
ylabel('logarithmic score');
set(gca,'ytick',[])
xlabel('time (days)')