# Confirmatory factor analysis

In statistics, confirmatory factor analysis (CFA) is a special form of factor analysis, most commonly used in social research.[1] It is used to test whether measures of a construct are consistent with a researcher's understanding of the nature of that construct (or factor). As such, the objective of confirmatory factor analysis is to test whether the data fit a hypothesized measurement model. This hypothesized model is based on theory and/or previous analytic research.[2] CFA was first developed by Jöreskog[3] and has built upon and replaced older methods of analyzing construct validity such as the MTMM Matrix as described in Campbell & Fiske (1959).[4]

In confirmatory factor analysis, the researcher first develops a hypothesis about what factors he believes are underlying the measures he has used (e.g., "Depression" being the factor underlying the Beck Depression Inventory and the Hamilton Rating Scale for Depression) and may impose constraints on the model based on these a priori hypotheses. By imposing these constraints, the researcher is forcing the model to be consistent with his theory. For example, if it is posited that there are two factors accounting for the covariance in the measures, and that these factors are unrelated to one another, the researcher can create a model where the correlation between factor A and factor B is constrained to zero. Model fit measures could then be obtained to assess how well the proposed model captured the covariance between all the items or measures in the model. If the constraints the researcher has imposed on the model are inconsistent with the sample data, then the results of statistical tests of model fit will indicate a poor fit, and the model will be rejected. If the fit is poor, it may be due to some items measuring multiple factors. It might also be that some items within a factor are more related to each other than others.

For some applications, the requirement of "zero loadings" (for indicators not supposed to load on a certain factor) has been regarded as too strict. A newly developed analysis method, "exploratory structural equation modeling", specifies hypotheses about the relation between observed indicators and their supposed primary latent factors while allowing for estimation of loadings with other latent factors as well.[5]

## Statistical model

In confirmatory factor analysis researchers are typically interested in studying the degree to which responses on a p x 1 vector of observable random variables can be used to assign a value one or more unobserved variable(s) η they are designed to tap into. The investigation is largely accomplished by estimating and evaluating the loadings of each item used to tap aspects of the unobserved latent variable. That is, y[i] is the vector of observed responses predicted by the unobserved latent variable ${\displaystyle \xi }$, which is defined as:

${\displaystyle Y=\Lambda \xi +\epsilon }$,

where ${\displaystyle Y}$ is the p x 1 vector of observed random variables, ${\displaystyle \xi }$ is the unobserved latent variables, or variables in the multidimensional case, and ${\displaystyle \Lambda }$ is a p x k matrix with k equal to the number of latent variables.[6] Since, ${\displaystyle Y}$ are imperfect measures of ${\displaystyle \xi }$, the model also consists of error, ${\displaystyle \epsilon }$. Estimates in the maximum likelihood (ML) case generated by iteratively minimizing the fit function,

${\displaystyle F_{ML}=ln|\Lambda \Omega \Lambda {'}+I-diag(\Lambda \Omega \Lambda {'})|+tr(R(\Lambda \Omega \Lambda {'}+I-diag(\Lambda \Omega \Lambda {'})^{-1})-ln(R)-p}$

where ${\displaystyle \Lambda \Omega \Lambda {'}+I-diag(\Lambda \Omega \Lambda {'})}$ is the variance-covariance matrix implied by the proposed factor analysis model and ${\displaystyle R}$ is the observed variance-covariance matrix.[6] That is, values are found for freed model parameters that minimizes the difference between the mode-implied variance-covariance matrix and observed variance-covariance matrix.

## Alternative estimation strategies

Although numerous algorithms have been used to estimate CFA models, ML remains the primary estimation procedure.[7] That being said, CFA models are often applied to data conditions that deviate from the normal theory requirements for valid ML estimation. For example, social scientists often estimate CFA models with non-normal data and indicators scaled using discrete ordered categories.[8] Accordingly, alternative algorithms have been developed that attend to the diverse data conditions applied researchers encounter. The alternative estimators have been characterized into two general type: (1) robust and (2) limited information estimator.[9]

When ML is implemented with data that deviates away from the assumptions of normal theory, CFA models may produce biased parameter estimates and misleading conclusions.[10] Robust estimation typically attempts to correct the problem by adjusting the normal theory model χ2 and standard errors.[9] For example, Satorra and Bentler (1994) recommended using ML estimation in the usual way and subsequently dividing the model χ2 by a measure of the degree of multivariate kurtosis.[11] An added advantage of robust ML estimators is their availability in common SEM software (e.g., LAVAAN).[12]

Unfortunately, robust ML estimators can become untenable under common data conditions. In particular, when indicators are scaled using few response categories (e.g., disagree, neutral, agree) robust ML estimators tend to perform poorly.[10] Limited information estimators, such as weighted least squares (WLS), are likely a better choice when manifest indicators take on an ordinal form.[13] Broadly, limited information estimators attend to the ordinal indicators by using polychoric correlations to fit CFA models.[14] Polychoric correlations capture the covariance between two latent variables when only their categorized form is observed, which is achieved largely through the estimation of threshold parameters.[15]

## Exploratory factor analysis

Both exploratory factor analysis (EFA) and confirmatory factor analysis (CFA) are employed to understand shared variance of measured variables that is believed to be attributable to a factor or latent construct. Despite this similarity, however, EFA and CFA are conceptually and statistically distinct analyses.

The goal of EFA is to identify factors based on data and to maximize the amount of variance explained.[16] The researcher is not required to have any specific hypotheses about how many factors will emerge, and what items or variables these factors will comprise. If these hypotheses exist, they are not incorporated into and do not affect the results of the statistical analyses. By contrast, CFA evaluates a priori hypotheses and is largely driven by theory. CFA analyses require the researcher to hypothesize, in advance, the number of factors, whether or not these factors are correlated, and which items/measures load onto and reflect which factors.[17] As such, in contrast to exploratory factor analysis, where all loadings are free to vary, CFA allows for the explicit constraint of certain loadings to be zero.

EFA is sometimes reported in research when CFA would be a better statistical approach.[18] It has been argued that CFA can be restrictive and inappropriate when used in an exploratory fashion.[19] However, the idea that CFA is solely a “confirmatory” analysis may sometimes be misleading, as modification indices used in CFA are somewhat exploratory in nature. Modification indices show the improvement in model fit if a particular coefficient were to become unconstrained.[20] Likewise, EFA and CFA do not have to be mutually exclusive analyses; EFA has been argued to be a reasonable follow up to a poor-fitting CFA model.[21]

## Structural equation modeling

Structural equation modeling software is typically used for performing confirmatory factor analysis. LISREL,[22] EQS,[23] AMOS,[24] Mplus[25] and laavan package in R[26] are popular software programs. CFA is also frequently used as a first step to assess the proposed measurement model in a structural equation model. Many of the rules of interpretation regarding assessment of model fit and model modification in structural equation modeling apply equally to CFA. CFA is distinguished from structural equation modeling by the fact that in CFA, there are no directed arrows between latent factors. In other words, while in CFA factors are not presumed to directly cause one another, SEM often does specify particular factors and variables to be causal in nature. In the context of SEM, the CFA is often called 'the measurement model', while the relations between the latent variables (with directed arrows) are called 'the structural model'.

## Evaluating model fit

Most statistical methods only require one statistical test to determine the significance of the analyses. However, in CFA, several statistical tests are used to determine how well the model fits to the data.[16] Note that a good fit between the model and the data does not mean that the model is “correct”, or even that it explains a large proportion of the covariance. A “good model fit” only indicates that the model is plausible.[27] When reporting the results of a confirmatory factor analysis, one is urged to report: a) the proposed models, b) any modifications made, c) which measures identify each latent variable, d) correlations between latent variables, e) any other pertinent information, such as whether constraints are used.[28] With regard to selecting model fit statistics to report, one should not simply report the statistics that estimate the best fit, though this may be tempting. Though several varying opinions exist, Kline (2010) recommends reporting the Chi-squared test, the Root mean square error of approximation (RMSEA), the comparative fit index (CFI), and the standardised root mean square residual (SRMR).[1]

### Absolute fit indices

Absolute fit indices determine how well the a priori model fits, or reproduces the data.[29] Absolute fit indices include, but are not limited to, the Chi-Squared test, RMSEA, GFI, AGFI, RMR, and SRMR.[30]

#### Chi-squared test

The chi-squared test indicates the difference between observed and expected covariance matrices. Values closer to zero indicate a better fit; smaller difference between expected and observed covariance matrices.[20] Chi-squared statistics can also be used to directly compare the fit of nested models to the data. One difficulty with the chi-squared test of model fit, however, is that researchers may fail to reject an inappropriate model in small sample sizes and reject an appropriate model in large sample sizes.[20] As a result, other measures of fit have been developed.

#### Root mean square error of approximation

The root mean square error of approximation (RMSEA) avoids issues of sample size by analyzing the discrepancy between the hypothesized model, with optimally chosen parameter estimates, and the population covariance matrix.[30] The RMSEA ranges from 0 to 1, with smaller values indicating better model fit. A value of .06 or less is indicative of acceptable model fit.[31]

#### Root mean square residual and standardized root mean square residual

The root mean square residual (RMR) and standardized root mean square residual (SRMR) are the square root of the discrepancy between the sample covariance matrix and the model covariance matrix.[30] The RMR may be somewhat difficult to interpret, however, as its range is based on the scales of the indicators in the model (this becomes tricky when you have multiple indicators with varying scales; e.g., two questionnaires, one on a 0–10 scale, the other on a 1–3 scale).[1] The standardized root mean square residual removes this difficulty in interpretation, and ranges from 0 to 1, with a value of .08 or less being indicative of an acceptable model.[31]

#### Goodness of fit index and adjusted goodness of fit index

The goodness of fit index (GFI) is a measure of fit between the hypothesized model and the observed covariance matrix. The adjusted goodness of fit index (AGFI) corrects the GFI, which is affected by the number of indicators of each latent variable. The GFI and AGFI range between 0 and 1, with a value of over .9 generally indicating acceptable model fit.[32]

### Relative fit indices

Relative fit indices (also called “incremental fit indices”[33] and “comparative fit indices”[34]) compare the chi-square for the hypothesized model to one from a “null”, or “baseline” model.[29] This null model almost always contains a model in which all of the variables are uncorrelated, and as a result, has a very large chi-square (indicating poor fit).[30] Relative fit indices include the normed fit index and comparative fit index.

#### Normed fit index and non-normed fit index

The normed fit index (NFI) analyzes the discrepancy between the chi-squared value of the hypothesized model and the chi-squared value of the null model.[35] However, NFI tends to be negatively biased.[36] The non-normed fit index (NNFI; also known as the Tucker-Lewis index, as it was built on an index formed by Tucker and Lewis, in 1973[37]) resolves some of the issues of negative bias, though NNFI values may sometimes fall beyond the 0 to 1 range.[34] Values for both the NFI and NNFI should range between 0 and 1, with a cutoff of .95 or greater indicating a good model fit.[31]

#### Comparative fit index

The comparative fit index (CFI) analyzes the model fit by examining the discrepancy between the data and the hypothesized model, while adjusting for the issues of sample size inherent in the chi-squared test of model fit,[20] and the normed fit index.[34] CFI values range from 0 to 1, with larger values indicating better fit. Previously, a CFI value of .90 or larger was considered to indicate acceptable model fit.[31] However, recent studies have indicated that a value greater than .90 is needed to ensure that misspecified models are not deemed acceptable (Hu & Bentler, 1999). Thus, a CFI value of .95 or higher is presently accepted as an indicator of good fit (Hu & Bentler, 1999).

## Identification and underidentification

To estimate the parameters of a model, the model must be properly identified. That is, the number of estimated (unknown) parameters (q) must be less than or equal to the number of unique variances and covariances among the measured variables; p(p + 1)/2. This equation is known as the "t rule". If there is too little information available on which to base the parameter estimates, then the model is said to be underidentified, and model parameters cannot be estimated appropriately.[38]

## References

1. ^ a b c Kline, R. B. (2010). Principles and practice of structural equation modeling (3rd ed.). New York, New York: Guilford Press.
2. ^ Preedy, V. R., & Watson, R. R. (2009) Handbook of Disease Burdens and Quality of Life Measures. New York: Springer.
3. ^ Jöreskog, K. G. (1969). A general approach to confirmatory maximum likelihood factor analysis. Psychometrika, 34(2), 183-202.
4. ^ Campbell, D. T. & Fisk, D. W. (1959). Convergent and discriminant validation by the multitrait-multimethod matrix. Psychological Bulletin, 56, 81-105.
5. ^ Asparouhov, T. & Muthén, B. (2009). Exploratory structural equation modeling. Structural Equation Modeling, 16, 397-438
6. ^ a b Yang-Wallentin, Fan; Jöreskog, Karl G.; Luo, Hao (2010-07-13). "Confirmatory Factor Analysis of Ordinal Variables With Misspecified Models". Structural Equation Modeling: A Multidisciplinary Journal. 17 (3): 392–423. doi:10.1080/10705511.2010.489003. ISSN 1070-5511.
7. ^ Flora, David B.; Curran, Patrick J. "An Empirical Evaluation of Alternative Methods of Estimation for Confirmatory Factor Analysis With Ordinal Data.". Psychological Methods. 9 (4): 466–491. doi:10.1037/1082-989x.9.4.466.
8. ^ Millsap, Roger E.; Yun-Tein, Jenn (2004-07-01). "Assessing Factorial Invariance in Ordered-Categorical Measures". Multivariate Behavioral Research. 39 (3): 479–515. doi:10.1207/S15327906MBR3903_4. ISSN 0027-3171.
9. ^ a b Bandalos, Deborah L. (2014-01-02). "Relative Performance of Categorical Diagonally Weighted Least Squares and Robust Maximum Likelihood Estimation". Structural Equation Modeling: A Multidisciplinary Journal. 21 (1): 102–116. doi:10.1080/10705511.2014.859510. ISSN 1070-5511.
10. ^ a b Li, Cheng-Hsien (2015-07-15). "Confirmatory factor analysis with ordinal data: Comparing robust maximum likelihood and diagonally weighted least squares". Behavior Research Methods. 48 (3): 936–949. doi:10.3758/s13428-015-0619-7. ISSN 1554-3528.
11. ^ Bryant, Fred B.; Satorra, Albert (2012-07-20). "Principles and Practice of Scaled Difference Chi-Square Testing". Structural Equation Modeling: A Multidisciplinary Journal. 19 (3): 372–398. doi:10.1080/10705511.2012.687671. ISSN 1070-5511.
12. ^
13. ^ Rhemtulla, Mijke; Brosseau-Liard, Patricia É.; Savalei, Victoria. "When can categorical variables be treated as continuous? A comparison of robust continuous and categorical SEM estimation methods under suboptimal conditions.". Psychological Methods. 17 (3): 354–373. doi:10.1037/a0029315.
14. ^ Yang-Wallentin, Fan; Jöreskog, Karl G.; Luo, Hao (2010-07-13). "Confirmatory Factor Analysis of Ordinal Variables With Misspecified Models". Structural Equation Modeling: A Multidisciplinary Journal. 17 (3): 392–423. doi:10.1080/10705511.2010.489003. ISSN 1070-5511.
15. ^ Olsson, Ulf. "Maximum likelihood estimation of the polychoric correlation coefficient". Psychometrika. 44 (4): 443–460. doi:10.1007/BF02296207. ISSN 0033-3123.
16. ^ a b Suhr, D. D. (2006) - “Exploratory or confirmatory factor analysis?” in Statistics and Data Analysis, 31, Retrieved April 20, 2012, from http://www2.sas.com/proceedings/sugi31/200-31.pdf
17. ^ Thompson, B. (2004). Exploratory and confirmatory factor analysis: Understanding concepts and applications. Washington, DC, US: American Psychological Association.
18. ^ Levine, T. R. (2005). Confirmatory factor analysis and scale validation in communication research. Communication Research Reports, 22(4), 335-338.
19. ^ Browne, M. W. (2001). An overview of analytic rotation in exploratory factor analysis. Multivariate Behavioral Research, 36, 111-150.
20. ^ a b c d Gatignon, H. (2010). Confirmatory Factor Analysis in Statistical analysis of management data. DOI: 10.1007/978-1-4419-1270-1_4
21. ^ Schmitt, T. A. (2011). Current methodological considerations in exploratory and confirmatory factor analysis. Journal of Psychoeducational Assessment, 29(4), 304-321.
22. ^ CFA with LISREL
23. ^ Byrne, B. M. (2006). Structural equation modeling with EQS: Basic concepts, application, and programming. New Jersey: Lawrence Elbaum Associates.
24. ^ CFA using AMOS
25. ^ Mplus homepage
26. ^ [1]
27. ^ Schermelleh-Engel, K.,Moosbrugger, H., & Müller, H. (2003). Evaluating the fit of structural equation models: Tests of significance and descriptive goodness-of-fit measures, Methods of Psychological Research Online, 8(2), 23-74
28. ^ Jackson, D. L., Gillaspy, J. A., & Purc-Stephenson, R. (2009). Reporting practices in confirmatory factor analysis: An overview and some recommendations. Psychological Methods, 14(1), 6-23.
29. ^ a b McDonald, R. P., & Ho, M. H. R. (2002). Principles and practice in reporting statistical equation analyses. Psychological Methods, 7(1), 64-82
30. ^ a b c d Hooper, D., Coughlan, J., & Mullen, M.R. (2008). Structural equation modelling: Guidelines for determining model fit. Journal of Business Research Methods, 6, 53–60
31. ^ a b c d Hu, L., & Bentler, P. M. (1999). Cutoff criteria for fit indexes in covariance structure analysis: Conventional criteria versus new alternatives. Structural Equation Modeling, 6(1), 1-55.
32. ^ Baumgartner, H., & Hombur, C. (1996). Applications of structural equation modeling in marketing and consumer research: A review. International Journal of Research in Marketing, 13, 139-161.
33. ^ Tanaka, J. S. (1993). Multifaceted conceptions of fit in structure equation models. In K. A. Bollen & J.S. Long (Eds.), Testing structural equation models (pp. 136-162). Newbury Park, CA: Sage.
34. ^ a b c Bentler, P. M. (1990). Comparative fit indexes in structural models. Psychological Bulletin, 107(2), 238-46.
35. ^ Bentler, P. M., & Bonett, D. G. (1980). Significance tests and goodness of fit in the analysis of covariance structures. Psychological Bulletin, 88, 588-606.
36. ^ . Bentler, P. M. (1990). Comparative fit indexes in structural models. Psychological Bulletin, 107(2), 238-46.
37. ^ Tucker, L. R., & Lewis, C. (1973). A reliability coefficient for maximum likelihood factor analysis. Psychometrika, 38, 1-10.
38. ^ Babyak, M. A., & Green, S. B. (2010). Confirmatory factor analysis: An introduction for psychosomatic medicine researchers. Psychosomatic Medicine, 72, 587-597.