ssgraph-package | R Documentation |
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.
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
Reza Mohammadi a.mohammadi@uva.nl
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
## 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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.