Bayesian analysis of contingency tables: a simulation and graphics-based approach

In this paper we present a simulation and graphics-based model checking and model comparison methodology for the Bayesian analysis of contingency tables. We illustrate the approach by testing the hypotheses of independence and symmetry on complete and incomplete simulated tables.     

Models for two-way contingency tables

We derive a theoretical and parametrized form for a two-way contingency table by using matrix calculus and Kronecker products.     

Generalised quasi‐symmetry models for ordinal contingency tables

This paper extends the ordinary quasi‐symmetry (QS) model for square contingency tables with commensurable classification variables. The proposed generalised QS model is defined in terms of odds ratios that apply to ordinal variables. In particular, we present QS models based on global, cumulative and continuation odds ratios and discuss their properties. Finally, the conditional generalised QS model is introduced for local and global odds ratios. These models are illustrated through the analysis of two data sets.     

Minimum chi-square estimator and test for incomplete multidimensional contingency tables

A modified transformed chi-square statistic is defined for testing hypotheses of quasi-independence in the incomplete multi-dimensional contingency table and a simple method for determining degrees of freedom is given. A modified transformed chi-squareestimator of the expected cell frequencies is given in closed form for a general class of exact linear constraints. The co-variance matrix of estimated cell frequencies is derived under the assumption of a conditional Poisson distribution.     

Bayesian inference for contingency tables with given marginals

Summary In this paper we introduce a class of prior distributions for contingency tables with given marginals. We are interested in the structrre of concordance/discordance of such tables. There is actually a minor limitation in that the marginals are required to assume only rational values. We do argue, though, that this is not a serious drawback for all applicatory purposes. The posterior and predictive distributions given anM-sample are computed. Examples of Bayesian estimates of some classical indices of concordance are also given. Moreover, we show how to use simulation in order to overcome some difficulties which arise in the computation of the posterior distribution.     

Multidimensional contingency tables with missing data

This paper is concerned wim ine maximum likelihood estimation and the likelihood ratio test for hierarchical loglinear models of multidimensional contingency tables with missing data. The problems of estimation and test for a high dimensional contingency table can be reduced into those for a class of low dimensional tables. In some cases, the incomplete data in the high dimensional table can become complete in the low dimensional tables through the reduction can indicate how much the incomplete data contribute to the estimation and the test.     

Bayesian analysis for incomplete multi-way contingency tables with nonignorable nonresponse

We propose Bayesian methods with five types of priors to estimate cell probabilities in an incomplete multi-way contingency table under nonignorable nonresponse. In this situation, the maximum likelihood (ML) estimates often fall in the boundary solution, causing the ML estimates to become unstable. To deal with such a multi-way table, we present an EM algorithm which generalizes the previous algorithm used for incomplete one-way tables. Three of the five types of priors were previously introduced while the other two are newly proposed to reflect different response patterns between respondents and nonrespondents. Data analysis and simulation studies show that Bayesian estimates based on the old three priors can be worse than the ML regardless of occurrence of boundary solution, contrary to previous studies. The Bayesian estimates from the two new priors are most preferable when a boundary solution occurs. We provide an illustrating example using data for a study of the relationship between a mother's smoking and her newborn's weight.     

A comparison of tests of homogeneity for sparse contingency tables

Five tests of homogeneity for a 2x(k+l) contingency table are compared using Monte Carlo techniques. For these studiesit is assumed that k becomes large in such a way that thecontingency table is sparse for 2xk of the cells, but the sample size in two of the cells remains large. The test statistics studied are: the chi-square approximation to the Pearson test statistic, the chi-square approximation to the likelihood ratio statistic, the normal approximation to Zelterman's (1984)the normal approximation to Pearson's chi-square, and the normal approximation to the likelihood ratio statistic. For the range of parameters studied the chi-square approximation to Pearson's statistic performs consistently well with regard to its size and power.     

Bayesian methods for contingency tables using Gibbs sampling

Cell counts in contingency tables can be smoothed using loglinear models. Recently, sampling-based methods such as Markov chain Monte Carlo (MCMC) have been introduced, making it possible to sample from posterior distributions. The novelty of the approach presented here is that all conditional distributions can be specified directly, so that straight-forward Gibbs sampling is possible. Thus, the model is constructed in a way that makes burn-in and checking convergence a relatively minor issue. The emphasis of this paper is on smoothing cell counts in contingency tables, and not so much on estimation of regression parameters. Therefore, the prior distribution consists of two stages. We rely on a normal nonconjugate prior at the first stage, and a vague prior for hyperparameters at the second stage. The smoothed counts tend to compromise between the observed data and a log-linear model. The methods are demonstrated with a sparse data table taken from a multi-center clinical trial. The research for the first author was supported by Brain Pool program of the Korean Federation of Science and Technology Societies. The research for the second author was partially supported by KOSEF through Statistical Research Center for Complex Systems at Seoul National University.     

The stability of several measures of association in small contingency tables

Measures of association are often used to describe the relationship between row and column variables in two—dimensional contingency tables. It is not uncommon in biomedical research to categorize continuous variables to obtain a two—dimensional table. In these situations it is desirable that the measure of association not be too sensitive to changes in the number of categories or to the choice of cut points. To accomplish this objective we attempt to find a measure of association that closely approximates the corresponding measure of association for the underlying distribution.Measures that are close to the underlying measure for various table sizes andcutpoints are called stable measures.     

Bayesian Analysis of Contingency Tables

ABSTRACT

The display of the data by means of contingency tables is used in different approaches to statistical inference, for example, to broach the test of homogeneity of independent multinomial distributions. We develop a Bayesian procedure to test simple null hypotheses versus bilateral alternatives in contingency tables. Given independent samples of two binomial distributions and taking a mixed prior distribution, we calculate the posterior probability that the proportion of successes in the first population is the same as in the second. This posterior probability is compared with the p-value of the classical method, obtaining a reconciliation between both results, classical and Bayesian. The obtained results are generalized for r × s tables.     

Bayesian estimation of proportions with a cross-entropy prior

This paper suggests estimators of the frequencies (N₈) or proportions {N₈/N) of N distinguishable objects contained in S categories; given various types of information, We consider information in the form of exact constraints on the N₈, sample frequencies, and frequencies of related data, The analysis uses Bayesian methods, where the prior distribution is assumed to be a function of the cross-entropy between the N₈ and a reference distribution, We show the relationship between our estimator and the log-linear and logit models and also present a sampling experiment to compare our proposed estimator with the iterated proportional fitting estimator.     

Comparison of dependent two-way contingency tables

Comparison of two-way contingency tables using measures of association is considered. Multiple comparison procedures for dependent tables are proposed, enabling us to compare tables that are faces from larger multl-dimensional tables. An example

is given to Illustrate the analysis of two 2 × 2-tables formed

from a 2⁴-table.     

Bayesian analysis of quasi-life tables

Quasi-life tables, in which the data arise from many concurrent, independent, discrete-time renewal processes, were defined by Baxter (1994, Biometrika 81:567–577), who outlined some methods for estimation. The processes are not observed individually; only the total numbers of renewals at each time point are observed. Crowder and Stephens (2003, Lifetime Data Anal 9:345–355) implemented a formal estimating-equation approach that invokes large-sample theory. However, these asymptotic methods fail to yield sensible estimates for smaller samples. In this paper, we implement a Bayesian analysis based on MCMC computation that works equally well for large and small sample sizes. We give three simulated examples, studying the Bayesian results, the impact of changing prior specification, and empirical properties of the Bayesian estimators of the lifetime distribution parameters. We also study the Baxter (1994, Biometrika 81:567–577) data, and uncover structure that has not been commented upon previously.     

Outlier detection in contingency tables using decomposable graphical models

For high-dimensional data, it is a tedious task to determine anomalies such as outliers. We present a novel outlier detection method for high-dimensional contingency tables. We use the class of decomposable graphical models to model the relationship among the variables of interest, which can be depicted by an undirected graph called the interaction graph. Given an interaction graph, we derive a closed-form expression of the likelihood ratio test (LRT) statistic and an exact distribution for efficient simulation of the test statistic. An observation is declared an outlier if it deviates significantly from the approximated distribution of the test statistic under the null hypothesis. We demonstrate the use of the LRT outlier detection framework on genetic data modeled by Chow–Liu trees.     

A nonparametric procedure for testing partially ranked data

In consumer preference studies, it is common to seek a complete ranking of a variety of, say N, alternatives or treatments. Unfortunately, as N increases, it becomes progressively more confusing and undesirable for respondents to rank all N alternatives simultaneously. Moreover, the investigators may only be interested in consumers’ top few choices. Therefore, it is desirable to accommodate the setting where each survey respondent ranks only her/his most preferred k (k?N) alternatives. In this paper, we propose a simple procedure to test the independence of N alternatives and the top-k ranks, such that the value of k can be predetermined before securing a set of partially ranked data or be at the discretion of the investigator in the presence of complete ranking data. The asymptotic distribution of the proposed test under root-n local alternatives is established. We demonstrate our procedure with two real data sets.     

A model for comparisons with ordered categorical data

A model developed by Andrich for ordered categorical data is extended to develop tests for treatment effects with paired or matched samples. In particular, this includes analysis for pre-post studies and crossover designs. Some advantages of this model are that it allows for misclassification of subjects, yields reasonable conditional requirements for exact analysis, a normal approximation is good for all but the smallest of sample sizes, and it is relatively simple mathematically. Furthermore, the form of the tests derived are logical extensions of tests for unordered categories.