## File Name: pmle_eval_prior.R
## File Version: 0.03
pmle_eval_prior <- function(prior, pars, eps=1e-100)
{
NP <- length(pars)
priorval <- rep(NA, NP)
eps <- 1e-100
for (pp in 1:NP){
prior_arg_pp <- prior[[pp]][[2]]
prior_arg_pp[[1]] <- pars[pp]
priorval[pp] <- log( do.call( prior[[pp]][[1]], prior_arg_pp ) + eps )
}
return(priorval)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.