Abstract:To understand psychological data, it is crucial to examine the structure and dimensions of variables. In this study, we examined alternative estimation algorithms to the conventional GLASSO-based exploratory graph analysis (EGA) in network psychometric models to assess the dimensionality structure of the data. The study applied Bayesian conjugate or Jeffreys' priors to estimate the graphical structure and then used the Louvain community detection algorithm to partition and identify groups of nodes, which allowed the detection of the multi- and unidimensional factor structures. Monte Carlo simulations suggested that the two alternative Bayesian estimation algorithms had comparable or better performance when compared with the GLASSO-based EGA and conventional parallel analysis (PA). When estimating the multidimensional factor structure, the analytically based method (i.e., EGA.analytical) showed the best balance between accuracy and mean biased/absolute errors, with the highest accuracy tied with EGA but with the smallest errors. The sampling-based approach (EGA.sampling) yielded higher accuracy and smaller errors than PA; lower accuracy but also lower errors than EGA. Techniques from the two algorithms had more stable performance than EGA and PA across different data conditions. When estimating the unidimensional structure, the PA technique performed the best, followed closely by EGA, and then EGA.analytical and EGA.sampling. Furthermore, the study explored four full Bayesian techniques to assess dimensionality in network psychometrics. The results demonstrated superior performance when using Bayesian hypothesis testing or deriving posterior samples of graph structures under small sample sizes. The study recommends using the EGA.analytical technique as an alternative tool for assessing dimensionality and advocates for the usefulness of the EGA.sampling method as a valuable alternate technique. The findings also indicated encouraging results for extending the regularization-based network modeling EGA method to the Bayesian framework and discussed future directions in this line of work. The study illustrated the practical application of the techniques to two empirical examples in R.

Dimension constraints improve hypothesis testing for large-scale, graph-associated, brain-image data

Hypothesis Testing of Matrix Graph Model with Application to Brain Connectivity Analysis.

Multiple Comparison Procedures for Neuroimaging Genomewide Association Studies

Bayesian Inference and Testing of Group Differences in Brain Networks

Simultaneous Inference in Multiple Matrix-Variate Graphs for High-Dimensional Neural Recordings

Large-Scale Simultaneous Testing Using Kernel Density Estimation

Optimal False Discovery Rate Control for Large Scale Multiple Testing with Auxiliary Information

Graph Combinatorics Based Group-Level Network Inference for Brain Connectome Analysis

Nonparametric Bayes multiresolution testing for high-dimensional rare events

Whole-brain Prediction Analysis with GraphNet

Multiscale statistical testing for connectome-wide association studies in fMRI

Valid Two-Sample Graph Testing via Optimal Transport Procrustes and Multiscale Graph Correlation with Applications in Connectomics

Large-scale Multiple Testing of Cross-covariance Functions with Applications to Functional Network Models

MATRIX GRAPH HYPOTHESIS TESTING AND APPLICATION IN BRAIN CONNECTIVITY ALTERNATION DETECTION

Exploring Estimation Procedures for Reducing Dimensionality in Psychological Network Modeling

A Hypothesis Testing for Large Weighted Networks with Applications to Functional Neuroimaging Data.

Simultaneous Cluster Structure Learning and Estimation of Heterogeneous Graphs for Matrix‐variate Fmri Data

A semiparametric two-sample hypothesis testing problem for random dot product graphs

Large-scale directed network inference with multivariate transfer entropy and hierarchical statistical testing

Large-scale adaptive multiple testing for sequential data controlling false discovery and nondiscovery rates

False Discovery Rate Control for Lesion‐Symptom Mapping With Heterogeneous Data via Weighted p‐Values