ssgraph-package {ssgraph}

Bayesian Graphical Estimation using Spike-and-Slab Priors

Description

The R package ssgraph is for Bayesian estimation of graphical models by using spike-and-slab priors. The package is implemented the recent improvements in the Bayesian graphical models' literature, including Wang (2015). To speed up the computations, the computationally intensive tasks of the package are implemented in C++ in parallel using OpenMP.

How to cite this package

To cite ssgraph in publications use:

Mohammadi R. (2020). ssgraph: Bayesian Graphical Estimation using 
Spike-and-Slab Priors, R package version 1.11, 
https://cran.r-project.org/package=ssgraph

Author(s)

Reza Mohammadi [e-mail]

References

Wang, H. (2015). Scaling it up: Stochastic search structure learning in graphical models, Bayesian Analysis, 10(2):351-377

George, E. I. and McCulloch, R. E. (1993). Variable selection via Gibbs sampling. Journal of the American Statistical Association, 88(423):881-889

Griffin, J. E. and Brown, P. J. (2010) Inference with normal-gamma prior distributions in regression problems. Bayesian Analysis, 5(1):171-188

Mohammadi, A. et al (2017). Bayesian modelling of Dupuytren disease by using Gaussian copula graphical models, Journal of the Royal Statistical Society: Series C, 66(3):629-645

Mohammadi, R. and Wit, E. C. (2019). BDgraph: An R Package for Bayesian Structure Learning in Graphical Models, Journal of Statistical Software, 89(3):1-30

Mohammadi, A. and Wit, E. C. (2015). Bayesian Structure Learning in Sparse Gaussian Graphical Models, Bayesian Analysis, 10(1):109-138

Examples

## Not run: 

library( ssgraph )

# Generating multivariate normal data from a 'random' graph
data.sim <- bdgraph.sim( n = 100, p = 8, size = 10, vis = TRUE )

# Running algorithm based on GGMs
ssgraph.obj <- ssgraph( data = data.sim, iter = 5000, save = TRUE )

summary( ssgraph.obj )

# To compare the result with true graph

compare( pred = ssgraph.obj, actual = data.sim, 
         main = c( "Target", "ssgraph" ), vis = TRUE )
         
plotroc( pred = sgraph.obj, actual = data.sim )

## End(Not run)