Abstract
The correlated trait–correlated method (CTCM) model for the analysis of multitrait–multimethod (MTMM) data is known to suffer convergence and admissibility (C&A) problems. We describe a little known and seldom applied reparameterized version of this model (CTCM-R) based on Rindskopf’s reparameterization of the simpler confirmatory factor analysis model. In a Monte Carlo study, we compare the CTCM, CTCM-R, and the correlated trait–correlated uniqueness (CTCU) models in terms of C&A, model fit, and parameter estimation bias. The CTCM-R model largely avoided C&A problems associated with the more traditional CTCM model, producing C&A solutions nearly as often as the CTCU model, but also avoiding parameter estimation biases known to plague the CTCU model. As such, the CTCM-R model is an attractive alternative for the analysis of MTMM data.
The choice of a particular confirmatory factor analytic (CFA) model for the analysis of multitrait–multimethod (MTMM) data remains controversial. Two of the most popularly adopted models, namely, the correlated trait–correlated method model (CTCM) and the correlated trait–correlated uniqueness model (CTCU), have been shown to suffer from empirical estimation (Conway, Lievens, Scullen, & Lance, 2004; Marsh & Bailey, 1991) and theoretical problems (Lance, Noble, & Scullen, 2002), respectively. More than 30 years ago, Rindskopf (1983) presented a reparameterization of the CFA model that effectively implements nonnegative restrictions on uniqueness estimates and thus avoids Heywood cases, but this model has received very little attention in the literature.
In the current study, we (a) show how Rindskopf’s reparameterized CFA model is extended to MTMM data and (b) compare the CTCM-R model to the other two frequently used CFA models (i.e., the CTCM and CTCU models) in a comprehensive Monte Carlo simulation study. The three models were empirically evaluated based on model convergence and admissibility (C&A), model goodness of fit, and accuracy of the parameter estimates. In the following sections, after a brief review of relevant literature on the strength and limitations of the CTCM and CTCU models, we describe the CTCM-R model in more detail, and then present descriptions of the simulation study design, our results, and a brief discussion of the importance of our findings.
Background
One of the main sources of controversy in the MTMM literature has concerned the optimal analytic approach for MTMM data, and consensus has evolved over time that the most appropriate analytic approach is some form of structural equation modeling (SEM; Eid, Lischetzke, & Nussbeck, 2006; Kenny & Kashy, 1992; Schmitt & Stults, 1986). Two of the most frequently adopted SEMs are the CTCM (or one or more of its special cases, see Widaman, 1985) and the CTCU models (Kenny, 1976; Marsh, 1989).
The CTCM model can be written as
where MTMM is the p×p MTMM correlation matrix where, typically (but not necessarily), p = T*M (where T refers to the number of Traits and M refers to the number of Methods).
Marsh (1989) and Marsh and Bailey (1991) elaborated on the CTCU model developed by Kenny (1976) as an alternative to the CTCM model that largely avoids empirical estimation problems associated with it. The CTCU model can be written as
where again
Several studies have compared the CTCM model and the CTCU model with simulated data (Conway et al., 2004; Marsh & Bailey, 1991; Zhang, Jin, Leite, & Algina, 2014) or previous published data (Lance et al., 2002). This literature indicates that while the CTCU model largely solves C&A problems associated with the CTCM model, it suffers from a number of conceptual limitations (Lance et al., 2002) and yields upwardly biased estimates of Trait factor loadings and correlations when Method factors are correlated (Conway et al., 2004), and Lance, Dawson, Birklebach, and Hoffman’s (2010) and Lance, Fan, Siminovsky, Morgan, and Shaikh’s (2014) reviews show that this is routinely the case. The sparse literature on the CTCM-R model (Dillon, Kumar, & Mulani, 1987; Lance & Fan, 2016; Marsh, 1989) suggests that it may retain the positive features of the CTCM model and yet yield C&A rates that are comparable to the CTCU model.
CTCM-R Model
In the traditional CFA model (
Here,
The CFA-R model has been implemented successfully in only a few CFA (La Du & Tanaka, 1989; Nagy, Trautwein, & Lüdtke, 2010) and SEM studies (Hashoul-Andary et al., 2016; Walker & Weber, 1987) where it effectively resolved inadmissibly issues. However, it appears that the CFA-R model has been studied systematically only once before using simulated data. Dillon et al. (1987) compared (a) the CFA-R model, (b) an alternative parameterization presented by Bentler (1976), and (c) simply fixing offending estimates (i.e., Heywood cases) to zero on an ad hoc basis in a simple CFA model containing five manifest and two latent variables. Based on their findings, Dillon et al. (1987) had little to recommend for the first two approaches, saying “though these two approaches represent theoretically elegant ways of handling this problem, they have limited practical usefulness compared with the simpler approach in which the offending parameter estimate is fixed at zero” (p. 134). We disagree with respect to the CFA-R model: fixing offending estimates to zero is an ad hoc symptom-based approach, while the CFA-R model represents a more general and preemptive model-based strategy.
The CFA-R model may be extended to the CTCM-R model, written as
where
The CTCM-R model has also been applied in only a few MTMM studies (e.g., Kinicki, McKee-Ryan, Schriesheim, & Carson, 2002; Lance et al., 2010; Marsh, 1989), and Lance and Fan (2016) and Lance et al. (2014) found that it performed effectively in the reanalysis of previously published MTMM data. However, we know of no systematic attempt to study the performance of the CTCM-R model relative to the CTCM and CTCU models in the analysis of simulated data with known population parameters.
Study Purpose
The goal of this study was to systematically compare the CTCM-R model relative to the CTCM and CTCU models using Monte Carlo simulations under a variety of study conditions that are representative of characteristics of MTMM matrices present in the literature. We chose the CTCM and CTCU models as comparison models because (a) of their popularity in the literature and (b) their relative advantages and disadvantages are already well known. Our general expectation was that the CTCM-R model would demonstrate superior C&A rates as compared to the CTCM model and simultaneously avoid parameter estimation bias that is known to plague the CTCU model. That is, we expected that the CTCM-R model would preserve the strengths and avoid the weaknesses of both the CTCM and CTCU models.
Method
This simulation study had four major steps: (a) generating model-implied population covariance matrices based on predefined population values; (b) analyzing sample data with the CTCM, CTCM-R, and CTCU models; (c) obtaining model convergence and admissibility rates, goodness-of-fit indices, and parameter estimates; and (d) examining and comparing the three models’ performance in terms of model C&A rates, model–data fit, and bias in parameter estimates (i.e., deviation of model estimates from true population values).
Population Characteristics
Population data were simulated to represent realistic MTMM data toward maximizing the generalizability of its findings to actual MTMM data found in the literature. As such, population values for the current simulation study came from the latest large-scale review of MTMM studies to date (Lance et al., 2014) and represented high, medium, and low values on the various characteristics of MTMM data that were manipulated (see below). Table 1 shows a summary of characteristics of MTMM matrices reviewed by Lance et al. (2014) and the population parameters and their levels chosen for the present study. Values chosen for this study represented a wide variety of subject areas (e.g., assessment centers, personality, attitudes, job performance, self-concept) and measurement methods (e.g., response formats, test forms, rater sources) and included values that others have used in MTMM simulation studies (Conway et al., 2004; Lance, Woehr, & Meade, 2007). As such, the simulation was a completely crossed factorial design with 4 (Sample Size: 100, 250, 500, 1,000) × 3 (Trait Factor Loading Size: .31, .50, .69) × 3 (Trait Correlation Size: .07, .36, .61) × 3 (Method Factor Loading Size: .17, .28, .48) × 3 (Method Correlation Size: .02, .29, .58) × 2 (Number of Traits Factors: 3, 5) × 2 (Number of Method Factors: 3, 5), = 1,296 unique conditions.
Population Values of Simulation Study.
Note. MTMM = multitrait–multimethod. k = 258 studies.
Procedure
We generated population correlation matrices using parameter values shown in Table 1. Specifically, population MTMM matrices were created using combinations of predefined population values shown in Table 1 with the CTCM model as the generating population model. Note that within a specified condition, all factor loadings and correlations were varied ±.10 so that their averages were equal to the specified population values. Table 2 shows one example of the
Example of Lambda and Phi Matrices for Generating a 5T3M Population Matrix.
Note. MTMM = multitrait–multimethod. The pattern matrix for 5T3M population MTMM matrix with averaged Method Loading at .28, Trait loading at .5, Method Correlation at .29, Trait correlation at .36.
We used LISREL 8 (Jöreskog & Sörbom, 1993) to analyze the generated sample data with each of the three targeted CFA-MTMM models (i.e., CTCM, CTCM-R, and CTCU). We harvested data from LISREL output files that contain model convergence and admissibility indicators, a number of model-fit indices (e.g., comparative fit index, Tucker–Lewis index, χ2, root mean square error of approximation), and parameter estimates that were all saved externally and analyzed for model performance comparisons.
Results
We conducted six 4 × 3 × 3 × 3 × 3 × 2 × 2 analyses of variance using the following dependent variables:
A binary variable coded 1 = the solution was convergent within 1,000 iterations and admissible (contained no standardized factor loadings or factor correlations >1.0 and no negative uniquenesses, where applicable) = 1, or 0 = the solution was nonconvergent and/or inadmissible (the solution did not converge within 1,000 iterations and/or contained at least one standardized factor loading or factor correlation >1.0 or at least one negative uniqueness)
A binary indicator of model goodness-of-fit coded 1 = the solution met Hu and Bentler’s (1999) cutoff criteria for Bentler’s (1990) comparative fit index and the Tucker–Lewis index ≥.95, the root mean square error of approximation ≤.06, and the standardized root mean square residual ≤.08, or 0 = the solution failed to satisfy one or more of these cutoff criteria
Bias in estimation of Trait and Method factor loadings and correlations measured as the deviation between the estimates and the true population values
Due to the extremely large sample size, every effect was statistically significant, even using an alpha rate as stringent as .0001. As such, we chose to report results associated with a minimally practically significant effect size indicated by η2≥ .02. 1 Using this criterion, no higher order effects beyond two-way interactions were indicated as being practically significant.
As Table 3 shows, and as expected, the CTCM model suffered C&A problems with an overall lower C&R rate (29.2%) than the other two models (η2 = .237), despite the fact that it was the generating model. The CTCM-R model, even though it is mathematically equivalent to the CTCM model, returned C&A solutions much more often (62.5%), approaching that of the CTCU model (87.5%). As such, even though the CTCM and CTCM-R are mathematically equivalent, the CTCM-R parameterization returned C&A solutions more than twice as often as the more traditional CTCM model.
Results: Model Convergence and Admissibility (C&A).
Note. C&A = model convergence and admissibility; CTCM = correlated trait–correlated method; CTCU = correlated trait–correlated uniqueness.
Table 3 also shows an effect of sample size on C&A (η2 = .077): not surprisingly, sample estimates became more stable with larger sample sizes. C&A rates also increased with larger population Trait (η2 = .043) and Method (η2 = .050) factor loadings. Since both Trait and Method factors represent sources of systematic variance, this suggests that, in general, C&A solutions are achieved more often with more reliable data.
Table 4 shows results for two small effects for analytic model and Trait loading level. Here, the best model fit is for the CTCM model and in one sense this comes as no surprise because it is the generating model. However, even though it is the generating model the CTCM model reached C&A solutions relatively infrequently. Thus, when a C&A solution is achieved for the CTCM model it is usually a very well fitting solution. The CTCM-R model is more forgiving in the sense that it achieves C&A solutions far more often than the CTCM model but sometimes at the cost of obtaining a reproduced matrix whose fit to the data does not satisfy all of the Hu and Bentler (1999) cutoff criteria for the goodness-of-fit indices we analyzed. Of course, the CTCU model was the most forgiving of all in terms of achieving C&A but had the lowest percentage of solutions yielding good model fit because, of course, it is the wrong model (i.e., it was not the generating model). Still, model fit was good for the CTCU model almost 92% of the time. This is consistent with Lance et al.’s (2007) findings that good model fit is not necessarily informative as to the correctness of the model fit to the data. Table 4 also shows that higher population Trait loadings also led to more frequent C&A solutions and better model fit, again likely due to variables’ higher reliabilities in these conditions.
Results: Good Model Fit.
Note. C&A = model convergence and admissibility; CTCM = correlated trait–correlated method; CTCU = correlated trait-correlated uniqueness.
Table 5 shows results for parameter estimation bias. Results for the CTCM and CTCM-R models indicated negligible bias for Trait factor loadings and slightly negative (conservative) bias for Trait correlations. Conversely, estimates for both Trait factor loadings and correlations were significantly positively biased for the CTCU model and this is consistent with previous research (Conway et al., 2004). We also found two interaction effects on Trait factor loading bias and these are shown in Figures 1 and 2. Figure 1 shows that Trait factor loading bias remains negligible across levels of Method factor loadings for the CTCM and CTCM-R models, but as Method loadings increase, positive bias in the estimation of Trait factor loadings increases dramatically for the CTCU model (η2 = .021). Similarly, Figure 2 shows that bias in Trait loading estimates for the CTCM and CTCM-R models remains negligible across variations in Method correlations but as Methods correlate more highly, substantial bias in Trait factor loading estimates also incurs for the CTCU model (η2 = .024). These findings are consistent with Conway et al.’s (2004) results: When Method factor loadings or correlations are low, minimal bias incurs for estimates of Trait factor loadings and correlations for the CTCU model, but when Method factor loadings and correlations are high, substantial bias incurs. This would not be a problem were it not for the fact that several large-scale reviews of MTMM-related research have estimated mean Method correlations at .34 (Doty & Glick, 1998), .37 (Buckley, Cote, & Comstock, 1990), .48 (Cote & Buckley, 1987), .52 (Lance et al., 2010), and .63 (Williams, Cote, & Buckley, 1989)—values that Figures 1 and 2 show introduce significant bias into CTCU estimates of Trait factor loadings and correlations.
Results: Parameter Estimation Bias.
Note. CTCM = correlated trait–correlated method; CTCU = correlated trait-correlated uniqueness.

Interaction effect of method loading level and analytic model on trait factor loading bias.

Interaction effect of method correlation level and analytic model on trait factor correlation bias.
Discussion
This study presented a comprehensive examination of the CTCM-R model in relation to two other widely applied CFA-MTMM models, the CTCM and CTCU models, in terms of C&A, model fit, and accuracy of parameter estimates. In general terms, the most important finding from this study was that the CTCM-R model at the same time preserved positive features of both the CTCM and CTCU models and avoided their shortcomings. The CTCM-R model is just as faithful to Campbell and Fiske’s (1959) original theoretical presentation of the MTMM methodology as is its mathematically identical CTCM model but largely avoids its C&A problems. On the other hand, the CTCM-R model maintained C&A rates that approached those of the CTCU model without suffering its conceptual limitations and inherent estimation biases.
The CTCM and CTCM-R models yielded slightly upwardly biased estimates for Method factor loadings (mean bias = .027 and .055, respectively) and Method correlations (mean bias = .051 and .055, respectively) and this was unexpected as the CTCM model was the generating model. But these biases pale in comparison to the large negative biases incurred by the CTCU model for Method loadings (mean bias = −.313) and correlations (mean bias = −.304) by virtue of the fact that these parameters are not estimated under the CTCU model. In effect, the omission of Method factors under the CTCU model creates an unmeasured variables problem (James, 1980) that, under routine application of the CTCU model, (a) renders Method variance unquantified in an amalgam of Method, unique and nonsystematic effort variance, and (b) leads to the biases in Trait factor loadings and correlations shown in Table 5. We do note, however, that Conway (1998) and Scullen (1999) have proposed two-step procedures by which Method variance can be estimated from uniquenesses’ covariances. However, to our knowledge, these procedures have been reported being used only four times previously (Li, Hughes, Kwok, & Hsu, 2012; Lievens & Conway, 2001; Siers & Christiansen, 2013; Van Iddekinge, Raymark, Eidson, & Attenweiler, 2004), and the accuracy of these procedures has not been investigated. This is an obvious direction for future work.
Cautions and Limitations
First, the CTCM-R model retains the positive features of the CTCM model and largely avoids its C&A problems, but it does not eliminate all possible sources of nonconvergent and inadmissible solutions. Model misspecifications can still lead to Trait and/or Method factor correlation estimates >|1.0|. Note that the CTCU model also does not solve this problem for Trait factors and skirts the issue entirely with respect to Method factors because they are not modeled.
Second, it could be argued that the CTCM-R model is unnecessary as most SEM software allows the user to intervene when faced with inadmissible estimates. For example, EQS automatically fences parameter estimates to boundary values (Bentler, 2006), and inequality constraints can be invoked in LISREL (Jöreskog & Sörbom, 1993) and Mplus (Muthén & Muthén, 1998-2012). However, these “fixes” may obscure important sources of model misspecification and are often undertaken in a post hoc, atheoretical fashion (e.g., Bowler & Woehr, 2006). Savalei and Kolenikov (2008) provide a discussion of the trade-offs between diagnosticity of unconstrained but improper solutions versus constrained and interpretable solutions. As a model-based (vs. symptom-based) solution, the CTCM-R model seeks a pre-emptive solution to inadmissible estimates and uninterpretable solutions.
Third, we note that there are many other mathematical models for MTMM data that we did not consider here. For example, Widaman (1985) presented a taxonomy of CFA models for MTMM data, the majority of which were special cases of the CTCM model. One of these, a correlated trait–uncorrelated method (CTUM) model is similar to the CTCU model in many respects 2 and has been shown to produce similarly biased model parameter estimates (Lance & Fan, 2016). Like the CTCU model, Eid’s (2000) CTC(M-1) model exhibits better C&A rates than the CTCM model but provides an incomplete representation of the MTMM factor space and for this reason we chose not to study it here. Also, we chose not to study Campbell and O’Connell’s (1967) direct product model as it is concerned with the estimation of interactive Trait × Method relationships whereas we investigated only linear effects. Finally, it should be noted that there is an emerging literature on generalizability theory (GT) approaches to the analysis of MTMM-like data (e.g., Jackson, Michaelides, Dewberry, & Kim, 2016; Putka & Hoffman, 2013; Woehr, Putka, & Bowler, 2012). GT offers certain advantages over CFA approaches to the analysis of MTMM data. For example, GT is typically capable of modeling (many) more data facets and can much more easily accommodate ill-structured data (Putka, Lance, Le, & McCloy, 2011) compared with CFA, but GT necessarily constrains measurement facets and elements within facets to be orthogonal and can be computationally intensive compared to CFA (Woehr et al., 2012).
Conclusion
The present study shows that the CTCM-R model simultaneously retains the attractive features of both the CTCM and CTCU models and avoids their limitations. Combined with evidence that the CTCM-R model performs well in the analysis of actual MTMM data (e.g., Lance et al., 2014) the present results with simulated data indicate that the CTCM-R model is a viable candidate analytic model for MTMM data.
Footnotes
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
