R/logpriorh2.R

logpriorh2 <- function(h2, prior_alpha, prior_beta)
{
    logp = prior_alpha * log(prior_beta) - lgamma(prior_alpha) - 1.0*(prior_alpha + 1.0)*log(h2) - prior_beta/h2	
    return(logp)
}

Try the bbefkr package in your browser

Any scripts or data that you put into this service are public.

bbefkr documentation built on May 2, 2019, 3:04 a.m.