errors with intrinsic prior

[merliseclyde]

REPREX

test_that("intrinsic prior for GLM", {
  data(Pima.tr, package = "MASS")
  expect_error( pima.BAS <- bas.glm(type ~ .,
                      data = Pima.tr, n.models = 2^7, method = "BAS",
                      betaprior = intrinsic(n = nrow(Pima.tr)), family = binomial(),
                      modelprior = uniform()
  ))
#  expect_equal(0, sum(pima.BAS$shrinkage > 1))
})

Error in bas.glm(type ~ ., data = Pima.tr, n.models = 2^7, method = "BAS", :
REAL() can only be applied to a 'numeric', not a 'integer'