Robustness Testing in SEM: Building a Robustness-Test Matrix
Robustness is not running many variations and selecting the most attractive result. A useful check targets a specific threat, retains the primary model as reference and discloses changing conclusions.
Prespecify the primary model, then construct a threat-check-evidence matrix covering distributions and estimators, outliers, missingness, specification, endogeneity, nonlinearity, heterogeneity and measurement. CB-SEM and PLS-SEM require different tools; procedures cannot be transferred unchanged.
Robustness checks start with vulnerable assumptions
- Main conclusion
Identify the conclusion to examine
- Assumptions
Choose assumptions to vary
- Alternatives
Use defensible alternative specifications
- Comparison
Report stability and changes
Conceptual diagram. Do not retain only checks that support a preferred result.
Source: Lê Minh · Conceptual illustration, not software output.1. What does robustness testing ask?
A conclusion is robust when its direction, magnitude and substantive meaning do not change severely under plausible analytical alternatives. Each variation must represent a credible design, measurement or estimation risk.
Robust standard errors, robustness testing and replication differ. Robust corrections address certain distributional violations; robustness testing compares specifications; replication tests new data. None replaces the others.
Vaithilingam et al. (2024) review PLS-SEM robustness practice and emphasize matching techniques to threats. This article extends threat-based checking to CB-SEM while identifying PLS-specific tools.
- Specify the primary model before inspecting alternatives.
- Link each check to a specific threat.
- Compare direction, magnitude and uncertainty.
- Do not report only favorable alternatives.
- Use tools appropriate to CB-SEM or PLS-SEM.
- Disclose weakened conclusions.
2. Fix the primary model before robustness testing
The primary specification follows theory, design and the analysis plan. It defines variables, relationships, controls, estimator, missing-data treatment, scales and exclusions. Changing it after every result turns robustness into specification searching. Kline (2023) connects specification and assessment with theory, diagnostics and assessable comparisons.
Create a version table first. Each row should identify the analysis, one main change, the threat and the result being compared. This clarifies every alternative’s purpose.
Robustness cannot rescue an unidentified causal model. Cross-sectional, self-reported, nonrandom data still limit causal claims despite stable coefficients. Antonakis et al. (2010) warn that endogeneity can invalidate causal interpretation even with good fit.
| ID | Change from the primary model | Threat | Results compared |
|---|---|---|---|
| M0 | No change | Primary model | Estimate, SE or CI, fit, R2 or Q2 |
| R1 | Alternative estimator or bootstrap | Nonnormality | Direction, CI and fit |
| R2 | Retain/remove influential cases | Outliers | Estimate change |
| R3 | Another plausible control specification | Omitted variable | Target estimate |
| R4 | Groups or segments | Heterogeneity | Path differences |
| R5 | Add a nonlinear component | Nonlinearity | Shape and predictive relevance |
3. Map threats before selecting techniques
A long test list does not automatically establish credibility. Ask which assumptions may fail, through what mechanism, which conclusions would be affected and how far the available data can investigate.
Not every study needs every check. Same-source cross-sectional studies should prioritize measurement quality, common method and causal limitations. Group comparisons need invariance. Prediction-oriented PLS-SEM also requires out-of-sample evaluation.
The following matrix selects checks by risk rather than available software.
| Threat | Signal | Possible check | Evidence to report |
|---|---|---|---|
| Nonnormality | Heavy tails, high kurtosis CR | Robust SEs, bootstrap, suitable estimator | Estimate, SE or CI, fit |
| Outliers/influence | Large Mahalanobis distances, unusual residuals | Case review and sensitivity | Results retaining and addressing cases |
| Missing | Missingness rates or patterns differ by group | FIML, multiple imputation, sensitivity | Assumed mechanism and differences |
| Misspecification | Residual patterns or theoretical conflict | Alternative plausible models | Fit and target parameters |
| Endogeneity | Reverse causality, omitted causes | Design, IV, Gaussian copula where suitable for PLS | Identification conditions |
| Heterogeneity | Observed groups or latent mixtures | Invariance, MGA, segmentation | Differences and group sizes |
| Nonlinearity | Residual patterns, theory | Justified quadratic or interaction terms | Effect shape |
| Common method | One source, same occasion | Procedural design and marker/method factors | Limitations and sensitivity |
4. Sensitivity to distribution and estimator
In CB-SEM, comparing conventional ML with robust corrections or bootstrap assesses sensitivity of SEs, CIs and fit to nonnormality. Curran, West and Finch (1996) show interactions with misspecification; Savalei (2014) explains that robust corrections adjust inference, not incorrect models. Yuan and Bentler (2000) discuss likelihood methods for nonnormal missing data. Explain which estimator best fits assumptions rather than selecting the favorable p-value.
For AMOS, Bollen and Stine (1992) provide a bootstrap goodness-of-fit approach under nonnormality, while bootstrap CIs support some parameter inference. Few-category ordinal data or scaled statistics may require software with suitable categorical or robust estimators.
In PLS-SEM, bootstrap is common for inference but does not check every threat. Investigate endogeneity, nonlinearity, heterogeneity and prediction separately when justified.
Bootstrap usually quantifies sampling uncertainty. It becomes part of a robustness plan when linked to a specific question and compared with the primary model.
5. Outliers, influential cases and missing data
Unusual observations can strongly influence covariances and paths. Verify data, assess whether cases belong to the population, then compare retention and treatment under prespecified criteria. Do not delete simply because a p-value changes.
Describe missingness rates, patterns and assumed mechanisms. Complete-case analysis can reduce power and bias results when missingness relates to study variables. Consider FIML in CB-SEM under appropriate assumptions, or multiple imputation preserving the analytical structure.
Focus sensitivity on target parameters. If a path changes sign or its CI widens substantially under another plausible missing-data treatment, reduce certainty.
| Version | N | Target estimate | 95% CI | Conclusion |
|---|---|---|---|---|
| Primary model | 420 | 0,24 | [0,10; 0,38] | Positive effect |
| Without 3 influential cases | 417 | 0,22 | [0,08; 0,36] | No substantive change |
| FIML instead of complete cases | 438 | 0,19 | [0,04; 0,34] | Weaker, same direction |
| Multiple imputation | 438 | 0,20 | [0,05; 0,35] | Consistent with FIML |
6. Specification sensitivity and endogeneity
Alternatives must be theoretically plausible, not every imaginable configuration. Examples include models with and without a prespecified control or competing causal orders the design cannot distinguish. Identify prespecified and post hoc variations.
Endogeneity can arise from omitted variables, measurement error, reverse causality or selection. SEM fit does not establish exogeneity. Responses require better design, temporal data, instruments with defensible identification or specialized techniques with explicit assumptions.
Gaussian copula may investigate a form of PLS-SEM endogeneity when the regressor is nonnormal and other conditions hold. It cannot replace identification arguments. IV or simultaneous CB-SEM models also need explicit specification and identification.
A coefficient stable across controls may remain biased if every specification shares the same omitted cause or measurement error.
7. Heterogeneity, invariance and nonlinearity
For theoretically meaningful groups, assess invariance before comparing paths. CB-SEM multigroup testing usually proceeds from configuration through loading, intercept and residual constraints according to objectives. PLS-SEM uses MICOM to assess composite invariance before MGA.
Unobserved heterogeneity can make average effects conceal opposing mechanisms. Segmentation or mixtures may help, but solutions need stability, sufficient group sizes and substantive interpretation. Do not choose group count merely for favorable p-values.
Investigate nonlinearity when theory or residuals suggest it. Quadratic or interaction terms need predicted-value plots and supported data ranges. A nonsignificant linear coefficient does not establish no relationship when the true shape is curved.
| Threat | CB-SEM | PLS-SEM | Shared requirement |
|---|---|---|---|
| Invariance | Multigroup CFA | MICOM | Comparable measurement |
| Path differences | Constrained multigroup test | MGA or permutation | Theoretically justified groups and adequate samples |
| Unobserved heterogeneity | Mixture SEM | Segmentation | Stable and meaningful solutions |
| Nonlinearity | Latent interaction or quadratic term | Quadratic or interaction term | Hypotheses and supported data ranges |
| Prediction | Suitable holdout or cross-validation | PLSpredict or CVPAT | Separate fit from prediction |
8. A reproducible robustness workflow
Use this as an analysis appendix. The number of variations depends on risks, but readers should see the primary model before its alternatives.
Fix the estimand and primary model
State the target parameter/conclusion, model, estimator and analysis data.
List threats
Use the design, measures, data sources and CB-SEM/PLS-SEM objectives.
Choose checks for priority threats
Do not run a technique merely because software provides it.
Assign variation IDs
Record one main change from the primary model and the comparison result.
Check measurement
Ensure admissible solutions in every version and invariance for group comparisons.
Use consistent reporting
Save estimates, uncertainty, fit or prediction for favorable and unfavorable results.
Compare directions, magnitudes and CIs
Go beyond p-values and identify substantively meaningful changes.
Explain inconsistencies
Connect changes with threats rather than selecting the best-looking version.
Moderate conclusions when needed
State sensitivity, scope and the need for new data.
Preserve code, output and comparisons
Enable reproduction of the primary model and every sensitivity check.
9. Hypothetical example: how strong is a conclusion?
This example is illustrative. Primary CB-SEM on N = 438 estimates Information Quality → Usage Intention at beta = .24, 95% CI [.10, .38]. A robust estimator yields beta = .23, 95% CI [.08, .38]. FIML instead of complete cases yields beta = .19, 95% CI [.04, .34].
Adding a prespecified control reduces beta to .14, 95% CI [-.01, .29]. Direction remains positive but uncertainty includes 0. Evidence is therefore sensitive to controls rather than fully robust.
An honest report states stability to estimator and missing-data treatment but weaker evidence under a plausible control specification. Do not report only the first two alternatives.
| Version | Threat | beta | 95% CI | Assessment |
|---|---|---|---|---|
| M0 Primary | Primary model | .24 | [.10, .38] | Positive evidence |
| R1 Robust estimator | Nonnormality | .23 | [.08, .38] | Stable |
| R2 FIML | Missing | .19 | [.04, .34] | Weaker, same conclusion |
| R3 Control set | Omitted variable | .14 | [-.01, .29] | Specification-sensitive |
| R4 Remove influential cases | Outliers | .22 | [.07, .37] | Stable |
10. APA 7 reporting template and checklist
This template uses hypothetical values. Include a robustness table in the article or appendix, including alternatives that weaken conclusions.
Reporting example in English
The primary model indicated an association of [X] with [Y], beta = .24, 95% CI [.10, .38]. Results were similar using [robust estimator], beta = .23, 95% CI [.08, .38], and FIML for missing data, beta = .19, 95% CI [.04, .34]. Adding prespecified controls reduced the estimate to beta = .14, 95% CI [-.01, .29]. Conclusions were therefore stable to estimator and missing-data treatment but sensitive to the control specification. Results were interpreted as conditional associations rather than causal confirmation.
- Prespecify the estimand and primary model before sensitivity checks.
- Link every check to a specific threat.
- Distinguish robust correction from robustness testing.
- Use tools appropriate for CB-SEM or PLS-SEM.
- Compare estimates and CIs, not p-values alone.
- Inspect original data before treating outliers.
- State missing-data assumptions.
- Consider misspecification and endogeneity in the design.
- Assess invariance before group comparisons.
- Report alternatives that weaken conclusions.
- Do not equate stability with causal evidence.
- Preserve code, output and version tables for reproduction.
- Identify hypothetical examples.
- Match citations and references.
Conclusion
Robustness testing is valuable as a risk-assessment design rather than post hoc operations. Retain the primary model as reference and justify every alternative before inspecting results.
A conclusion can be robust, conditionally robust or sensitive. All are useful when honestly reported because they clarify the scope of the evidence.
References
- Antonakis, J., Bendahan, S., Jacquart, P., & Lalive, R. (2010). On making causal claims: A review and recommendations. The Leadership Quarterly, 21(6), 1086-1120.https://doi.org/10.1016/j.leaqua.2010.10.010
- Bollen, K. A., & Stine, R. A. (1992). Bootstrapping goodness-of-fit measures in structural equation models. Sociological Methods & Research, 21(2), 205-229.https://doi.org/10.1177/0049124192021002004
- Curran, P. J., West, S. G., & Finch, J. F. (1996). The robustness of test statistics to nonnormality and specification error in confirmatory factor analysis. Psychological Methods, 1(1), 16-29.https://doi.org/10.1037/1082-989X.1.1.16
- Kline, R. B. (2023). Principles and practice of structural equation modeling (5th ed.). Guilford Press.https://www.guilford.com/books/Principles-and-Practice-of-Structural-Equation-Modeling/Rex-Kline/9781462551910
- Savalei, V. (2014). Understanding robust corrections in structural equation modeling. Structural Equation Modeling: A Multidisciplinary Journal, 21(1), 149-160.https://doi.org/10.1080/10705511.2013.824793
- Vaithilingam, S., Ong, C. S., Moisescu, O. I., & Nair, M. S. (2024). Robustness checks in PLS-SEM: A review of recent practices and recommendations for future applications in business research. Journal of Business Research, 173, 114465.https://doi.org/10.1016/j.jbusres.2023.114465
- Yuan, K. H., & Bentler, P. M. (2000). Three likelihood-based methods for mean and covariance structure analysis with nonnormal missing data. Sociological Methodology, 30(1), 165-200.https://doi.org/10.1111/0081-1750.00078
Methodological note: Example numbers are hypothetical. MICOM, MGA and Gaussian copula belong to PLS-SEM workflows and cannot be transferred unchanged to CB-SEM. Robustness analysis does not turn observational data into causal evidence.
