sampleKappa | R Documentation |
This is an implementation of the algorithm of Forbes and Mardia (2015) to
sample from the Bessel exponential distribution, which is the conditional
distribution of the concentration parameter of a von Mises distribution
given the mean mu
. The distribution is proportional to exp(-
η g κ)/I_0(κ)^η. Note that beta_0
in Forbes and
Mardia (2015) is renamed g
here.
sampleKappa(etag, eta)
etag |
Numeric; This is |
eta |
Integer; This is the posterior sample size, which is n + c where
c is the number of observations contained in the conjugate prior. For
uninformative, |
A sampled value kappa
from the Bessel exponential
distribution.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.