https://github.com/Data2Dynamics/d2d
Tip revision: d72d4c223fb8adc0e3e44a1d3d1eafce18b995ae authored by Andreas Raue on 06 March 2015, 00:17:07 UTC
Increase C version code, after Joep's submits today to prevent compatibility problems
Increase C version code, after Joep's submits today to prevent compatibility problems
Tip revision: d72d4c2
arPlotPriors.m
for i=1:50:length(ar.p)
figure;
if i+50 > length(ar.p)
xdata = linspace(i,length(ar.p),length(ar.p)-i+1);
types = ar.type(i:length(ar.p));
means = ar.mean(i:length(ar.p));
stds = ar.std(i:length(ar.p));
means(types==0) = nan;
stds(types==0) = nan;
errorbar(xdata,means,stds,'b.');
xlim([i length(ar.p)]);
ylim([-5.5 3.5]);
set(gca,'tickdir','out');
hold on
plot(xdata,ar.p(i:length(ar.p)),'r*');
set(gca,'XTick',xdata);
set(gca,'XTickLabel',ar.pLabel(i:length(ar.p)));
view(gca,[90 90]);
hold off
else
xdata = linspace(i,i+49,50);
types = ar.type(i:i+49);
means = ar.mean(i:i+49);
stds = ar.std(i:i+49);
means(types==0) = nan;
stds(types==0) = nan;
errorbar(xdata,means,stds,'b.');
xlim([i i+49]);
ylim([-5.5 3.5]);
set(gca,'tickdir','out');
hold on
plot(xdata,ar.p(i:i+49),'r*');
set(gca,'XTick',xdata);
set(gca,'XTickLabel',ar.pLabel(i:i+49));
view(gca,[90 90]);
hold off
end
end
% figure;
% hold on
% errorbar(linspace(1,length(ar.p),length(ar.p)),ar.mean,ar.std,'b.');
% xlim([0 length(ar.p)]);
% plot(linspace(1,length(ar.p),length(ar.p)),ar.p,'r*');
% hold off