Latent variable models with nonparametric interaction effects of latent variables

Longitudinal Modeling of Age-Dependent Latent Traits with Generalized Additive Latent and Mixed Models

We present generalized additive latent and mixed models (GALAMMs) for analysis of clustered data with responses and latent variables depending smoothly on observed variables. A scalable maximum likelihood estimation algorithm is proposed, utilizing the Laplace approximation, sparse matrix computation, and automatic differentiation. Mixed response types, heteroscedasticity, and crossed random effects are naturally incorporated into the framework. The models developed were motivated by applications in cognitive neuroscience, and two case studies are presented. First, we show how GALAMMs can jointly model the complex lifespan trajectories of episodic memory, working memory, and speed/executive function, measured by the California Verbal Learning Test (CVLT), digit span tests, and Stroop tests, respectively. Next, we study the effect of socioeconomic status on brain structure, using data on education and income together with hippocampal volumes estimated by magnetic resonance imaging. By combining semiparametric estimation with latent variable modeling, GALAMMs allow a more realistic representation of how brain and cognition vary across the lifespan, while simultaneously estimating latent traits from measured items. Simulation experiments suggest that model estimates are accurate even with moderate sample sizes.

Structure detection of semiparametric structural equation models with Bayesian adaptive group lasso

Structural equation models (SEMs) are widely recognized as the most important statistical tool for assessing the interrelationships among latent variables. This study develops a Bayesian adaptive group least absolute shrinkage and selection operator procedure to perform simultaneous model selection and estimation for semiparametric SEMs, wherein the structural equation is formulated using the additive nonparametric functions of observed and latent variables. We propose the use of basis expansions to approximate the unknown functions. By introducing adaptive penalties to the groups of basis expansions, the nonlinear, linear, or non-existent effects of observed and latent variables in the structural equation can be automatically detected. A simulation study demonstrates that the proposed method performs satisfactorily. This paper presents an application of revealing the observed and latent risk factors of diabetic kidney disease.

Semiparametric Bayesian inference on skew-normal joint modeling of multivariate longitudinal and survival data

We propose a semiparametric multivariate skew-normal joint model for multivariate longitudinal and multivariate survival data. One main feature of the posited model is that we relax the commonly used normality assumption for random effects and within-subject error by using a centered Dirichlet process prior to specify the random effects distribution and using a multivariate skew-normal distribution to specify the within-subject error distribution and model trajectory functions of longitudinal responses semiparametrically. A Bayesian approach is proposed to simultaneously obtain Bayesian estimates of unknown parameters, random effects and nonparametric functions by combining the Gibbs sampler and the Metropolis-Hastings algorithm. Particularly, a Bayesian local influence approach is developed to assess the effect of minor perturbations to within-subject measurement error and random effects. Several simulation studies and an example are presented to illustrate the proposed methodologies.