hold off
plot([M.t(:,1:2:end) M.r(:,1)],'-.');
hold on
plot(M.est);
co = get(gca,'ColorOrder');
for i=1:5
ind =1:10:90;
h=errorbar(ind,M.est(ind,i),...
M.est(ind,i)-M.lb(ind,i),M.est(ind,i)-M.ub(ind,i),'.');
set(h,'color',co(i,:));
end

set(gca,'YLim',[-0.7,2]);