The R
package bgms provides tools for Bayesian analysis of graphical models describing networks of variables. The package uses Markov chain Monte Carlo methods combined with a pseudolikelihood approach to estimate the posterior distribution of model parameters. Gibbs variable selection (George and McCulloch 1993) is used to model the underlying network structure of the graphical model. By imposing a discrete spike and slab prior on the pairwise interactions, it is possible to shrink the interactions to exactly zero. The Gibbs sampler embeds a Metropolis approach for mixtures of mutually singular distributions (Gottardo and Raftery 2008) to account for the discontinuity at zero. The goal is to provide these tools for Markov Random Field (MRF) models for a wide range of variable types in the bgms package, and it currently provides them for analyzing networks of binary and/or ordinal variables (Marsman and Haslbeck 2023).
Multivariate analysis using graphical models has received much attention in the recent psychological and psychometric literature (Robinaugh et al. 2020; Marsman and Rhemtulla 2022; Contreras et al. 2019). Most of these graphical models are Markov Random Field (MRF) models, whose graph structure reflects the conditional associations between variables (Kindermann and Snell 1980). In these models, a missing edge between two variables in the network implies that these variables are independent, given the remaining variables (Lauritzen 2004). In other words, the remaining variables of the network fully account for the potential association between the unconnected variables.
Testing the structure of the MRF requires us to determine the plausibility of the opposing hypotheses of conditional dependence and conditional independence. That is, how plausible is it that the observed data come from a network with a structure that includes the edge between two variables compared to a network structure that excludes that edge? Frequentist approaches are limited in this regard because they can only reject the conditional independence hypothesis, not support it (Wagenmakers et al. 2018; Wagenmakers 2007). This leads to the problem that if an edge is excluded, we do not know whether this is because the edge is absent in the population or because we lack the power to reject the null hypothesis of independence. To avoid this problem, we will use a Bayesian approach using Bayes factors (Kass and Raftery 1995)). The inclusion Bayes factor (Huth et al. 2023; Sekulovski et al. 2024) allows us to quantify how much the data support both conditional dependence -evidence of edge presence- or conditional independence -evidence of edge absence. It also allows us to conclude that there is limited support for either hypothesis (Dienes 2014)-an absence of evidence.
Contreras, A., I. Nieto, C. Valiente, R. Espinosa, and C. Vazquez. 2019. âThe Study of Psychopathology from the Network Analysis Perspective: A Systematic Review.â
Psychotherapy and Psychosomatics88 (2): 71â83.
https://doi.org/10.1159/000497425.
Dienes, Z. 2014. âUsing Bayes to Get the Most Out of Non-Significant Results.â
Frontiers in Psychology5 (781): 1â17.
https://doi.org/10.3389/fpsyg.2014.00781.
George, E. I., and R. E. McCulloch. 1993. âVariable Selection via Gibbs Sampling.â
Journal of the American Statistical Association88 (423): 881â89.
https://doi.org/10.1080/01621459.1993.10476353.
Gottardo, R., and A. E. Raftery. 2008. âMarkov Chain Monte Carlo with Mixtures of Mutually Singular Distributions.â
Journal of Computational and Graphical Statistics17 (4): 949â75.
https://doi.org/10.1198/106186008X386102.
Huth, K., J. de Ron, A. E. Goudriaan, K. Luigjes, R. Mohammadi, R. J. van Holst, E.-J. Wagenmakers, and M. Marsman. 2023. âBayesian Analysis of Cross-Sectional Networks: A Tutorial in R and JASP.â
Advances in Methods and Practices in Psychological Science6 (4): 1â18.
https://doi.org/10.1177/25152459231193334.
Kass, R. E., and A. E. Raftery. 1995. âBayes Factors.â
Journal of the American Statistical Association90 (430): 773â95.
https://doi.org/10.2307/2291091.
Kindermann, R., and J. L. Snell. 1980. Markov Random Fields and Their Applications. Vol. 1. Contemporary Mathematics. Providence: American Mathematical Society.
Lauritzen, S. L.. 2004. Graphical Models. Oxford: Oxford University Press.
Marsman, M., and J. M. B. Haslbeck. 2023. âBayesian Analysis of the Ordinal Markov Random Field.â
PsyArXiv.
https://doi.org/10.31234/osf.io/ukwrf.
Marsman, M., and M. Rhemtulla. 2022. âGuest Editorsâ Introduction to the Special Issue âNetwork Psychometrics in Actionâ: Methodological Innovations Inspired by Empirical Problems.â
Psychometrika87 (1): 1â11.
https://doi.org/10.1007/s11336-022-09861-x.
Robinaugh, D. J., R. H. A. Hoekstra, E. R. Toner, and D. Borsboom. 2020. âThe Network Approach to Psychopathology: A Review of the Literature 2008â2018 and an Agenda for Future Research.â
Psychological Medicine50: 353â66.
https://doi.org/10.1017/S0033291719003404.
Sekulovski N, Keetelaar S, Huth K, Wagenmakers E.-J, van Bork R, van den Bergh D, Marsman M (2024). âTesting conditional independence in psychometric networks: An analysis of three bayesian methods.â
Multivariate Behavioral Research, 1â21.
doi:10.1080/00273171.2024.2345915.
Wagenmakers, E.-J. 2007. âA Practical Solution to the Pervasive Problems of p Values.â
Psychonomic Bulletin & Review14: 779â804.
https://doi.org/10.3758/BF03194105.
Wagenmakers, E.-J., M. Marsman, T. Jamil, A. Ly, J. Verhagen, J. Love, R. Selker, et al. 2018. âBayesian Inference for Psychology. Part I: Theoretical Advantages and Practical Ramifications.â
Psychonomic Bulleting & Review25 (1): 58â76.
https://doi.org/10.3758/s13423-017-1343-3.
RetroSearch is an open source project built by @garambo | Open a GitHub Issue
Search and Browse the WWW like it's 1997 | Search results from DuckDuckGo
HTML:
3.2
| Encoding:
UTF-8
| Version:
0.7.4