supplementaries/Mode Jumping MCMC/supplementary/examples/BAS archive/BAS/R/modelpriors.R

In aliaksah/EMJMCMC2016: EMJMCMC

uniform = function() {
  structure(list(family="Uniform",hyper.parameters=.5), class="prior")
}

Bernoulli=function(probs=0.5) {
  if (length(probs) == 1) {
    if (probs == .5) structure(list(family="Uniform",hyper.parameters=.5), class="prior")
    else structure(list(family="Bernoulli", hyper.parameters=probs), class="prior")}
  else   structure(list(family="Bernoulli", hyper.parameters=probs), class="prior")
}

beta.binomial=function(alpha=1.0, beta=1.0) {
    structure(list(family="Beta-Binomial", hyper.parameters=c(alpha, beta)),
              class="prior")
 }