In this paper, we consider a class of models for two-way matrices with binary entries of 0 and 1. First, we consider Boolean matrix decomposition, conceptualize it as a latent response model (LRM) and, by making use of this conceptualization, generalize it to a larger class of matrix decomposition models. Second, probability matrix decomposition (PMD) models are introduced as a probabilistic version of this larger class of deterministic matrix decomposition models. Third, an algorithm for the computation of the maximum likelihood (ML) and the maximum a posteriori (MAP) estimates of the parameters of PMD models is presented. This algorithm is an EM-algorithm, and is a special case of a more general algorithm that can be used for the whole class of LRMs. And fourth, as an example, a PMD model is applied to data on decision making in psychiatric diagnosis.

A general approach to the analysis of subjective categorical data is considered, in which agreement matrices of two or more raters are directly expressed in terms of error and agreement parameters. The method provides focused analyses of ratings from several raters for whom ratings have measurement error distributions that may induce bias in the evaluation of substantive questions of interest. Each rater's judgment process is modeled as a mixture of two components: an error variable that is unique for the rater in question as well as an agreement variable that operationalizes the “true” values of the units of observation. The statistical problems of identification, estimation, and testing of such measurement models are discussed.

The general model is applied in several special cases. The most simple situation is that underlying Cohen's Kappa, where two raters place units into unordered categories. The model provides a generalization and systematization of the Kappa-idea to correct for agreement by chance. In applications with typical research designs, including a between-subjects design and a mixed within-subjects, between-subjects design, the model is shown to disentangle structural and measurement components of the observations, thereby controlling for possible confounding effects of systematic rater bias. Situations considered include the case of more than two raters as well as the case of ordered categories. The different analyses are illustrated by means of real data sets.

Spanning nearly sixty years of research, statistical network analysis has passed through (at least) two generations of researchers and models. Beginning in the late 1930's, the first generation of research dealt with the distribution of various network statistics, under a variety of null models. The second generation, beginning in the 1970's and continuing into the 1980's, concerned models, usually for probabilities of relational ties among very small subsets of actors, in which various simple substantive tendencies were parameterized. Much of this research, most of which utilized log linear models, first appeared in applied statistics publications.

But recent developments in social network analysis promise to bring us into a third generation. The Markov random graphs of Frank and Strauss (1986) and especially the estimation strategy for these models developed by Strauss and Ikeda (1990; described in brief in Strauss, 1992), are very recent and promising contributions to this field. Here we describe a large class of models that can be used to investigate structure in social networks. These models include several generalizations of stochastic blockmodels, as well as models parameterizing global tendencies towards clustering and centralization, and individual differences in such tendencies. Approximate model fits are obtained using Strauss and Ikeda's (1990) estimation strategy.

In this paper we describe and extend these models and demonstrate how they can be used to address a variety of substantive questions about structure in social networks.

Some historical background and preliminary technical information are first presented, and then a number of hidden, but important, methodological aspects of dual scaling are illustrated and discussed: normed versus projected weights, the amount of information accounted for by each solution, a perfect solution to the problem of multidimensional unfolding, multidimensional quantification space, graphical display, number-of-option problems, option standardization versus item standardization, and asymmetry of symmetric (dual) scaling. Contrary to the common perception that dual scaling and similar quantification methods are now mathematically transparent, the present study demonstrates how much more needs to be clarified for routine use of the method to arrive at valid conclusions. Data analysis must be carried out in such a way that common sense, intuition and sound logic will prevail.

We develop simple noniterative estimators of the polyserial correlation coefficient. A general relationship between the polyserial correlation and the point polyserial correlation is exploited to give extensions of Pearson's, Brogden's, and Lord's biserial estimators to the multicategory setting. The small sample and asymptotic properties of these estimators are studied in some detail. A comparison with maximum likelihood estimates shows that Lord's polyserial estimator is fairly efficient across three probability models.

For each Rasch (Masters) partial credit item, there exists a set of independent Rasch binary and indecomposable trinary items for which the sum of the scores and the partial credit score have identical probability density functions. If each indecomposable trinary item is further expressed as the sum of two binary items, then the binary items are positively dependent and cannot be both of the Rasch type.

Multitrait-Multimethod (MTMM) matrices are often analyzed by means of confirmatory factor analysis (CFA). However, fitting MTMM models often leads to improper solutions, or non-convergence. In an attempt to overcome these problems, various alternative CFA models have been proposed, but with none of these the problem of finding improper solutions was solved completely. In the present paper, an approach is proposed where improper solutions are ruled out altogether and convergence is guaranteed. The approach is based on constrained variants of components analysis (CA). Besides the fact that these methods do not give improper solutions, they have the advantage that they provide component scores which can later on be used to relate the components to external variables. The new methods are illustrated by means of simulated data, as well as empirical data sets.

This paper is concerned with the geometric properties of dissimilarity coefficients defined on finite sets and especially with their Euclidean nature. We present several particular transformations which preserve Euclideanarity and we complete, through the study of a one-parameter family, the current knowledge of the metric and Euclidean structure of coefficients based on binary data. These results are directly deduced from two theorems which prove the positive semi-definite status of some quadratic forms which play a large role in some definitions of dissimilarity commonly used.

Chang and Stout (1993) presented a derivation of the asymptotic posterior normality of the latent trait given examinee responses under nonrestrictive nonparametric assumptions for dichotomous IRT models. This paper presents an extention of their results to polytomous IRT models in a fairly straightforward manner. In addition, a global information function is defined, and the relationship between the global information function and the currently used information functions is discussed. An information index that combines both the global and local information is proposed for adaptive testing applications.

Item response theory models posit latent variables to account for regularities in students' performances on test items. Wilson's “Saltus” model extends the ideas of IRT to development that occurs in stages, where expected changes can be discontinuous, show different patterns for different types of items, or even exhibit reversals in probabilities of success on certain tasks. Examples include Piagetian stages of psychological development and Siegler's rule-based learning. This paper derives marginal maximum likelihood (MML) estimation equations for the structural parameters of the Saltus model and suggests a computing approximation based on the EM algorithm. For individual examinees, empirical Bayes probabilities of learning-stage are given, along with proficiency parameter estimates conditional on stage membership. The MML solution is illustrated with simulated data and an example from the domain of mixed number subtraction.

Loglinear unidimensional and multidimensional Rasch models are considered for the analysis of repeated observations of polytomous indicators with ordered response categories. Reparameterizations and parameter restrictions are provided which facilitate specification of a variety of hypotheses about latent processes of change. Models of purely quantitative change in latent traits are proposed as well as models including structural change. A conditional likelihood ratio test is presented for the comparison of unidimensional and multiple scales Rasch models. In the context of longitudinal research, this renders possible the statistical test of homogeneity of change against subject-specific change in latent traits. Applications to two empirical data sets illustrate the use of the models.

Joint correspondence analysis is a technique for constructing reduced-dimensional representations of pairwise relationships among categorical variables. The technique was proposed by Greenacre as an alternative to multiple correspondence analysis. Joint correspondence analysis differs from multiple correspondence analysis in that it focuses solely on between-variable relationships. Greenacre described one alternating least-squares algorithm for conducting joint correspondence analysis. Another alternating least-squares algorithm is described in this article. The algorithm is guaranteed to converge, and does so in fewer iterations than does the algorithm proposed by Greenacre. A modification of the algorithm for handling Heywood cases is described. The algorithm is illustrated on two data sets.

Generalized bilinear models are presented for the statistical analysis of two-way arrays. These models combine bilinear models and generalized linear modeling, and yield a family of models that includes many existing models, as well as suggest other potentially useful ones. This approach both unifies and extends models for two-way arrays, including the ability to treat response and explanatory variables differently in the models, and the incorporation of external information about the variables directly into the analysis. A unifying framework for the generalized bilinear models is provided by considering four particular cases which have been proposed and used in the existing statistical literature. A three-step procedure is proposed to analyze data sets by generalized bilinear models. Two data sets of different nature are analyzed.

The RC(M) association model (Goodman, 1979, 1985, 1986, 1991) is useful for analyzing the relationship between the variables of a 2-way cross-classification. The models presented here are generalizations of the RC(M) association model for 3-way tables. The family of models proposed here, “3-mode association” models, use Tucker's 3-mode components model (Tucker, 1964, 1966; Kroonenberg, 1983) to represent either the three factor interaction or the combined effects of two and three factor interactions. An example from a study in developmental psychology (Kramer & Gottman, 1992) is provided to illustrate the usefulness of the proposed models.

As a multivariate model of the number of events, Rasch's multiplicative Poisson model is extended such that the parameters for individuals in the prior gamma distribution have continuous covariates. The parameters for individuals are integrated out and the hyperparameters in the prior distribution are estimated by a numerical method separately from difficulty parameters that are treated as fixed parameters or random variables. In addition, a method is presented for estimating parameters in Rasch's model with missing values.

The purpose of this paper is to present a hypothesis testing and estimation procedure, Crossing SIBTEST, for detecting crossing DIF. Crossing DIF exists when the difference in the probabilities of a correct answer for the two examinee groups changes signs as ability level is varied. In item response theory terms, crossing DIF is indicated by two crossing item characteristic curves. Our new procedure, denoted as Crossing SIBTEST, first estimates the matching subtest score at which crossing occurs using least squares regression analysis. A Crossing SIBTEST statistic then is used to test the hypothesis of crossing DIF. The performance of Crossing SIBTEST is evaluated in this study.

This paper presents a stochastic multidimensional unfolding (MDU) procedure to spatially represent individual differences in phased or sequential decision processes. The specific application or scenario to be discussed involves the area of consumer psychology where consumers form judgments sequentially in their awareness, consideration, and choice set compositions in a phased or sequential manner as more information about the alternative brands in a designated product/service class are collected. A brief review of the consumer psychology literature on these nested congnitive sets as stages in phased decision making is provided. The technical details of the proposed model, maximum likelihood estimation framework, and algorithm are then discussed. A small scale Monte Carlo analysis is presented to demonstrate estimation proficiency and the appropriateness of the proposed model selection heuristic. An application of the methodology to capture awareness, consideration, and choice sets in graduate school applicants is presented. Finally, directions for future research and other potential applications are given.

This paper introduces dynamic latent-class models for the analysis and interpretation of stability and change in recurrent choice data. These latent-class models provide a nonparametric representation of individual taste differences. Changes in preferences are modeled by allowing for individual-level transitions from one latent class to another over time. The most general model facilitates a saturated representation of class membership changes. Several special cases are presented to obtain a parsimonious description of latent change mechanisms. An easy to implement EM algorithm is derived for parameter estimation. The approach is illustrated by a detailed analysis of a purchase incidence data set.

The main purpose of this paper is to investigate the sensitivity analysis of structural equation model when minor perturbation is introduced. Some influence measure that based on the general case weight perturbation is derived for the generalized least squares estimation. An influence measure that related to the Cook's distance is also developed for the special case deletion perturbation scheme. Using the proposed methodology, the influential observation in a data set can be detected. Moreover, the general theory can be applied to detect the influential parameters in a model. Finally, some illustrative artificial and real examples are presented.

In a broad class of item response theory (IRT) models for dichotomous items the unweighted total score has monotone likelihood ratio (MLR) in the latent trait θ. In this study, it is shown that for polytomous items MLR holds for the partial credit model and a trivial generalization of this model. MLR does not necessarily hold if the slopes of the item step response functions vary over items, item steps, or both. MLR holds neither for Samejima's graded response model, nor for nonparametric versions of these three polytomous models. These results are surprising in the context of Grayson's and Huynh's results on MLR for nonparametric dichotomous IRT models, and suggest that establishing stochastic ordering properties for nonparametric polytomous IRT models will be much harder.