Hi,
I found out that the Raftery and Lewis (1992) convergence test (called with “raftery_lewis_diagnostics” in the estimation command) is not implemented when “mh_nblocks=2”. Is there a reason for that? I know that when there are two chains, the Brooks and Gelman tests are possible, but why does it mean that the Raftery and Lewis test will not be implemented, since it was requested in the estimation command? It is all in the function ‘McMCDiagnostics.m’.
Yaakov