Fix the inconsistency between mh_drop and the number of mh draws discarded
when computing the Brooks and Gelman convergence diagnostics.time-shift
parent
1ae1b5ff20
commit
b762f010cd
|
@ -91,6 +91,7 @@ localVars.nblck = nblck;
|
||||||
localVars.NumberOfMcFilesPerBlock = NumberOfMcFilesPerBlock;
|
localVars.NumberOfMcFilesPerBlock = NumberOfMcFilesPerBlock;
|
||||||
localVars.Origin = Origin;
|
localVars.Origin = Origin;
|
||||||
localVars.StepSize = StepSize;
|
localVars.StepSize = StepSize;
|
||||||
|
localVars.mh_drop = options_.mh_drop;
|
||||||
localVars.NumberOfDraws = NumberOfDraws;
|
localVars.NumberOfDraws = NumberOfDraws;
|
||||||
localVars.NumberOfLines = NumberOfLines;
|
localVars.NumberOfLines = NumberOfLines;
|
||||||
localVars.time = time;
|
localVars.time = time;
|
||||||
|
|
|
@ -65,7 +65,7 @@ for j=fpar:npar,
|
||||||
ligne = 0;
|
ligne = 0;
|
||||||
for iter = Origin:StepSize:NumberOfDraws
|
for iter = Origin:StepSize:NumberOfDraws
|
||||||
ligne = ligne+1;
|
ligne = ligne+1;
|
||||||
linea = ceil(0.5*iter);
|
linea = ceil(mh_drop*iter);
|
||||||
n = iter-linea+1;
|
n = iter-linea+1;
|
||||||
cinf = round(n*ALPHA/2);
|
cinf = round(n*ALPHA/2);
|
||||||
csup = round(n*(1-ALPHA/2));
|
csup = round(n*(1-ALPHA/2));
|
||||||
|
|
Loading…
Reference in New Issue