Abstract:
|
Gaussian graphical models, where it is assumed that the variables of interest jointly follow multivariate normal distributions with sparse precision matrices, have been used to study intrinsic dependence among several variables, but the Gaussianity assumption may be restrictive in many applications. A nonparanormal graphical model is a nonparametric generalization of a Gaussian graphical model for continuous variables where it is assumed that the variables follow a Gaussian graphical model only after some unknown smooth monotone transformation. We consider a Bayesian approach in the nonparanormal graphical model by putting priors on the unknown transformations through a random series based on B-splines where the coefficients are ordered to induce monotonicity. A truncated normal prior leads to partial conjugacy in the model and is useful for posterior simulation using Gibbs sampling. On the underlying precision matrix of the transformed variables, we consider a continuous shrinkage prior on its Cholesky decomposition and use an efficient posterior Gibbs sampling scheme. We study numerical performance of the proposed method through a simulation study and apply it on a real dataset.
|