MR2 error

[kspych]

I am trying to run the MR2 function with 1 exposure and 4 responses (see code below).
MR2_output <- MR2(beta_Y, beta_X, EVgamma = c(1, 1),niter = 7500, burnin = 2500, thin = 5, monitor = 500, seed = 28061971)

However, I am receiving the following error:
Error in if (log(runif(1)) < logAccRatio) { :
missing value where TRUE/FALSE needed

The example 1 simulation ran just fine. Is this happening because I only have 1 exposure?

[lb664]

Yes, you are right. We never thought to have so few exposures. The problem is the BetaBinomial prior, and indeed you see an error message on the MH acceptance probability which includes the distribution of the BetaBinomial prior. Even with a parsimonious parametrisation (1,1), it may require no less than 4 exposures.

Anyway, your comment is interesting, and in the next release, we will add a (less flexible) Binomial prior with user-specific probability of inclusion. Then the number of exposures should not cause any problems.