The hardware and bandwidth for this mirror is donated by dogado GmbH, the Webhosting and Full Service-Cloud Provider. Check out our Wordpress Tutorial.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]dogado.de.

lavDiag

Diagnostics and Visualization for Latent Variable Models

Overview

lavDiag extends the lavaan ecosystem with a suite of diagnostic, visualization, and empirical-fit tools for latent variable models (CFA, SEM, and related frameworks). It provides fast, parallel-safe computation of factor scores, model-based predictions, and empirical versus model fit curves for both continuous and ordinal indicators.

All functions are designed to work seamlessly with single-group and multi-group models, returning tidy tibble outputs ready for plotting or downstream analysis.

Key Features

• 🔹 Fast parallel factor scores lavPredict_parallel() — a robust, ordinal-aware replacement for lavaan::lavPredict().

• 🔹 Augmentation and diagnostics augment() — attaches model predictions, residuals, SEs, and CIs to observed data.

• 🔹 Model-based grids prepare() — generates smooth latent grids and model-based item curves for continuous, ordinal, and mixed indicators.

• 🔹 Empirical vs. model item curves item_data() + item_plot() — compare model-implied and empirical (GAM-based) curves with metrics like R², RMSE, and MAE.

• 🔹 Residual diagnostics

  • resid_cor() — extract residual correlations
  • resid_corrplot() — visualize residual correlation matrices
  • resid_qq() — Q–Q plots of residual z-statistics
  • hopper_plot() — visualize residual “hopper” patterns

• 🔹 Model summaries and estimates model_info() and parameter_estimates() — consistent schema for metadata and estimates.

• 🔹 CFA visualization plot_cfa() — clean diagrams using visNetwork.

Installation

From GitHub

# install.packages("remotes")
remotes::install_github("reckak/lavDiag")

From CRAN (once accepted)

install.packages("lavDiag")

Quick Example

library(lavaan)
library(lavDiag)

# Example CFA
HS.model <- '
  visual  =~ x1 + x2 + x3
  textual =~ x4 + x5 + x6
  speed   =~ x7 + x8 + x9
'

fit <- cfa(HS.model, data = HolzingerSwineford1939, meanstructure = TRUE)

# Augment observed data with model predictions and residuals
aug <- augment(fit)

# Compute and visualize item-level empirical fit
it <- item_data(fit)
item_plot(it)

# Residual correlation plot
resid_corrplot(fit, type = "cor.bentler")

Parallelization

Most computationally intensive functions (e.g., lavPredict_parallel(), item_data(), prepare()) use future and furrr backends for safe parallelism. Backends are configurable via .set_future_plan().

Dependencies

Core dependencies include:

• lavaan, dplyr, tidyr, purrr, tibble
• future, furrr, mgcv, corrplot, ggplot2, visNetwork

All functions use consistent tidy-style output and rlang-safe programming.

Contributing

Issues, pull requests, and feedback are welcome!

• 📬 Report bugs: GitHub Issues
• 💡 Feature requests: open an issue or start a discussion

Citation

If you use lavDiag in your research, please cite it as:

Rečka, K. (2025). lavDiag: Diagnostics and Visualization for Latent Variable Models. GitHub: https://github.com/reckak/lavDiag

License

MIT License © 2025 Karel Rečka

These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.
Health stats visible at Monitor.