permits analysis of incomplete Metropolis run

time-shift
Michel Juillard 2011-10-13 17:12:54 +02:00
parent e92c54ae3a
commit 19a5ccb6f1
1 changed files with 7 additions and 4 deletions

View File

@ -38,10 +38,13 @@ logPost = [];
params = [];
blck = 1;
for i=1:record.LastFileNumber
o=load([M.fname '/metropolis/' M.fname '_mh' int2str(i) '_blck' ...
int2str(blck)]);
logPost = [logPost; o.logpo2];
params = [params; o.x2];
fname = [M.fname '/metropolis/' M.fname '_mh' int2str(i) '_blck' ...
int2str(blck) '.mat'];
if exist(fname,'file')
o=load(fname);
logPost = [logPost; o.logpo2];
params = [params; o.x2];
end
end
figure;