Testing a unidimensional,qualitative unfolding model for attitudinal or developmental data

Assuming that subject responses rank order stimuli by preference, statistical methods are presented for testing the hypothesis that responses conform to a unidimensional, qualitative unfolding model and to an a priori stimulus ordering. The model postulates that persons and stimulus variables are ordered along a single continuum and that subjects most prefer stimuli nearest their own position. The underlying continuum need not form an interval scale of the stimulus attribute. The general assumptions of the test for the unfolding model make it suitable for the analysis of structure in attitude responses, preference data, and developmental stage data.This research was supported by a grant from the U.S. Public Health Service (Grant No. 1-R01-MH27861-01) to the University of Minnesota. I wish to thank Sanford Weisberg for his helpful suggestions. I also wish to thank Karen Kitchener and Patricia King for letting me use their data.     

Acquiescent Responding in Balanced Multidimensional Scales and Exploratory Factor Analysis

Personality tests often consist of a set of dichotomous or Likert items. These response formats are known to be susceptible to an agreeing-response bias called acquiescence. The common assumption in balanced scales is that the sum of appropriately reversed responses should be reasonably free of acquiescence. However, inter-item correlation (or covariance) matrices can still be affected by the presence of variance due to acquiescence. To analyse these correlation matrices, we propose a method that is based on an unrestricted factor analysis and can be applied to multidimensional scales. This method obtains a factor solution in which acquiescence response variance is isolated in an independent factor. It is therefore possible, without the potentially confounding effect of acquiescence, to: (a) examine the dominant factors related to content latent variables; and (b) estimate participants’ factor scores on content latent variables. This method, which is illustrated by two empirical data examples, has proved to be useful for improving the simplicity of the factor structure. This research was partially supported by a grant from the Spanish Ministry of Science and Technology (SEJ2005-09170-C04-04/PSIC), and a grant from the Catalan Ministry of Universities, the Research and Information Society (2005SGR00017). The authors are obliged to the team of reviewers for helpful comments on an earlier version of this paper.     

Latent curve analysis

As a method for representing development, latent trait theory is presented in terms of a statistical model containing individual parameters and a structure on both the first and second moments of the random variables reflecting growth. Maximum likelihood parameter estimates and associated asymptotic tests follow directly. These procedures may be viewed as an alternative to standard repeated measures ANOVA and to first-order auto-regressive methods. As formulated, the model encompasses cohort sequential designs and allow for period or practice effects. A numerical illustration using data initially collected by Nesselroade and Baltes is presented.The authors wish to thank John Nesselroade for providing us the data for our illustration and Karen Paul and Connie Tilse for assisting in the data analysis. This research was supported by a grant (No. AG03164) from the National Institute on Aging to the senior author.     

Comparative evaluation of two superior stopping rules for hierarchical cluster analysis

A split-sample replication stopping rule for hierarchical cluster analysis is compared with the internal criterion previously found superior by Milligan and Cooper (1985) in their comparison of 30 different procedures. The number and extent of overlap of the latent population distributions was systematically varied in the present evaluation of stopping-rule validity. Equal and unequal population base rates were also considered. Both stopping rules correctly identified the actual number of populations when there was essentially no overlap and clusters occupied visually distinct regions of the measurement space. The replication criterion, which is evaluated by clustering of cluster means from preliminary analyses that are accomplished on random partitions of an original data set, was superior as the degree of overlap in population distributions increased. Neither method performed adequately when overlap obliterated visually discernible density nodes.This research was supported in part by NIMH grant 5R01 MH 32457 14.     

Components of reliability

As an alternative to the analysis of variance approach to reliability a multiple-factor analysis approach is illustrated. The one-factor and the multiple-factor models for reliability are compared. Tests on the latent roots associated with the principal components of intercorrelation matrices are used to determine the number of components to be retained. Conditions under which one or more of the principal components should be utilized are discussed.This investigation was completed at the Oregon Research Institute and supported by NIMH grant MH-0892-01 and NSF grant GB-257. Computations were performed at the Western Data Processing Center and the Health Sciences Computing Facility, University of California at Los Angeles.     

Linear structural equations with latent variables

An interdependent multivariate linear relations model based on manifest, measured variables as well as unmeasured and unmeasurable latent variables is developed. The latent variables include primary or residual common factors of any order as well as unique factors. The model has a simpler parametric structure than previous models, but it is designed to accommodate a wider range of applications via its structural equations, mean structure, covariance structure, and constraints on parameters. The parameters of the model may be estimated by gradient and quasi-Newton methods, or a Gauss-Newton algorithm that obtains least-squares, generalized least-squares, or maximum likelihood estimates. Large sample standard errors and goodness of fit tests are provided. The approach is illustrated by a test theory model and a longitudinal study of intelligence.This investigation was supported in part by a Research Scientist Development Award (KO2-DA00017) and a research grant (DA01070) from the U. S. Public Health Service.     

Neuroscientific Measures of Covert Behavior

In radical behaviorism, the difference between overt and covert responses does not depend on properties of the behavior but on the sensitivity of the measurement tools employed by the experimenter. Current neuroscientific research utilizes technologies that allow measurement of variables that are undetected by the tools typically used by behavior analysts. Data from a specific neuroscientific technique, event-related potential (ERP), suggest that emission of otherwise covert responses can be indexed and that such covert responses are sensitive to stimulus control and selection by consequences. The P3 ERP effect is proposed as indicative of emission. Moreover, ERP results in semantic priming experiments suggest that operants are sensitive to changes in stimulus control even when they are not emitted (latent responses). Changes in response strength of latent responses as a function of stimulus control can in fact be measured by reaction time data and an ERP dependent variable called the N400 effect. If the interpretations provided in this paper are accurate, an index of covertly emitted operants (P3 effect) constitutes experimental evidence suggesting the validity of a Skinnerian radical behaviorist perspective on behavior. Moreover, in a Skinnerian paradigm, measured fluctuations in the response strength of latent operants as a function of environmental changes (N400 effect) would validate Palmer's (2009) concept of the repertoire.     

Bayesian Semiparametric Structural Equation Models with Latent Variables

Structural equation models (SEMs) with latent variables are widely useful for sparse covariance structure modeling and for inferring relationships among latent variables. Bayesian SEMs are appealing in allowing for the incorporation of prior information and in providing exact posterior distributions of unknowns, including the latent variables. In this article, we propose a broad class of semiparametric Bayesian SEMs, which allow mixed categorical and continuous manifest variables while also allowing the latent variables to have unknown distributions. In order to include typical identifiability restrictions on the latent variable distributions, we rely on centered Dirichlet process (CDP) and CDP mixture (CDPM) models. The CDP will induce a latent class model with an unknown number of classes, while the CDPM will induce a latent trait model with unknown densities for the latent traits. A simple and efficient Markov chain Monte Carlo algorithm is developed for posterior computation, and the methods are illustrated using simulated examples, and several applications.     

A multidimensional item response model: Constrained latent class analysis using the gibbs sampler and posterior predictive checks

In this paper it will be shown that a certain class of constrained latent class models may be interpreted as a special case of nonparametric multidimensional item response models. The parameters of this latent class model will be estimated using an application of the Gibbs sampler. It will be illustrated that the Gibbs sampler is an excellent tool if inequality constraints have to be taken into consideration when making inferences. Model fit will be investigated using posterior predictive checks. Checks for manifest monotonicity, the agreement between the observed and expected conditional association structure, marginal local homogeneity, and the number of latent classes will be presented.This paper is supported by grant S40-645 of the Dutch Organization for Scientific Research (NWO).     

A Polychoric Instrumental Variable (PIV) Estimator for Structural Equation Models with Categorical Variables

This paper presents a new polychoric instrumental variable (PIV) estimator to use in structural equation models (SEMs) with categorical observed variables. The PIV estimator is a generalization of Bollen’s (Psychometrika 61:109–121, 1996) 2SLS/IV estimator for continuous variables to categorical endogenous variables. We derive the PIV estimator and its asymptotic standard errors for the regression coefficients in the latent variable and measurement models. We also provide an estimator of the variance and covariance parameters of the model, asymptotic standard errors for these, and test statistics of overall model fit. We examine this estimator via an empirical study and also via a small simulation study. Our results illustrate the greater robustness of the PIV estimator to structural misspecifications than the system-wide estimators that are commonly applied in SEMs. Kenneth Bollen gratefully acknowledges support from NSF SES 0617276, NIDA 1-RO1-DA13148-01, and DA013148-05A2. Albert Maydeu-Olivares was supported by the Department of Universities, Research and Information Society (DURSI) of the Catalan Government, and by grant BSO2003-08507 from the Spanish Ministry of Science and Technology. We thank Sharon Christ, John Hipp, and Shawn Bauldry for research assistance. The comments of the members of the Carolina Structural Equation Modeling (CSEM) group are greatly appreciated. An earlier version of this paper under a different title was presented by K. Bollen at the Psychometric Society Meetings, June, 2002, Chapel Hill, North Carolina.     

A general latent trait model for response processes

The purpose of the current paper is to propose a general multicomponent latent trait model (GLTM) for response processes. The proposed model combines the linear logistic latent trait (LLTM) with the multicomponent latent trait model (MLTM). As with both LLTM and MLTM, the general multicomponent latent trait model can be used to (1) test hypotheses about the theoretical variables that underlie response difficulty and (2) estimate parameters that describe test items by basic substantive properties. However, GLTM contains both component outcomes and complexity factors in a single model and may be applied to data that neither LLTM nor MLTM can handle. Joint maximum likelihood estimators are presented for the parameters of GLTM and an application to cognitive test items is described.This research was partially supported by the National Institute of Education grant number NIE-6-7-0156 to Susan Embretson (Whitely), principal investigator. However the optinions expressed herein do not necessarily reflect the position or policy of the National Institute of Education, and no official endorsement by the National Institute of Education should be inferred.     

Potentiation of latent inhibition

Rats were given exposure either to an odor (almond) or a compound of odor plus taste (almond plus saline), prior to training in which the odor served as the conditioned stimulus. It was found, for both appetitive and aversive procedures, that conditioning was retarded by preexposure (a latent inhibition effect), and the extent of the retardation was greater in rats preexposed to the compound (i.e., latent inhibition to the odor was potentiated by the presence of the taste). In contrast, the presence of the taste during conditioning itself overshadowed learning about the odor. We argue that the presence of the salient taste in compound with the odor enhances the rate of associative learning, producing a rapid loss in the associability of the odor. This loss of associability will generate both overshadowing and the potentiation of latent inhibition that is observed after preexposure to the compound.     

Asymptotic identifiability of nonparametric item response models

The identifiability of item response models with nonparametrically specified item characteristic curves is considered. Strict identifiability is achieved, with a fixed latent trait distribution, when only a single set of item characteristic curves can possibly generate the manifest distribution of the item responses. When item characteristic curves belong to a very general class, this property cannot be achieved. However, for assessments with many items, it is shown that all models for the manifest distribution have item characteristic curves that are very near one another and pointwise differences between them converge to zero at all values of the latent trait as the number of items increases. An upper bound for the rate at which this convergence takes place is given. The main result provides theoretical support to the practice of nonparametric item response modeling, by showing that models for long assessments have the property of asymptotic identifiability. The research was partially supported by the National Institute of Health grant R01 CA81068-01.     

Local influence analysis of nonlinear structural equation models

By regarding the latent random vectors as hypothetical missing data and based on the conditional expectation of the complete-data log-likelihood function in the EM algorithm, we investigate assessment of local influence of various perturbation schemes in a nonlinear structural equation model. The basic building blocks of local influence analysis are computed via observations of the latent variables generated by the Metropolis-Hastings algorithm, while the diagnostic measures are obtained via the conformal normal curvature. Seven perturbation schemes, including some perturbation schemes on latent vectors, are investigated. The proposed procedure is illustrated by a simulation study and a real example. Acknowledgment: This research is fully supported by a grant (CUHK 4243/02H) from the Research Grant Council of the Hong Kong Special Administration Region. The authors are indebted to ICPSR and the relevant funding agency for allowing use of their data, and to the Editor and reviewers for their valuable comments for improving the paper.     

Multilevel and Latent Variable Modeling with Composite Links and Exploded Likelihoods

Composite links and exploded likelihoods are powerful yet simple tools for specifying a wide range of latent variable models. Applications considered include survival or duration models, models for rankings, small area estimation with census information, models for ordinal responses, item response models with guessing, randomized response models, unfolding models, latent class models with random effects, multilevel latent class models, models with log-normal latent variables, and zero-inflated Poisson models with random effects. Some of the ideas are illustrated by estimating an unfolding model for attitudes to female work participation. We wish to thank The Research Council of Norway for a grant supporting our collaboration.     

A note on some modifications of latent roots and vectors

When an arbitrary positive scalar matrix is added to a correlation matrix the latent roots of the sum are equal to the corresponding roots of the correlation matrix plus an amount equal to the scalar number of the scalar matrix. The latent vectors of the sum are identical with those of the correlation matrix. An approximation to these relationships is suggested for the case in which the sum is of a correlation matrix and of a positive semidefinite diagonal matrix. The approximation is used to allow the solution of a characteristic problem for a correlation matrix with unities in the main diagonal to provide a family of solutions for the same correlation matrix.This research has been supported by a grant from the National Institute of Mental Health, MH 7864-01.     

On the misuse of manifest variables in the detection of measurement bias

Measurement invariance (lack of bias) of a manifest variableY with respect to a latent variableW is defined as invariance of the conditional distribution ofY givenW over selected subpopulations. Invariance is commonly assessed by studying subpopulation differences in the conditional distribution ofY given a manifest variableZ, chosen to substitute forW. A unified treatment of conditions that may allow the detection of measurement bias using statistical procedures involving only observed or manifest variables is presented. Theorems are provided that give conditions for measurement invariance, and for invariance of the conditional distribution ofY givenZ. Additional theorems and examples explore the Bayes sufficiency ofZ, stochastic ordering inW, local independence ofY andZ, exponential families, and the reliability ofZ. It is shown that when Bayes sufficiency ofZ fails, the two forms of invariance will often not be equivalent in practice. Bayes sufficiency holds under Rasch model assumptions, and in long tests under certain conditions. It is concluded that bias detection procedures that rely strictly on observed variables are not in general diagnostic of measurement bias, or the lack of bias.Preparation of this article was supported in part by PSC-CUNY grant #661282 to Roger E. Millsap.     

Building an identifiable latent class model with covariate effects on underlying and measured variables

In recent years, latent class models have proven useful for analyzing relationships between measured multiple indicators and covariates of interest. Such models summarize shared features of the multiple indicators as an underlying categorical variable, and the indicators' substantive associations with predictors are built directly and indirectly in unique model parameters. In this paper, we provide a detailed study on the theory and application of building models that allow mediated relationships between primary predictors and latent class membership, but that also allow direct effects of secondary covariates on the indicators themselves. Theory for model identification is developed. We detail an Expectation-Maximization algorithm for parameter estimation, standard error calculation, and convergent properties. Comparison of the proposed model with models underlying existing latent class modeling software is provided. A detailed analysis of how visual impairments affect older persons' functioning requiring distance vision is used for illustration.This work was supported by National Institute on Aging (NIA) Program Project P01-AG-10184-03 and National Institutes of Mental Health grant R01-MH-56639-01A1. Dr. Bandeen-Roche is a Brookdale National Fellow. The authors wish to thank Drs. Gary Rubin and Sheila West for kindly making the Salisbury Eye Evaluation data available. We also thank the Editor, the Associate Editor, and three referees for their valuable comments.     

Comparing Latent Structures of the Grade of Membership,Rasch, and Latent Class Models

This paper focuses on model interpretation issues and employs a geometric approach to compare the potential value of using the Grade of Membership (GoM) model in representing population heterogeneity. We consider population heterogeneity manifolds generated by letting subject specific parameters vary over their natural range, while keeping other population parameters fixed, in the marginal space (based on marginal probabilities) and in the full parameter space (based on cell probabilities). The case of a 2 × 2 contingency table is discussed in detail, and a generalization to 2^J tables with J ≥ 3 is sketched. Our approach highlights the main distinction between the GoM model and the probabilistic mixture of classes by demonstrating geometrically the difference between the concepts of partial and probabilistic memberships. By using the geometric approach we show that, in special cases, the GoM model can be thought of as being similar to an item response theory (IRT) model in representing population heterogeneity. Finally, we show that the GoM item parameters can provide quantities analogous to more general logistic IRT item parameters. As a latent structure model, the GoM model might be considered a useful alternative for a data analysis when both classes of extreme responses, and additional heterogeneity that cannot be captured by those latent classes, are expected in the population. This work was supported by Award #1R03 AG18986-01 from the National Institute on Aging and NIH grant #1R01 CA94212-01. The presentation of the ideas in this paper owes much to discussions with Stephen Fienberg and Brian Junker, Carnegie Mellon University. The author thanks Jim Ramsay and two anonymous reviewers for their valuable comments on earlier drafts of this paper.