--- title: Partial Least Squares Structural Equation Modeling url: https://www.emergentmind.com/topics/partial-least-squares-structural-equation-modeling-pls-sem-0a1d5ba7-5bbc-4706-b00e-b54ba9783ab3 type: topic --- # Partial Least Squares Structural Equation Modeling Partial Least Squares Structural Equation Modeling (PLS-SEM) is a variance-based structural equation modeling methodology that enables simultaneous estimation of relationships among latent constructs and their manifest indicators, while maximizing the explained variance of endogenous variables. Distinct from covariance-based SEM, PLS-SEM is particularly suited for complex models, small-to-moderate sample sizes, non-normal indicator distributions, and theory development or prediction-focused research contexts. ## 1. Theoretical Foundations and Model Structure PLS-SEM decomposes models into two subcomponents: - **Outer (Measurement) Model**: Specifies how manifest variables (indicators) relate to latent variables (LVs). Reflective models assume the LV causes the observed measures, specified as $X_{pk} = \alpha_{pk} + \lambda_{pk}\xi_k + \varepsilon_{pk}$. Formative models define the LV as a linear combination of its measurements, as $\xi_k = \delta_{0k} + \sum_{p=1}^{P_k} \delta_{pk}X_{pk} + \zeta_k$ [2210.12520]. - **Inner (Structural) Model**: Specifies hypothesized causal relationships among LVs, typically in the form $\xi_j = \beta_{j0} + \sum_{k\in \mathrm{Pred}(j)}\beta_{jk}\xi_k + \eta_j$, where $\beta_{jk}$ are path coefficients [2210.12520, 2209.04900]. PLS-SEM maximizes the explained variance ($R^2$) of endogenous LVs, supports both reflective and formative blocks, and operates under minimal distributional assumptions—no requirement for multivariate normality or large samples [2307.11895]. ## 2. Estimation Algorithms and Computational Workflow The canonical PLS-SEM estimation follows Wold’s path-modeling iteration: - **Initialization**: Standardize indicators, and assign unit or correlation-based weights to indicators within each LV block [2209.04900, 2307.11895]. - **Outer Approximation (Measurement Model)**: Compute weighted LV scores as $ŷ_k = \sum w_{kj}x_j$. Update weights iteratively, typically as $w^{(t+1)}_{kj} \propto \operatorname{cov}(x_j, ŷ_k)/\operatorname{var}(x_j)$ [2209.04900]. - **Inner Approximation (Structural Model)**: Generate “inner” latent variable estimates as linear combinations of neighboring LV scores, per the model’s path structure (e.g., $u_j = \sum_k b_{jk}t_k$) [2307.11895]. Estimate path coefficients via ordinary least squares regression on the LV scores. - **Weight Updating and Convergence**: Iterate the above steps until the change in weights is less than a pre-specified tolerance. This two-stage algorithm enables estimation of both model parameters and LV scores simultaneously, enabling subsequent resampling-based inference and prediction [2311.03636]. ## 3. Model Evaluation: Measurement and Structural Components Rigorous model assessment encompasses both measurement and structural aspects: - **Indicator Reliability**: Outer loadings $\lambda_{ij} \geq 0.70$ are preferred, but values $0.50-0.70$ may be retained if overall construct validity is established [2209.04900]. All variance inflation factors (VIFs) should be $<5$ to rule out multicollinearity. - **Convergent Validity**: Average Variance Extracted (AVE)—$AVE_k = (\sum \lambda_{ki}^2)/m$—must exceed $0.50$ [2307.11895]. - **Construct Reliability**: Cronbach's $\alpha$ $(\alpha > 0.70)$ and composite reliability measures (e.g., $CR_k$) evaluate internal consistency [2311.03636]. - **Discriminant Validity**: The Fornell–Larcker criterion ($\sqrt{AVE_j} > \max_{k \neq j} \mathrm{Corr}(\eta_j, \eta_k)$), cross-loadings, and the Heterotrait–Monotrait Ratio (HTMT; $<$0.85 or 0.90) confirm construct distinctiveness [2307.11895, 2209.04900]. - **Structural Model Assessment**: Path coefficients ($\beta$), associated $t$-values (based on $5000$ bootstrap samples), and $p$-values are reported for hypothesis significance [2209.04900, 2311.03636]. Explanatory power is indicated by $R^2$ statistics and effect sizes $f^2$ [2307.11895]. - **Predictive Relevance**: Out-of-sample forecast metrics such as $Q^2$ (via blindfolding or PLSPredict), root mean square error (RMSE), and mean absolute error (MAE) against a linear model benchmark are reported [2209.04900]. ## 4. Special Topics: Segmentation, Ordinal Data, and Cyclical Causality ### Simultaneous Clustering and PLS-SEM PLS-SEM-KM integrates K-means clustering with the PLS-SEM algorithm, optimizing both cluster assignments and SEM parameters jointly. Unlike sequential workflows (PLS→clustering), the simultaneous approach produces clusters homogeneous with respect to structural relationships, enhancing segment validity—confirmed empirically by ARI improvements and simulation benchmarks against FIMIX-PLS [1810.07677]. ### Ordinal Partial Least Squares (OPLS) Traditional PLS-SEM is suboptimal for ordinal data with few categories (e.g., $I=4,5$). OPLS addresses this by employing a polychoric correlation substitution for the Pearson covariance matrix in all algorithmic computations. This adjustment substantially reduces negative bias in path coefficient estimation in small-category ordinal settings. For $I \geq 7$, OPLS converges to standard PLS-SEM estimates [1212.5049]. ### Modeling Cyclic (Reciprocal) Effects Standard SEM prohibits cyclic paths. A two-step approach with PLS-SEM enables modeling reciprocal causality in cross-sectional data: first estimate the acyclic model, then re-specify models with feedback paths using LV scores from Step 1 as proxies for 'lagged' values. Bootstrap-based parametric tests compare the strength of forward and cyclic effects [2210.12520]. This technique enabled demonstration that internet usage intensity both results from and reinforces digital skills and physical access. ## 5. Power Analysis and Sample Size Planning Statistical power is critical for PLS-SEM study design. The inverse square root method is the current standard: - **Required Sample Size**: $N = (p_\alpha / \beta_{\min})^2,$ where $p_\alpha$ corresponds to the chosen significance threshold (e.g., $p_\alpha=2.486$ at $\alpha=0.05$ for 80% power), and $\beta_{\min}$ is the smallest effect size of interest. For $\beta_{\min}=0.5$, $N=25$ [2511.14546]. - **Minimum Detectable Effect Size (MDES)**: $\beta_{\min}=p_\alpha / \sqrt{N}$. For $N=68$, $\beta_{\min} \approx 0.30$. Use of the "PLS-SEM-power" R package and Shiny application operationalizes this process for both a priori sample size determination and post hoc sensitivity analysis [2511.14546]. The method assumes one path at a time and requires $N>10$ for valid application. ## 6. Empirical Applications and Practical Workflow Empirical PLS-SEM studies typically adhere to the following protocol [2209.04900, 2307.11895, 2311.03636]: 1. Develop and operationalize a model with constructs and indicators (using Likert or similar scales). 2. Screen data and specify reflective/formative measurement blocks. 3. Apply the PLS algorithm (SmartPLS, R’s seminr, or equivalent). 4. Evaluate measurement validity (AVE, $\alpha$, CR, HTMT) and structural relationships ($\beta$, $R^2$, $f^2$). 5. Test hypotheses through bootstrapping. 6. Conduct out-of-sample predictive validation (PLS-Predict, $Q^2$, RMSE, MAE). 7. Interpret significant and non-significant effects, with attention to the largest effect sizes and practical relevance. 8. Report statistical power analysis, indicating both a priori planning and actual sensitivity achieved [2511.14546]. Empirical studies illustrate these steps in diverse domains, including technology adoption [2209.04900], digital transformation [2307.11895], and educational technology [2311.03636]. For ordinal data, OPLS should be selected for improved bias properties [1212.5049]. For population heterogeneity or market segmentation, joint PLS-SEM-KM estimation is preferable [1810.07677]. For causal feedback, iterative two-step PLS-SEM estimation is needed [2210.12520]. ## 7. Methodological Innovations and Adaptations Recent methodological advances include: - Introduction of simultaneous clustering within PLS-SEM (PLS-SEM-KM) to address unobserved heterogeneity [1810.07677]. - Ordinal PLS adaptation for manifest measures with few ordered categories, utilizing polychoric correlations and latent threshold modeling [1212.5049]. - Two-stage modeling of cyclic effects with cross-sectional data, enabling quantification of feedback mechanisms in socio-technical systems [2210.12520]. - Shiny/R tools for power analysis, automating computation of required $N$ and MDES via the inverse square root method [2511.14546]. - Application of advanced validation diagnostics—including $Q^2$ and out-of-sample RMSE/MAE—for robust predictive assessment [2307.11895, 2209.04900]. These innovations are anchored in systematic empirical testing and are integrated in current best-practice PLS-SEM workflows. --- For technical implementation and further detail, consult workflow exemplars and software scripts as reported in domain applications [2209.04900, 2307.11895, 2311.03636]. Source: https://www.emergentmind.com/topics/partial-least-squares-structural-equation-modeling-pls-sem-0a1d5ba7-5bbc-4706-b00e-b54ba9783ab3