Next Article in Journal
How Anxious and Stressed Are Emerging Adults? The Role of Mindfulness and Intolerance of Uncertainty
Previous Article in Journal
Editorial to the Special Issue “Feature Papers in Psychometrics and Educational Measurement”
 
 
Please note that, as of 22 March 2024, Psych has been renamed to Psychology International and is now published here.
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Comparison of Methods for Determining the Number of Factors to Retain with Exploratory Factor Analysis of Dichotomous Data

Department of Educational Psychology, Ball State University, Muncie, IN47306, USA
Psych 2023, 5(3), 1004-1018; https://doi.org/10.3390/psych5030067
Submission received: 9 August 2023 / Revised: 5 September 2023 / Accepted: 7 September 2023 / Published: 13 September 2023
(This article belongs to the Section Psychometrics and Educational Measurement)

Abstract

:
Exploratory factor analysis (EFA) is a very widely used statistical procedure in the social and behavioral sciences. This technique features in validity studies, as well as investigations of latent structure underlying observed measurements. A primary aspect of using EFA is determining the number of factors to retain. In addition to theoretical considerations, a variety of statistical tools have been developed and recommended for use in assisting researchers with respect to factor retention. Some research has been conducted to investigate the accuracy of these methods in the case of continuous factor indicators. The purpose of the current simulation study was to extend this earlier work to situations in which the indicator variables are dichotomous, as with questionnaire or test items. Results of this study revealed that an approach based on the combined results of the empirical Kaiser criterion, comparative data, and Hull methods, as well as Gorsuch’s CNG scree plot test by itself, all yielded accurate results with respect to the number of factors to retain. Implications for practice are discussed.

1. Introduction

Exploratory factor analysis (EFA) is a widely used tool in the psychological sciences. For example, in a recent review of the Psycinfo search engine for the 10 years between 31 July 2013 and 31 July 2023, using the keyword “exploratory factor analysis”, a total of 560,312 separate journal articles were identified. Given its ubiquity and importance in assessing the viability of psychological theories, the appropriate application of EFA by researchers is of key importance in order to ensure that statistical results are as accurate as possible. Perhaps the central decision to be made by researchers using EFA is the number of factors to retain [1]. Many methods exist for determining the number of factors to retain when using EFA, and a number of studies have compared these methods in the presence of normally distributed observed indicator variables; e.g., [2]. However, in many psychological and educational applications, the observed indicator variables associated with the use of EFA are dichotomous in nature e.g., [3,4]. For example, popular scales of mood [5] and anxiety [6] rely on items that force respondents to select whether particular symptoms are present or absent. Likewise, personality inventories [7], academic admissions tests [8], and numerous tests of cognitive ability [9] consist of items with dichotomous responses. Therefore, it is of key importance for researchers and data analysts to have insights regarding optimal techniques for determining the number of factors to retain in the case of categorical, and in particular dichotomous, indicator variables. Although research has been conducted to examine the performance of techniques for determining the number of factors to retain in EFA, these have focused on continuous indicator variables [2]. Given the aforementioned ubiquity of dichotomous items in in social science research, it is worthwhile to investigate the performance of factor determination methods for such scales.
The purpose of this study was to build on earlier work [2,10] by examining a wide range of statistics for determining the number of factors to retain in the context of dichotomous observed indicator variables. The prior research in this area has focused primarily on approaches for determining the number of factors to retain when the indicator variables are continuous [2]. While certainly informative, this earlier work does not pertain directly to the case where indicators consist of dichotomous variables, such as item responses, which are common in the social sciences, as discussed above. There has been some strong recent work examining the performance of parallel analysis and revised parallel analysis in the context of EFA with dichotomous items [10] However, this work has not included the full range of statistical tools available to determine the number of factors to retain, beyond the parallel analysis framework. The current study extends on this earlier work by comparing the full array of methods heretofore used with continuous indicators with dichotomous items.
The remainder of the manuscript is structured as follows: First is a brief review of the one-factor model, followed by a description of the various factor retention methods examined in this study. Next is a review of prior research focused on factor retention techniques with continuous and categorical indicators. The simulation methodology used in this study is then described, followed by the results, and finally a discussion of the findings from the current study in light of prior research.

2. Exploratory Factor Analysis

The EFA model can be expressed as:
Y = Λ ξ + Ψ
where
  • Y = matrix of observed indicator variables
  • ξ = matrix of factor(s)
  • Λ = matrix of factor loadings relating indicators to factor(s)
  • Ψ = matrix of unique random errors associated with the observed indicators
Of particular interest in the context of EFA is the factor loading matrix, Λ , which links the observed indicator variables to the latent factors. Indicator factor combinations with non-trivial loadings, e.g., >0.3 [11], are said to be associated with one another. Based on the indicator groupings, researchers can make inferences regarding the nature of the latent traits underlying the observed indicators. A number of methods are available for estimating the model parameters in Equation (1), a step often referred to as factor extraction. Perhaps the most widely used of these methods are maximum likelihood estimation (MLE) and principal axis factoring (PAF). After the initial factor loadings are extracted, they are typically transformed (i.e., rotated) in order to improve interpretability of the results through the mathematical encouragement of simple structure in Λ ; i.e., each indicator will be associated with a single factor.

3. Determining the Number of Factors to Retain

EFA is, by definition, an exploratory statistical technique. As such, researchers using it do not explicitly link individual indicator variables to specific latent factors. In addition, typical EFA practice involves an exploration of multiple possible factor solutions, with the optimal such approach being selected as the “winner”. A key aspect of selecting the number of factors to retain is matching the factor loading pattern with theory. The ultimate decision regarding the EFA results to retain should be conceptually coherent and theoretically sound with respect to how the indicators group together with the factors [12]. In addition, researchers can make use of several statistically based techniques for ascertaining the number of factors to retain. These methods do not refer to theory, but rather provide insights regarding the number of factors to retain based upon statistical criteria. Although these purely statistical approaches cannot be in isolation to determine the number of factors to retain, they are very useful, in conjunction with theoretical considerations, to help researchers ascertain the number of factors to retain. Following is a description of several such techniques that have been shown in prior research to be useful tools for determining the number of factors to retain.

3.1. Chi-Square Test

One of the more widely used approach for determining the number of factors to retain is the chi-square statistic associated with MLE. MLE utilizes a minimization process in which model parameter estimates (i.e., factor loadings, variances, and error variances) that yield a predicted covariance matrix ( Σ ) of the indicator variables that is very close in value to the observed covariance matrix (S). The value of the function that results from this minimization process can be converted to a chi-square statistic that can be used to test the null hypothesis that Σ = S ; i.e., the model predicted covariance matrix is equal to the actual covariance matrix. This null hypothesis states that the fit of the model is perfect, which is actually rarely achieved in practice [13,14]. Therefore, many close-fitting models that prove useful may be rejected by the chi-square goodness of fit test, thereby leading a researcher to conclude that the model does not provide a good fit to the data, when in fact the proposed model does provide reasonably good data model fit. For this reason, this test may not be particularly useful for assessing the fit of an individual model [14,15].
Despite the limited utility of the MLE chi-square test for assessing the feasibility of a single model, it has been shown to be an effective tool for comparing EFA models with differing numbers of factors. This approach involves obtaining the chi-square statistic for each of several EFA results (e.g., 1 factor, 2 factors, 3 factors) and then calculating the difference in values between adjacent factor models. These differences follow the chi-square distribution if the observed indicators come from a multivariate normal distribution [16]. Thus, the chi-square difference statistic can be used to test the null hypothesis of equivalent model fit between a pair of EFA models. As an example, the researcher might fit EFA models with 1 to 4 factors. Then, the differences in chi-square fit statistic values between 1 and 2, 2 and 3, and 3 and 4 factors are calculated. If the 1 vs. 2, and 2 vs. 3 chi-square differences values are statistically significant, but the 3 vs. 4 are not, then the researcher would retain 3 factors.
This approach, commonly referred to as the sequential model test (SMT), has been shown to be an effective tool when the assumption of multivariate normality is met, and particularly when factors are relatively highly correlated [2]. However, when data do not follow a multivariate normal distribution, the population model is mis-specified, and/or minor factors are present in the population, SMT has tendency to overextract; i.e., recommend the retention of too many factors [2,17]. SMT was included in the current study in order to serve as a baseline approach, given its popularity in applied settings. However, given that the data used here were dichotomous, thus violating the assumption of multivariate normality, it was not anticipated that this approach would be particularly effective in determining the number of factors to retain.

3.2. Parallel Analysis

Another inferential approach for determining the number of factors to retain is parallel analysis (PA). This method was first described by Horn [18], and involves the generation of synthetic data that has the same marginal properties (i.e., means and variances) as the actual observed data, but which has no underlying latent structure; i.e., 0 factors. Most often, the synthetic data are generated nonparametrically through the random mixing of indicator variable values within variables across observations. Typically, a large number (e.g., 1000) of these synthetic datasets is generated, for each of which EFA is conducted. The eigenvalues from these analyses, which reflect the explained variance associated with each factor, are then used to create distributions of eigenvalues that would be expected if no factor structure is present. The eigenvalues obtained using the observed data are then compared to these distributions in order to determine the number of factors to retain. A factor is retained if its observed eigenvalue is larger than the 95th percentile of the distribution of null factor eigenvalues generated from the random data. For example, if the observed first, second, and third eigenvalues were 4.5, 2.1, and 0.9 and the 95th percentiles of the first, second, and third eigenvalues were 1.4, 1.2, and 1, the researcher would retain two factors. Researchers [17] extended the use of PA based on the tetrachoric correlation matrix of dichotomous item response data and principal components analysis. Given that the data used in this study were dichotomous item responses, this approach based on the tetrachoric correlation was used in the current study. The tetrachoric correlation has been shown to be most appropriate for estimating correlation coefficients with dichotomous variables [19].

3.3. Revised Parallel Analysis

Researchers have pointed out that standard PA may have some limitations with respect to assessing the magnitude of all but the largest eigenvalues produced by EFA [20]. The rationale behind this argument is that it is only for the assessment of the first eigenvalue that the underlying latent structure is assumed to include no factors. If the first eigenvalue for the observed data exceeds the 95th percentile of the reference distribution, the presence of one factor can be inferred. Thus, the second observed eigenvalue should be compared with a null distribution assuming a single latent trait. Likewise, the third observed eigenvalue should be compared with a null distribution assuming two latent traits, and so on. In order to address this issue, [20] proposed the use of what they called revised PA (RPA), in which the comparison data are generated assuming that k − 1 latent variables underlie the observed data, rather than 0 factors as is the case with traditional PA. Thus, when testing for 4 factors, RPA would generate data from a 3-factor model, holding the marginal distributions of the variables to be equal to those in the observed sample. In other respects, the methodology for RPA is the same as for standard PA.

3.4. Comparative Data Method

Concurrent with the work by [20,21] described the comparative data method (CD) for determining the number of factors to retain. CD generates a sample of 10,000 cases based on a correlation matrix associated with k − 1 latent variables. Next, 500 random samples of n observations (where n is equal to the observed data sample size) are drawn from the 10,000 cases, and principal components analysis (PCA) is conducted in order to obtain eigenvalues. Inference regarding the number of factors to retain is then carried out as with PA and RPA. For PA, RPA, and CD in this study, the tetrachoric correlation was used.

3.5. Minimum Average Partial

Velicer [22] proposed a method for determining the number of factors to retain that is based upon an examination of the average squared partial correlations among the observed indicators, accounting for the influence of the latent variables. This method, called minimum average partial (MAP), is carried out using a multiple-step procedure. In the first step, the correlations among the observed variables are calculated, squared, and then the squares are averaged. In the second step, the squared correlations among the indicators are again calculated and averaged, in this case after partialing out the first latent variable obtained using EFA. In the third step, the average squared correlation among the observed variables is again calculated, this time partialing out the first two latent variables. These steps are repeated for the first p − 1 factors, where p is the number of observed indicators. The researcher then would retain the number of factors corresponding to the minimum average squared partial correlation, as this corresponds to the point at which the maximum amount of systematic variance in the observed indicators is accounted for by the latent variables. Thus, if the average squared partial correlation values were 0.2, 0.15, 0.09, 0.03, 0.07, 0.11, and 0.12, the researcher would select the 4-factor solution. Simulation research has consistently demonstrated that MAP is one of the more accurate methods for determining the number of factors to retain [21,23,24,25].

3.6. Empirical Scree Plot Methods

A plot of the eigenvalues by the factor number, known as the scree plot, was proposed by [26] determining the optimal EFA solution. The purpose of this plot is to examine the relationship between the number of factors and the amount of variance explained in the observed variables (as measured by the eigenvalues), with an eye toward identifying for which factor number the amount of explained variance declines sharply. In practice, this is done using a scatterplot with the eigenvalue on the y-axis, the factor number on the x-axis, and a line connecting points in the plot. As noted above, the eigenvalues decrease in value from the first through the last factor, which is reflected in the scatterplot. The researcher using this approach examines the plot, looking for the point where the line connecting the eigenvalues begins to flatten out in its rate of decline. This point will correspond to the number of factors that should be retained.
In practice, the scree plot has been found to be relatively inaccurate with respect to the number of factors to retain, in large part perhaps due to the highly subjective interpretive criterion [27]. Thus, multiple more objective approaches based on the scree plot have been discussed in the literature. Gorsuch’s [1] CNG scree test involves the calculation of the slope linking the first three eigenvalues, then the calculation of the slope linking eigenvalues 2, 3, and 4, then the slope linking eigenvalues 3, 4, and 5, and so on. The researcher then compares these slopes with one another, and selects the number of factors where the difference between the slopes is greatest. Thus, for example, if the largest difference between slope values lies between the line for points 2, 3, and 4 versus the line for points 3, 4, and 5, we would retain 4 factors. Zoski and Jurs [28] suggested a variant (NMREG) of the Gorsuch approach in which pairs of regression equations are estimated using all of the data points, rather than just sets of 3 at a time. Thus, for p indicator variables, the following pairs of equations would be considered:
Line 1 (eigenvalues 1, 2, and 3)     Line 2 (eigenvalues 4 through p)
Line 3 (eigenvalues 1, 2, 3, and 4)     Line 4 (eigenvalues 5 through p)
Line 5 (eigenvalues 1, 2, 3, 4, and 5)     Line 6 (eigenvalues 6 through p)
The slopes for the lines in each pair (e.g., line 1 versus line 2) are then compared using a t-test, and the number of factors to be retained is associated with the maximum t value. As an example, if the maximum t statistic is associated with the comparison between lines 3 and 4, then 4 factors (corresponding the largest factor number in line 3) would be retained.

3.7. Hull Method

Lorenzo-Seva, et al. [29] described an approach for determining the number of factors to retain from an EFA that is related to the scree plot approaches described above. The Hull technique is based on model goodness of fit test results using the following steps.
  • Fit EFA models with varying numbers of factors (e.g., 1, 2, 3, etc.).
  • For each model from step 1, calculate a goodness of fit statistic ( G m and model degrees of freedom d f m ), such as the CFI.
  • Compare the values of G m across factor models. A model is rejected if its value of G m is lower than that of an adjacent model with fewer factors.
  • Models from those remaining after step 3 are rejected as viable if their G m value lies below a line connecting points in a plot of G m and d f m .
  • Repeat step 4 until no nonviable solutions are remaining.
  • Select the number of factors, m, for which the following statistic is maximized from among the set of potentially viable solutions:
G m G m 1 / d f m d f m 1 G m + 1 G m / d f m + 1 d f m
Prior research has shown that the Hull method using the CFI as G m is an effective tool for identifying the number of factors with continuous indicators [2,29], particularly when it is used in combination with other statistics. This is the approach that was employed in the current study.

3.8. Empirical Kaiser Criterion

Braeken and van Assen [30] described a method for determining the number of factors based upon the eigenvalues, thereby distinguishing it from the Hull approach, and placing it in the broad family including the objective scree plot techniques and PA/RPA. This empirical Kaiser criterion (EKC) is an outgrowth of the familiar Kaiser criterion, in which all factors with eigenvalues greater than 1 are retained. A problem with the original Kaiser approach is that it does not account for sampling variability inherent in analyses that do not involve the entire population [27]. The EKC method is based upon an assumption that the distribution of eigenvalues follows the Marcenko-Pastur distribution [31]. The upper bound of this distribution is calculated as:
λ 1 , r e f = 1 + p n 2
where
  • p = number of indicators
  • n = sample size
The largest eigenvalue from the EFA (corresponding to the first factor) is compared to this maximum value. The remaining eigenvalues are compared to the variance corrected reference eigenvalues as calculated below:
λ m , r e f = m a x p i = 0 m 1 λ m p m + 1 1 + p n 2 , 1
The number of factors to be retained corresponds to the number of eigenvalues that are greater than the reference value in Equation (4). For example, consider the case where the first 4 observed eigenvalues are 3.5, 2.4, 1.7, and 1.1. If there are 20 indicators and a sample size of 200, the upper bound from Equation (3) is 1.73. Thus, the first factor would be retained. If the subsequent reference eigenvalues calculated from Equation (4) were 1.5, 1.4, and 1.25, the researcher would retain a total of 3 factors.

3.9. Combination Approach

In addition to the individual statistics described above, a combination approach to determining the number of factors to retain was also investigated in the current study. This combination technique included a vote-counting approach to EKC, Hull, and CD, per Auerswald and Moshagen [2]. Specifically, the number of factors to retain corresponded to the largest vote count from among the methods. If multiple numbers of factors received the same number of votes, the most parsimonious model from the vote winners was selected. It is important to note that the prior authors also included SMT in their combination approach. However, given that the dichotomous data used in the current study do not conform to the multivariate normal assumption underlying MLE, the SMT approach was excluded from the combination method utilized here.

3.10. Limited Information Item Factor Analysis

McDonald [32] introduced the Normal-Ogive Harmonic Analysis Robust Method (NOHARM) approach to fitting factor analysis models to dichotomous item response data. The NOHARM model takes the form:
P U j i = 1 | θ = N β j 0 + β j 1 θ 1 i + + β j k θ k i
where
  • U j i = response to item j for respondent i; 1 = item endorsement
  • θ = vector of latent traits
  • θ k i = level on latent trait k for respondent i
  • β j 0 = intercept for item j
  • β j k = coefficient relating item j to trait k
Parameter estimates for the model in Equation (5) can be obtained using a limited information approach based on unweighted least squares (ULS).
Fit of the NOHARM model to the data can be assessed using a chi-square statistic, as described by Gessaroli and De Champlain [33]:
χ G D 2 = N 3 l = 2 J j = 1 J 1 z j l ( r )
where
  • N = total sample size
  • z j l ( r ) = Fisher’s z-transformed residual correlation, r, for item pair jl
In order to determine the number of factors using χ G D 2 , a sequence of models is fit to the data, differing based upon the number of factors, much as with SMT. For each such model, χ G D 2 is calculated and the difference between this statistic for models with k and k − 1 number of factors is calculated. This Δ χ G D 2 statistic follows an approximate χ 2 distribution [33] with the degrees of freedom equal to the difference in degrees of freedom for the two models. If the Δ χ G D 2 is statistically significant, the researcher concludes that the EFA model with k factors provides better fit than the k − 1 factor. The k factor solution is then compared to the k + 1 solution, and the procedure is repeated until a non-statistically significant conclusion is achieved, at which point the more parsimonious solution in the comparison (e.g., k vs. k + 1) is selected. For example, if the Δ χ G D 2 values comparing 1 vs. 2 and 2 vs. 3 factors were statistically significant, but the value comparing 3 vs. 4 factors was not, the researcher would retain 3 factors.

3.11. Full Information Item Factor Analysis

An alternative to modeling dichotomous item response data with more than one latent trait comes in the form of the full information multidimensional IRT (MIRT) model. The 2-parameter logistic (2PL) MIRT model (which corresponds to the standard factor model) for dichotomous data takes the form:
P U j i = 1 | θ = e α j 1 θ 1 j + + α j k θ k i + γ j 1 + e α j 1 θ 1 j + + α j k θ k i + γ j
where
  • α j k = discrimination parameter for item j on latent trait k
  • γ j = difficulty for item j
Maximum likelihood can be used to obtain the parameters in Equation (7) based on the EM algorithm [34]. Alternative estimation approaches such as either a Metropolis–Hastings Robbins–Monro or Bayes can also be used in conjunction with the MIRT model. The EM algorithm was used in the current study. As with the standard EFA and NOHARM models, an exploratory implementation of the MIRT model can be used when the researcher is unsure as to the number of underlying latent traits for a given dataset. Further details regarding the implementation of the exploratory MIRT model in the R software environment can be found in Chalmers [26].
Determination of the number of factors to retain when using the full information MIRT model involves use of the likelihood ratio test, in much the same way as the SMT and Δ χ G D 2 described above. For each prospective model, the exploratory MIRT model is fit to the data and the chi-square goodness of fit statistic is obtained. As with SMT and NOHARM, the difference between these statistics is compared to the chi-square distribution with degrees of freedom equal to the difference in the degrees of freedom for the two models. An approach similar to that described above for NOHARM and the SMT can be used with exploratory MIRT in order to determine the optimal number of factors to retain. As an example, if the tests comparing 1 vs. 2 factors and 2 vs. 3 factors were statistically significant, but the test comparing 3 vs. 4 factors was not, the researcher would retain 3 factors from the EFA.

3.12. Prior Research Investigating Factor Retention Methods

There is a relatively large body of literature investigating the performance of multiple methods for identifying the number of factors to retain in the context of EFA. This work has been particularly focused on the case of continuous (and primarily multivariate normal) indicator variables. For example, research has demonstrated that PA and its variants are among the most effective methods for identifying the number of factors to retain in the context of EFA [35,36,37,38,39]. Early simulation results consistently demonstrated that PA tended to identify the correct number of factors underlying a set of observed indicators more frequently than did most other alternative approaches, such as the scree plot, the eigenvalue greater than 1 rule, and proportion of variance explained by the factors. More recent work has shown that RPA and CDM may yield somewhat more accurate results than standard PA, and that RPA is perhaps the best performer most consistently across a variety of conditions [40]. Auerswald and Moshagen [2] found that PA based on PCA (rather than the common factor model) yielded more accurate results than did RPA based on the common factor model. Similarly, Guo and Choi [10] found that PA and RPA based on PCA was able to accurately identify the number of factors to retain when the indicators were dichotomous in nature, provided that the sample is 500 or larger. Lim and Jahng [38] conducted a simulation study comparing various versions of PA and found that standard PA yielded accurate results to within 1 factor (plus or minus) of the actual number of latent traits underlying the observed data.
With respect to the objective methods associated with the scree plot, simulation research has shown that the Zoski and Jurs approach and the CNG scree test are both very effective at determining the number of factors to retain, assuming that there are at least 3 latent variables present in the data ([27]. This work was limited to the case of multivariate normal indicator variables. Thus, the current research was designed to extend upon these earlier results by applying the objective scree plot techniques to the case of dichotomous indicator variables.
Auerswald and Moshagen [2] compared the performance of several criteria for determining the number of factors to retain when using EFA. They concluded that using results from a combination of statistics, rather than a single value, might prove to be the most effective. Based upon their simulation study, Auerswald and Moshagen recommended that when researchers fit EFA models with continuous indicator variables, the results of the SMT and either RPA, Hull, or EKC agree on the number of factors to retain, they settle on this value. However, if the SMT and the other technique do not suggest the same number of factors, then the researcher should rely primarily on CD, EKC, or RPA solely. Finally, their results demonstrated that using a combination of approaches typically yields more accurate results than does any single approach.
With respect to χ G D 2 , Finch and Habing [41,42] showed that a variant of Δ χ G D 2 performed well, particularly with more factors present in the population. This approach was more effective than alternatives when the underlying latent trait was normally distributed, whereas when the trait was skewed, nonparametric alternatives were better able to control the Type I error rate [41]. Svetina and Levy [43] extended this work and found that the Δ χ G D 2 statistic was more accurate than nonparametric alternatives, particularly for data when some indicators were associated with multiple factors, the test was longer, and the factors were more highly correlated.

3.13. Study Goals

The primary goal of this simulation study was to compare several methods for determining the number of factors to retain in the context of EFA with dichotomous indicator variable data. This study was designed to expand on earlier work, particularly that of Guo and Choi [10] and Auerswald and Moshagen [2]. As discussed above, prior work with continuous indicator variables indicated that using a combination of methods may be the optimal approach for this purpose [2]. In addition, Guo and Choi showed that using PA based on the tetrachoric correlation matrix and PCA can provide accurate results for dichotomous item response data. With respect to specific techniques, CD, EKC, or RPA may yield the most accurate results with respect to the number of factors to retain, based on prior research. The current study builds upon this earlier work by investigating the factor retention problem when indicators are dichotomous in nature using a wide array of possible methods for this purpose. Several of the methods examined in earlier research are included in this study, as are approaches designed specifically for categorical indicator variables, including NOHARM and the MIRT model. In addition, methods that have been found to be effective for a larger number of indicators than is present in the current work (i.e., [10]) were also included in this study. Thus, the current study was designed to extend two recent studies on EFA and factor retention with a focus on dichotomous indicators and a large array of methods. This study provides direct comparisons among a wide range of methods that have not been compared in this way before. The techniques are compared with respect to the proportion of instances in which the correct number of factors is recommended by each method, as well as the mean number of factors suggested for retention.

4. Methods

In order to address the goals of this study, a Monte Carlo simulation study was conducted. For each combination of conditions (described below), 1000 replications were used. The simulations were carried out with the R software system, version 4.1.2 [44]. Dichotomous item response data were simulated based on a structural equation model with either 1 or 3 latent traits, using the MonteCarloSEM R package [45]. The item thresholds used in data generation were drawn from a standardized mathematics test, with a mean of 0.04 and ranging between −2.3 and 2.7. The tetrachoric correlation matrix was estimated for the set of dichotomous indicators and then served as the dataset for each of the factor retention methods included in the study.
In addition to the methods used to determine the number of factors to retain, other manipulated factors included sample size, the number of underlying latent traits, interfactor correlation, the number of indicators per factor, and the factor loading value. The sample sizes used in this study were 200, 500, 1000, and 2000. These values were selected based on prior research in the area of factor analysis [2,10,46]. Data were generated either from a 1 or a 3-factor model, allowing for assessment of the factor retention methods for both unidimensional and multidimensional data. There were 5, 10, or 20 indicator variables per factor, for each of the factor number conditions, reflecting what might be seen with short subscales (5 items) to relatively long subscales (20 items). Results for the 10 and 20 indicator conditions were found to be virtually identical and therefore, in order to simplify the presentation of results somewhat, only those for 5 and 10 indicators are reported below. The factor loadings were set at either 0.6 or 0.8, which, when combined with the number of indicators, resulted in scale reliability values between 0.74 and 0.97. In addition, these loading values were selected in order to examine performance of the factor retention methods when more than 50% of the variance in each indicator was associated with the latent trait (loading = 0.8) and when less than 50% of the indicator variance was explained by the factor (loading = 0.6). Three conditions were included for interfactor correlations, including 0.25, 0.50, and 0.75. These values were selected so as to reflect small, medium, and large relationships among the factors [47].
The methods used to determine the number of factors to retain included EKC, Hull, SMT, CD, Combined, MAP, NMREG, CNG, PA, RPA, MIRT, and NOHARM. These methods were applied using the R packages ‘nFactors’ [48], ‘EFA.dimensions’ [49], ‘mirt’ [50], and ‘sirt’ [51]. The outcomes of interest were the proportion of replications for which each method correctly identified the number of factors to retain and the mean number of factors retained across the 1000 replications for each combination of conditions. In order to identify which combinations of the manipulated conditions were associated with the correct proportion of factors retained, analysis of variance (ANOVA) was used. For each combination of study conditions, the proportion of correct cases was calculated and served as the dependent variable in the ANOVA. In addition to statistical significance, the eta-squared effect size was also calculated for each term in the model. ANOVA models were fit separately for the 1 and 3-factor cases, due to the difference in the terms available for each.

5. Results

5.1. Three Factors

5.1.1. Factor Retention Accuracy Rate

For the case when 3 factors were present in the population, ANOVA identified the interactions of the factor determination method, number of indicators, and interfactor correlation ( F 22,286 = 7.03 ,   η 2 = 0.79 ), and the factor determination method, number of indicators, and sample size ( F 66,286 = 2.18 ,   η 2 = 0.33 ) as being statistically significantly related to the proportion of replications with the correct number of factors being retained. Table 1 includes the proportion of replications with the correct number of factors by method, number of indicators, and interfactor correlations, in the 3-factors case. When there were 10 indicators per factor, EKC and Combined exhibited the highest proportion of correct number of factors retained, across interfactor correlation levels. In addition, CD, CNG and MAP had rates of correct identification of 1.00 for an interfactor correlation of 0.25 and 0.50, but rates below 0.90 when factors had a correlation of 0.75. Of these latter techniques, CD had the highest accuracy rate for the highest interfactor correlation.
When each factor had 5 indicators, the accuracy rates were generally lower for all methods across interfactor correlation conditions. The only exception to this pattern was CNG at correlations of 0.25 or 0.50, for which the accuracy rate was 1.00. In addition, EKC, CD and Combined had accuracy rates of 0.99 for a correlation of 0.25. As the interfactor correlation increased in value, the accuracy rates for all of the methods studied here declined in the 5-indicator condition. At the highest interfactor correlation value, MIRT, NOHARM, and CNG had the highest accuracy rates in the 5-indicators case. The SMT had the lowest rates across all conditions in Table 1.
The proportion of replications with the correct number of factors identified by number of indicators, sample size, and method appear in Table 2. Across methods, the rate of correct factor number identification was higher for larger samples and more indicators. For 5 indicators, the highest accuracy rates were associated with CNG, except for when N = 2000. In this latter case, PA had an accuracy rate of 1.00, with CD being the second most accurate, followed by CNG. When there were 10 indicators per factor, EKC and the Combined method had the highest accuracy rates. Across conditions, SMT had the lowest rates of accurately determining the number of factors to retain. Finally, when there were 3 factors, each method studied here had higher accuracy rates for retaining the number of factors when the loadings were 0.8 as opposed to 0.6 (Table 3). This result was most marked for MAP and CD, and least so for MIRT and NOHARM.

5.1.2. Mean Number of Factors Retained

The number of factors retained by method, number of indicators, and interfactor correlation in the 3-factor condition appear in Table 4. Across conditions, SMT, PA, and RPA all overestimated the number of factors to retain. In addition, NMREG also overestimated the number of factors present when the interfactor correlation was 0.5 or 0.75. Among the other methods, the mean number of factors retained was closest to the actual value of 3 for EKC, CD, Combined, CNG, MAP, MIRT, and NOHARM in the 10-indicator case, particularly with correlations below 0.75. At this highest correlation value, EKC, Combined, MIRT, and NOHARM were closest to the actual data generating factor number. With 5 indicators, the interfactor correlation had a greater impact on the performance of the methods studied here. When the correlation was 0.25 or 0.50, CNG suggested retaining 3 factors on average, whereas EKC, CD, and Combined all had a mean number of factors retained of 2.8 or higher. At the highest interfactor correlation, EKC, CD, and Combined all suggested that fewer factors be retained, with means of 2 or fewer. In contrast, CNG had a mean of 3.4 factors with 5 indicators and an interfactor correlation of 0.75, suggesting that it tended to overfactor when the latent traits were most strongly related to one another. Finally, in the 5-indicators condition, MIRT tended to overfactor and NOHARM tended to underfactor.
The results in Table 5 show the mean number of factors retained for each method by the number of indicators and the sample size. Many of these patterns are similar to those in Table 4 and thus will not be repeated here. Of particular import in this table is that the mean number of factors retained was closest to the population value of 3 with larger samples for EKC, CD, and Combined. In addition, the mean number of factors for PA was approximately 3 for sample size of 2000, regardless of the number of indicators. The method with the mean number of factors consistently closest to the population value of 3 across conditions was CNG. The mean number of latent traits retained by method and factor loadings for the 3-factor condition appear in Table 6. As noted above, the mean number of factors to retain was closest to the population value of 3 for CNG across factor loading values, with MIRT having the next closest mean number of factors.

5.2. 1-Factor

5.2.1. Factor Retention Accuracy Rate

For the 1-factor condition, ANOVA identified the interaction of factor retention method by the number of indicators per factor and sample size ( F 66,286 = 2.91 ,     η 2 = 0.54 ), and method by factor loading magnitude ( F 11,286 = 3.01 ,   η 2 = 0.17 ) as being statistically significantly related to the factor retention accuracy rate. A number of methods exhibited accurate rates of 1.00 across sample size and number of indicators, including EKC, CD, Combined, MAP, and CNG. In contrast, SMT, NMREG, and RPA had the lowest accuracy rates across number of indicators and sample sizes. Accuracy rates were generally higher in the 10-indicator condition, particularly for MIRT and NOHARM. The accuracy rates by factor loading for the 1-factor case appears in Table 7. For the most accurate methods, EKC, CD, MAP, CNG, and Combined, the loading values were not associated with the proportion of cases for which 1 factor was correctly identified. Among the other methods, accuracy for PA and NOHARM was most strongly impacted by loading value, where accuracy was markedly higher for the 0.80 condition.

5.2.2. Mean Number of Factors Retained

The mean number of factors retained by factor loading in the 1-factor case appears in Table 8. EKC, CD, MAP, CNG, and Combined had mean number of factors retained of 1, matching the population value of 1. In addition, the Hull technique had a mean just above 1. The mean number of factors retained by method, sample size, and method in the 1-factor case appear in Table 8. In keeping with the patterns described above, EKC, CD, MAP, CNG, and Combined had a mean number of factors retained of 1 across conditions, with the mean for Hull also being very close to 1. NOHARM and MIRT both performed better (i.e., had a mean closer to 1) for more indicators, as well as for larger samples in the case of MIRT.

6. Discussion

The goal of this study was to investigate a wide array of statistical methods for determining the number of factors to retain when using EFA with dichotomous indicators. With respect to the statistical methods under investigation in this work, the results presented above suggest that researchers should consider using the Combined approach based on EKC, Hull, and CD. This result was in keeping with Auerswald and Moshagen [2], who also found that a combination approach was effective in ascertaining the number of factors to retain when the indicators were normally distributed. In addition, the CNG method alone also proved to be very effective in revealing the number of factors to retain. Indeed, in the most challenging case with 5 indicators per factor and an interfactor correlation of 0.75, CNG yielded the most accurate results across methods. However, when 10 indicators were present, CNG was not the most accurate method in the 0.75 interfactor correlation condition. Thus, it cannot be uniformly recommended for use in all cases.
When 10 indicators per factor were found to be present, multiple methods yielded highly accurate results for determining the number of factors to retain. Specifically, both the EKC and Combination techniques were highly accurate regardless of the interfactor correlation. In addition, MAP, CD, and CNG yielded perfectly accurate results for the 0.25 and 0.50 interfactor correlation conditions. However, when the correlations were 0.75, they were somewhat less accurate than EKC and Combination.
When a single factor was present in the population, a number of methods studied here yielded highly accurate results, including EKC, CD, Combination, MAP, and CNG. Thus, when a single latent trait underlies the observed indicators, researchers can have confidence that any of these methods is likely to yield accurate results regarding the number of factors to retain, assuming conditions similar to those used in this study. Of course, in reality the data analyst will not know how many factors are actually present. Nonetheless, it does appear that if theory suggests a single factor is present and the methods mentioned above as being most accurate are used and yield a 1-factor solution, the researcher may be fairly confident in this result.
Several methods did not perform particularly well and may not be optimal for use in conditions similar to those used in this simulation study. In particular, PA and RPA, both of which performed well in previous studies [10], did not do so well under most conditions used in this study. Both methods tended to overfactor the data, leading to the low accuracy rates described above. It should be noted that in the previous work focusing on PA/RPA, there were more indicators per factor than were used in the current study. Thus, it seems that when the number of indicators per factor is 5, 10, or 20, neither PA nor RPA provides as accurate results as some other methods studied here, or as they do themselves with 30 or 60 indicators [10]. Similarly, the chi-square test based upon NOHARM also yielded lower accuracy rates than a number of the other techniques used in this study, particularly with 5 indicators. However, these results were due to a tendency to underfactor in the 3 latent trait condition, rather than overfactor, as was the case for PA/RPA. The MIRT model approach also yielded somewhat lower accuracy rates than a number of other methods in the 5 indicators 3-factors condition, but with a tendency to overfactor. When a single factor was present, PA, RPA, MIRT, and NOHARM all yielded overfactoring results. Taken together, these findings suggest that when statistical methods did not recommend retaining the correct number of factors, there were multiple ways in which they reached this incorrect decision. Sometimes overfactoring was the cause, whereas in other cases it was due to underfactoring. In any case, these methods may not be optimal for use with dichotomous indicators, particularly in light of the fact that other approaches (e.g., CNG, Combination and EKC, CD, or MAP in some cases) do appear to yield accurate results.

7. Implications for Practice

As described above, the results presented here provide some guidelines for researchers and data analysts. First, researchers and data analysts should consider the combined approach based on EKC, Hull, and CD for determining the number of factors to retain. The CNG technique would be a reasonable second choice for this purpose. Second, having more indicators per factor (presuming loadings in the range simulated here) will yield more accurate factor retention results. Thus, all other things being equal, researchers should seek to use longer scales. When scales are relatively short (e.g., 5 indicators per factor) and the interfactor correlation is high (e.g., 0.75), researchers may find it difficult to identify the correct number of factors to retain in an EFA, no matter the methodology used. Finally, larger samples were associated with more accurate results, in general. Therefore, researchers using EFA should attempt to obtain samples of 500 or more, particularly if they are working with a relatively small number (e.g., 5) of indicators per factor.

8. Limitations of This Study

As with all studies, there are limitations to the current work that future research should address. First, not all possible conditions that might be seen in practice were included in the current study. The goal of this work was to focus on dichotomous indicators such as might be encountered in many survey or testing environments. However, it is certainly the case that many scales rely on polytomous items that yield ordinal, rather than dichotomous, data. Thus, a limitation of this study was the fact that only items with two categories were included. Future work should include a focus on ordinal items. In addition, this study focused on cases with either 1 or 3 latent traits. However, many situations in practice involve scales with more possible latent traits, including assessments of cognitive ability, mood, or executive functioning. Thus, future work should include more latent traits while also examining a relatively small number of indicators, as was the case in this study. The current study included only latent traits that followed the standard normal distribution, which will not be present in all research contexts. For this reason, future work should consider cases where the latent trait is not normally distributed but rather might be skewed and/or highly kurtotic. Another limitation of the current study was with respect to the number of indicators per factor. In the current study, there were either 5 or 10 indicators per factor. In some instances, scales may be comprised of more indicators than this, e.g., 20, 30, or even more. Thus, future work should examine the performance of these methods with more indicators per factor.

9. Conclusions

The current study was designed to extend upon prior studies comparing methods for determining the number of factors to retain in the context of EFA. This earlier work focused on a wide array of methods in the context of normally distributed indicators [2] or PA/RPA with dichotomous indicators [10]. This study included a wide array of techniques in the context of dichotomous indicators with a small number of indicators per factor. The results suggest that data analysts and researchers might be best served to use CNG when there are 5 indicators per factor, and CNG, EKC, MAP, or the combination techniques for 10 indicators per factor. In such cases, one can anticipate accurate results regarding the number of factors to retain, particularly when the interfactor correlations are 0.5 or lower. Of course, regardless of what the statistical methods might suggest, researchers should always ensure that their ultimate findings are grounded in theory. If the optimal factor retention solution based on statistical results does not yield a result that is theoretically sound, the researcher should reconsider their analyses. However, when used in conjunction with theory, several methods studied here can help researchers identify the optimal number of factors to retain from an EFA.

Funding

This research received no external funding.

Data Availability Statement

Simulation code available upon request.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Gorsuch, R.L. Factor Analysis, 2nd ed.; Lawrence Erlbaum Associates, Publishers: Hillsdale, NJ, USA, 1983. [Google Scholar]
  2. Auerswald, M.; Moshagen, M. How to determine the number of factors to retain in exploratory factor analysis: A comparison of extraction methods under realistic conditions. Psychol. Methods 2019, 24, 468–491. [Google Scholar] [CrossRef] [PubMed]
  3. Boltzmann, M.; Schmidt, S.B.; Gutenbrunner, C.; Krauss, J.K.; Hoglinger, G.U.; Weimar, C.; Rollnik, J.D. Validity of the Early Functional Ability scale (EFA) among critically ill patients undergoing early neurological rehabilitation. BMC Neurol. 2022, 22, 333. [Google Scholar] [CrossRef]
  4. Selau, T.; da Silva, M.A.; de Mendonca Filho, E.J.; Bandeira, D.R. Evidence of validity and reliability of the adaptive functioning scale for intellectual disability (EFA-DI). Psicol. Reflex. E Crit. 2020, 33, 26. [Google Scholar] [CrossRef] [PubMed]
  5. Beck, A.T.; Steer, R.A.; Garbin, M.G. Psychometric properties of the Beck Depression Inventory: Twenty-five years of evaluation. Clin. Psychol. Rev. 1988, 8, 77–100. [Google Scholar] [CrossRef]
  6. Beck, A.T.; Epstein, N.; Brown, G.; Steer, R.A. An inventory for measuring clinical anxiety: Psychometric properties. J. Consult. Clin. Psychol. 1988, 56, 893–897. [Google Scholar] [CrossRef]
  7. Selbom, M. The MMPI-2-Restructured Form (MMPI-2-RF): Assessment of personality and psychopathology in the Twenty-First Century. Annu. Rev. Clin. Psychol. 2019, 15, 149–177. [Google Scholar] [CrossRef]
  8. Kuncel, N.R.; Wee, S.; Serafin, L.; Hezlett, S.A. The validity of the Graduate Record Examination for Master’s and Doctoral programs: A meta-analytic investigation. Educ. Psychol. Meas. 2010, 70, 340–352. [Google Scholar] [CrossRef]
  9. Wechsler, D. Wechsler Adult Intelligence Scale-Fourth Edition (WAIS-IV); American Psychological Association: Washington, DC, USA, 2008. [Google Scholar]
  10. Guo, W.; Choi, Y.-J. Assessing dimensionality of IRT models using traditional revised parallel analyses. Educ. Psychol. Meas. 2023, 83, 609–629. [Google Scholar] [CrossRef]
  11. Tabachnick, B.G.; Fidell, L.S. Using Multivariate Statistics; Pearson: New York, NY, USA, 2019. [Google Scholar]
  12. Finch, W.H. Exploratory Factor Analysis; Sage: Thousand Oaks, CA, USA, 2019. [Google Scholar]
  13. Bollen, K.A. Structural Equations with Latent Variables; John Wiley & Sons: New York, NY, USA, 1989. [Google Scholar]
  14. Tong, X.; Bentler, P.M. Evaluation of a New Mean Scaled and Moment Adjusted Test Statistic for SEM. Struct. Equ. Model. 2013, 20, 148–156. [Google Scholar] [CrossRef]
  15. Kim, J.-O.; Mueller, C.W. Factor Analysis: Statistical Methods and Practical Issues; Sage: Thousand Oaks, CA, USA, 1978. [Google Scholar]
  16. Brown, T.A. Confirmatory Factor Analysis for Applied Research; The Guilford Press: New York, NY, USA, 2015. [Google Scholar]
  17. Hayashi, K.; Bentler, P.M.; Yuan, K.-H. On the likelihood ratio test for the number of factors in exploratory factor analysis. Struct. Equ. Model. 2007, 14, 505–526. [Google Scholar] [CrossRef]
  18. Horn, J.L. A Rationale and Test for the Number of Factors in Factor Analysis. Psychometrika 1965, 30, 179–185. [Google Scholar] [CrossRef]
  19. Agresti, A. Categorical Data Analysis; Jon Wiley & Sons: New York, NY, USA, 2013. [Google Scholar]
  20. Green, S.B.; Levy, R.; Thompson, M.S.; Lu, M.; Lo, W.-J. A Proposed Solution to the Problem with using Completely Random Data to Assess the Number of Factors with Parallel Analysis. Educ. Psychol. Meas. 2012, 72, 357–374. [Google Scholar] [CrossRef]
  21. Ruscio, J.; Roche, B. Determining the Number of Factors to Retain in an Exploratory Factor Analysis using Comparison Data of Known Factorial Structure. Psychol. Assess. 2012, 24, 282–292. [Google Scholar] [CrossRef]
  22. Velicer, W.F. Determining the Number of Components from the Matrix of Partial Correlations. Psychometrika 1976, 41, 321–327. [Google Scholar] [CrossRef]
  23. Caron, P.-O. Minimum Average Partial Correlation and Parallel Analysis: The Influence of Oblique Structures. Commun. Stat. -Simul. Comput. 2019, 40, 2110–2117. [Google Scholar] [CrossRef]
  24. Garrido, L.E.; Abad, F.J.; Ponsoda, V. Performance of Velicer’s Minimum Average Partial Factor Retention Method with Categorical Variables. Educ. Psychol. Meas. 2011, 71, 551–570. [Google Scholar] [CrossRef]
  25. Zwick, W.R.; Velicer, W.F. Comparison of Five Rules for Determining the Number of Components to Retain. Psychol. Bull. 1986, 99, 432–442. [Google Scholar] [CrossRef]
  26. Cattell, R.B. The Scree Test for the Number of Factors. Multivar. Behav. Res. 1966, 2, 245–276. [Google Scholar] [CrossRef]
  27. Raiche, G.; Walls, T.A.; Magis, D.; Riopel, M.; Blais, J.-G. Non-Graphical Solutions for Cattell’s Scree Test. Methodology 2012, 9, 23–29. [Google Scholar] [CrossRef]
  28. Zoski, K.W.; Jurs, S. Using Multiple Regression to Determine the Number of Factors to Retain in Factor Analysis. Mult. Linear Regres. Viewp. 1993, 20, 5–9. [Google Scholar]
  29. Lorenzo-Seva, U.; Timmerman, M.E.; Kiers, H.A. The Hull method for selecting the number of common factors. Multivar. Behav. Res. 2011, 46, 340–364. [Google Scholar] [CrossRef]
  30. Braeken, J.; van Assen, M.A. An empirical Kaiser criterion. Psychol. Methods 2017, 22, 450–466. [Google Scholar] [CrossRef] [PubMed]
  31. Marcenko, V.A.; Pastur, L.A. Distribution of eigenvalues for some sets of random matrices. Math. USSR-Shornik 1967, 1, 457–483. [Google Scholar] [CrossRef]
  32. McDonald, R.P. Nonlinear Factor Analysis; Psychometric Monographs, No. 15; Psychometric Society: Richmond, VA, USA, 1967. [Google Scholar]
  33. Gessaroli, M.E.; De Champlain, A.F. Using an approximate Chi-square statistic to test the number of dimensions underlying the responses to a set of items. J. Educ. Meas. 1996, 33, 157–179. [Google Scholar] [CrossRef]
  34. Bock, R.D.; Aitkin, M. Marginal maximum likelihood estimation of item parameters: Application of an EM algorithm. Psychom. 1981, 46, 443–459. [Google Scholar] [CrossRef]
  35. Fabrigar, L.R.; Wegener, D.T. Exploratory Factor Analysis; Oxford University Press: Oxford, UK, 2011. [Google Scholar]
  36. Green, S.; Xu, Y.; Thompson, M.S. Relative accuracy of two modified parallel analysis methods that use the proper reference distribution. Educ. Psychol. Meas. 2018, 78, 589–604. [Google Scholar] [CrossRef]
  37. Green, S.B.; Thompson, M.S.; Levy, R.; Lo, W.-J. Type I and II Error Rates and Overall Accuracy of the Revised Parallel Analysis Method for Determining the Number of Factors. Educ. Psychol. Meas. 2015, 75, 428–457. [Google Scholar] [CrossRef]
  38. Lim, S.; Jahng, S. Determining the number of factors using parallel analysis and its recent variants. Psychol. Methods 2019, 24, 452–467. [Google Scholar] [CrossRef]
  39. Preacher, K.J.; MacCallum, R.C. Repairing Tom Swift’s Electric Factor Analysis Machine. Underst. Stat. 2003, 2, 13–43. [Google Scholar] [CrossRef]
  40. Green, S.B.; Redell, N.; Thompson, M.S.; Levy, R. Accuracy of Revised and Traditional Parallel Analyses for Assessing Dimensionality with Binary Data. Educ. Psychol. Meas. 2016, 76, 5–21. [Google Scholar] [CrossRef]
  41. Finch, W.H.; Habing, B. Performance of DIMTEST and NOHARM based statistics for testing unidimensionality. Appl. Psychol. Meas. 2007, 31, 292–307. [Google Scholar] [CrossRef]
  42. Finch, W.H.; Habing, B. Comparison of NOHARM and DETECT: Counting dimensions and allocating items. J. Educ. Meas. 2005, 42, 149–170. [Google Scholar] [CrossRef]
  43. Svetina, D.; Levy, R. Dimensionality in compensatory MIRT when complex structure exists: Evaluation of DETECT and NOHARM. J. Exp. Educ. 2016, 84, 398–420. [Google Scholar] [CrossRef]
  44. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2023. [Google Scholar]
  45. Orcan, F. MonteCarloSEM: An R package to simulate data for SEM. Int. J. Assess. Tools Educ. 2021, 8, 704–713. [Google Scholar] [CrossRef]
  46. Worthington, R.L.; Whittaker, T.A. Scale development research: A content analysis and recommendations for best practices. Couns. Psychol. 2006, 34, 806–838. [Google Scholar] [CrossRef]
  47. Cohen, J. Statistical Power Analysis for the Behavioral Sciences, 2nd ed.; Lawrence Erlbaum Associates, Publishers: Hillsdale, NJ, USA, 1988. [Google Scholar]
  48. Raiche, G.; Magis, D. nFactors: Parallel Analysis and Other Non-Graphical Solutions to the Cattell Scree Test. An R Software Library, R package version 2.4.1.1. 2022.
  49. O’Connor, B.P. EFA.dimensions: Exploratory Factor Analysis Functions for Assessing Dimensionality. A R Software Library, R package version 0.1.7.7. 2023.
  50. Chalmers, R.P. MIRT: A multidimensional item response theory package for the R environment. J. Stat. Softw. 2012, 48, 1–29. [Google Scholar] [CrossRef]
  51. Robitzsch, A. SIRT: Supplementary Item Response Theory Models. An R Software Library, R package version 3.13-228. 2022.
Table 1. Proportion of replications with correct number of factors, number of indicators, and interfactor correlation: 3 factors.
Table 1. Proportion of replications with correct number of factors, number of indicators, and interfactor correlation: 3 factors.
I *CEKCHullSMTCDCombMAPNMREGCNGPARPAMIRTNOHARM
50.250.99510.63620.00000.99110.99510.52300.41101.00000.69860.58050.54270.0408
0.500.84610.54950.00000.89110.84610.43140.50881.00000.67590.56480.59760.0334
0.750.00000.02230.00000.47390.00000.17310.07350.62720.66010.47380.51190.0443
100.251.00000.83640.00001.00001.00001.00000.78401.00000.79780.70000.87650.8234
0.501.00000.80920.00001.00001.00001.00000.81381.00000.77530.65030.81250.7294
0.750.94190.74090.00000.81160.94560.59660.00290.74320.74640.60430.69910.6402
* I = Number of indicators per factor; C = Interfactor correlation.
Table 2. Proportion of replications with correct number of factors, number of indicators, and sample size: 3 factors.
Table 2. Proportion of replications with correct number of factors, number of indicators, and sample size: 3 factors.
I *NEKCHullSMTCDCombMAPNMREGCNGPARPAMIRTNOHARM
52000.50730.33870.00000.57470.50730.51090.08280.91810.50000.40000.51260.0440
5000.64260.39650.00000.80870.64260.43870.33630.91170.55450.50040.55090.0473
10000.66480.44710.00000.88730.66480.36570.46340.90430.76920.71220.56000.0409
20000.66670.45410.00000.95600.66670.00000.55270.93261.00000.91240.57950.0257
102000.95670.66780.00000.87500.96170.86830.43430.85230.64170.58340.74570.6408
5000.96780.72550.00000.88600.96780.87670.44760.86560.76670.62470.78450.6780
10000.99800.89130.00000.98780.99800.88000.59230.87580.89180.80000.81560.7632
20001.00000.89740.00001.00001.00000.83720.66010.90930.99250.94130.84810.8155
* I = Number of indicators per factor; N = Sample size.
Table 3. Proportion of replications with correct number of factors, by number of number of factors and factor loading value.
Table 3. Proportion of replications with correct number of factors, by number of number of factors and factor loading value.
F *LEKCHullSMTCDCombMAPNMREGCNGPARPAMIRTNOHARM
10.601.00000.98740.00001.00001.00001.00000.13941.00000.87290.77750.59330.2946
0.801.00000.99430.00001.00001.00001.00000.09551.00000.90110.83830.65060.6423
30.600.78150.55330.00000.78780.78280.39390.39790.85610.72490.66700.64230.3473
0.800.84130.67940.00000.95610.84130.91480.48610.92670.83340.80000.69150.3792
* F = Number of factors; L = Factor loading value.
Table 4. Number of factors retained by number of indicators, and interfactor correlation: 3 factors.
Table 4. Number of factors retained by number of indicators, and interfactor correlation: 3 factors.
I *CEKCHullSMTCDCombMAPNMREGCNGPARPAMIRTNOHARM
50.252.99512.47917.66642.99422.99512.33862.45693.00004.08364.74933.44351.4236
0.502.84162.15147.72852.85832.84161.93374.13053.00003.8774.60053.39491.3478
0.751.17131.06017.56712.04381.17131.44116.26103.41673.46214.44753.48351.3869
100.253.00002.755520.34503.00003.00003.00002.79773.00003.77583.61033.11863.1470
0.503.00002.592120.30863.00003.00003.00005.62863.00003.46003.58813.15393.1386
0.752.92812.504120.06462.64562.94192.350214.72733.27543.20633.39142.96992.8878
* I = Number of indicators per factor; C = Interfactor correlation.
Table 5. Number of factors retained by number of indicators, and sample size: 3 factors.
Table 5. Number of factors retained by number of indicators, and sample size: 3 factors.
I *NEKCHullSMTCDCombMAPNMREGCNGPARPAMIRTNOHARM
52002.18721.83837.91882.24192.18722.11545.10663.09054.13244.52963.45741.3956
5002.35651.89127.69642.67992.35652.02244.19023.09563.55733.88803.44091.4396
10002.42901.94307.67822.83452.42901.90413.84383.10153.28253.42983.42421.3893
20002.40661.93316.99122.91212.40661.24753.69843.06743.00753.08853.44001.3198
102002.93832.412220.55982.75002.95672.79509.44933.17013.89664.02952.93232.8738
5002.96782.475520.65002.78362.96782.80999.31353.16513.31553.64453.10633.0665
10002.99802.786320.27422.99392.99802.79887.20153.12613.24503.32563.15193.1230
20003.00002.794819.47353.00003.00002.72994.90723.09443.00003.01103.18443.2781
* I = Number of indicators per factor; N = Sample size.
Table 6. Number of factors retained by number of number of factors and factor loading value.
Table 6. Number of factors retained by number of number of factors and factor loading value.
F *LEKCHullSMTCDCombMAPNMREGCNGPARPAMIRTNOHARM
10.601.00001.01267.14641.00001.00001.00005.86781.00001.81111.96481.43532.1216
0.801.00001.00578.69571.00001.00001.00007.17581.00001.31671.47111.47431.4117
30.602.63322.168213.39742.62752.63781.94797.06293.15643.63793.89413.30001.9109
0.802.72772.410015.47342.92282.72772.85935.03143.08243.34603.58033.24882.3896
* F = Number of factors; L = Factor loading value.
Table 7. Proportion of replications with correct number of factors, number of indicators, and sample size: 1 factor.
Table 7. Proportion of replications with correct number of factors, number of indicators, and sample size: 1 factor.
I *NEKCHullSMTCDCombMAPNMREGCNGPARPAMIRTNOHARM
52001.00000.98480.00001.00001.00001.00000.11041.00000.60320.57120.44690.1762
5001.00000.98720.00001.00001.00001.00000.17081.00000.67190.50040.47770.1985
10001.00000.99310.00001.00001.00001.00000.27041.00000.80030.71220.40200.3272
20001.00001.00000.00001.00001.00001.00000.28021.00001.00000.91240.48860.4463
102001.00000.97740.00001.00001.00001.00000.00001.00000.68900.62930.85070.6884
5001.00000.99380.00001.00001.00001.00000.00001.00000.81140.79000.81620.7029
10001.00000.99950.00001.00001.00001.00000.03851.00000.91710.86570.78700.6860
20001.00001.00000.00001.00001.00001.00000.04731.00001.0000.97030.67810.6960
* I = Number of indicators per factor; N = Sample size.
Table 8. Number of factors retained by number of indicators, and sample size: 1 factor.
Table 8. Number of factors retained by number of indicators, and sample size: 1 factor.
I *NEKCHullSMTCDCombMAPNMREGCNGPARPAMIRTNOHARM
52001.00001.01524.35201.00001.00001.00003.27681.00001.84382.36511.65192.2753
5001.00001.01284.32821.00001.00001.00002.89221.00001.45131.81921.65952.2958
10001.00001.00004.20281.00001.00001.00002.58451.00001.15161.39491.76232.0383
20001.00001.00694.24461.00001.00001.00002.58821.00001.00001.36291.66161.7961
102001.00001.006212.22841.00001.00001.000011.63581.00001.49381.62961.35751.3419
5001.00001.006212.06171.00001.00001.000011.32811.00001.27281.51421.21841.3425
10001.00001.022611.63531.00001.00001.00009.79621.00001.09961.29071.18381.3406
20001.00001.000011.09031.00001.00001.00008.72661.00001.00001.07121.16281.3475
* I = Number of indicators per factor; N = Sample size.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Finch, W.H. A Comparison of Methods for Determining the Number of Factors to Retain with Exploratory Factor Analysis of Dichotomous Data. Psych 2023, 5, 1004-1018. https://doi.org/10.3390/psych5030067

AMA Style

Finch WH. A Comparison of Methods for Determining the Number of Factors to Retain with Exploratory Factor Analysis of Dichotomous Data. Psych. 2023; 5(3):1004-1018. https://doi.org/10.3390/psych5030067

Chicago/Turabian Style

Finch, W. Holmes. 2023. "A Comparison of Methods for Determining the Number of Factors to Retain with Exploratory Factor Analysis of Dichotomous Data" Psych 5, no. 3: 1004-1018. https://doi.org/10.3390/psych5030067

APA Style

Finch, W. H. (2023). A Comparison of Methods for Determining the Number of Factors to Retain with Exploratory Factor Analysis of Dichotomous Data. Psych, 5(3), 1004-1018. https://doi.org/10.3390/psych5030067

Article Metrics

Back to TopTop