Forming bootstrap confidence intervals and examining bootstrap distributions of standardized coefficients in structural equation modelling: A simplified workflow using the R package semboottools.

Standardized coefficients - including factor loadings, correlations, and indirect effects - are fundamental to interpreting structural equation modeling (SEM) results in psychology. However, they often exhibit skewed sampling distributions in finite samples, which are not captured by conventional symmetric confidence intervals (CIs). Methods such as bootstrap CI that do not impose symmetry are more appropriate for these coefficients. Despite its popularity, the widely used R package lavaan (version 0.6-19 or earlier) provides limited bootstrap support for standardized coefficients. Specifically, its function standardizedSolution() uses the delta method for CIs and lacks bootstrap p values. It provides a flexible and powerful function, bootstrapLavaan(), for bootstrapping, and it can be used to form bootstrap CIs for the standardized coefficients. However, this function requires a certain level of R coding skills. Moreover, no built-in functions are available to inspect bootstrap distributions, which are recommended for assessing the stability of the bootstrap estimates. To address these limitations, we developed the semboottools R package, which provides a simple workflow in SEM to form bootstrap confidence intervals for unstandardized and standardized estimates of model and user-defined parameters. It allows researchers to generate percentile or bias-corrected bootstrap CIs, standard errors, asymmetric p values, compare the bootstrap CIs with other CI methods (e.g., delta method), and visualize the distributions of bootstrap estimates - with minimal coding effort. We believe the tool can facilitate researchers in easily forming bootstrap CIs, comparing different CI methods to assess the need for bootstrapping, and examining the distribution of bootstrap estimates to assess their stability.