Network Analysis in Psychopathology: Practical Guide in R (qgraph)

Network analysis applied to psychopathology has changed the way many researchers think about mental disorders. In contrast to the traditional model, which assumes that symptoms are indicators of an underlying latent disease (just as fever is an indicator of an infection), the network approach proposes that the symptoms themselves and their interactions are the disorder. There is no hidden entity called "depression" that causes sadness, insomnia, and fatigue; rather, these symptoms activate each other, forming a network that can become trapped in a pathological state.

This perspective has roots in dynamical systems theory and emerged strongly in the early 2010s, driven especially by the work of Denny Borsboom and his collaborators at the University of Amsterdam. The idea is not entirely new (functional and behavioral approaches have been questioning the notion of disorder as an entity for decades), but network analysis provided for the first time concrete statistical tools to model these interactions between symptoms quantitatively.

Anhedonia hub Sadness Fatigue Insomnia Generalized anxiety Rumination Concen- tration Irrita- bility Depressive symptom network EBIC-glasso · k = 8 symptoms · edge width = partial-correlation magnitude |r| > 0.40 |r| 0.20 to 0.40 |r| < 0.20
Symptom network estimated with EBIC-glasso. Anhedonia operates as a hub: it concentrates the strongest connections with sadness, fatigue, and insomnia.

How a symptom network is constructed

In practice, a psychopathological network is constructed from questionnaire or clinical interview data in which multiple symptoms are assessed. Each symptom is represented as a node, and the connections between nodes (the edges) reflect statistical associations between symptoms after controlling for all others. The most widely used method for continuous data is the regularized partial correlation network, where a LASSO-type penalty (specifically, the graphical LASSO or gLASSO) is applied to eliminate spurious connections and retain only the most robust associations. For ordinal or binary data, Ising models or polychoric correlations are used.

The visual result is a graph where symptoms with strong associations appear connected by thick lines and close to each other, while symptoms with little relationship appear distant and without direct connection. This visualization, although attractive, can be misleading if not interpreted carefully. Node positioning algorithms (such as Fruchterman-Reingold) optimize the visual layout, but the spatial distance between nodes does not always have a clear substantive interpretation. What matters are the edges and their weights, not the absolute positions in the graph.

Centrality measures and their interpretation

One of the most widely used (and most debated) concepts in psychopathological network analysis is centrality. Centrality measures attempt to identify which symptoms are most important or influential within the network. The three most common are strength (the sum of the absolute weights of all connections of a node), betweenness (how many shortest paths between pairs of nodes pass through that node), and closeness (the inverse of the sum of the shortest distances to all other nodes).

Of these three, strength is the most stable and the one that replicates best across studies. Betweenness and closeness have been shown to be highly unstable in typical psychopathological networks, to the point that several authors advise against their use in this context. Hallquist et al. (2019) and others have argued convincingly that these measures, originally developed for social networks with thousands of nodes, do not perform well in small networks with 15 or 20 symptoms, which is the usual scenario in psychopathology.

Beyond classic centrality measures, there has been growing interest in bridge centrality, which identifies symptoms that connect different communities or clusters within the network. The clinical idea is appealing: these bridge symptoms could be responsible for comorbidity between disorders. If insomnia connects the depressive symptom cluster with the anxious symptom cluster, intervening on insomnia could reduce both. It is an interesting hypothesis, but for now the empirical evidence supporting it is limited.

Limitations and current state of the field

Network analysis in psychopathology has received important criticisms that are worth knowing. The most fundamental is that networks estimated from cross-sectional data (which constitute the vast majority) do not allow causal or temporal inferences. The fact that insomnia and fatigue are strongly connected in a cross-sectional network does not mean that one causes the other, nor that intervening on one will improve the other. Forbes et al. (2017) demonstrated that the structure of cross-sectional networks can differ substantially from the structure of temporal networks estimated with longitudinal data, which calls into question the clinical utility of cross-sectional networks.

Another important limitation is replicability. Psychopathological networks are sensitive to sample composition, to the instruments used, and to the researcher's analytical decisions (what type of regularization to use, what hyperparameters to select, how to treat ordinal data). Some studies have found that the global network structure replicates reasonably across samples, but the centrality measures of specific nodes are much less stable. This is problematic if the promise of the approach is precisely to identify key symptoms for intervention.

Despite these limitations, network analysis remains a valuable tool in psychopathological research. The most promising developments include temporal networks estimated with ecological momentary assessment (EMA) data, Bayesian network models that allow the incorporation of prior knowledge, and longitudinal panel networks that combine between-subjects and within-subjects information. The psychonetrics package has also consolidated confirmatory and panel network models, which let you test a theoretical network structure rather than merely estimating one exploratorily. The field is maturing, and with that maturation comes a more realistic understanding of what these techniques can and cannot offer. If you work with clinical data and are interested in exploring this approach, the bootnet package in R is a good starting point, as it includes stability and accuracy procedures that will help you evaluate the robustness of your results.

Before you submit: network papers get picked apart on exactly the points above. Centrality reported without any stability check, a cross-sectional network described in causal language, betweenness and closeness listed as if they were as trustworthy as strength. If your manuscript is already written, run it through the AI Paper Reviewer: a free Reviewer 2 style pre-review that tells you, in under a minute, which single problem is most likely to sink the submission, while you can still fix it instead of reading it in the editor's letter three months from now.

If your thesis or paper uses a network model and you want to be sure centrality and stability are correctly estimated and reported, my statistical consulting fits the network in R (qgraph, bootnet) and prepares the report for your committee or reviewer. And if your interest is identifying patient profiles rather than symptom-to-symptom interaction, latent class analysis may fit you better.

Keep reading

All blog articles