Log penalized complexity precision prior like INLA's internal — log_prec_prior • ktools

Log penalized complexity precision prior like INLA's internal

log_prec_prior(log_prec, u, alpha)

Arguments

log_prec: log precision
u: u
alpha: alpha

Examples

sd <- seq(0.001, 2, 0.01)
prec <- sd2prec(sd)
lprec <- log(prec)
log_prec_prior(lprec, 0.33, 0.01) %>%
  exp() %>%
  plot(sd, .)
log_prec_prior(lprec, 1, 0.01) %>%
  exp() %>%
  lines(sd, .)