Psychometric Network Models from Time-Series and Panel Data

Sacha Epskamp

doi:10.1007/s11336-020-09697-3

Psychometric Network Models from Time-Series and Panel Data

Published online by Cambridge University Press: 01 January 2025

Sacha Epskamp

Show author details

Sacha Epskamp*: Affiliation:
University of Amsterdam
*: Correspondence should be made to Sacha Epskamp, Department of Psychology: Psychological Methods Groups, University of Amsterdam, PO Box 15906, 1001 NK Amsterdam, The Netherlands. Email: sacha.epskamp@gmail.com

Article contents

Abstract
Introduction
Preliminary Topics
Time-Series Data: The ts-lvgvar Model
Panel Data: The panel-lvgvar Model
Simulation Studies
Ergodicity and Generalizability
Discussion
Footnotes
References

Rights & Permissions

Abstract

Researchers in the field of network psychometrics often focus on the estimation of Gaussian graphical models (GGMs)—an undirected network model of partial correlations—between observed variables of cross-sectional data or single-subject time-series data. This assumes that all variables are measured without measurement error, which may be implausible. In addition, cross-sectional data cannot distinguish between within-subject and between-subject effects. This paper provides a general framework that extends GGM modeling with latent variables, including relationships over time. These relationships can be estimated from time-series data or panel data featuring at least three waves of measurement. The model takes the form of a graphical vector-autoregression model between latent variables and is termed the ts-lvgvar when estimated from time-series data and the panel-lvgvar when estimated from panel data. These methods have been implemented in the software package psychonetrics, which is exemplified in two empirical examples, one using time-series data and one using panel data, and evaluated in two large-scale simulation studies. The paper concludes with a discussion on ergodicity and generalizability. Although within-subject effects may in principle be separated from between-subject effects, the interpretation of these results rests on the intensity and the time interval of measurement and on the plausibility of the assumption of stationarity.

Keywords

network psychometrics Gaussian graphical model structural equation modeling dynamics time-series data panel data

Information

Type: Theory and Methods
Information: Psychometrika , Volume 85 , Issue 1 , March 2020 , pp. 206 - 231

DOI: https://doi.org/10.1007/s11336-020-09697-3 [Opens in a new window]
Creative Commons: This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Copyright: Copyright © 2020 The Author(s)

1. Introduction

Researchers in the field of network psychometrics (Marsman et al., Reference Marsman, Borsboom, Kruis, Epskamp, van Bork, Waldorp and Marsman2018) study the estimation of multivariate statistical models in attempts to map out the complex interplay of interactions between variables. This field emerged from the network perspective on psychology (Borsboom, Reference Borsboom2017; Cramer, Waldorp, van der Maas, & Borsboom, Reference Cramer, Waldorp, van der Maas and Borsboom2010)—departing from the latent variable model and instead conceptualizing observed variables (e.g., attitudes, symptoms, and moods) as causal agents in a complex interplay of psychological (and other) components. In recent years, network models have become useful additions to the psychometric toolbox. For example, network models provide useful visualization tools to check the fit of latent variable models (Epskamp, Cramer, Waldrop, Schmittmann, & Borsboom, Reference Epskamp, Cramer, Waldrop, Schmittmann and Borsboom2012), are closely tied and often equivalent to latent variable models (Marsman et al., Reference Marsman, Borsboom, Kruis, Epskamp, van Bork, Waldorp and Marsman2018), are often uniquely identified (Epskamp, Waldorp, Mõttus, & Borsboom, Reference Epskamp, Waldorp, Mõttus and Borsboom2018), and may provide exploratory insight into the underlying factor structure by investigating its clustering (Golino & Epskamp, Reference Golino and Epskamp2017).

In a psychometric network model,Footnote ¹ variables are represented by nodes that are connected by edges, which are weighted according to some statistic. In this paper, I focus on two particular network models now routinely used in the analysis of continuous data: the Gaussian graphical model (GGM; Epskamp, Waldorp, et al., Reference Epskamp, Waldorp, Mõttus and Borsboom2018; Lauritzen, Reference Lauritzen1996) and the graphical vector-autoregression model (GVAR; Epskamp, Waldorp, et al., Reference Epskamp, Waldorp, Mõttus and Borsboom2018; Wild et al., Reference Wild, Eichler, Friederich, Hartmann, Zipfel and Herzog2010). The GGM forms an undirected network model, in which edges represent partial correlation coefficients. The GGM is closely tied, but not fully equivalent, to the directed structures typically used in structural equation modeling (SEM; Epskamp, Rhemtulla, & Borsboom, Reference Epskamp, Rhemtulla and Borsboom2017) and is often estimated from cross-sectional data (Fried et al., Reference Fried, van Borkulo, Cramer, Boschloo, Schoevers and Borsboom2017), in which many subjects are measured only once. The (lag-1) GVAR model takes the form of a generalization of the GGM in single-subject time series (Epskamp, Waldorp, et al., Reference Epskamp, Waldorp, Mõttus and Borsboom2018). In the GVAR, temporal dependencies are modeled via a regression on the previous measurement occasion, which leads to a matrix of regression coefficients that can also be used to draw a directed network model (Bringmann et al., Reference Bringmann, Vissers, Wichers, Geschwind, Kuppens, Peeters and Tuerlinckx2013)—often termed the temporal network because it encodes predictive effects over time. The remaining variances and covariances (i.e., the covariance structure after controlling for the previous measurement occasion) can be modeled as a GGM, which is also termed the contemporaneous network. When time series of multiple subjects are available, a third GGM can be formed on the between-subject effects (relationships between stable means)—also termed the between-subject network.

Current practices in network psychometrics feature two pressing limitations. (1) Although network modeling has been proposed as an alternative to latent variable modeling (i.e., covariation is caused by one or more unobserved common causes), the complete omission of latent variables may be one step too far (Bringmann & Eronen, Reference Bringmann and Eronen2018; Fried & Cramer, Reference Fried and Cramer2017; Guyon, Falissard, & Kop, Reference Guyon, Falissard and Kop2017). Network models for observed variables rely on the assumption that all causally interacting variables are observed without error. However, a certain level of measurement error (Schuurman, Houtveen, & Hamaker, Reference Schuurman, Houtveen and Hamaker2015) should be assumed in psychological data. (2) Network models are now often estimated from cross-sectional data. However, such results are not reflective of within-subject dynamics over time (Bos et al., Reference Bos, Snippe, de Vos, Hartmann, Simons, van der Krieke and Wichers2017; Molenaar, Reference Molenaar2004). Cross-sectional analysis cannot distinguish between within- and between-subject variances (Hamaker, Reference Hamaker, Mehl and Conner2012) or indeed investigate temporal effects over time. In principle, panel data, in which many subjects are measured a few times, can be used to distinguish between within- and between-subject effects (Hamaker, Kuiper, & Grasman, Reference Hamaker, Kuiper and Grasman2015), but have only been sparingly discussed in network psychometrics (e.g., Rhemtulla, Van Bork, & Cramer, Reference Rhemtulla, Van Bork and Cramer2017). Prior literature has tackled these topics separately. For example, Epskamp, Rhemtulla, and Borsboom (Reference Epskamp, Rhemtulla and Borsboom2017) addressed Limitation 1 by proposing to form a psychometric framework to incorporate GGMs in latent variable models, both at the latent level (latent network models) and at the residual level (residual network models). The multi-level GVAR model (Epskamp, Waldorp, et al., Reference Epskamp, Waldorp, Mõttus and Borsboom2018) can overcome Limitation 2 but does not include latent variables and requires intensive repeated measures of many subjects.

The goal in this paper is to combine the above-described solutions into a general framework. The framework emerges by combining a general factor model and a GVAR model and is referred to as the lvgvar model.Footnote ² I will discuss estimation in two very distinct settings:

• Time-series data of a single subject. This setting will be termed ts-lvgvar.
• Panel data of many subjects measured on a few (at least three) measurement occasions. This setting will be termed panel-lvgvar.

These settings differ crucially from one another. ts-lvgvar concerns a fixed subject (fixed subject p), which involves variation over measurement occasions (random time T). Meanwhile, panel-lvgvar concerns variation over subjects (random subject P) at a few fixed time points $t, t + 1, t + 2, \dots$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$t, t+1, t+2, \ldots $$\end{document} . The panel-lvgvar is a multi-level model with random effects on the mean structure. I showcase both methods using two empirical examples and assess the performance of model search strategies through two large-scale simulation studies. Finally, I discuss the generalizability of results between these settings in detail. The lvgvar framework is implemented in the open-source software package psychonetrics.Footnote ³

2. Preliminary Topics

2.1. Notation

This paper uses a similar notation as Epskamp, Waldorp, et al. (Reference Epskamp, Waldorp, Mõttus and Borsboom2018). Roman letters indicate observed variables, and Greek letters indicate parameters or latent variables. Bold-faced lowercase letters indicate vectors, and bold-faced uppercase letters indicate matrices. Normal-faced lowercase letters indicate fixed values, and normal-faced uppercase letters indicate random variables. These notations are also used in the subscript of a vector. For example, $y_{p, t}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {y}_{p,t}$$\end{document} indicates a fixed response pattern $y$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {y}$$\end{document} for subject p on fixed measurement occasion t; $y_{P, t}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {y}_{P,t}$$\end{document} indicates the response vector of a random subject P at fixed measurement occasion t; and $y_{p, T}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {y}_{p,T}$$\end{document} indicates the response vector of a fixed subject p at a random measurement occasion T.

2.2. Estimation

Suppose n cases (people or time points in this paper) are measured on $n_{z}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n_z$$\end{document} variables, with $z_{c}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {z}_c$$\end{document} denoting the response vector of case c, which is the cth row of a data matrix $Z$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {Z}$$\end{document} . The vector $z_{C}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {z}_C$$\end{document} will be assumed normally distributed:

\begin{matrix} z_{C} \sim N (μ, Σ), \end{matrix}

in which $μ$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\mu }$$\end{document} represents a mean vector and $Σ$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\Sigma }$$\end{document} represents a variance–covariance matrix.Footnote ⁴ Let $\bar{z}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\bar{\pmb {z}}$$\end{document} represent the sample means:

\begin{matrix} \bar{z} = \frac{1}{n} \sum_{c = 1}^{n} z_{c} . \end{matrix}

Furthermore, let $\bar{S}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\bar{\pmb {S}}$$\end{document} represent the sample variance–covariance matrix:

\begin{matrix} \bar{S} = \frac{1}{n} \sum_{c = 1}^{n} (z_{c} - \bar{z}) {(z_{c} - \bar{z})}^{⊤} . \end{matrix}

Using these summary statistics and assuming multivariate normality, we can estimate parameters by minimizing the following fit function:

(1)

\begin{matrix} F_{ML} = trace (S, Σ^{- 1}) + {(\bar{z} - μ)}^{⊤} Σ^{- 1} (\bar{z} - μ) - ln | Σ^{- 1} |, \end{matrix}

which is proportional to $- 2 / n$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$-2/n$$\end{document} times the log-likelihood of the data. Full information maximum likelihood estimation (FIML) can be used when data are missing. To this end, the data can be subdivided in subsets of data that have the same missingness patterns. The FIML fit function to be minimized takes the following form:

(2)

\begin{matrix} F_{FIML} = \frac{1}{n} \sum_{i} n_{i} (trace (S_{i}, Σ_{i}^{- 1}) + {(\bar{z_{i}} - μ_{i})}^{⊤} Σ_{i}^{- 1} (\bar{z_{i}} - μ_{i}) - ln | Σ_{i}^{- 1} |), \end{matrix}

in which $n_{i}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n_i$$\end{document} represents the sample size of subset i; $S_{i}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {S}_i$$\end{document} is the sample variance–covariance matrix of subset i (note, $S_{i} = O$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {S}_i = \pmb {O}$$\end{document} if $n_{i} = 1$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n_i = 1$$\end{document} ); $\bar{z_{i}}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\bar{\pmb {z}_i}$$\end{document} conotes the sample means of subset i (note, it is the same as the observed score if $n_{i} = 1$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n_i = 1$$\end{document} ); $Σ_{i}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\Sigma }_i$$\end{document} is a subset of $Σ$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\Sigma }$$\end{document} that only contains elements of observed data in subset i; and $μ_{i}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\mu }_i$$\end{document} is a subset of $μ$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\mu }$$\end{document} that only contains elements of observed data in subset i.

2.3. Factor Model

Let $η_{p, t}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\eta }_{p,t}$$\end{document} indicate a length $n_{η}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n_\eta $$\end{document} vector of variables of interest at time point t for subject p. If $η$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\eta }$$\end{document} is not observed, it may be assumed that $η$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\eta }$$\end{document} linearly causes observed indicators in a length $n_{y}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n_y$$\end{document} vector $y$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {y}$$\end{document} according to a linear factor model:

(3)

\begin{matrix} y_{p, t} = τ + Λ η_{p, t} + ε_{p, t}, \end{matrix}

in which $τ$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\tau }$$\end{document} is a general intercept vector; $Λ$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\Lambda }$$\end{document} is a matrix of factor loadings; and $ε_{p, t}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\varepsilon }_{p,t}$$\end{document} is a vector of residuals. Assume that the residuals, denoting measurement error or nuisance variables, have mean $0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {0}$$\end{document} across people but not necessarily across time:

\begin{matrix} E (ε_{P, t}) & = 0 No general bias across people. \\ E (ε_{p, T}) & = μ_{p}^{(ε)} Subject p may respond consistently biased . \end{matrix}

It is possible, without loss of information, to regard any of the above variable vectors as composites of the mean of subject p, denoted with $μ_{p}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\mu }_p$$\end{document} , and deviation from that mean in measurement t, denoted with $ξ_{p, t}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\xi }_{p,t}$$\end{document} :

\begin{matrix} y_{p, t} & = μ_{p}^{(y)} + ξ_{p, t}^{(y)} & Observed scores are composed of a mean plus deviation. \\ η_{p, t} & = μ_{p}^{(η)} + ξ_{p, t}^{(η)} & Latent variables are composed of a mean plus deviation. \\ ε_{p, t} & = μ_{p}^{(ε)} + ξ_{p, t}^{(ε)} & Residuals are composed of general bias plus time-specific bias, \end{matrix}

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \pmb {y}_{p,t}&= \pmb {\mu }^{(\pmb {y})}_{p} + \pmb {\xi }^{(\pmb {y})}_{p,t}&\text {Observed scores are composed of a mean plus deviation.}\\ \pmb {\eta }_{p,t}&= \pmb {\mu }^{(\pmb {\eta })}_{p} + \pmb {\xi }^{(\pmb {\eta })}_{p,t}&\text {Latent variables are composed of a mean plus deviation.} \\ \pmb {\varepsilon }_{p,t}&= \pmb {\mu }^{(\pmb {\varepsilon })}_{p} + \pmb {\xi }^{(\pmb {\varepsilon })}_{p,t}&\text {Residuals are composed of general bias plus time-specific bias}, \end{aligned}$$\end{document}

with the expected value of all $ξ$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\xi }$$\end{document} vectors concerning either random people or random measurement occasions set to equal zero and no assumed temporal dependency between the residual deviations. It is then possible to formulate the measurement model of Eq. (3) in a within-subject and a between-subject part:

(4)

\begin{matrix} μ_{p}^{(y)} & = τ + Λ μ_{p}^{(η)} + μ_{p}^{(ε)} & Between-subject measurement model stable over time. \\ ξ_{p, t}^{(y)} & = Λ ξ_{p, t}^{(η)} + ξ_{p, t}^{(ε)} & Within-subject measurement model varying over time. \end{matrix}

Note, the lack of subscripts on $Λ$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\Lambda }$$\end{document} and $τ$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\tau }$$\end{document} indicates an assumption of measurement invariance across subjects and time (Adolf, Schuurman, Borkenau, Borsboom, & Dolan, Reference Adolf, Schuurman, Borkenau, Borsboom and Dolan2014) and across within- and between-subject models.Footnote ⁵ Because of the assumption of measurement invariance across subjects in $τ$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\tau }$$\end{document} , no intercept is included in the within-subject measurement model.

2.4. Network Models

2.4.1. The Gaussian Graphical Model

The GGM can be formed as a model for a variance–covariance matrix $Σ$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\Sigma }$$\end{document} using the following notation (Epskamp, Rhemtulla, & Borsboom, Reference Epskamp, Rhemtulla and Borsboom2017):

(5)

\begin{matrix} Σ = Δ {(I - Ω)}^{- 1} Δ . \end{matrix}

In this expression, $Δ$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\Delta }$$\end{document} is a diagonal scaling matrix that controls the variances, and $Ω$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\Omega }$$\end{document} is a square symmetrical model matrix with zeroes on the diagonal and partial correlation coefficients on the off diagonal. Note, the expression above is identical to inverting $Σ$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\Sigma }$$\end{document} , standardizing the result, and multiplying all off-diagonal elements by $- 1$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$-1$$\end{document} . The sparsity (zeroes) in off-diagonal elements of $Ω$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\Omega }$$\end{document} will equal the sparsity in the precision matrix $Σ^{- 1}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\Sigma }^{-1}$$\end{document} , which can equivalently be modeled to obtain a GGM. However, the expression above has some benefits, namely in that the sign of the obtained parameters is in line with the interpretation and that the expression allows users to constrain partial correlations equally across groups while allowing the scale to vary freely (Kan, van der Maas, & Levine, Reference Kan, van der Maas and Levine2019).

2.4.2. The Graphical Vector-Autoregression Model

The GVAR model utilizes two network structures to extend the GGM when observations are temporally dependent: the temporal network and the contemporaneous network. Assuming stationarity in all parameters over time, the model can be written as a regression on the previous time point. The typical formation for observed variables takes the following form (ignoring a subscript p for subject):

\begin{matrix} y_{t} = μ + B (y_{t - 1} - μ) + ζ_{t}, ζ_{T} \sim N (0, Σ^{(ζ)}), \end{matrix}

in which $ζ_{t}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\zeta }_{t}$$\end{document} represents a vector of normally distributed innovations. Equivalently, the model can be written in terms of a conditional normal distribution:

\begin{matrix} y_{T} ∣ y_{T - 1} = y_{t - 1} \sim N (μ + B (y_{t - 1} - μ), Σ^{(ζ)}) . \end{matrix}

The matrix $B$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {B}$$\end{document} encodes temporal within-subject effects, and its transpose is typically used to display a personalized weighted directed network of temporal effects (Bringmann et al., Reference Bringmann, Vissers, Wichers, Geschwind, Kuppens, Peeters and Tuerlinckx2013). The matrix can also be standardized to partial directed correlations (Wild et al., Reference Wild, Eichler, Friederich, Hartmann, Zipfel and Herzog2010). The matrix $Σ^{(ζ)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\Sigma }^{(\pmb {\zeta })}$$\end{document} encodes within-subject contemporaneous effects—associations between variables in the same measurement occasion after taking temporal effects into account—and can be used to obtain a personalized undirected network structure with a GGM (Epskamp, Van Borkulo, et al., Reference Epskamp, Van Borkulo, Van Der Veen, Servaas, Isvoranu, Riese and Cramer2018) using Expression (5). The assumption of stationarity can further be used to obtain expressions for the stationary variance–covariance structure as well as the lag-k variance–covariance structure, as depicted below. Finally, the mean structure itself may vary across subjects, which could be used to form a GGM on the between-subject level: the between-subject network.

3. Time-Series Data: The ts-lvgvar Model

The ts-lvgvar concerns a fixed subject p and takes measurement T as random. Assuming multivariate normality for all variables, the following may be formulated:

\begin{matrix} y_{p, T} & \sim N (μ_{p}^{(y)}, Σ_{p, 0}^{(y)}) & Within-subject distribution of observed variables. \\ η_{p, T} & \sim N (μ_{p}^{(η)}, Σ_{p, 0}^{(η)}) & Within-subject distribution of latent variables. \end{matrix}

More generally, the within-subject lag-k covariances can be defined as:

\begin{matrix} cov (y_{p, T + k}, y_{p, T}) & = Σ_{p, k}^{(y)} Observed variables within-subject lag- k covariances. \\ cov (η_{p, T + k}, η_{p, T}) & = Σ_{p, k}^{(η)} Latent variables within-subject lag- k covariances, \end{matrix}

with the superscript indicating the variable of interest. As seen in Eq. (4), the expected mean vector $μ_{p}^{(y)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\mu }^{(\pmb {y})}_{p}$$\end{document} is a composite of the intercept, the latent means, and the systematic bias in subject p. Because $μ_{p}^{(ε)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\mu }^{(\pmb {\varepsilon })}_{p}$$\end{document} has the same length as $y_{p}^{(ε)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {y}^{(\pmb {\varepsilon })}_{p}$$\end{document} , it is not possible to estimate this bias from time-series data of a single subject, or even from an analysis with a few subjects, even when equality constraints are placed on the intercepts and factor-loading structure (because the latent means may vary, which may be the topic of interest). Therefore, an identifying assumption $μ_{p}^{(ε)} = 0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\mu }^{(\pmb {\varepsilon })}_{p} = \pmb {0}$$\end{document} is required in the ts-lvgvar setting. Note, any within-subject variation is due to variations in deviations from the mean, which is denoted here with $ξ$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\xi }$$\end{document} . As a consequence, the within-subject mean structure takes the following form:

\begin{matrix} μ_{p}^{(y)} = τ + Λ μ_{p}^{(η)} Stationary within-subject means, \end{matrix}

and the variance–covariance structure takes the following form:

(6)

This matrix captures within-subject variation, with k indicating the lag, which may differ person to person.

The relationships between the latent variables may be modeled as a lag-1 GVAR. Assuming stationarity, the structural model then becomes:

\begin{matrix} vec (Σ_{p, 0}^{(η)}) & = {(I \otimes I - B_{p} \otimes B_{p})}^{- 1} vec (Σ_{p, 0}^{(ζ)}) Stationary latent variation. \\ Σ_{p, k}^{(η)} & = B_{p} Σ_{p, k - 1}^{(η)} Stationary lag- k covariances between latent variables . \end{matrix}

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathrm {vec}\left( \pmb {\Sigma }^{(\pmb {\eta })}_{p, 0}\right)&= \left( \pmb {I} \otimes \pmb {I} - \pmb {B}_p \otimes \pmb {B}_p \right) ^{-1} \mathrm {vec}\left( \pmb {\Sigma }^{(\pmb {\zeta })}_{p,0}\right) \quad \text {Stationary latent variation.}\\ \pmb {\Sigma }^{(\pmb {\eta })}_{p, k}&= \pmb {B}_p \pmb {\Sigma }^{(\pmb {\eta })}_{p, k-1} \quad \hbox { Stationary lag-}k \hbox { covariances between latent variables}. \end{aligned}$$\end{document}

The innovation variance–covariance matrix can further be modeled as a GGM to obtain a latent contemporaneous network:

\begin{matrix} Σ_{p}^{(ζ)} & = Δ_{p}^{(ζ)} {(I - Ω_{p}^{(ζ)})}^{- 1} Δ_{p}^{(ζ)} . \end{matrix}

This leads to two network structures: a temporal network modeled with $B_{p}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {B}_p $$\end{document} , and a contemporaneous network modeled with $Ω_{p}^{(ζ)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\Omega }^{(\pmb {\zeta })}_p$$\end{document} . Typical SEM identifying constraints are needed to make the model discernible, such as having nonnegative degrees of freedom, placing $μ_{p}^{(η)} = 0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\mu }^{(\pmb {\eta })}_{p} = \pmb {0}$$\end{document} in a single-subject setting, and constraining either the first factor loadings in $Λ$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\Lambda }$$\end{document} or the diagonal elements of $Δ_{p}^{(ζ)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\Delta }^{(\pmb {\zeta })}_p$$\end{document} to one.

3.1. Estimation

With the exception of modeling contemporaneous relationships as a GGM, the ts-lvgvar takes the form of a dynamic factor model (Molenaar, Reference Molenaar1985, Reference Molenaar2017), which is often estimated by using a Toeplitz matrix (Hamaker, Dolan, & Molenaar, Reference Hamaker, Dolan and Molenaar2002). Let $z_{t}^{⊤} = (\begin{matrix} y_{p, t}^{⊤} & y_{p, t + 1}^{⊤} \end{matrix})$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {z}_{t}^{\top } = \begin{bmatrix}\pmb {y}_{p,t}^{\top }&\pmb {y}_{p,t+1}^{\top } \end{bmatrix}$$\end{document} represent a pair of consecutive responses, and let $Σ = var (z_{T})$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\Sigma } = \mathrm {var}(\pmb {z}_T)$$\end{document} . It follows that $Σ$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\Sigma }$$\end{document} will take the form of a Toeplitz matrix:

\begin{matrix} Σ = (\begin{matrix} Σ^{*} & Σ_{p, 1}^{(y) ⊤} \\ Σ_{p, 1}^{(y)} & Σ_{p, 0}^{(y)} \end{matrix}) . \end{matrix}

It is evident that $Σ^{*} = Σ_{p, 0}^{(y)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\Sigma }^{*} = \pmb {\Sigma }_{p,0}^{(\pmb {y})}$$\end{document} . However, modeling these blocks with the same parameters may lead to a false number of degrees of freedom, as the data will be copied to obtain these blocks. To this end, I keep the model for $Σ^{*}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\Sigma }^{*}$$\end{document} saturated with a unique set of parameters. The same holds for the mean structure:

\begin{matrix} μ = (\begin{matrix} μ^{*} \\ μ^{(y)} \end{matrix}) . \end{matrix}

Before fitting such a model, the data need to be structured in a certain way. In general, an augmented data matrix can be created by copying the data, shifting it by one row, and appending it to the former data matrix. However, it may be warranted to remove certain pairs (e.g., removing effects that occur over night). If someone is measured three times per day for three days, indicating that Day 1 consists of measurements 1, 2, 3; Day 2 of 4, 5, 6; and Day 3 of 7, 8, 9. The data may then be structured as follows:

(7)

\begin{matrix} Z = (\begin{matrix} . & y_{p, 1}^{⊤} \\ y_{p, 1}^{⊤} & y_{p, 2}^{⊤} \\ y_{p, 2}^{⊤} & y_{p, 3}^{⊤} \\ . & y_{p, 4}^{⊤} \\ y_{p, 4}^{⊤} & y_{p, 5}^{⊤} \\ y_{p, 5}^{⊤} & y_{p, 6}^{⊤} \\ . & y_{p, 7}^{⊤} \\ y_{p, 7}^{⊤} & y_{p, 8}^{⊤} \\ y_{p, 8}^{⊤} & y_{p, 9}^{⊤} \end{matrix}) . \end{matrix}

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \pmb {Z} = \begin{bmatrix} . &{}\quad \pmb {y}_{p,1}^{\top } \\ \pmb {y}_{p,1}^{\top } &{}\quad \pmb {y}_{p,2}^{\top } \\ \pmb {y}_{p,2}^{\top } &{}\quad \pmb {y}_{p,3}^{\top } \\ . &{}\quad \pmb {y}_{p,4}^{\top } \\ \pmb {y}_{p,4}^{\top } &{}\quad \pmb {y}_{p,5}^{\top } \\ \pmb {y}_{p,5}^{\top } &{}\quad \pmb {y}_{p,6}^{\top } \\ . &{}\quad \pmb {y}_{p,7}^{\top } \\ \pmb {y}_{p,7}^{\top } &{}\quad \pmb {y}_{p,8}^{\top } \\ \pmb {y}_{p,8}^{\top } &{}\quad \pmb {y}_{p,9}^{\top } \end{bmatrix}. \end{aligned}$$\end{document}

The parameters may now be estimated by optimizing the FIML fit function (2). It should be noted that fitting such a time-series model to the Toeplitz matrix may not be the best course of action because the resulting fit function does not represent the true likelihood of the data. A different method would be to construct one large variance–covariance matrix for the entire vectorized dataset and subsequently evaluate the likelihood by treating this as a single observation (Ciraki, Reference Ciraki2007, p. 90). Ciraki, (Reference Ciraki2007) provides derivatives for these (more) general dynamic SEMs—with the exception that the contemporaneous effects are modeled at the variance–covariance level rather than at the GGM level. This estimation method, however, is computationally very expensive and has therefore not been implemented in the psychonetrics package.

3.2. Empirical Example

To exemplify the ts-lvgvar, I analyzed time-series data previously studied by Wichers, Groot, Psychosystems, Group, and Others (Reference Wichers, Groot, Psychosystems and Lenin2016) and made public by Kossakowski, Groot, Haslbeck, Borsboom, and Wichers (Reference Kossakowski, Groot, Haslbeck, Borsboom and Wichers2017). This dataset concerns a 57-year-old man with a history of major depression, who was measured over the course of 239 days through the experience sampling method (ESM) using a digital device with a touch screen. During the study, the participant reduced the intake of antidepressants and relapsed into a clinical depression (Wichers et al., Reference Wichers, Groot, Psychosystems and Lenin2016). Because of the assumed stationarity in the ts-lvgvar, I selected only the previously unstudied post-assessment phase (days 156 to 239) during which medication levels were no longer changed. Furthermore, I selected 28 items that were continuous and varied substantively across the sample. These items were designed to measure mood states, pathological symptoms, self-esteem, and physical condition. Before analyzing the data, I tested each individual item for a linear trend by regressing the item scores on the time variable and replaced scores for each significant ( $α = 0.05$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha = 0.05$$\end{document} ) regression with the residuals.

Measurement Model Formation Because the original study did not measure items according to a predefined measurement model, I set out to find a factor structure that leads to a limited number of indicators per factor. To this end, I first performed a parallel analysis on the data, which led me to retain five factors. Next, I performed an exploratory factor analysis (EFA) on the dataset using an oblimin rotation. In order to get a comparable number of indicators for each factor and to reduce the number of indicators to a number the software could handle (given decent computation speed), I determined the three indicators with the strongest absolute factor loadings for each of the five factors and discarded all other indicators. This led to the retention of 14 indicatorsFootnote ⁶: “irritated,” “satisfied,” “lonely,” “anxious,” “enthusiastic,” “guilty,” “strong,” “restless,” “agitated,” “worry,” “ashamed,” “tired,” “headache,” and “sleepy,” which were all measured using 7-point scales. Finally, I formed the “confirmatory” modelFootnote ⁷ by making all factor loadings that had a stronger absolute value than 0.25 in the EFA solution free parameters and by constraining all other factor loadings to zero. Both the parallel analysis and EFA were performed using version 1.8.12 of the psych package for R (Revelle, Reference Revelle2018). Figure 2a visualizes the standardized factor loadings of the final estimated ts-lvgvar model (explained below), with the factors reordered to improve readability. The factors can roughly be interpreted as positive affect or positive activation (F1), self-consciousness (F2), anxiety (F3), irritability or negative activation (F3), somatization (F4), and anxiety (F5).

Model Estimation I fitted the ts-lvgvar model using version 0.4 of the psychonetrics package (code available in supplementary materials at https://osf.io/z5hbs/). The augmented data in the final model, as shown in Eq. (7), contained 486 rows of observations with $21.4 %$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$21.4\%$$\end{document} missingness. Using FIML estimation, I first estimated a model in which the latent network structures (temporal and contemporaneous) were fully connected. The residual variances were estimated using a Cholesky decomposition, Footnote ⁸ which ensured that all residual variances were nonnegative. Although exact fit was rejected, $χ^{2} (234) = 447.07, p < 0.001$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\chi ^2(234) = 447.07, p < 0.001$$\end{document} , the model showed adequate close fit. The root-mean-square error of approximation (RMSEA) was 0.043 ( $95 % CI$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$95\%\,\mathrm {CI}$$\end{document} 0.037–0.049), and most incremental fit indices were acceptable ( $NFI = 0.91$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {NFI} = 0.91$$\end{document} , $PNFI = 0.74$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {PNFI} = 0.74$$\end{document} , $TLI = 0.95$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {TLI} = 0.95$$\end{document} , $NNFI = 0.95$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {NNFI} = 0.95$$\end{document} , $RFI = 0.89$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {RFI} = 0.89$$\end{document} , $IFI = 0.96$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {IFI} = 0.96$$\end{document} , $RNI = 0.96$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {RNI} = 0.96$$\end{document} , $CFI = 0.96$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {CFI} = 0.96$$\end{document} ). Next, I fixed to zero all edges from the contemporaneous and temporal network that were not significant at $α = 0.01$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha = 0.01$$\end{document} , and refit the model. Finally, I performed stepwise model search to find a model with optimal Bayesian information criterion (BIC). The details of this algorithm are further explained in Fig. 1. This pruned model did not fit significantly worse than the original model, $Δ χ^{2} (25) = 36.80, p = 0.06$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \chi ^2(25) = 36.80, p = 0.06$$\end{document} , featured a lower AIC ( $Δ AIC = 13.2$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \mathrm {AIC} = 13.2$$\end{document} ) and a lower BIC ( $Δ BIC = 117.85$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \mathrm {BIC} = 117.85$$\end{document} ). The pruned model featured an acceptable RMSEA of 0.042 ( $95 % CI$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$95\%\,\mathrm {CI}$$\end{document} 0.036–0.048), with mostly acceptable incremental fit indices ( $NFI = 0.91$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {NFI} = 0.91$$\end{document} , $PNFI = 0.82$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {PNFI} = 0.82$$\end{document} , $TLI = 0.95$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {TLI} = 0.95$$\end{document} , $NNFI = 0.95$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {NNFI} = 0.95$$\end{document} , $RFI = 0.90$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {RFI} = 0.90$$\end{document} , $IFI = 0.95$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {IFI} = 0.95$$\end{document} , $RNI = 0.95$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {RNI} = 0.95$$\end{document} , $CFI = 0.95$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {CFI} = 0.95$$\end{document} ).

Figure 1. Model search strategy used for given $α$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document} level (here, $α = 0.01$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha = 0.01$$\end{document} was used) in the ts-lvgvar and panel-lvgvar example analyses. For the ts-lvgvar, model selection is performed on the temporal network $B_{p}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {B}_p$$\end{document} and the contemporaneous network $Ω_{p}^{(ζ)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\Omega }^{(\pmb {\zeta })}_p$$\end{document} , and for the panel-lvgvar model selection is performed on the temporal network $B_{*}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \pmb {B}_*$$\end{document} , the contemporaneous network $\underset{within}{Ω^{(ζ)}}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\underset{\mathrm {within}}{\pmb {\Omega }^{(\pmb {\zeta })}} $$\end{document} , and the between-subject network $\underset{between}{Ω^{(ζ)}}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\underset{\mathrm {between}}{\pmb {\Omega }^{(\pmb {\zeta })}}$$\end{document} . This algorithm has been implemented in the modelsearch function in the psychonetrics package.

Figure 2 shows the temporal effects, standardized to partial directed correlations (Wild et al., Reference Wild, Eichler, Friederich, Hartmann, Zipfel and Herzog2010), and the estimated contemporaneous partial correlation network. Table 1 shows the numeric estimates of these networks and the model-implied marginal contemporaneous correlations. In the estimated model, each of the five factors featured strong positive autoregressions, also termed inertia. At the temporal level, deviations in the second factor (self-consciousness) predicted lower levels of the first factor (positive activation) over time. At the contemporaneous level, the first factor played a central role: High levels of positive activation seemed to be uniquely associated with lower levels in all of the other four negative factors.

Figure 2. Results of the ts-lvgvar analysis of the post-assessment data from Kossakowski et al. (Reference Kossakowski, Groot, Haslbeck, Borsboom and Wichers2017). Blue edges indicate positive effects, and ref edges indicate negative effects. The estimated within-subject latent network structures (b and d) were estimated together with the measurement model in a. Panels c and e show how often each edge was included in a 25% case-drop bootstrap (Color figure online).

Table 1. Numeric results of the ts-lvgvar analysis of the post-assessment data from Kossakowski et al. (Reference Kossakowski, Groot, Haslbeck, Borsboom and Wichers2017).

Bootstrap Results To assess the stability of the estimation algorithm, I performed 1000 case-drop bootstraps (Epskamp, Borsboom, & Fried, Reference Epskamp, Borsboom and Fried2017) in which I dropped $25 %$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$25\%$$\end{document} of the data and reestimated the model structure using the search strategy of Fig. 1. The $25 %$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$25\%$$\end{document} rate was chosen in line with previous simulation studies on the use of case-drop bootstraps in cross-sectional networks, which suggests dropping a minimum of $25 %$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$25\%$$\end{document} cases may be needed to assess stability of the network structure. To take temporal dependency into account, I dropped a block of data rather than individual rows from the dataset in each bootstrap sample, starting from a measurement occasion chosen at random. Figures 2e and 3e visualize the number of times that each signed edge would be included in the estimated model, and Table 2 shows the raw number of times each parameter was included (ignoring sign). These results indicate a high level of stability in the inertia parameters (included in 749 to 1000 bootstrap samples) as well as the contemporaneous edges (included in 945 to 1000 bootstrap samples). At the temporal level, it is noteworthy that the temporal edge from F2 to F1 was only included in 290 bootstrap samples, and a temporal edge in the reverse direction (F1 to F2) was included 444 times. Even though both edges could be included in principle, none of the bootstrap samples featured both edges simultaneously included, indicating that in 734 of the bootstrap samples one temporal edge of any direction was included between F1 and F2. In the original sample, models with both edges or only $F 1 \to F 2$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F1 \rightarrow F2$$\end{document} fitted only slightly worse than the model with only $F 2 \to F 1$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F2 \rightarrow F1$$\end{document} in both AIC ( $Δ AIC = - 1.16$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \mathrm {AIC} = -1.16$$\end{document} and $Δ AIC = - 10.08$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \mathrm {AIC} = -10.08$$\end{document} , respectively) and BIC ( $Δ BIC = - 5.35$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \mathrm {BIC} = -5.35$$\end{document} and $Δ BIC = - 10.08$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \mathrm {BIC} = -10.08$$\end{document} , respectively). Thus, there is some evidence that the relationship between F1 and F2 may be reciprocal. Noteworthy at the contemporaneous level is that some negative partial correlations were included in bootstrap samples where one would expect positive effects: negative relationships between F2 and F3 were included 333 times and negative relationships between F2 and F5 were included 110 times. These negative edges are noteworthy because these would be expected when F1 is a common effect of the other factors (Epskamp, Waldorp, et al., Reference Epskamp, Waldorp, Mõttus and Borsboom2018). However, because these were only included sparingly, they cannot be interpreted substantively.

Table 2. The number of times each parameter was included in the case-drop bootstrap ts-lvgvar analysis of the post-assessment data from Kossakowski et al. (Reference Kossakowski, Groot, Haslbeck, Borsboom and Wichers2017).

Each replication (100 in total) was based on a 75% subsample of the original dataset. Bold-faced values indicate parameters that were included in the original analysis.

4. Panel Data: The panel-lvgvar Model

The panel-lvgvar concerns the variation across random people in a few fixed measurement occasions $t, t + 1, \dots$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$t, t+1,\ldots $$\end{document} (i.e., the distributions over people).

\begin{matrix} y_{P, t} & \sim N (μ_{t}^{(y)}, Σ_{t}^{(y)}) & Cross-sectional distribution of observed variables. \\ η_{P, t} & \sim N (μ_{t}^{(η)}, Σ_{t}^{(η)}) & Cross-sectional distribution of latent variables, \end{matrix}

Stationarity is assumed in time, and therefore, it is assumed that these structures do not change over time. Hence, subscript t in the mean vectors and in the variance–covariance matrices above can be replaced with the symbol $†$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\dagger $$\end{document} :

\begin{matrix} μ_{t}^{(\dots)} = μ_{†}^{(\dots)}, Σ_{t}^{(\dots)} = Σ_{†}^{(\dots)} \forall_{t} Equal means and (co)variances for all time points. \end{matrix}

From the measurement model in Eq. (4), the expected value becomes:

\begin{matrix} μ_{†}^{(y)} = τ + Λ μ_{†}^{(η)} Expected cross-sectional observed scores . \end{matrix}

As for any fixed time point t, $cov (μ_{P}^{(y)}, ξ_{P, t}^{(y)}) = O$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm {cov}\left( \pmb {\mu }^{(\pmb {y})}_{P}, \pmb {\xi }^{(\pmb {y})}_{P,t}\right) = \pmb {O}$$\end{document} by definition, the lag-k variance–covariance structure can be derived as:

(8)

\begin{matrix} Σ_{k}^{(y)} = var (μ_{P}^{(y)}) + cov (ξ_{P, † + k}^{(y)}, ξ_{P, †}^{(y)}) . \end{matrix}

Note that $var (μ_{P}^{(y)})$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm {var}\left( \pmb {\mu }^{(\pmb {y})}_{P}\right) $$\end{document} relates to purely between-subject variation. The between-subject variance–covariance structure can therefore be defined as follows:

(9)

\begin{matrix} \underset{between}{Σ^{(η)}} & = var (μ_{P}^{(η)}) & between-subject latent (co)variances, \\ \underset{between}{Σ^{(ε)}} & = var (μ_{P}^{(ε)}) & between-subject residual (co)variances, \\ \underset{between}{Σ^{(y)}} & = var (μ_{P}^{(y)}) = Λ \underset{between}{Σ^{(η)}} Λ^{⊤} + \underset{between}{Σ^{(ε)}} & between-subject (co)variances . \end{matrix}

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \underset{\mathrm {between}}{\pmb {\Sigma }^{(\pmb {\eta })}}&= \mathrm {var}\left( \pmb {\mu }^{(\pmb {\eta })}_{P}\right)&\text {between-subject latent (co)variances,} \nonumber \\ \underset{\mathrm {between}}{\pmb {\Sigma }^{(\pmb {\varepsilon })}}&= \mathrm {var}\left( \pmb {\mu }^{(\pmb {\varepsilon })}_{P}\right)&\text {between-subject residual (co)variances,} \nonumber \\ \underset{\mathrm {between}}{\pmb {\Sigma }^{(\pmb {y})}}&= \mathrm {var}\left( \pmb {\mu }^{(\pmb {y})}_{P}\right) = \pmb {\Lambda } \underset{\mathrm {between}}{\pmb {\Sigma }^{(\pmb {\eta })}} \pmb {\Lambda }^{\top } + \underset{\mathrm {between}}{\pmb {\Sigma }^{(\pmb {\varepsilon })}}&\text {between-subject (co)variances}. \end{aligned}$$\end{document}

This matrix equals the first part of Eq. (8). As for the second part, by virtue of the law of total expectation, the fixed-effect latent structure takes the following form:

(10)

\begin{matrix} \underset{within}{Σ_{*, k}^{(y)}} = E_{P} ({cov}_{T}, (ξ_{P, T + k}^{(y)}, ξ_{P, T}^{(y)} ∣ P)) = E_{T} ({cov}_{P}, (ξ_{P, T + k}^{(y)}, ξ_{P, T}^{(y)} ∣ T)) = {cov}_{P} (ξ_{P, † + k}^{(y)}, ξ_{P, †}^{(y)}), \end{matrix}

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \underset{\mathrm {within}}{\pmb {\Sigma }_{*,k}^{(\pmb {y})}} = \mathcal {E}_P\left( \mathrm {cov}_T\left( \pmb {\xi }^{(\pmb {y})}_{P,T+k}, \pmb {\xi }^{(\pmb {y})}_{P,T} \mid P \right) \right) = \mathcal {E}_T\left( \mathrm {cov}_P\left( \pmb {\xi }^{(\pmb {y})}_{P,T+k}, \pmb {\xi }^{(\pmb {y})}_{P,T} \mid T \right) \right) = \mathrm {cov}_P\left( \pmb {\xi }^{(\pmb {y})}_{P,\dagger +k}, \pmb {\xi }^{(\pmb {y})}_{P,\dagger }\right) , \end{aligned}$$\end{document}

in which subscripts T and P are added to the expectation and covariance operators to make clear with respect to which aspect the expectation or covariance is taken. From left to right, $E_{P} ({cov}_{T}, (ξ_{P, T + k}^{(y)}, ξ_{P, T}^{(y)} ∣ P))$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {E}_P\left( \mathrm {cov}_T\left( \pmb {\xi }^{(\pmb {y})}_{P,T+k}, \pmb {\xi }^{(\pmb {y})}_{P,T} \mid P \right) \right) $$\end{document} indicates the fixed-effect structure of the within-subject deviations—the within-subject structure of the average subject. The next expression, $E_{T} ({cov}_{P}, (ξ_{P, T + k}^{(y)}, ξ_{P, T}^{(y)} ∣ T))$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {E}_T\left( \mathrm {cov}_P\left( \pmb {\xi }^{(\pmb {y})}_{P,T+k}, \pmb {\xi }^{(\pmb {y})}_{P,T} \mid T \right) \right) $$\end{document} indicates the expected between-subject variation over all time points, which due to the assumption of stationarity reduces to the final part of the expression: ${cov}_{P} (ξ_{P, † + k}^{(y)}, ξ_{P, †}^{(y)})$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm {cov}_P\left( \pmb {\xi }^{(\pmb {y})}_{P,\dagger +k}, \pmb {\xi }^{(\pmb {y})}_{P,\dagger }\right) $$\end{document} , which is also present in Eq. (8). This indicates that under the assumption of homogeneity—every subject follows the same within-subject process—the second part of Eq. (8) will equal the within-subject variance–covariance structure for every subject. If homogeneity is violated, the structure can instead be interpreted as the within-subject variance–covariance structure of the average subject.

Equations (9) and (10) may be used to write (8) as follows:

(11)

which forms the basis for estimating the panel-lvgvar model. The within-subject part can further be modeled as a latent variable model using the fixed-effect variance–covariance matrices for the latent variables and residuals:

\begin{matrix} \underset{within}{Σ_{*, k}^{(y)}} = (\begin{matrix} Λ \underset{within}{Σ_{*, k}^{(η)}} Λ^{⊤} + \underset{within}{Σ_{*, k}^{(ε)}} & if k = 0 \\ Λ \underset{within}{Σ_{*, k}^{(η)}} Λ^{⊤} \end{matrix}) . \end{matrix}

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \underset{\mathrm {within}}{\pmb {\Sigma }_{*,k}^{(\pmb {y})}} = {\left\{ \begin{array}{ll} \pmb {\Lambda } \underset{\mathrm {within}}{\pmb {\Sigma }_{*,k}^{(\pmb {\eta })}} \pmb {\Lambda }^{\top } + \underset{\mathrm {within}}{\pmb {\Sigma }_{*,k}^{(\pmb {\varepsilon })}} &{}\text{ if } k = 0 \nonumber \\ \pmb {\Lambda } \underset{\mathrm {within}}{\pmb {\Sigma }_{*,k}^{(\pmb {\eta })}} \pmb {\Lambda }^{\top } &{} \end{array}\right. }. \end{aligned}$$\end{document}

The structural model becomes:

\begin{matrix} vec (\underset{within}{Σ_{*, 0}^{(η)}}) & = {(I \otimes I - B_{*} \otimes B_{*})}^{- 1} vec (\underset{within}{Σ_{*, 0}^{(ζ)}}), \\ \underset{within}{Σ_{*, k}^{(η)}} & = B_{*} \underset{within}{Σ_{*, k - 1}^{(η)}}, \end{matrix}

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathrm {vec}\left( \underset{\mathrm {within}}{\pmb {\Sigma }_{*,0}^{(\pmb {\eta })}}\right)&= \left( \pmb {I} \otimes \pmb {I} - \pmb {B}_* \otimes \pmb {B}_* \right) ^{-1} \mathrm {vec}\left( \underset{\mathrm {within}}{\pmb {\Sigma }_{*,0}^{(\pmb {\zeta })}}\right) , \\ \underset{\mathrm {within}}{\pmb {\Sigma }_{*,k}^{(\pmb {\eta })}}&= \pmb {B}_* \underset{\mathrm {within}}{\pmb {\Sigma }_{*,k-1}^{(\pmb {\eta })}}, \end{aligned}$$\end{document}

in which $B_{*}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {B}_*$$\end{document} represents a matrix of fixed-effect temporal effects and $Σ_{*, k}^{(ζ)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\Sigma }_{*,k }^{(\pmb {\zeta })}$$\end{document} the lag-k fixed-effect contemporaneous variance–covariance structure. Finally, the two latent variance–covariance structures in the model could further be modeled as GGMs:

\begin{matrix} \underset{within}{Σ_{*}^{(ζ)}} & = \underset{within}{Δ^{(ζ)}} {(I - \underset{within}{Ω^{(ζ)}})}^{- 1} \underset{within}{Δ^{(ζ)}} & Within-subject latent contemporaneous network . \\ \underset{between}{Σ^{(η)}} & = \underset{between}{Δ^{(η)}} {(I - \underset{between}{Ω^{(η)}})}^{- 1} \underset{between}{Δ^{(η)}} & Between-subject latent network . \end{matrix}

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \underset{\mathrm {within}}{\pmb {\Sigma }_{*}^{(\pmb {\zeta })}}&= \underset{\mathrm {within}}{\pmb {\Delta }^{(\pmb {\zeta })}} \left( \pmb {I} - \underset{\mathrm {within}}{\pmb {\Omega }^{(\pmb {\zeta })}} \right) ^{-1} \underset{\mathrm {within}}{\pmb {\Delta }^{(\pmb {\zeta })}}&\text {Within-subject latent contemporaneous network}.\\ \underset{\mathrm {between}}{\pmb {\Sigma }^{(\pmb {\eta })}}&= \underset{\mathrm {between}}{\pmb {\Delta }^{(\pmb {\eta })}} \left( \pmb {I} - \underset{\mathrm {between}}{\pmb {\Omega }^{(\pmb {\eta })}} \right) ^{-1} \underset{\mathrm {between}}{\pmb {\Delta }^{(\pmb {\eta })}}&\text {Between-subject latent network}. \end{aligned}$$\end{document}

This leads to three network structures: (1) the (fixed-effect) within-subject temporal network modeled with $B_{*}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {B}_*$$\end{document} , (2) the (fixed-effect) within-subject latent contemporaneous network modeled with $\underset{within}{Ω^{(ζ)}}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\underset{\mathrm {within}}{\pmb {\Omega }^{(\pmb {\zeta })}}$$\end{document} , and (3) the between-subject latent network modeled with $\underset{between}{Ω^{(η)}}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\underset{\mathrm {between}}{\pmb {\Omega }^{(\pmb {\eta })}}$$\end{document} .

4.1. Estimation

Estimation of the panel-lvgvar in panel data with $n_{t}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n_t$$\end{document} measurements can be done by forming the vector:

\begin{matrix} z_{p} = (\begin{matrix} y_{p, 1} \\ y_{p, 2} \\ ⋮ \\ y_{p, n_{t}} \end{matrix}), \end{matrix}

such that the mean structure repeats the stationary mean:

\begin{matrix} μ = (\begin{matrix} μ^{(y)} \\ μ^{(y)} \\ ⋮ \\ μ^{(y)} \end{matrix}) . \end{matrix}

The variance–covariance structure then takes the form of a block Toeplitz matrix:

\begin{matrix} Σ = (\begin{matrix} Σ_{†, 0}^{(y)} & Σ_{†, 1}^{(y) ⊤} & \dots & Σ_{†, n_{t} - 1}^{(y) ⊤} \\ Σ_{†, 1}^{(y)} & Σ_{†, 0}^{(y)} & \dots & Σ_{†, n_{t} - 2}^{(y) ⊤} \\ ⋮ & ⋮ & ⋱ & ⋮ \\ Σ_{†, n_{t} - 1}^{(y)} & Σ_{†, n_{t} - 2}^{(y)} & \dots & Σ_{†, 0}^{(y)} . \end{matrix}) . \end{matrix}

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \pmb {\Sigma } = \begin{bmatrix} \pmb {\Sigma }_{\dagger ,0}^{(\pmb {y})} &{} \pmb {\Sigma }_{\dagger ,1}^{(\pmb {y})\top } &{} \cdots &{} \pmb {\Sigma }_{\dagger ,n_t - 1}^{(\pmb {y})\top } \\ \pmb {\Sigma }_{\dagger ,1}^{(\pmb {y})} &{} \pmb {\Sigma }_{\dagger ,0}^{(\pmb {y})} &{} \cdots &{} \pmb {\Sigma }_{\dagger ,n_t - 2}^{(\pmb {y})\top } \\ \vdots &{} \vdots &{} \ddots &{} \vdots \\ \pmb {\Sigma }_{\dagger ,n_t - 1}^{(\pmb {y})} &{} \pmb {\Sigma }_{\dagger ,n_t - 2}^{(\pmb {y})} &{} \cdots &{} \pmb {\Sigma }_{\dagger ,0}^{(\pmb {y})}. \end{bmatrix}. \end{aligned}$$\end{document}

The model can then be estimated by minimizing Eq. (1) or Eq. (2). When not all variables are observed at all waves (or entire waves of data are missing), elements in $μ$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\mu }$$\end{document} and rows and columns in $Σ$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\Sigma }$$\end{document} corresponding to missing variables can be cut out, or FIML estimation can be used when data are treated as missing.

4.2. Empirical Example

To exemplify the panel-lvgvar, I investigated a subset from the Longitudinal Internet Studies for the Social Sciences (LISS; Scherpenzeel & Das, Reference Scherpenzeel, Das, Das, Ester and Kaczmirek2010) panel administered by CentERdata (Tilburg University, The Netherlands). The LISS panel is a representative sample of Dutch individuals who participate in a yearly survey on several factors.Footnote ⁹ To minimize the amount of missing data and to exemplify the usage of the method on as little as three waves of data and the usage of the method in waves with unequal time differences between consecutive waves, I used three of the waves from the LISS Core Study on personality: 2013, 2014, and 2017. Table 3 shows the items I selected to measure five factors: self-esteem, pessimism, optimism, life satisfaction, positive affect, and negative affect. Pessimism and optimism were assessed with the LOT-R scale (Carver, Scheier, & Segerstrom, Reference Carver, Scheier and Segerstrom2010), and life satisfaction was assessed with the Satisfaction With Life scale (Diener, Emmons, Larsen, & Griffin, Reference Diener, Emmons, Larsen and Griffin1985). Because the current software does not handle more than about 30 variables, I used shorter scales than typical to assess self-esteem, positive affect, and negative affect. Self-esteem was assessed with a three-item scale accredited to “Radboud University Nijmegen, Netherlands” in the codebook, rather than the longer Rosenberg scale also present in the data (Rosenberg, Reference Rosenberg1965). Positive and negative affects were assessed with the PANAS scale (Watson, Clark, & Tellegen, Reference Watson, Clark and Tellegen1988), which was shortened to four indicators per factor. To this end, I first fit panel-lvgvar models for positive and negative affects separately and retained the four items that featured the highest absolute factor loadings for both factors. Finally, to make the analysis fully reproducible and to exemplify the estimation of the panel-lvgvar from summary statistics, I only retained cases with no missing data. Thus, the final sample included 2, 998 cases administered three times on 22 items per wave. The variance–covariance matrix and mean vector of the data, as well as the code to reproduce the main analysis, are available in the supplementary.

Table 3. Measurement model for LISS data example.

Model Estimation I fit the panel-lvgvar model using version 0.4 of the psychonetrics package (code available in supplementary materials at https://osf.io/z5hbs/). First, I fit a model in which all edges were included in the temporal, contemporaneous, and between-subject networks and estimated all residual structures using a Cholesky decomposition. The model featured a numeric estimate of approximately zero for the between-subject residual variance for the item SE2.Footnote ¹⁰ I fixed this parameter to zero (note, the within-subject residual was not constrained to zero, so the item did not become a perfect indicator as a result) and reestimated the model. Although exact fit was rejected, $χ^{2} (2118) = 5686.50, p < 0.001$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\chi ^2(2118) = 5686.50, p < 0.001$$\end{document} , the model showed excellent close fit with an RMSEA of 0.024 ( $95 % CI$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$95\%\,\mathrm {CI}$$\end{document} 0.023–0.024) and good incremental fit ( $NFI = 0.95$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {NFI} = 0.95$$\end{document} , $PNFI = 0.91$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {PNFI} = 0.91$$\end{document} , $TLI = 0.97$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {TLI} = 0.97$$\end{document} , $NNFI = 0.97$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {NNFI} = 0.97$$\end{document} , $RFI = 0.95$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {RFI} = 0.95$$\end{document} , $IFI = 0.97$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {IFI} = 0.97$$\end{document} , $RNI = 0.97$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {RNI} = 0.97$$\end{document} , $CFI = 0.97$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {CFI} = 0.97$$\end{document} ). Next, I performed the same model search strategy used above in the ts-lvgvar example and further detailed in Fig. 1. While the pruned model fitted significantly worse than the original model, $Δ χ^{2} (36) = 64.15, p < 0.01$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \chi ^2(36) = 64.15, p < 0.01$$\end{document} , it was better in terms of AIC ( $Δ A I C = 7.8$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta AIC = 7.8$$\end{document} ) and BIC ( $Δ B I C = 224.05$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta BIC = 224.05$$\end{document} ). The pruned model also showed excellent fit, with an RMSEA of 0.024 ( $95 % CI$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$95\%\,\mathrm {CI}$$\end{document} 0.023–0.024) and good incremental fit ( $NFI = 0.95$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {NFI} = 0.95$$\end{document} , $PNFI = 0.93$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {PNFI} = 0.93$$\end{document} , $TLI = 0.97$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {TLI} = 0.97$$\end{document} , $NNFI = 0.97$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {NNFI} = 0.97$$\end{document} , $RFI = 0.95$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {RFI} = 0.95$$\end{document} , $IFI = 0.97$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {IFI} = 0.97$$\end{document} , $RNI = 0.97$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {RNI} = 0.97$$\end{document} , $CFI = 0.97$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {CFI} = 0.97$$\end{document} ).

Table 4. Numeric results of the panel-lvgvar analysis of the LISS Core Study on personality.

Figure 3 portrays the estimated factor loadings and network structures in the pruned model, and Table 4 shows the numeric estimates of the standardized network parameters. The temporal network mainly features pathways from life satisfaction to self-esteem and optimism, and subsequently to pessimism and negative affect. At the contemporaneous and between-subject levels, optimism seems to play a central role in the networks—being connected to all other variables in both networks with consistently strong edges. Most noteworthy in both the contemporaneous and between-subject networks are the rather strong positive edges between positive affect and negative affect. Table 4 shows that the model-implied correlations in these two levels are rather weak and also positive at the contemporaneous level ( $r = 0.05$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r = 0.05$$\end{document} at the contemporaneous level and $r = - 0.11$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r = -0.11$$\end{document} at the between-subject level). These weak marginal correlations are in line with previous literature (e.g., Tuccitto, Giacobbi Jr, & Leite, Reference Tuccitto, Giacobbi and Leite2010), suggesting that the correlation between positive and negative affects might not be as strongly negative as expected and that the factors may be (near) orthogonal instead. At the partial correlation levels, these weak marginal correlations lead to stronger partial correlations of an unexpected sign because the correlation is weaker than can be expected due to the links between positive/negative affect and third variables, such as self-esteem and optimism. It may even be that these variables act as a common effect between positive and negative affects (De Ron, Fried, & Epskamp, Reference De Ron, Fried and Epskamp2019; Epskamp, Waldorp, et al., Reference Epskamp, Waldorp, Mõttus and Borsboom2018), in which case a strong partial correlation of an unexpected sign may also be expected. To check whether the results were influenced by the subset of indicators chosen to assess positive affect and negative affect, I also estimated a panel-lvgvar model for only positive and negative affects using all indicators. These showed similarly weak marginal correlations ( $r = 0.003$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r = 0.003$$\end{document} at the contemporaneous level and $r = - 0.056$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r = -0.056$$\end{document} at the between-subject level).

Bootstrap Results To assess the stability of the estimation algorithm, I again analyzed a 1000 case-drop bootstrap samples on the final panel-lvgvar model. In each of these samples, $25 %$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$25\%$$\end{document} of the participants were dropped at random from the original sample, and panel-lvgvar models were estimated using the algorithm described in Fig. 1. The number of times each edge was included is visually shown in Fig. 3 and numerically shown in Table 5. Across all three networks, edges that were included in the original analyses were also likely to be included in the case-drop bootstrapped analyses, and mostly edges that were not included in the original analysis were also less likely to be included in the case-drop bootstrapped analyses. This indicated a high level of stability in the estimated network structures, in particular in the between-subject network. In the contemporaneous network, the edge between negative affect and optimism was estimated only 543 times and an extra edge between positive affect and life satisfaction was included 99 times. The temporal network featured lower levels of stability: several edges that were included in the original sample were included less than 400 times in the bootstrap samples, and some additional edges (most noteworthy self-esteem to negative affect, negative affect to pessimism, self-esteem to pessimism, and self-loops on pessimism and optimism) were included over 200 times. As such, there is some evidence that the temporal network is less sparse than portrayed in the original analysis.

Figure 3. Results of the panel-lvgvar analysis of the LISS Core Study on personality. Blue edges indicate positive effects, and ref edges indicate negative effects. The estimated fixed-effect within-subject latent networks are shown in b and c and the estimated between-subject network in d. These are estimated jointly with the measurement model parameters represented in a. Panels e, f and g show the inclusion proportion of each edge in a 25% case-drop bootstrap (Color figure online).

Table 5. The number of times each parameter was included in the case-drop bootstrap panel-lvgvar analysis of the LISS Core Study on personality.

Each replication (1000 in total) was based on a 75% subsample of the original dataset. Bold-faced values indicate parameters that were included in the original analysis.

5. Simulation Studies

To assess the performance of model search strategies in ts-lvgvar and panel-lvgvar, I performed a two simulation studies by simulating new datasets using the parameters from the models shown in Figs. 2 and 3, respectively. I varied the number of cases (time points) between 50, 100, 250, 500, and 1000 for the ts-lvgvar simulation study and (subjects) between 500, 1000, 2500, 5000, and 10, 000 for the panel-lvgvar simulations. After simulating data, I estimated the model structure using the algorithm presented in Fig. 1 and two variants: a variant in which no model search is performed after the initial pruning step 2, and a variant in which the model search from step 3 onwards was replaced with a stepwise up model search strategy in which the edge with the strongest modification index was added to the model until BIC was no longer improved (implemented in the stepup function in psychonetrics). In addition, I tested four variants of step 2 Fig. 1: in addition to initially removing edges at an $α$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document} level, I also removed edges in step 2 using Bonferroni (Bland & Altman, Reference Bland and Altman1995), Holm (Holm, Reference Holm1979), and false discovery rate (FDR; Benjamini & Hochberg, Reference Benjamini and Hochberg1995) adjustments. Finally, I varied $α$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document} between 0.01 and 0.05 in the entire algorithm. This led to 24 different estimation algorithms. Each condition was replicated 100 times, leading to 12, 000 simulated datasets in both the ts-lvgvar and panel-lvgvar simulation studies. In these simulations, I investigated the success rate of estimation. (Some estimations failed due to numerical optimization issues such as non-positive definite matrices.) In the successfully estimated models, I investigated the sensitivity, specificity, and the correlation between original (true in the simulation) and estimated parameters (Epskamp & Fried, Reference Epskamp and Fried2018).

Figure 4 shows the results of the simulation studies. All estimation algorithms were conservative at every sample size (high specificity), and increased in sensitivity and correlation as sample size increased. Generally, results using $α = 0.05$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha = 0.05$$\end{document} were a bit less conservative and a bit more sensitive than results using $α = 0.01$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha = 0.01$$\end{document} . While very conservative, estimators that only used adjusted significance level in a pruning step without subsequent stepwise model search performed poor in terms of sensitivity, while not performing much better than estimators followed by stepwise model search strategies in terms of specificity. The estimator used in empirical examples above (sig_0.01 -> modelsearch) is highlighted in Fig. 4 and performed well on all measures. It should be noted, however, that the estimators using only significance pruning were much faster than the “stepup” estimator and even more faster than the “modelsearch” estimators: in the ts-lvgvar: estimating a model with extensive model search on average took 10.19 minutes in ts-lvgvar simulations and 35.15 minutes in panel-lvgvar simulations. Finally, the relatively large number of failed estimations is noteworthy, especially at lower sample sizes. Often, estimation failed due to numerical reasons leading to non-positive definite information matrices (e.g., due to non-positive definite variance–covariance structures). Future improvements to the psychonetrics package and other software methods implementing the ts-lvgvar and panel-lvgvar frameworks could improve computation speed of the model (e.g., by relying more on compiled programming languages) as well as the rate of successful model computations (e.g., by using more sophisticated starting values and/or more constrains on the parameter space).

Figure 4. Simulation results for ts-lvgvar and panel-lvgvar model estimation algorithms implemented in psychonetrics. Data were generated under parameters from Fig. 2 for the ts-lvgvar simulations and under parameters from Fig. 3 for the panel-lvgvar simulations, and each condition was replicated 100 times. The lines display the average values over these replications. “Sensitivity” denotes the proportion of true edges that were also included in the estimated model, “specificity” denotes the proportion of true missing edges that were also not included in the estimated model, “correlation” denotes the Pearson correlation between true and estimated edge weights, and “success” denotes the proportion of models estimations that ran without errors. The thick highlighted line represents the algorithm used in the empirical examples in this paper.

6. Ergodicity and Generalizability

The derivations in this paper can be used to show some important relationships between within- and between-subject effects as well as the correspondences between models in which time is random (time-series data) and models in which the subject is random (cross-sectional and panel designs). In particular, the expression in Eq. (11) has several important implications. First, it shows that the variance–covariance matrix of a cross-sectional set equates to a blend of within- and between-subject variances (Hamaker, Reference Hamaker, Mehl and Conner2012). As such, from a single cross-sectional study, it is impossible to determine whether two variables covary due to trait-level covariation between the means (between-subject) or due to state-level covariation between deviations from the mean (within-subject). Second, when investigating covariation over time between consecutive measurement occasions, nonzero covariances are not indicative per se of within-subject Granger causality (prediction over time) because nonzero covariances over time may also result due to between-subject effects. This is the core argument of Hamaker et al. (Reference Hamaker, Kuiper and Grasman2015), who warn against a within-subject interpretation of the cross-lagged panel model unless a random intercept (here modeled with $μ_{p}^{(η)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\mu }^{(\pmb {\eta })}_{p}$$\end{document} and $μ_{p}^{(ε)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\mu }^{(\pmb {\varepsilon })}_{p}$$\end{document} ) is included in the model.

Assuming full homogeneity in mean and variance–covariance structure—every subject has the same mean, and deviations in every subject from that mean have the same variance–covariance structure—Eq. (9) demonstrates that the between-subject covariance structure would become a matrix of zeroes (there are no differences in means), whereas Eq. (10) demonstrates that the within-subject fixed-effect structure would equal the within-subject variance–covariance structure for every subject. Now, when adding to this the assumption of stationarity, such that Eq. (11) is valid for any fixed measurement occasion t, the cross-sectional variance–covariance structure may equal the within-subject variance–covariance structure for every subject. When all variables are assumed normally distributed, such that no higher-order moments are present, these conditions then satisfy the assumption of ergodicity (Molenaar, Reference Molenaar2004), which is frequently cited to be a requirement for generalizing cross-sectional studies to the within-subject level (e.g., Fisher, Medaglia, & Jeronimus, Reference Fisher, Medaglia and Jeronimus2018).

Although the terms within-subject and between-subject—or alternatively within-person and between-persons—are commonly used in the literature and in this paper, they may lead to a faulty conclusion that between-subject effects necessarily do not occur within a person. The separation of within- and between-subject effects in Eq. (11) hinges on the assumption of stationarity, which is necessarily false in humans with a fixed life span. As such, the assumption may be better phrased as local stationarity, indicating that stationarity is assumed in the period of measurement only. Consider the effect between urbanicity and cannabis use (Isvoranu, Borsboom, van Os, & Guloksuz, Reference Isvoranu, Borsboom, van Os and Guloksuz2016). If these two variables are measured on three consecutive days, then urbanicity is not likely to vary within subject, leading to a between-subject interpretation of the correlation between the two variables. However, if these two variables are instead measured once a year for three consecutive years, within-subject covariation may be found (e.g., when a person moves to a city in which cannabis is more readily available, that person may use more cannabis compared to his or her average). As such, there is some gray area among within- and between-subject effects: When measurements are too close together, $\underset{between}{Σ^{(η)}}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\underset{\mathrm {between}}{\pmb {\Sigma }^{(\pmb {\eta })}}$$\end{document} may itself be interpreted as a mixture between slow within-subject effects and true between-subject effects. Likewise, if a particular questionnaire contained trait-level questions (e.g., “in the last year, did you smoke Cannabis”), $\underset{within}{Σ_{*, k}^{(η)}}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\underset{\mathrm {within}}{\pmb {\Sigma }_{*,k}^{(\pmb {\eta })}}$$\end{document} may not accurately reflect valid daily fluctuations in the variables of interest and may be better interpreted as autocorrelated measurement error.

Note that in this paper the factor-loading structure does not differ between the within-subject and the between-subject levels. Typical multi-level software will allow for the specification of different factor-loading structures on both levels. However, when the factor-loading structure is not equal between the two levels, it is highly questionable whether the same factor is measured. For example, a goal of a study might be to obtain factor scores for a set of subjects at a given time point. If the within- and between-subject factor-loading structure is not the same, then these cannot be combined into factor scores, and only estimates of the means per subject can be obtained. The within-subject factor structure then reduces to correlated residuals over time. To be clear, it is not required that factors vary on both within- and between-subject levels: A factor-loading structure can easily be constrained equally over both levels, whereas a factor itself does not vary within one of the levels (e.g., genetic information). In principle, measurement invariance tests can be performed across people to test whether the factor structure is the same (Adolf et al., Reference Adolf, Schuurman, Borkenau, Borsboom and Dolan2014).

7. Discussion

I presented a general framework in which undirected network models can be used to model relationships between latent variables in both time-series data of a single subject (ts-lvgvar) and in panel data of many subjects measured on at least three occasions (panel-lvgvar). This presents a unification of network psychometrics, such that (dynamic) SEMs that include relationships over time at a latent level are combined with undirected network models. Two empirical examples, one using time-series data and one using panel data, as well as two large-scale simulation studies showed that network structures can reliably be estimated from data. These methods have been implemented in the software package psychonetrics, which is exemplified in the analysis code provided in the supplementary materials. I conclude the paper with a discussion on ergodicity and generalizability, discussing that although within-subject effects may in principle be separated from between-subject effects, the interpretation of these results rests on the intensity and intervals between measurements as well as the plausibility of the assumption of stationarity.

Limitations Several limitations that are often raised for discrete time-series models (e.g., Epskamp, Waldorp, et al., Reference Epskamp, Waldorp, Mõttus and Borsboom2018) may also be raised for the current work: (1) both the ts-lvgvar and the panel-lvgvar crucially rely on the assumptions of (local) stationarity and normality, which are unlikely to hold in practice. In principle, the assumption of stationarity is checked by assessing the model fit of the panel-lvgvar data (because it imposes unique constraints across all waves of measurement) but not by assessing the model fit of the ts-lvgvar when the data are fit to a lag-1 Toeplitz matrix as discussed in this paper. (2) Although methods exist to check for (multivariate) normality, it is not yet clear how to proceed when normality is violated because normality can be violated in many ways, and violations of normality have not yet been investigated in detail in the estimation of GGMs. To this end, it is vital to supplement any analyses with data-driven bootstrap methods, as done here, for example to check for the stability of the results by taking subsamples of the data. (3) The choice of lag interval and the duration of the study is not trivial. This choice will crucially impact the interpretation of the results. For example, although the provided empirical examples from both the ts-lvgvar and panel-lvgvar featured temporal networks, the interpretation of these networks is vastly different because the ts-lvgvar example concerned a lag interval of a few hours, whereas the panel-lvgvar example concerned the lag interval of one year. A further limitation is that only lag-1 relationships were used in the modeling framework, which is the minimum lag required to handle temporal dependency in the data.

In addition, the introduction of latent variables in these network models raises questions of measurement. In the empirical examples of the current paper, the measurement models used are far from confirmatory. The aim of this paper is to present methods for confirmatory fit and exploratory estimation of network structures at a latent level, while the measurement model is assumed known. To exemplify the method, however, I used a highly exploratory routine to obtain a measurement mode in the ts-lvgvar example. Other exploratory tools to assess the measurement structure could have been used, such as exploratory graph analysis (Golino & Epskamp, Reference Golino and Epskamp2017), which is based on clustering in estimated GGMs, or P-factor technique (Zevon & Tellegen, Reference Zevon and Tellegen1982). Such methods may well lead to alternative factor structures and as a result different latent network structures. Furthermore, certain choices made in the panel-lvgvar model, such as the number of indicators for positive and negative affects, may have impacted the results. Ideally, these methods would be used on data that feature clear and previously studied measurement models (Flake & Fried, Reference Flake and Fried2019).

Important to note is that while this paper extensively discussed issues pertaining to ergodicity and the ability to draw within-person conclusions from designs based on models in which the subject is random (cross-sectional and panel designs), the presented methods do not solve these issues. While I show that the temporal and contemporaneous structures from the panel-lvgvar directly relate to fixed-effect temporal and contemporaneous structures that may be obtained from time-series analysis, it may be stressed that such fixed effects can only say something about average effects in a population. Only when there is full homogeneity in within-person variance–covariance structures (an unlikely assumption) do these fixed-effect structures truly relate exactly to the within-person structures of every individual subject.

Related Work The presented modeling framework is closely related to several other frameworks: (1) The ts-lvgvar reduces to a general GVAR model (Epskamp, Waldorp, et al., Reference Epskamp, Waldorp, Mõttus and Borsboom2018) when all variables in the network are observed—the factor-loading matrix then reduces to an identity matrix and the residual variances are set to zero. In this case, the model reduces further to a general vector-autoregression model if the contemporaneous level is modeled as a variance–covariance matrix rather than a GGM. (2) When latent variables are used in the ts-lvgvar, but the contemporaneous level is not modeled using a GGM, the model reduces to a type of dynamic factor model (Molenaar, Reference Molenaar1985). (3) The panel-lvgvar can be considered a multi-level GVAR model (Epskamp, Waldorp, et al., Reference Epskamp, Waldorp, Mõttus and Borsboom2018) if all variables in the networks are observed without measurement error—with the exception that all network parameters are fixed (only intercepts are random). Extending the panel-lvgvar framework to include random effects on the network parameters will likely not be trivial, but may form an important direction of future research. (4) When the contemporaneous and between-subject effects of the panel-lvgvar are modeled as variance–covariance matrices and all variables in the networks are observed, the model reduces to a random intercept cross-lagged panel model (Hamaker et al., Reference Hamaker, Kuiper and Grasman2015), with the exception that stationarity is used to avoid estimating the variance–covariance structure of the first measurement. (5) When only cross-sectional data are used, the panel-lvgvar reduces to a latent network model (Epskamp, Rhemtulla, & Borsboom, Reference Epskamp, Rhemtulla and Borsboom2017), after identifying all within-subject variance–covariance structures to be zero. Subsequently, when in cross-sectional data all variables in the network are assumed observed without measurement error, the model reduces to a standard GGM (Epskamp & Fried, Reference Epskamp and Fried2018). (6) When the contemporaneous and between-subject levels of the panel-lvgvar are modeled using variance–covariance matrices, the model can be seen as a special case of general dynamic SEM models (Ciraki, Reference Ciraki2007). When the data contain many persons and time series, it is in principle possible to obtain both fixed-effect and person-wise estimates (random effects), but this is beyond the scope of the current paper (e.g., Asparouhov, Hamaker, & Muthén, Reference Asparouhov, Hamaker and Muthén2018; Epskamp, Waldorp, et al., Reference Epskamp, Waldorp, Mõttus and Borsboom2018).

Future Directions The current paper may motivate several lines of research: (1) The current formulations of the ts-lvgvar and panel-lvgvar only model contemporaneous relationships through undirected networks. It may instead be of interest to model these through directed relationships, as is typical in SEM and in structural vector-autoregression (Epskamp, Waldorp, et al., Reference Epskamp, Waldorp, Mõttus and Borsboom2018; Gates & Molenaar, Reference Gates and Molenaar2012; Gates, Molenaar, Hillary, Ram, & Rovine, Reference Gates, Molenaar, Hillary, Ram and Rovine2010). Furthermore, both frameworks may be combined into a unified framework that allows for the estimation of mixed directional and undirectional models. (2) Multi-group models for both the ts-lvgvar (in which case a “group” would be a single subject) and panel-lvgvar may be investigated. The psychonetrics package allows for such multi-group analyses as well as testing for measurement invariance and homogeneity in network structures and mean structures. Future researchers could investigate the performance of such tests and provide detailed guidelines on how to assess measurement invariance and homogeneity in these models. (3) The psychonetrics package also allows for the modeling of residual structures as a GGM, which may be used to model residual network models as further discussed by Epskamp, Rhemtulla, and Borsboom (Reference Epskamp, Rhemtulla and Borsboom2017); however, this was also beyond the scope of this paper. (4) More detailed investigation on departures from normality in the ts-lvgvar and panel-lvgvar (or network models in general) may be investigated. For example, through the use of robust maximum likelihood estimation (Satorra & Bentler, Reference Satorra and Bentler1994), threshold models (Muthén, Reference Muthén1984) or network models for binary data (Marsman et al., Reference Marsman, Borsboom, Kruis, Epskamp, van Bork, Waldorp and Marsman2018), or mixed categorical and continuous data (Haslbeck & Waldorp, Reference Haslbeck and Waldorp2015). (5) While reported simulation studies investigated several different model search strategies, many more model search strategies could be envisioned which may perform better than the strategies used in this paper. Future research could therefore investigate alternative estimation procedures. Of note, regularization techniques are often used now in the estimation of psychological network models for both time-series and cross-sectional data (Epskamp, Waldorp, et al., Reference Epskamp, Waldorp, Mõttus and Borsboom2018; Epskamp & Fried, Reference Epskamp and Fried2018). These techniques have, however, also been criticized for poorer performance than unregularized estimation, especially in large samples (Williams & Rast, Reference Williams and Rast2018), and in recent literature unregularized GGM estimation grew more common (e.g., Isvoranu et al., Reference Isvoranu, Guloksuz, Epskamp, van Os, Borsboom and Investigators2019). Promising lines of research on regularization involving latent variables (e.g., Chandrasekaran, Parrilo, & Willsky, Reference Chandrasekaran, Parrilo and Willsky2012; Yuan, Reference Yuan2012; Jacobucci, Grimm, & McArdle, Reference Jacobucci, Grimm and McArdle2016) exist and may lead to regularized network estimation methods for the ts-lvgvar and panel-lvgvar.

Conclusion The ts-lvgvar and panel-lvgvar frameworks extend the current toolbox of network psychometrics by (1) including latent variable models, (2) separating within- and between-subject variances, and (3) allowing for the estimation of within-subject temporal network models. Note that although this paper only contains examples with latent variables, the models can readily be used to estimate network models between observed variables by setting factor loadings to an identity matrix and setting all residual variances to zero. As such, the panel-lvgvar in particular offers a powerful new method for estimating fixed-effect temporal and contemporaneous networks from large sample panel data rather than intensive time-series data.

Acknowledgements

I would like to thank Katharina Jorgensen for detailed comments on an earlier draft of this paper. This work was carried out on the Dutch national e-infrastructure with the support of SURF Cooperative: all analyses were performed on the SurfSara HPC Cloud and all simulations were performed on the SurfSara Lisa Cluster. The LISS panel data were collected by CentERdata (Tilburg University, The Netherlands) through its MESS project funded by the Netherlands Organization for Scientific Research.

Footnotes

Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

¹ Note: the network models discussed here include nodes that represent random variables. The term network model is also often used to describe models that feature random edges and/or the observed variables (e.g., social networks). For an extensive discussion on the differences between these kinds of models, see Epskamp, Fried, et al. (Reference Epskamp, Fried, van Borkulo, Robinaugh, Marsman, Dalege and Cramer2018).

² An alternative perspective is that the lvgvar is a combination of a dynamic factor model (Molenaar, Reference Molenaar1985, Reference Molenaar2017) with a GGM on the autoregression’s residuals.

³ http://psychonetrics.org.

⁴ Henceforth, superscripts on μ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\mu }$$\end{document} and Σ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb {\Sigma }$$\end{document} will clarify the variable vector over which the means/variances are defined. For notational simplicity, I drop superscript ( z ) \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(\pmb {z})$$\end{document} .

⁵ In the ts-lvgvar, the model is estimated for an individual subject, so measurement invariance assumptions across subjects and across within- and between-subject models are not needed per se there for estimation. Furthermore, measurement invariance across subjects can be tested in the ts-lvgvar through multi-group models already implemented in psychonetrics, but this was beyond the scope of the current paper.

⁶ The item “lonely” featured in the top three indicators of two factors.

⁷ It should be noted that there is nothing “confirmatory” about this measurement model because it was obtained using a highly exploratory method.

⁸ Because the residual variance–covariance matrix is diagonal, its Cholesky decomposition is simply a diagonal matrix with the square root of the residual variances on the diagonal.

⁹ Detailed information and access to the data is available at www.lissdata.nl.

¹⁰ Further inspection revealed that this was a result of using the Cholesky decomposition; estimating the residual matrices using a variance–covariance matrix led to a Heywood case with a negative residual variance. These are hard to interpret, so I continued using a Cholesky decomposition and instead fixed this residual variance to zero.

References

Adolf, J.,Schuurman, N.K., Borkenau, P.,Borsboom, D., & Dolan, C.V., (2014). Measurement invariance within and between individuals: A distinct problem in testing the equivalence of intra-and inter-individual model structures. Frontiers in Psychology, 5, 883CrossRef Google Scholar PubMed

Asparouhov, T., Hamaker, E.L., & Muthén, B., (2018). Dynamic structural equation models. Structural Equation Modeling: A Multidisciplinary Journal, 25(3), 359–388CrossRef Google Scholar

Benjamini, Y., & Hochberg, Y., (1995). Controlling the false discovery rate: A practical and powerful approach to multiple testing. Journal of the Royal statistical society: series B (Methodological), 57(1), 289–300CrossRef Google Scholar

Bland, J.M., & Altman, D.G., (1995). Multiple significance tests: The Bonferroni method. BMJ, 310(6973), 170CrossRef Google Scholar PubMed

Borsboom, D.,(2017). A network theory of mental disorders. World Psychiatry, 16(1), 5–13CrossRef Google Scholar PubMed

Bos, F.M., Snippe, E., de Vos, S., Hartmann, J.A., Simons, C.J., van der Krieke, L., Wichers, M.,(2017). Can we jump from cross-sectional to dynamic interpretations of networks implications for the network perspective in psychiatry. Psychotherapy and Psychosomatics, 3, 175–177CrossRef Google Scholar

Bringmann, L.F., & Eronen, M.I.,(2018). Don’t blame the model: Reconsidering the network approach to psychopathology. Psychological Review, 125(4), 606–615CrossRef Google Scholar PubMed

Bringmann, L.F., Vissers, N., Wichers, M., Geschwind, N., Kuppens, P., Peeters, F., Tuerlinckx, F.,(2013). A network approach to psychopathology: New insights into clinical longitudinal data. PLoS ONE, 8(4), e60188CrossRef Google Scholar PubMed

Carver, C.S., Scheier, M.F., & Segerstrom, S.C., Optimism. Clinical Psychology Review 2010 30(7), 879–889CrossRef Google Scholar

Chandrasekaran, V., Parrilo, P.A., & Willsky, A.S.,(2012). Latent variable graphical model selection via convex optimization (with discussion) The Annals of Statistics, 40(4), 1935–1967Google Scholar

Ciraki, D., (2007). Dynamic structural equation models: Estimation and interference (Unpublished doctoral dissertation). London School of Economics and Political Science (United Kingdom).Google Scholar

Cramer, AOJWaldorp, L.J., van der Maas, HLJ & Borsboom, D.,(2010). Comorbidity: A network perspective. Behavioral and Brain Sciences, 33 2–3137–150CrossRef Google Scholar PubMed

De Ron, J., Fried, E., I., & Epskamp, S., (2019). Psychological networks in clinical populations: A tutorial on the consequences of Berkson’s bias. https://doi.org/10.31234/osf.io/5t8zw.CrossRef Google Scholar

Diener, E., Emmons, R.A., Larsen, R.J., & Griffin, S.,(1985). The satisfaction with life scale. Journal of Personality Assessment, 49(1), 71–75CrossRef Google Scholar PubMed

Epskamp, S., Borsboom, D., & Fried, E.I.,(2017). Estimating psychological networks and their accuracy: A tutorial paper. Behavior Research Methods,. https://doi.org/10.3758/s13428-017-0862-1.Google Scholar

Epskamp, S., Cramer, AOJWaldrop, L.J., Schmittmann, V.D., & Borsboom, D.,(2012). qgraph. Network visualizations of relationships in psychometric data. Journal of Statistical Software, 48(4), 1–18CrossRef Google Scholar

Epskamp, S., & Fried, E.I.,(2018). A tutorial on regularized partial correlation networks. Psychological Methods, 23(4), 617–634. https://doi.org/10.1037/met0000167.CrossRef Google Scholar PubMed

Epskamp, S., Fried, E.I., van Borkulo, C.D., Robinaugh, D.J., Marsman, M., Dalege, J., Cramer, AOJ (2018). Investigating the utility of fixed-margin sampling in network psychometrics. Multivariate Behavioral Research,. https://doi.org/10.1080/00273171.2018.1489771.Google Scholar PubMed

Epskamp, S., Rhemtulla, M., & Borsboom, D.,(2017). Generalized network pschometrics: Combining network and latent variable models. Psychometrika, 82(4), 904–927. https://doi.org/10.1007/s11336-017-9557-x.CrossRef Google Scholar PubMed

Epskamp, S., Van Borkulo, C., Van Der Veen, D., Servaas, M., Isvoranu, A-MRiese, H., Cramer, A.,(2018). Personalized network modeling in psychopathology: The importance of contemporaneous and temporal connections. Clinical Psychological Science, 6(4), 416–427. https://doi.org/10.17605/OSF.IO/JNPRZ.CrossRef Google Scholar PubMed

Epskamp, S., Waldorp, L.J., Mõttus, R., & Borsboom, D.,(2018). The Gaussian graphical model in cross-sectional and time-series data. Multivariate Behavioral Research, 53(4), 453–480. https://doi.org/10.1080/00273171.2018.1454823.CrossRef Google Scholar PubMed

Fisher, A.J., Medaglia, J.D., & Jeronimus, B.F.,(2018). Lack of group-to-individual generalizability is a threat to human subjects research. Proceedings of the National Academy of Sciences, 115(27), E6106–E6115CrossRef Google Scholar PubMed

Flake, J. K., & Fried, E.I., (2019). Measurement schmeasurement: Questionable measurement practices and how to avoid them. https://doi.org/10.31234/osf.io/hs7wm.CrossRef Google Scholar

Fried, E.I., & Cramer, A.O.,(2017). Moving forward: Challenges and directions for psychopathological network theory and methodology. Perspectives on Psychological Science, 12(6), 999–1020CrossRef Google Scholar PubMed

Fried, E.I., van Borkulo, C.D., Cramer, AOJBoschloo, L., Schoevers, R.A., & Borsboom, D.,(2017). Mental disorders as networks of problems: A review of recent insights. Social Psychiatry and Psychiatric Epidemiology, 52(1), 1–10. https://doi.org/10.1007/s00127-016-1319-z.CrossRef Google Scholar PubMed

Gates, K.M., & Molenaar, P.C.,(2012). Group search algorithm recovers effective connectivity maps for individuals in homogeneous and heterogeneous samples. NeuroImage, 63(1), 310–319. https://doi.org/10.1016/j.neuroimage.2012.06.026.CrossRef Google Scholar PubMed

Gates, K.M., Molenaar, P.C., Hillary, F.G., Ram, N., & Rovine, M.J.,(2010). Automatic search for fMRI connectivity mapping: An alternative to Granger causality testing using formal equivalences among SEM path modeling, VAR, and unified SEM. NeuroImage, 50(3), 1118–1125. https://doi.org/10.1016/j.neuroimage.2009.12.117.CrossRef Google Scholar PubMed

Golino, H.F., & Epskamp, S.,(2017). Exploratory graph analysis: A new approach for estimating the number of dimensions in psychological research. PLoS ONE, 12(6), e0174035. https://doi.org/10.1371/journal.pone.0174035.CrossRef Google Scholar PubMed

Guyon, H., Falissard, B., & Kop, J-L(2017). Modeling psychological attributes in psychology—An epistemological discussion: Network analysis vs. latent variables. Frontiers in Psychology, 8, 798CrossRef Google Scholar

Hamaker, E. L.,(2012). Why researchers should think “within-person”: A paradigmatic rationale. In Mehl, M.R. & Conner, T.S. (Eds.), Handbook of research methods for studying daily life (pp. 43-61).Google Scholar

Hamaker, E.L., Dolan, C.V., & Molenaar, P.C.,(2002). On the nature of sem estimates of arma parameters. Structural Equation Modeling, 9(3), 347–368CrossRef Google Scholar

Hamaker, E.L., Kuiper, R.M., & Grasman, R.P.,(2015). A critique of the cross-lagged panel model. Psychological Methods, 20(1), 102CrossRef Google Scholar PubMed

Haslbeck, J.M.B., , & Waldorp, L.J., (2015). mgm: Structure estimation for time-varying mixed graphical models in high-dimensional data. arxiv:1510.06871.Google Scholar

Holm, S.,(1979). A simple sequentially rejective multiple test procedure. Scandinavian Journal of Statistics, 6(2), 65–70Google Scholar

Isvoranu, A.M., Borsboom, D., van Os, J., & Guloksuz, S.,(2016). A network approach to environmental impact in psychotic disorders: Brief theoretical framework. Schizophrenia Bulletin, 42(4), 870–873CrossRef Google Scholar PubMed

Isvoranu, A.-M., Guloksuz, S., Epskamp, S.,van Os, J., Borsboom, D., Investigators, G., et al. (2019). Toward incorporating genetic risk scores into symptom networks of psychosis. Psychological Medicine. https://doi.org/10.1017/S003329171900045X.CrossRef Google Scholar

Jacobucci, R., Grimm, K.J., & McArdle, J.J.,(2016). Regularized structural equation modeling. Structural Equation Modeling: A Multidisciplinary Journal, 23(4), 555–566CrossRef Google Scholar PubMed

Kan, K-Jvan der Maas, H.L.,& Levine, S.Z.,(2019). Extending psychometric network analysis: Empirical evidence against g in favor of mutualism? Intelligence, 73, 52–62CrossRef Google Scholar

Kossakowski, J.J., Groot, P.C., Haslbeck, J.M.B.Borsboom, D., & Wichers, M.,(2017). Data from ‘critical slowing down as a personalized early warning signal for depression’ Journal of Open Psychology Data, 5(1), 1. https://doi.org/10.5334/jopd.29.CrossRef Google Scholar

Lauritzen, S.L.,(1996). Graphical models,Oxford:Clarendon PressCrossRef Google Scholar

Marsman, M., Borsboom, D., Kruis, J., Epskamp, S., van Bork, R., Waldorp, L., Marsman, M.,(2018). An introduction to network psychometrics: Relating ising network models to item response theory models. Multivariate Behavioral Research, 53(1), 15–35. https://doi.org/10.1080/00273171.2017.1379379.CrossRef Google Scholar PubMed

Molenaar, P.C.,(1985). A dynamic factor model for the analysis of multivariate time series. Psychometrika, 50(2), 181–202CrossRef Google Scholar

Molenaar, P.C.,(2017). Equivalent dynamic models. Multivariate Behavioral Research, 52(2), 242–258CrossRef Google Scholar PubMed

Molenaar, P.C.M.(2004). A manifesto on psychology as idiographic science: Bringing the person back into scientific psychology, this time forever. Measurement: Interdisciplinary Research & Perspective, 2(4), 201–218. https://doi.org/10.1207/s15366359mea0204_1.Google Scholar

Muthén, B.,(1984). A general structural equation model with dichotomous, ordered categorical, and continuous latent variable indicators. Psychometrika, 49(1), 115–132CrossRef Google Scholar

Revelle, W., (2019). psych: Procedures for personality and psychological research. Evanston, IL: Northwestern University. https://CRAN.R-project.org/package=psych (version 1.9.12).Google Scholar

Rhemtulla, M.,Van Bork, R., & Cramer, A.O.J., (2017). Cross-lagged network models. Multivariate Behavioral Research. Preprint from https://osf.io/r24q6/.Google Scholar

Rosenberg, M.,(1965). Rosenberg self-esteem scale (rse) Acceptance and Commitment Therapy. Measures Package, 61, 52Google Scholar

Satorra, A., & Bentler, P.M., (1994). Corrections to test statistics and standard errors in covariance structure analysis. In A. von Eye & C. C. Clogg (Eds.), Latent variables analysis: Applications for developmental research (pp. 399–419). Sage Publications,Inc.,Google Scholar

Scherpenzeel, A.C., & Das, M., (2010). True” longitudinal and probability-based internet panels: Evidence from the netherlands. In Das, M.Ester, P., & Kaczmirek, L., (Eds.), Social and behavioral research and the internet: Advances in applied methods and research strategies (pp. 77–104).Google Scholar

Schuurman, N.K., Houtveen, J.H., Hamaker, E.L.,(2015). Incorporating measurement error in n = 1 psychological autoregressive modeling. Frontiers in Psychology, 6, 1038. https://doi.org/10.3389/fpsyg.2015.01038.CrossRef Google Scholar

Tuccitto, D.E., Giacobbi, P.R. Jr & Leite, W.L.,(2010). The internal structure of positive and negative affect: A confirmatory factor analysis of the panas. Educational and Psychological Measurement, 70(1), 125–141CrossRef Google Scholar

Watson, D., Clark, L.A., & Tellegen, A.,(1988). Development and validation of brief measures of positive and negative affect: The panas scales. Journal of Personality and Social Psychology, 54(6), 1063CrossRef Google Scholar PubMed

Wichers, M., Groot, P.C., Psychosystems, E.S.M.Lenin, E.W.S.etal(2016). Critical slowing down as a personalized early warning signal for depression. Psychotherapy and Psychosomatics, 85(2), 114–116CrossRef Google Scholar PubMed

Wild, B., Eichler, M., Friederich, H-CHartmann, M., Zipfel, S., & Herzog, W.,(2010). A graphical vector autoregressive modelling approach to the analysis of electronic diary data. BMC Medical Research Methodology, 10(1), 28CrossRef Google Scholar

Williams, D., R., & Rast, P., (2018). Back to the basics: Rethinking partial correlation network methodology. British Journal of Mathematical and Statistical Psychology. https://doi.org/10.1111/bmsp.12173.CrossRef Google Scholar

Yuan, M.,(2012). Discussion: Latent variable graphical model selection via convex optimization. The Annals of Statistics, 40, 1968–1972CrossRef Google Scholar

Zevon, M.A., & Tellegen, A.,(1982). The structure of mood change: An idiographic/nomothetic analysis. Journal of Personality and Social Psychology, 43(1), 111CrossRef Google Scholar

Figure 1. Model search strategy used for given α\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\alpha $$\end{document} level (here, α=0.01\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\alpha = 0.01$$\end{document} was used) in the ts-lvgvar and panel-lvgvar example analyses. For the ts-lvgvar, model selection is performed on the temporal network Bp\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\pmb {B}_p$$\end{document} and the contemporaneous network Ωp(ζ)\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\pmb {\Omega }^{(\pmb {\zeta })}_p$$\end{document}, and for the panel-lvgvar model selection is performed on the temporal network B∗\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$ \pmb {B}_*$$\end{document}, the contemporaneous network Ω(ζ)within\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\underset{\mathrm {within}}{\pmb {\Omega }^{(\pmb {\zeta })}} $$\end{document}, and the between-subject network Ω(ζ)between\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\underset{\mathrm {between}}{\pmb {\Omega }^{(\pmb {\zeta })}}$$\end{document}. This algorithm has been implemented in the modelsearch function in the psychonetrics package.

Figure 2. Results of the ts-lvgvar analysis of the post-assessment data from Kossakowski et al.(2017). Blue edges indicate positive effects, and ref edges indicate negative effects. The estimated within-subject latent network structures (b and d) were estimated together with the measurement model in a. Panels c and e show how often each edge was included in a 25% case-drop bootstrap (Color figure online).

Table 1. Numeric results of the ts-lvgvar analysis of the post-assessment data from Kossakowski et al.(2017).

Table 2. The number of times each parameter was included in the case-drop bootstrap ts-lvgvar analysis of the post-assessment data from Kossakowski et al.(2017).

Table 3. Measurement model for LISS data example.

Table 4. Numeric results of the panel-lvgvar analysis of the LISS Core Study on personality.

Figure 3. Results of the panel-lvgvar analysis of the LISS Core Study on personality. Blue edges indicate positive effects, and ref edges indicate negative effects. The estimated fixed-effect within-subject latent networks are shown in b and c and the estimated between-subject network in d. These are estimated jointly with the measurement model parameters represented in a. Panels e, f and g show the inclusion proportion of each edge in a 25% case-drop bootstrap (Color figure online).

Table 5. The number of times each parameter was included in the case-drop bootstrap panel-lvgvar analysis of the LISS Core Study on personality.

Figure 4. Simulation results for ts-lvgvar and panel-lvgvar model estimation algorithms implemented in psychonetrics. Data were generated under parameters from Fig. 2 for the ts-lvgvar simulations and under parameters from Fig. 3 for the panel-lvgvar simulations, and each condition was replicated 100 times. The lines display the average values over these replications. “Sensitivity” denotes the proportion of true edges that were also included in the estimated model, “specificity” denotes the proportion of true missing edges that were also not included in the estimated model, “correlation” denotes the Pearson correlation between true and estimated edge weights, and “success” denotes the proportion of models estimations that ran without errors. The thick highlighted line represents the algorithm used in the empirical examples in this paper.

Article contents

Psychometric Network Models from Time-Series and Panel Data

Abstract

Keywords

Information

1. Introduction

2. Preliminary Topics

2.1. Notation

2.2. Estimation

2.3. Factor Model

2.4. Network Models

2.4.1. The Gaussian Graphical Model

2.4.2. The Graphical Vector-Autoregression Model

3. Time-Series Data: The ts-lvgvar Model

3.1. Estimation

3.2. Empirical Example

4. Panel Data: The panel-lvgvar Model

4.1. Estimation

4.2. Empirical Example

5. Simulation Studies

6. Ergodicity and Generalizability

7. Discussion

Acknowledgements

Footnotes

References

Save article to Kindle

Save article to Dropbox

Save article to Google Drive

Reply to: Submit a response

Your details

You have entered the maximum number of contributors

Conflicting interests