search

merliseclyde / BAS

BAS R package https://merliseclyde.github.io/BAS/
https://merliseclyde.github.io/BAS/
GNU General Public License v3.0
41 stars 16 forks source link

Jeffreys prior with MCMC returns marginal prob > 1 #33

Closed merliseclyde closed 5 years ago

merliseclyde commented 5 years ago

rexpr: in test-bas-glm.R

test_that("Jeffreys & MCMC", {
 data(Pima.tr, package="MASS")
 pima_BAS <-  bas.glm(type ~ .,
                      data = Pima.tr, method = "MCMC",
                      betaprior = Jeffreys(),
                      family = binomial(),
                      modelprior = tr.beta.binomial(1, 1, 4))
 expect_equal(0, sum(pima_BAS$probne0 > 1))