Derivative Computations and Robust Standard Errors for Linear Mixed Effects Models in lme4

While robust standard errors and related facilities are available in R for many types of statistical models, the facilities are notably lacking for models estimated via lme4. This is because the necessary statistical output, including the Hessian and casewise gradient of random effect parameters, is not immediately available from lme4 and is not trivial to obtain. In this article, we supply and describe two new functions to obtain this output from Gaussian mixed models: estfun.lmerMod() and vcov.full.lmerMod(). We discuss the theoretical results implemented in the code, focusing on calculation of robust standard errors via package sandwich. We also use the Sleepstudy data to illustrate the code and compare it to a benchmark from package lavaan.Comment: Accepted at Journal of Statistical Softwar

Generalized Measurement Invariance Tests with Application to Factor Analysis

The issue of measurement invariance commonly arises in factor-analytic contexts, with methods for assessment including likelihood ratio tests, Lagrange multiplier tests, and Wald tests. These tests all require advance definition of the number of groups, group membership, and offending model parameters. In this paper, we construct tests of measurement invariance based on stochastic processes of casewise derivatives of the likelihood function. These tests can be viewed as generalizations of the Lagrange multiplier test, and they are especially useful for: (1) isolating specific parameters affected by measurement invariance violations, and (2) identifying subgroups of individuals that violated measurement invariance based on a continuous auxiliary variable. The tests are presented and illustrated in detail, along with simulations examining the tests' abilities in controlled conditions.measurement invariance, parameter stability, factor analysis, structural equation models

Testing non-nested structural equation models

In this paper, we apply Vuong's (1989) likelihood ratio tests of non-nested models to the comparison of non-nested structural equation models. Similar tests have been previously applied in SEM contexts (especially to mixture models), though the non-standard output required to conduct the tests has limited their previous use and study. We review the theory underlying the tests and show how they can be used to construct interval estimates for differences in non-nested information criteria. Through both simulation and application, we then study the tests' performance in non-mixture SEMs and describe their general implementation via free R packages. The tests offer researchers a useful tool for non-nested SEM comparison, with barriers to test implementation now removed.Comment: 24 pages, 6 figure

blavaan : Bayesian structural equation models via parameter expansion

This article describes blavaan, an R package for estimating Bayesian structural equation models (SEMs) via JAGS and for summarizing the results. It also describes a novel parameter expansion approach for estimating specific types of models with residual covariances, which facilitates estimation of these models in JAGS. The methodology and software are intended to provide users with a general means of estimating Bayesian SEMs, both classical and novel, in a straightforward fashion. Users can estimate Bayesian versions of classical SEMs with lavaan syntax, they can obtain state-of-the-art Bayesian fit measures associated with the models, and they can export JAGS code to modify the SEMs as desired. These features and more are illustrated by example, and the parameter expansion approach is explained in detail

Using Sunflower Plots and Classification Trees to Study Typeface Legibility

This article describes the application of sunflower plots and classification trees to the study of onscreen typeface legibility. The two methods are useful for describing high-dimensional data in an intuitive manner, which is crucial for interacting with both the typographers who design the typefaces and the practitioners who must make decisions about which typeface to use for specific applications. Furthermore, classification trees help us make specific recommendations for how much of a character attribute is “enough” to make it legible. We present examples of sunflower plots and classification trees using data from a recent typeface legibility experiment, and we present R code for replicating our analyses. Some familiarity with classification trees and logistic regression will be helpful to the reader

Computation and application of generalized linear mixed model derivatives using lme4

Maximum likelihood estimation of generalized linear mixed models(GLMMs) is difficult due to marginalization of the random effects. Computing derivatives of a fitted GLMM's likelihood (with respect to model parameters) is also difficult, especially because the derivatives are not by-products of popular estimation algorithms. In this paper, we describe GLMM derivatives along with a quadrature method to efficiently compute them, focusing on lme4 models with a single clustering variable. We describe how psychometric results related to IRT are helpful for obtaining these derivatives, as well as for verifying the derivatives' accuracies. After describing the derivative computation methods, we illustrate the many possible uses of these derivatives, including robust standard errors, score tests of fixed effect parameters, and likelihood ratio tests of non-nested models. The derivative computation methods and applications described in the paper are all available in easily-obtained R packages

Score‐based measurement invariance checks for Bayesian maximum‐a‐posteriori estimates in item response theory

A family of score-based tests has been proposed in recent years for assessing the invariance of model parameters in several models of item response theory (IRT). These tests were originally developed in a maximum likelihood framework. This study discusses analogous tests for Bayesian maximum-a-posteriori estimates and multiple-group IRT models. We propose two families of statistical tests, which are based on an approximation using a pooled variance method, or on a simulation approach based on asymptotic results. The resulting tests were evaluated by a simulation study, which investigated their sensitivity against differential item functioning with respect to a categorical or continuous person covariate in the two- and three-parametric logistic models. Whereas the method based on pooled variance was found to be useful in practice with maximum likelihood as well as maximum-a-posteriori estimates, the simulation-based approach was found to require large sample sizes to lead to satisfactory results

Efficient Bayesian Structural Equation Modeling in Stan

Structural equation models comprise a large class of popular statistical models, including factor analysis models, certain mixed models, and extensions thereof. Model estimation is complicated by the fact that we typically have multiple interdependent response variables and multiple latent variables (which may also be called random effects or hidden variables), often leading to slow and inefficient MCMC samples. In this paper, we describe and illustrate a general, efficient approach to Bayesian SEM estimation in Stan, contrasting it with previous implementations in R package blavaan (Merkle & Rosseel, 2018). After describing the approaches in detail, we conduct a practical comparison under multiple scenarios. The comparisons show that the new approach is clearly better. We also discuss ways that the approach may be extended to other models that are of interest to psychometricians.Comment: 21 pages, 5 figure