Description Usage Arguments Details Value References See Also Examples
These are the density and random generation functions for the HuangWand prior distribution for a covariance matrix.
1 2 3 4  dhuangwand(x, nu=2, a, A, log=FALSE)
dhuangwandc(x, nu=2, a, A, log=FALSE)
rhuangwand(nu=2, a, A)
rhuangwandc(nu=2, a, A)

x 
This is a k x k positivedefinite
covariance matrix Sigma for 
nu 
This is a scalar degrees of freedom parameter
nu. The default is 
a 
This is a positiveonly vector of scale parameters a of length k. 
A 
This is a positiveonly vector of of scale hyperparameters
A of length k. Larger values result in a more
uninformative prior. A default, uninformative prior is

log 
Logical. If 
Application: Continuous Multivariate
Density: p(theta) = W^(1)[nu+k1](2*nu*diag(1/a)) G^(1)(1/2, 1/A^2)
Inventor: Huang and Wand (2013)
Notation 1: theta ~ HW[nu](a, A)
Notation 2: p(theta) ~ HW[nu](theta  a, A)
Parameter 1: degrees of freedom nu
Parameter 2: scale a > 0
Parameter 3: scale A > 0
Mean:
Variance:
Mode:
Huang and Wand (2013) proposed a prior distribution for a covariance matrix that uses a hierarchical inverse Wishart. This is a more flexible alternative to the inverse Wishart distribution, and the HuangWand prior retains conjugacy. The Cholesky parameterization is also provided here.
The HuangWand prior distribution alleviates two main limitations of an inverse Wishart distribution. First, the uncertainty in the diagonal variances of a covariance matrix that is inverse Wishart distributed is represented with only one degrees of freedom parameter, which may be too restrictive. The HuangWand prior overcomes this limitation. Second, the inverse Wishart distribution imposes a dependency between variance and correlation. The HuangWand prior lessens, but does not fully remove, this dependency.
The standard deviations of a HuangWand distributed covariance matrix are halft distributed, as HT(nu, A). This is in accord with modern assumptions about distributions of scale parameters, and is also useful for sparse covariance matrices.
The rhuangwand
function allows either a
or A
to be
missing. When a
is missing, the covariance matrix is generated
from the hyperparameters. When A
is missing, the covariance
matrix is generated from the parameters.
dhuangwand
and dhuangwandc
give the density, and
rhuangwand
and rhuangwandc
generate random deviates.
Huang, A., Wand, M., et al. (2013), "Simple Marginally Noninformative Prior Distributions for Covariance Matrices". Bayesian Analysis, 8, p. 439–452.
dhalft
and
dinvwishart
1 2 3 4  library(LaplacesDemon)
dhuangwand(diag(3), nu=2, a=runif(3), A=rep(1e6,3), log=TRUE)
rhuangwand(nu=2, A=rep(1e6, 3)) #Missing a
rhuangwand(nu=2, a=runif(3)) #Missing A

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.
Please suggest features or report bugs in the GitHub issue tracker.
All documentation is copyright its authors; we didn't write any of that.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.