R/logpriorh2.R

Defines functions logpriorh2

Documented in logpriorh2

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

Try the bbemkr package in your browser

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

bbemkr documentation built on May 1, 2019, 10:53 p.m.