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.

A battery of thirty-six tests was given to a group of high-school seniors. The factorial analysis reveals essentially the same primary factors that were found in previous studies. The test battery reveals a simple structure.

This article devises a Bayesian multivariate formulation for analysis of ordinal data that records teacher classroom performance along multiple dimensions to assess aspects characterizing good instruction. Study designs for scoring teachers seek to measure instructional performance over multiple classroom measurement event sessions at varied occasions using disjoint intervals within each session and employment of multiple ratings on intervals scored by different raters; a design which instantiates a nesting structure with each level contributing a source of variation in recorded scores. We generally possess little a priori knowledge of the existence or form of a sparse generating structure for the multivariate dimensions at any level in the nesting that would permit collapsing over dimensions as is done under univariate modeling. Our approach composes a Bayesian data augmentation scheme that introduces a latent continuous multivariate response linked to the observed ordinal scores with the latent response mean constructed as an additive multivariate decomposition of nested level means that permits the extraction of de-noised continuous teacher-level scores and the associated correlation matrix. A semi-parametric extension facilitates inference for teacher-level dependence among the dimensions of classroom performance under multi-modality induced by sub-groupings of rater perspectives. We next replace an inverse Wishart prior specified for the teacher covariance matrix over dimensions of instruction with a factor analytic structure to allow the simultaneous assessment of an underlying sparse generating structure. Our formulation for Bayesian factor analysis employs parameter expansion with an accompanying post-processing sign re-labeling step of factor loadings that together reduce posterior correlations among sampled parameters to improve parameter mixing in our Markov chain Monte Carlo (MCMC) scheme. We evaluate the performance of our formulation on simulated data and make an application for the assessment of the teacher covariance structure with a dataset derived from a study of middle and high school algebra teachers.

Considerations of factor score estimates have concentrated on internal characteristics. This report considers external characteristics of four methods for determining factor score estimates; that is, relations of these estimates to measures on attributes not entered into the factor analysis. These external characteristics are important for many uses of factor score estimates. Findings are that different ones of the methods are appropriate for different uses.

The Maximum-likelihood estimator dominates the estimation of general structural equation models. Noniterative, equation-by-equation estimators for factor analysis have received some attention, but little has been done on such estimators for latent variable equations. I propose an alternative 2SLS estimator of the parameters in LISREL type models and contrast it with the existing ones. The new 2SLS estimator allows observed and latent variables to originate from nonnormal distributions, is consistent, has a known asymptotic covariance matrix, and is estimable with standard statistical software. Diagnostics for evaluating instrumental variables are described. An empirical example illustrates the estimator.

Diagnostic classification models are confirmatory in the sense that the relationship between the latent attributes and responses to items is specified or parameterized. Such models are readily interpretable with each component of the model usually having a practical meaning. However, parameterized diagnostic classification models are sometimes too simple to capture all the data patterns, resulting in significant model lack of fit. In this paper, we attempt to obtain a compromise between interpretability and goodness of fit by regularizing a latent class model. Our approach starts with minimal assumptions on the data structure, followed by suitable regularization to reduce complexity, so that readily interpretable, yet flexible model is obtained. An expectation–maximization-type algorithm is developed for efficient computation. It is shown that the proposed approach enjoys good theoretical properties. Results from simulation studies and a real application are presented.

A maximum likelihood procedure for combining standardized mean differences based on a noncentratt-distribution is proposed. With a proper data augmentation technique, an EM-algorithm is developed. Information and likelihood ratio statistics are discussed in detail for reliable inference. Simulation results favor the proposed procedure over both the existing normal theory maximum likelihood procedure and the commonly used generalized least squares procedure.

The recent “Every Student Succeed Act" encourages schools to use an innovative assessment to provide feedback about students’ mastery level of grade-level content standards. Mastery of a skill requires the ability to complete the task with not only accuracy but also fluency. This paper offers a new sight on using both response times and response accuracy to measure fluency with cognitive diagnosis model framework. Defining fluency as the highest level of a categorical latent attribute, a polytomous response accuracy model and two forms of response time models are proposed to infer fluency jointly. A Bayesian estimation approach is developed to calibrate the newly proposed models. These models were applied to analyze data collected from a spatial rotation test. Results demonstrate that compared with the traditional CDM that using response accuracy only, the proposed joint models were able to reveal more information regarding test takers’ spatial skills. A set of simulation studies were conducted to evaluate the accuracy of model estimation algorithm and illustrate the various degrees of model complexities.

The concept of sequential estimation is introduced in multidimensional scaling (MDS). The sequential estimation method developed in this paper refers to continually updating estimates of a configuration as new observations are added. This method has a number of advantages, such as a locally optimal design of the experiment can be easily constructed, and dynamic experimentation is made possible. Using artificial data, the performance of our sequential method is illustrated.

This paper describes the conjunctive counterpart of De Boeck and Rosenberg's hierarchical classes model. Both the original model and its conjunctive counterpart represent the set-theoretical structure of a two-way two-mode binary matrix. However, unlike the original model, the new model represents the row-column association as a conjunctive function of a set of hypothetical binary variables. The conjunctive nature of the new model further implies that it may represent some conjunctive higher order dependencies among rows and columns. The substantive significance of the conjunctive model is illustrated with empirical applications. Finally, it is shown how conjunctive and disjunctive hierarchical classes models relate to Galois lattices, and how hierarchical classes analysis can be useful to construct lattice models of empirical data.

About 27 years ago, a small group of students met with Professor Thurstone in Chicago to discuss methods of encouraging quantitative work in psychology. The initial group that was concerned about the slow rate of development of quantitative work in psychology included Jack Dunlap, Al Kurtz, Marion Richardson, John Stalnaker, G. Frederic Kuder, and Paul Horst. They had discussed the problem, had been helped a bit by Donald Paterson, and had decided that possibly if a magazine were set up to publish quantitative psychological material this would facilitate the development of the field. Persons who did good quantitative work, either theoretical or experimental, would thus have a forum where it would be accepted because it was high quality quantitative work, rather than being rejected because it was quantitative and hence “not of too great interest” to the readers.

In studies involving judgments of similarity or dissimilarity, a variety of other variables may also be measured. Examples might be direct ratings of the stimuli, pairwise preference judgments, and physical measurements of the stimuli with respect to various properties. In such cases, there are important advantages to joint analyses of the dissimilarity and collateral variables. A variety of models are described for relating these and algorithms described for fitting these to data. A number of hypothesis tests are developed and an example offered.

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.

In calculations of the discriminating-power parameter of the normal ogive model, Bock and Lieberman compared estimates derived from their maximum-likelihood solution with those derived from the heuristic solution. The two sets of estimates were in excellent agreement provided the heuristic solution used accurate tetrachoric correlation coefficients. Three computer methods for the calculation of the tetrachoric correlation were examined for accuracy and speed. The routine by Saunders was identified as an acceptably accurate method for calculating the tetrachoric correlation coefficient.

Latent class models in the social and behavioral sciences have remained structurally simple. One reason for this is that inference in statistical models can be computationally difficult. Methods for approximate inference, known as variational approximations, which have been developed in the machine learning, graphical modeling and statistical physics literatures, can be used to alleviate the computational difficulties of inference for latent variable models. The aim of the present article is to set these methods alongside some social and behavioral science literature to which they are relevant, and in particular to consider their potential for “categorical causal modeling”, using latent class analysis. We have collated a number of popular categorical-data models with latent variables and causal structure, typically incorporating a Markovian structure. The efficacy of the approximation methods has been demonstrated through simulations related to an important behavioral science model.