diff --git a/matlab/GetAllPosteriorDraws.m b/matlab/GetAllPosteriorDraws.m index c17de8ece2fc7c3c16b7bd44c1fcde9eb795e7da..856268d93383427f203b09943bb778bd3d7e1867 100644 --- a/matlab/GetAllPosteriorDraws.m +++ b/matlab/GetAllPosteriorDraws.m @@ -9,20 +9,21 @@ DirectoryName = CheckPath('metropolis'); Draws = zeros(NumberOfDraws*nblck,1); logpo = zeros(NumberOfDraws*nblck,1); ipost=0; -if column<0, +if column<=0, column=1; ipost=1; end - -for file = FirstMhFile:TotalNumberOfMhFile - for blck = 1:nblck +iline0=iline; +for blck = 1:nblck + iline=iline0; + for file = FirstMhFile:TotalNumberOfMhFile load([DirectoryName '/' M_.fname '_mh' int2str(file) '_blck' int2str(blck)],'x2','logpo2') NumberOfLines = size(x2(iline:end,:),1); Draws(linee:linee+NumberOfLines-1) = x2(iline:end,column); logpo(linee:linee+NumberOfLines-1) = logpo2(iline:end); linee = linee+NumberOfLines; + iline = 1; end - iline = 1; end if ipost,