R/00-prior-definitions.R

In awstringer1/casecrossover: Approximate Bayesian Inference for Case Crossover Models

# Script containing prior definitions
#

# PC PRECISION pc.prec
# Returns log of pc precision prior for theta, the log-precision
# Works for multidim theta
pcprec <- function(theta,u,alpha) {
  lambda <- -log(alpha)/u
  sum(log(lambda/2) - lambda * exp(-theta/2) - theta/2)
}