R/proposal_d.R

In aashen12/BARS: Bayesian Adaptive Regression Splines

proposal_d <- function(k) {
  d <- 1/3 #probability of death
  b <- 1/3 #probability of birth
  # NOTE: P(change) = 1/3
  num <- log(d) - log(k)
  den <- log(b) - log(2)
  return(den - num)
} #proposal in RJMCMC ratio for a death step