diagnostic plot error for bas.glm objects

[merliseclyde]

When using MCMC, the length of the MCMC estimates of the posterior inclusion probabilities is the number of models. This leads to an error in the diagnostic plot function.

REPREX:

data(Pima.tr, package = "MASS")
pima.MCMC <- bas.glm(type ~ .,
    data = Pima.tr, MCMC.iterations = 2^10,
    method = "MCMC", betaprior = aic.prior(),
    family = binomial(),
    modelprior = tr.poisson(ncol(Pima.tr) - 1, 7)
  )
  expect_null(diagnostics(pima.MCMC, type = "pip"))

In tests/testthat/test-bas_glm.R