Independence in Contingency Tables Using Simplicial Geometry

Frequently, contingency tables are generated in a multinomial sampling. Multinomial probabilities are then organized in a table assigning probabilities to each cell. A probability table can be viewed as an element in the simplex. The Aitchison geometry of the simplex identifies independent probability tables as a linear subspace. An important consequence is that, given a probability table, the nearest independent table is obtained by orthogonal projection onto the independent subspace. The nearest independent table is identified as that obtained by the product of geometric marginals, which do not coincide with the standard marginals, except in the independent case. The original probability table is decomposed into orthogonal tables, the independent and the interaction tables. The underlying model is log-linear, and a procedure to test independence of a contingency table, based on a multinomial simulation, is developed. Its performance is studied on an illustrative example.     

A versatile bivariate distribution on a bounded domain: another look at the product moment correlation

The Farlie-Gumbel-Morgenstern (FGM) family has been investigated in detail for various continuous marginals such as Cauchy, normal, exponential, gamma, Weibull, lognormal and others. It has been a popular model for the bivariate distribution with mild dependence. However, bivariate FGMs with continuous marginals on a bounded support discussed in the literature are only those with uniform or power marginals. In this paper we study the bivariate FGM family with marginals given by the recently proposed two-sided power (TSP) distribution. Since this family of bounded continuous distributions is very flexible, the properties of the FGM family with TSP marginals could serve as an indication of the structure of the FGM distribution with arbitrary marginals defined on a compact set. A remarkable stability of the correlation between the marginals has been observed.     

Exact Tests for 2 × 2 Contingency Tables

Fisher's exact test, difference in proportions, log odds ratio, Pearson's chi-squared, and likelihood ratio are compared as test statistics for testing independence of two dichotomous factors when the associated p values are computed by using the conditional distribution given the marginals. The statistics listed above that can be used for a one-sided alternative give identical p values. For a two-sided alternative, many of the above statistics lead to different p values. The p values are shown to differ only by which tables in the opposite tail from the observed table are considered more extreme than the observed table.     

Bivariate inverse Weibull distribution

ABSTRACT

Recently it is observed that the inverse Weibull (IW) distribution can be used quite effectively to analyse lifetime data in one dimension. The main aim of this paper is to define a bivariate inverse Weibull (BIW) distribution so that the marginals have IW distributions. It is observed that the joint probability density function and the joint cumulative distribution function can be expressed in compact forms. Several properties of this distribution such as marginals, conditional distributions and product moments have been discussed. We obtained the maximum likelihood estimates for the unknown parameters of this distribution and their approximate variance– covariance matrix. We perform some simulations to see the performances of the maximum likelihood estimators. One data set has been re-analysed and it is observed that the bivariate IW distribution provides a better fit than the bivariate exponential distribution.     

Copulas,degenerate distributions and quantile tests in competing risk problems

It is well known that, given observable data for a competing risk problem, there is always an independent model consistent with the data. It has been pointed out, however, that this independent model does not necessarily have to be one with proper marginals. One purpose of this paper is to explore the extent to which one might try to use the non-parametric assumption that the marginals are proper in order to test whether or not independence holds. This will lead us naturally to a closely related estimation problem—how we estimate the marginals given that a certain quantile of one variable is reached at the same time as a given quantile of the other variable. The problem will be considered using the copula-based approach of Zheng and Klein. Two different methods are discussed. One is a non-parametric maximum likelihood method. The other is a consistent estimator, called the bilinear adjustment estimator, that can be computed quickly and which thus lends itself more readily to simulation methods such as bootstrapping.The ultimate objective of the work in this paper is to provide an additional analytical tool to understand the effectiveness of preventive maintenance. Preventive maintenance censors component failure data and thus provides an important example of competing risk within the reliability area.     

A generalised NGINAR(1) process with inflated‐parameter geometric counting series

In this paper we propose a new stationary first‐order non‐negative integer valued autoregressive process with geometric marginals based on a generalised version of the negative binomial thinning operator. In this manner we obtain another process that we refer to as a generalised stationary integer‐valued autoregressive process of the first order with geometric marginals. This new process will enable one to tackle the problem of overdispersion inherent in the analysis of integer‐valued time series data, and contains the new geometric process as a particular case. In addition various properties of the new process, such as conditional distribution, autocorrelation structure and innovation structure, are derived. We discuss conditional maximum likelihood estimation of the model parameters. We evaluate the performance of the conditional maximum likelihood estimators by a Monte Carlo study. The proposed process is fitted to time series of number of weekly sales (economics) and weekly number of syphilis cases (medicine) illustrating its capabilities in challenging cases of highly overdispersed count data.     

Continuity of f-projections and applications to the iterative proportional fitting procedure*

This paper proves continuity of f-projections and the continuous dependence of the limit matrix of the iterative proportional fitting procedure (IPF procedure) on the given matrix as well as the given marginals under certain regularity constraints. For finite spaces, the concept of f-projections of finite measures on a compact and convex set is introduced and continuity of f-projections is proven. This result is applied to the IPF procedure. Given a nonnegative matrix as well as row and column marginals the IPF procedure generates a sequence of matrices, called the IPF sequence, by alternately fitting rows and columns to match their respective marginals. The procedure is equivalent to cyclic f-projections. If the IPF sequence converges, the application of the continuity of f-projections yields the continuous dependence of the limit matrix on the given matrix. By generalized convex programming and under some constraints, it is shown that the limit matrix of the IPF sequence continuously depends not only on the given matrix but also on the marginals.     

Inequalities for bivariate distribution with x ≤ y and marginals given

Upper and lower bounds on the joint bivariate distribution are found when the marginals are given and under the additional condition that X ≤ Y with probability one. The upper bound is the same as for the unrestricted bivariate distribution with marginals given, For the lower bound a simple inequality is derived which is exact, that is, achievable, in many cases including normal and exponential marginals.     

On the minimum expected quantity for the validity of the chi-squared test in 2 2 2 tables

A 2 2 2 contingency table can often be analysed in an exact fashion by using Fisher's exact test and in an approximate fashion by using the chi-squared test with Yates' continuity correction, and it is traditionally held that the approximation is valid when the minimum expected quantity E is E S 5. Unfortunately, little research has been carried out into this belief, other than that it is necessary to establish a bound E>E*, that the condition E S 5 may not be the most appropriate (Martín Andrés et al., 1992) and that E* is not a constant, but usually increasing with the growth of the sample size (Martín Andrés & Herranz Tejedor, 1997). In this paper, the authors conduct a theoretical experimental study from which they ascertain that E* value (which is very variable and frequently quite a lot greater than 5) is strongly related to the magnitude of the skewness of the underlying hypergeometric distribution, and that bounding the skewness is equivalent to bounding E (which is the best control procedure). The study enables estimating the expression for the above-mentioned E* (which in turn depends on the number of tails in the test, the alpha error used, the total sample size, and the minimum marginal imbalance) to be estimated. Also the authors show that E* increases generally with the sample size and with the marginal imbalance, although it does reach a maximum. Some general and very conservative validity conditions are E S 35.53 (one-tailed test) and E S 7.45 (two-tailed test) for alpha nominal errors in 1% h f h 10%. The traditional condition E S 5 is only valid when the samples are small and one of the marginals is very balanced; alternatively, the condition E S 5.5 is valid for small samples or a very balanced marginal. Finally, it is proved that the chi-squared test is always valid in tables where both marginals are balanced, and that the maximum skewness permitted is related to the maximum value of the bound E*, to its value for tables with at least one balanced marginal and to the minimum value that those marginals must have (in non-balanced tables) for the chi-squared test to be valid.     

Approximate theory of multiway block designs

Optimality properties of multiway block designs are deduced from the general results of J. Kiefer's approximate-design theory. In the model with additive effects these optimality properties solely depend on the two-dimensional marginals of the designs. Uniform designs, and designs whose two-dimensional marginals are products of the one-dimensional marginals, are shown to be optimal. Approximate Youden designs are introduced for the case when the support sets of the two-dimensional marginals are prescribed in advance. They are optimal in a relatively small class of competing designs only.     

Maximum Likelihood Estimation in the Bivariate Binomial (0,1) Distribution:Application TO 2×2 Tables

We consider a 2×2 contingency table, with dichotomized qualitative characters (A,A) and (B,B), as a sample of size n drawn from a bivariate binomial (0,1) distribution. Maximum likelihood estimates p?₁p?₂ and p? are derived for the parameters of the two marginals p₁p₂ and the coefficient of correlation. It is found that p? is identical to Pearson's (1904)?=(χ²/n)½, where ?² is Pearson's usual chi-square for the 2×2 table. The asymptotic variance-covariance matrix of p?_lp?₂and p is also derived.     

Power-normal distribution

Recently, Gupta and Gupta [Analyzing skewed data by power-normal model, Test 17 (2008), pp. 197–210] proposed the power-normal distribution for which normal distribution is a special case. The power-normal distribution is a skewed distribution, whose support is the whole real line. Our main aim of this paper is to consider bivariate power-normal distribution, whose marginals are power-normal distributions. We obtain the proposed bivariate power-normal distribution from Clayton copula, and by making a suitable transformation in both the marginals. Lindley–Singpurwalla distribution also can be used to obtain the same distribution. Different properties of this new distribution have been investigated in detail. Two different estimators are proposed. One data analysis has been performed for illustrative purposes. Finally, we propose some generalizations to multivariate case also along the same line and discuss some of its properties.     

A family of bivariate distributions with non-elliptical contours

In this paper, a family of copulas with two parameters is proposed and its dependence analysis is performed. The corresponding family of bivariate distributions with specified marginals is constructed. For normal marginals, the new distributions are non-elliptical and can be applied in data analysis. They provide various alternative hypotheses for testing normality. Finally, an example is given.     

Truncation invariant dependence structures

In this paper, a special class of m-dimensional distribution functions which can be uniquely determined in terms of their 2-dimensional marginals is studied. The members of the class can be characterized as having truncation invariant dependence structure. The representation given in this paper provides a physical meaning to the multivariate Cook-Johnson distribution, and introduces a systematic way of generating higher dimensional distributions by using rich 2-dimensional distributions provided that the 2-dimensional marginals are compatible. A class of 3-dimensional multivariate normal distribution has been generated and bounds in terms of lower dimensional marginals are provided.     

Dependence function for continuous bivariate densities

The notion of cross-product ratio for discrete two-way contingency table is extended to the case of continuous bivariate densities. This results in the “local dependence function” that measues the margin-free dependence between bivariate random variables. Properties and examples of the dependence function are discussed. The bivariate normal density plays a special role since it has constant dependence. Continuous bivariate densities can be constructed by specifying the dependence function along with two marginals in analogy to the construction of two-way contingency tables given marginals and patterns of interaction. The dependence function provides a partial ordering on bivariate dependence.     

On a class of bivariate mixed Sarmanov distributions

Multivariate distributions are more and more used to model the dependence encountered in many fields. However, classical multivariate distributions can be restrictive by their nature, while Sarmanov's multivariate distribution, by joining different marginals in a flexible and tractable dependence structure, often provides a valuable alternative. In this paper, we introduce some bivariate mixed Sarmanov distributions with the purpose to extend the class of bivariate Sarmanov distributions and to obtain new dependency structures. Special attention is paid to the bivariate mixed Sarmanov distribution with Poisson marginals and, in particular, to the resulting bivariate Sarmanov distributions with negative binomial and with Poisson‐inverse Gaussian marginals; these particular types of mixed distributions have possible applications in, for example modelling bivariate count data. The extension to higher dimensions is also discussed. Moreover, concerning the dependency structure, we also present some correlation formulas.     

Dependent models for observations which include angular ones

This paper discusses some stochastic models for dependence of observations which include angular ones. First, we provide a theorem which constructs four-dimensional distributions with specified bivariate marginals on certain manifolds such as two tori, cylinders or discs. Some properties of the submodel of the proposed models are investigated. The theorem is also applicable to the construction of a related Markov process, models for incomplete observations, and distributions with specified marginals on the disc. Second, two maximum entropy distributions on the cylinder are discussed. The circular marginal of each model is distributed as the generalized von Mises distribution which represents a symmetric or asymmetric, unimodal or bimodal shape. The proposed cylindrical model is applied to two data sets.     

A note on dependence modeling for Bernoulli variables

Understanding and modeling multivariate dependence structures depending upon the direction are challenging but an interest of theoretical and applied researchers. In this paper, we propose a characterization of tables generated by Bernoulli variables through the uniformization of the marginals and refer to them as Q-type tables. The idea is similar to the copulas. This approach helps to see the dependence structure clearly by eliminating the effect of the marginals that have nothing to do with the dependence structure. We define and study conditional and unconditional Q-type tables and provide various applications for them. The limitations of existing approaches such as Cochran-Mantel-Haenszel pooled odds ratio are discussed, and a new one that stems naturally from our approach is introduced.     

Parameter estimation for the dirichlet-multinomial distribution using supplementary beta-binomial data

We develop estimates for the parameters of the Dirichlet-multinomial distribution (DMD) when there is insufficient data to obtain maximum likelihood or method of moment estimates known in the literature. We do, however, have supplemetary beta-binomial data pertaining to the marginals of the DMD, and use these data when estimating the DMD parameters. A real situation and data set are given where our estimates are applicable.     

A Convergent Iterative Procedure for Constructing Bivariate Distributions

For many years there has been interest in families of bivariate distributions with the marginals as parameters. Questions of this kind arise if one is to build a stochastic model in a situation where one has some idea about the dependence structure and marginal distributions. In this article, among all bivariate distributions which satisfy the constraints imposed by the known marginals and/or dependence structure, one that has the maximum entropy is obtained by using iterative procedure, and its convergence is proved.