Estimating Multilevel Structural Equation Models with Random Slopes with Laplace and Variational Approximations

Bibliographic Data

ID21641776
AuthorsSteffen Nestler (0000-0001-9724-2441, University of Münster, corresponding author)
Year2026
Volume33
Issue3
Pages335-345
Publication date2026-05-04
Peer ReviewedYes
Open AccessYes
TypeARTICLE
VenueStructural Equation Modeling: A Multidisciplinary Journal (JOURNAL)
Journal identifiersISSN: 1070-5511 • E-ISSN: 1532-8007
PublisherInforma UK Limited (PUBLISHER • GB)
DOI10.1080/10705511.2026.2616824
OpenAlexW7128630210
LanguageEN
References cited28

Multilevel structural equation models (MSEMs) are an important statistical approach to analyze hierarchically nested data. While Bayesian methods are commonly used to estimate the MSEM parameters, maximum likelihood (ML) approaches are less often employed because of the computational challenges in the required numerical integration. Building on Rockwood’s reformulation of the MSEM, we investigate two computationally efficient approximation methods to the likelihood function—the Laplace approximation (LA) and the extended variational approximation (EVA)—by comparing their performance against Gauss-Hermite (GH) quadrature and a Bayesian approach in two simulation studies. Results demonstrate that LA and EVA provide accurate parameter estimates with substantially shorter computation times than GH, especially for a more complex model. Furthermore, LA and EVA showed almost no bias, good convergence rates, and appropriate coverage, particularly with larger sample sizes. Altogether, these findings suggest that LA and EVA are promising approaches for ML estimation of MSEMs with random slopes

Approximations of π · Laplace transform · Laplace's equation · Laplace's method · Spatial and Panel Data Analysis · Statistical Methods and Bayesian Inference · Statistical Methods and Inference

Tools

Open DOIOpen Access