summary & gelmanDiagnostics for uni-variate case : Error in if (plot == T & !is.na(diag$mpsrf)) do.call(gP, as.list(match.call())) :

[vincentchu-fr]

Hello,

I observed the following error when using the functions summary and gelmanDiagnostics for uni-variate problems: Error in if (plot == T & !is.na(diag$mpsrf)) do.call(gP, as.list(match.call())) :

In the following, you'll find the scripts which reproduce the issue. I don't know whether this is an issue of BayesianTools itself or it comes coda

library("BayesianTools") library("stats") ############

Create bayesian setup

bayesianSetup <- createBayesianSetup(likelihood = testDensityNormal, prior = createUniformPrior(lower = -10, upper = 10))

running MCMC

out = runMCMC(bayesianSetup = bayesianSetup)

summary(out)

plot(out)

gelmanDiagnostics(out, plot = T)

############

bayesianSetup <- createBayesianSetup(likelihood = testDensityMultiNormal, prior = createUniformPrior(lower = rep(-10, 3), upper = rep(10, 3)))

running MCMC

out = runMCMC(bayesianSetup = bayesianSetup)

summary(out)

plot(out)

gelmanDiagnostics(out, plot = T)