Extending the Debate Between Spearman and Wilson 1929: When do Single Variables Optimally Reproduce the Common Part of the Observed Covariances?

The covariances of observed variables reproduced from conventional factor score predictors are generally not the same as the covariances reproduced from the common factors. We sought to find a factor score predictor that optimally reproduces the common part of the observed covariances. It was found algebraically that—under some conditions—the single observed variable with highest loading on a factor reproduces the non-diagonal elements of the observed covariance matrix more exactly than the conventional factor score predictors. This finding is linked to Spearman's and Wilson's 1929 debate on the use of single variables as factor score predictors. A population-based and a sample-based simulation study confirmed the algebraic result that taking a single variable can outperform conventional factor score predictors in reproducing the non-diagonal covariances when the nonzero loading size and the number of nonzero loadings per factor are small. The results indicated that a weighted aggregation of variables does not necessarily lead to an improvement of the score over the variable with the highest loading.     

Model‐related factor score predictors for confirmatory factor analysis

The present paper introduces model‐related (MR) factor score predictors, which reflect specific aspects of confirmatory factor models. The development is mainly based on Schönemann and Steiger's regression score components, but it can also be applied to the factor score coefficients. It is shown that the rotation of factor score predictors has no impact on the covariance matrix reproduced from the corresponding regression component patterns. Thus, regression score components or factor score coefficients can be rotated in order to obtain the required properties. This idea is the basis for MR factor score predictors, which are computed by means of a partial Procrustes rotation towards a target pattern representing the interesting properties of a confirmatory factor model. Two examples demonstrate the construction of MR factor score predictors reflecting specific constraints of a factor model.     

Normal Versus Noncentral Chi-square Asymptotics of Misspecified Models

In this article, we present a Bayesian spatial factor analysis model. We extend previous work on confirmatory factor analysis by including geographically distributed latent variables and accounting for heterogeneity and spatial autocorrelation. The simulation study shows excellent recovery of the model parameters and demonstrates the consequences of ignoring spatial dependence. Specifically, we find inefficiency in the estimates of the factor score means and bias and inefficiency in the estimates of the corresponding covariance matrix. We apply the model to Schwartz value priority data obtained from 5 European countries. We show that the Schwartz motivational types of values, such as Conformity, Tradition, Benevolence, and Hedonism, possess high spatial autocorrelation. We identify several spatial patterns—specifically, Conformity and Hedonism have a country-specific structure, Tradition has a North–South gradient that cuts across national borders, and Benevolence has South–North cross-national gradient. Finally, we show that conventional factor analysis may lead to a loss of valuable insights compared with the proposed approach.     

Regression among factor scores

Structural equation models with latent variables are sometimes estimated using an intuitive three-step approach, here denoted factor score regression. Consider a structural equation model composed of an explanatory latent variable and a response latent variable related by a structural parameter of scientific interest. In this simple example estimation of the structural parameter proceeds as follows: First, common factor models areseparately estimated for each latent variable. Second, factor scores areseparately assigned to each latent variable, based on the estimates. Third, ordinary linear regression analysis is performed among the factor scores producing an estimate for the structural parameter. We investigate the asymptotic and finite sample performance of different factor score regression methods for structural equation models with latent variables. It is demonstrated that the conventional approach to factor score regression performs very badly. Revised factor score regression, using Regression factor scores for the explanatory latent variables and Bartlett scores for the response latent variables, produces consistent estimators for all parameters.     

An estimate of the covariance between variables which are not jointly observed

Situations sometimes arise in which variables collected in a study are not jointly observed. This typically occurs because of study design. An example is an equating study where distinct groups of subjects are administered different sections of a test. In the normal maximum likelihood function to estimate the covariance matrix among all variables, elements corresponding to those that are not jointly observed are unidentified. If a factor analysis model holds for the variables, however, then all sections of the matrix can be accurately estimated, using the fact that the covariances are a function of the factor loadings. Standard errors of the estimated covariances can be obtained by the delta method. In addition to estimating the covariance matrix in this design, the method can be applied to other problems such as regression factor analysis. Two examples are presented to illustrate the method. This research was partially supported by NIMH grant MH5-4576     

Properties of the maximum likelihood solution in factor analysis regression

Algebraic properties of the normal theory maximum likelihood solution in factor analysis regression are investigated. Two commonly employed measures of the within sample predictive accuracy of the factor analysis regression function are considered: the variance of the regression residuals and the squared correlation coefficient between the criterion variable and the regression function. It is shown that this within sample residual variance and within sample squared correlation may be obtained directly from the factor loading and unique variance estimates, without use of the original observations or the sample covariance matrix.     

Some rao-guttman relationships

An examination of the determinantal equation associated with Rao's canonical factors suggests that Guttman's best lower bound for the number of common factors corresponds to the number of positive canonical correlations when squared multiple correlations are used as the initial estimates of communality. When these initial communality estimates are used, solving Rao's determinantal equation (at the first stage) permits expressing several matrices as functions of factors that differ only in the scale of their columns; these matrices include the correlation matrix with units in the diagonal, the correlation matrix with squared multiple correlations as communality estimates, Guttman's image covariance matrix, and Guttman's anti-image covariance matrix. Further, the factor scores associated with these factors can be shown to be either identical or simply related by a scale change. Implications for practice are discussed, and a computing scheme which would lead to an exhaustive analysis of the data with several optional outputs is outlined.     

Structural equation modeling with heavy tailed distributions

Data in social and behavioral sciences typically possess heavy tails. Structural equation modeling is commonly used in analyzing interrelations among variables of such data. Classical methods for structural equation modeling fit a proposed model to the sample covariance matrix, which can lead to very inefficient parameter estimates. By fitting a structural model to a robust covariance matrix for data with heavy tails, one generally gets more efficient parameter estimates. Because many robust procedures are available, we propose using the empirical efficiency of a set of invariant parameter estimates in identifying an optimal robust procedure. Within the class of elliptical distributions, analytical results show that the robust procedure leading to the most efficient parameter estimates also yields a most powerful test statistic. Examples illustrate the merit of the proposed procedure. The relevance of this procedure to data analysis in a broader context is noted. The authors thank the editor, an associate editor and four referees for their constructive comments, which led to an improved version of the paper.     

On equivariance and invariance of standard errors in three exploratory factor models

Current practice in factor analysis typically involves analysis of correlation rather than covariance matrices. We study whether the standardz-statistic that evaluates whether a factor loading is statistically necessary is correctly applied in such situations and more generally when the variables being analyzed are arbitrarily rescaled. Effects of rescaling on estimated standard errors of factor loading estimates, and the consequent effect onz-statistics, are studied in three variants of the classical exploratory factor model under canonical, raw varimax, and normal varimax solutions. For models with analytical solutions we find that some of the standard errors as well as their estimates are scale equivariant, while others are invariant. For a model in which an analytical solution does not exist, we use an example to illustrate that neither the factor loading estimates nor the standard error estimates possess scale equivariance or invariance, implying that different conclusions could be obtained with different scalings. Together with the prior findings on parameter estimates, these results provide new guidance for a key statistical aspect of factor analysis.We gratefully acknowledge the help of the Associate Editor and three referees whose constructive comments lead to an improved version of the paper. This work was supported by National Institute on Drug Abuse Grants DA01070 and DA00017 and by the University of North Texas Faculty Research Grant Program.     

Descriptive axioms for common factor theory,image theory and component theory

Through an extension of work by Guttman, common factor theory, image theory, and component theory are derived from distinct minimum subsets of assumptions chosen out of a set of five possible assumptions. It is thence shown that the problem of indeterminacy of factor scores in the common factor model is precisely reflected in the problem of the non-orthogonality of anti-images. Indeed, image scores are determinate for the same reason that the usual estimates of factor scores are determinate, and image scores cannot be used as though they were factor scores for the same reason that factor score estimates cannot be used as though they were factor scores.     

时距估计中的锚定效应

采用两种时距估计方法,检验时距锚定值对时距估计的影响,并探讨了时距信息的心理表征方式。63名在校大学生参加了本次实验。实验1采用口头报告法,表明较大的时距锚定值（5s,5000ms）条件下,被试估对时距的估计值较大,而较小时距锚定值（1s,1000ms）条件下被试估计的时距值较小;语义相同但表述方式不同的锚定值（1s与1000ms,5s与5000ms）条件下的时距估计值没有显著差异。实验2采用产生法,进一步表明时距表述方式对产生时距没有显著影响。以上结果表明,时距估计受时距锚定值的影响,时距信息可能以语义形式进行表征,而不是简单的数字加单位的表层表征形式     

Using the expectation maximization algorithm to estimate coefficient alpha for scales with item-level missing data

A 2-step approach for obtaining internal consistency reliability estimates with item-level missing data is outlined. In the 1st step, a covariance matrix and mean vector are obtained using the expectation maximization (EM) algorithm. In the 2nd step, reliability analyses are carried out in the usual fashion using the EM covariance matrix as input. A Monte Carlo simulation examined the impact of 6 variables (scale length, response categories, item correlations, sample size, missing data, and missing data technique) on 3 different outcomes: estimation bias, mean errors, and confidence interval coverage. The 2-step approach using EM consistently yielded the most accurate reliability estimates and produced coverage rates close to the advertised 95% rate. An easy method of implementing the procedure is outlined.     

Pseudo maximum likelihood estimation and a test for misspecification in mean and covariance structure models

Using the theory of pseudo maximum likelihood estimation the asymptotic covariance matrix of maximum likelihood estimates for mean and covariance structure models is given for the case where the variables are not multivariate normal. This asymptotic covariance matrix is consistently estimated without the computation of the empirical fourth order moment matrix. Using quasi-maximum likelihood theory a Hausman misspecification test is developed. This test is sensitive to misspecification caused by errors that are correlated with the independent variables. This misspecification cannot be detected by the test statistics currently used in covariance structure analysis.For helpful comments on a previous draft of the paper we are indebted to Kenneth A. Bollen, Ulrich L. Küsters, Michael E. Sobel and the anonymous reviewers of Psychometrika. For partial research support, the first author wishes to thank the Department of Sociology at the University of Arizona, where he was a visiting professor during the fall semester 1987.     

A Class of Population Covariance Matrices in the Bootstrap Approach to Covariance Structure Analysis

Model evaluation in covariance structure analysis is critical before the results can be trusted. Due to finite sample sizes and unknown distributions of real data, existing conclusions regarding a particular statistic may not be applicable in practice. The bootstrap procedure automatically takes care of the unknown distribution and, for a given sample size, also provides more accurate results than those based on standard asymptotics. But the procedure needs a matrix to play the role of the population covariance matrix. The closer the matrix is to the true population covariance matrix, the more valid the bootstrap inference is. The current paper proposes a class of covariance matrices by combining theory and data. Thus, a proper matrix from this class is closer to the true population covariance matrix than those constructed by any existing methods. Each of the covariance matrices is easy to generate and also satisfies several desired properties. An example with nine cognitive variables and a confirmatory factor model illustrates the details for creating population covariance matrices with different misspecifications. When evaluating the substantive model, bootstrap or simulation procedures based on these matrices will lead to more accurate conclusion than that based on artificial covariance matrices.     

Structural Equation Modeling of Multivariate Time Series

The covariance structure of a vector autoregressive process with moving average residuals (VARMA) is derived. It differs from other available expressions for the covariance function of a stationary VARMA process and is compatible with current structural equation methodology. Structural equation modeling programs, such as LISREL, may therefore be employed to fit the model.

Particular attention is given to assumptions concerning the process before the first observation. An application to a repeated time series is used to demonstrate the effect of these assumptions on the structure of the reproduced covariance matrix.     

Best univocal estimates of orthogonal common factors

A general formulation is presented for obtaining conditionally unbiased, univocal common-factor score estimates that have maximum validity for the true orthogonal factor scores. We note that although this expression is formally different from both Bartlett's formulation and Heermann's approximate expression, all three, while developed from very different rationales, yield identical results given that the common-factor model holds for the data. Although the true factor score validities can be raised by a different non-orthogonal transformation of orthogonalized regression estimates—as described by Mulaik—the resulting estimates lose their univocality.     

Information Matrices and Standard Errors for MLEs of Item Parameters in IRT

The paper clarifies the relationship among several information matrices for the maximum likelihood estimates (MLEs) of item parameters. It shows that the process of calculating the observed information matrix also generates a related matrix that is the middle piece of a sandwich-type covariance matrix. Monte Carlo results indicate that standard errors (SEs) based on the observed information matrix are robust to many, but not all, conditions of model/distribution misspecifications. SEs based on the sandwich-type covariance matrix perform most consistently across conditions. Results also suggest that SEs based on other matrices are either not consistent or perform not as robust as those based on the sandwich-type covariance matrix or the observed information matrix.     

More efficient parameter estimates for factor analysis of ordinal variables by ridge generalized least squares

Data in psychology are often collected using Likert‐type scales, and it has been shown that factor analysis of Likert‐type data is better performed on the polychoric correlation matrix than on the product‐moment covariance matrix, especially when the distributions of the observed variables are skewed. In theory, factor analysis of the polychoric correlation matrix is best conducted using generalized least squares with an asymptotically correct weight matrix (AGLS). However, simulation studies showed that both least squares (LS) and diagonally weighted least squares (DWLS) perform better than AGLS, and thus LS or DWLS is routinely used in practice. In either LS or DWLS, the associations among the polychoric correlation coefficients are completely ignored. To mend such a gap between statistical theory and empirical work, this paper proposes new methods, called ridge GLS, for factor analysis of ordinal data. Monte Carlo results show that, for a wide range of sample sizes, ridge GLS methods yield uniformly more accurate parameter estimates than existing methods (LS, DWLS, AGLS). A real‐data example indicates that estimates by ridge GLS are 9–20% more efficient than those by existing methods. Rescaled and adjusted test statistics as well as sandwich‐type standard errors following the ridge GLS methods also perform reasonably well.     

Multilevel models for multiple-baseline data: modeling across-participant variation in autocorrelation and residual variance

Multilevel models (MLM) have been used as a method for analyzing multiple-baseline single-case data. However, some concerns can be raised because the models that have been used assume that the Level-1 error covariance matrix is the same for all participants. The purpose of this study was to extend the application of MLM of single-case data in order to accommodate across-participant variation in the Level-1 residual variance and autocorrelation. This more general model was then used in the analysis of single-case data sets to illustrate the method, to estimate the degree to which the autocorrelation and residual variances differed across participants, and to examine whether inferences about treatment effects were sensitive to whether or not the Level-1 error covariance matrix was allowed to vary across participants. The results from the analyses of five published studies showed that when the Level-1 error covariance matrix was allowed to vary across participants, some relatively large differences in autocorrelation estimates and error variance estimates emerged. The changes in modeling the variance structure did not change the conclusions about which fixed effects were statistically significant in most of the studies, but there was one exception. The fit indices did not consistently support selecting either the more complex covariance structure, which allowed the covariance parameters to vary across participants, or the simpler covariance structure. Given the uncertainty in model specification that may arise when modeling single-case data, researchers should consider conducting sensitivity analyses to examine the degree to which their conclusions are sensitive to modeling choices.     

"TO SEE OURSELVES AS OTHERS SEE US!": RELATIONS AMONG SELF-PERCEPTIONS,PEER-PERCEPTIONS,AND EXPECTED PEER-PERCEPTIONS OF PERSONALITY

The need for convergent validity in the measurement of personality variables and past failures to demonstrate such properties are briefly reviewed. Two apppoaches to the problem using several types of data obtained from various subject populations and employing alternative factor analytic methods and factor score estimation procedulres are described. The results indicate that convergent factor patterns and substantial convergence and discriminant validities among factor score estimates are obtainable when sufficient care is taken in the design and development of measure variables .of each type and when analysis methods appropriate for the problem are used.