Simulation for the complex Bingham distribution

The complex Bingham distribution is relevant for the shape analysis of landmark data in two dimensions. In this paper it is shown that the problem of simulating from this distribution reduces to simulation from a truncated multivariate exponential distribution. Several simulation methods are described and their efficiencies are compared.     

ESTIMATION OF THE CONCENTRATION PARAMETERS OF THE FISHER MATRIX DISTRIBUTION ON 50(3) AND THE BINGHAM DISTRIBUTION ON Sq, q≥ 2

Two procedures are considered for estimating the concentration parameters of the Fisher matrix distribution for rotations or orientations in three dimensions. The first is maximum likelihood. The use of a convenient 1-dimensional integral representation of the normalising constant, which greatly simplifies the computation, is suggested. The second approach exploits the equivalence of the Fisher distribution for rotations in three dimensions, and the Bingham distribution for axes in four dimensions. We describe a pseudo likelihood procedure which works for the Bingham distribution in any dimension. This alternative approach does not require numerical integration. Results on the asymptotic efficiency of the pseudo likelihood estimator relative to the maximum likelihood estimator are given, and the two estimators are compared in the analysis of a well-known vectorcardiography dataset.     

TESTING THE LARGEST OF A SET OF CORRELATION COEFFICIENTS

A previous paper which studied the distribution of the smallest distance between N independent random points on the surface of a sphere is generalised to higher dimensions in order to study the distribution of the largest sample correlation coefficient between a set of independent normally distributed variables. Inclusion-exclusion arguments provide accurate bounds for the tail of this distribution, and by another argument more exact bounds are also found, one of which is an improvement on the result in the previous paper. Bounds are also found for the power of the test against the alternative hypothesis that one only of the population correlation coefficients is non-zero. The test is also shown to be the likelihood ratio test against the latter alternative.     

The centered normal conditionals distribution

In this paper, we discuss some aspects of the distribution theory associated with the centered bivariate normal conditionals distribution including discussion of its marginal distributions. We calculate the maximum likelihood and pseudolikelihood estimators. We propose a simplified moment based method of estimation. Finally, we discuss generalizations to higher dimensions.     

Spatial Clustering Tests Based on the Domination Number of a New Random Digraph Family

We use the domination number of a parametrized random digraph family called proportional-edge proximity catch digraphs (PCDs) for testing multivariate spatial point patterns. This digraph family is based on relative positions of data points from various classes. We extend the results on the distribution of the domination number of proportional-edge PCDs, and use the domination number as a statistic for testing segregation and association against complete spatial randomness. We demonstrate that the domination number of the PCD has binomial distribution when size of one class is fixed while the size of the other (whose points constitute the vertices of the digraph) tends to infinity and has asymptotic normality when sizes of both classes tend to infinity. We evaluate the finite sample performance of the test by Monte Carlo simulations and prove the consistency of the test under the alternatives. We find the optimal parameters for testing each of the segregation and association alternatives. Furthermore, the methodology discussed in this article is valid for data in higher dimensions also.     

Nonparametric Bayes methods for directional data

A model for directional data in q dimensions is studied. The data are assumed to arise from a distribution with a density on a sphere of q — 1 dimensions. The density is unimodal and rotationally symmetric, but otherwise of unknown form. The posterior distribution of the unknown mode (mean direction) is derived, and small-sample posterior inference is discussed. The posterior mean of the density is also given. A numerical method for evaluating posterior quantities based on sampling a Markov chain is introduced. This method is generally applicable to problems involving unknown monotone functions.     

Sampling from compositional and directional distributions

This paper describes a method for sampling from a non-standard distribution which is important in both population genetics and directional statistics. Current approaches rely on complicated procedures which do not work well, if at all, in high dimensions and usual parameter set-ups. We use a Gibbs sampler which seems necessary in practical situations of high dimensions.     

Two-way analysis of variance for data from a concentrated bipolar Watson distribution

The bipolar Watson distribution is frequently used for modeling axial data. We extend the one-way analysis of variance based on this distribution to a two-way layout. We illustrate the method with directional data in three dimensions     

Bivariate extension of (dynamic) cumulative past entropy

Recently, the concept of cumulative residual entropy (CRE) has been studied by many researchers in higher dimensions. In this article, we extend the definition of (dynamic) cumulative past entropy (DCPE), a dual measure of (dynamic) CRE, to bivariate setup and obtain some of its properties including bounds. We also look into the problem of extending DCPE for conditionally specified models. Several properties, including monotonicity, and bounds of DCPE are obtained for conditional distributions. It is shown that the proposed measure uniquely determines the distribution function. Moreover, we also propose a stochastic order based on this measure.     

Modified first-difference estimator in a panel data model with unobservable factors both in the errors and the regressors when the time dimension is small

Panel data models with factor structures in both the errors and the regressors have received considerable attention recently. In these models, the errors and the regressors are correlated and the standard estimators are inconsistent. This paper shows that, for such models, a modified first-difference estimator (in which the time and the cross-sectional dimensions are interchanged) is consistent as the cross-sectional dimension grows but the time dimension is small. Although the estimator has a non standard asymptotic distribution, t and F tests have standard asymptotic distribution under the null hypothesis.     

A class of multivariate distribution-free tests of independence based on graphs

A class of distribution-free tests is proposed for the independence of two subsets of response coordinates. The tests are based on the pairwise distances across subjects within each subset of the response. A complete graph is induced by each subset of response coordinates, with the sample points as nodes and the pairwise distances as the edge weights. The proposed test statistic depends only on the rank order of edges in these complete graphs. The response vector may be of any dimensions. In particular, the number of samples may be smaller than the dimensions of the response. The test statistic is shown to have a normal limiting distribution with known expectation and variance under the null hypothesis of independence. The exact distribution free null distribution of the test statistic is given for a sample of size 14, and its Monte-Carlo approximation is considered for larger sample sizes. We demonstrate in simulations that this new class of tests has good power properties for very general alternatives.     

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.     

Alternative definitions of bivariate equilibrium distributions

The concept of equilibrium distribution plays an important role in survival analysis, reliability and insurance studies. If we consider the problem of extending this concept to higher dimensions, we do not have a unique solution. In this paper, alternative definitions of bivariate equilibrium distributions are studied and proposed. The Navarro et al. (2006) proposal is considered and some new results are given. We continue with the Gupta and Sankaran's (1998) definition. Necessary and sufficient conditions for its existence are stated and a characterization theorem is given. As a third alternative, a new definition based on conditional specification is introduced and several results are obtained. Reliability properties of the different versions are studied.     

Approximations to the multivariate normal integral

We propose a simple and efficient way to approximate multivariate normal probabilities using univariate and bivariate probabilities. The approximation is computationally tested for the trivariate and quadrivariate normal probabilities. A few problems of higher dimensions were also tested.     

Comparing methods for estimating the abilities for the multidimensional models of mixed item types

The maximum likelihood (MLE), the weighted maximum likelihood (WMLE), and the maximum a posteriori (MAP or BMLE) have been widely used to estimate ability parameters in item response theory (IRT), and their precisions and biases have been studied and compared. Multidimensional IRT (MIRT) has been shown to provide better subscore estimates in both paper-and-pencil and computer adaptive tests; thus, it is very important to have an accurate score estimate for the MIRT model. The purpose of this article is to compare the performances of the three estimation methods in the MIRT framework for tests of mixed item types that have both dichotomous and polytomously scored items, and for tests of mixed structured items (simple structured and complex structured). It is found that all three methods perform well for all conditions. For all models studied (one-, two-, three-, and four- dimensional model), WMLE has smaller BIAS and higher reliabilities, but larger RMSE and SE. WMLE and MLE are closer to each other than to BMLE. However, for higher dimensions, BMLE is recommended, especially when there are correlations between the dimensions.     

The Folded Logistic Distribution

ABSTRACT

Physical measurements like dimensions, including time, and angles in scientific experiments are frequently recorded without their algebraic sign. The directions of those physical quantities measured with respect to a frame of reference in most practical applications are considered to be unimportant and are ignored. As a consequence, the underlying distribution of measurements is replaced by a distribution of absolute measurements. When the underlying distribution is logistic, the resulting distribution is called the “folded logistic distribution”. Here, the properties of the folded logistic distribution will be presented and the techniques for estimating parameters will be given. The advantages of using this folded logistic distribution over the folded normal distribution will be discussed and some examples will be cited.     

Semiparametric Stein estimators

Stein [Stein, C. (1956). Inadmissibility of the usual estimator for the mean of a multivariate normal distribution. In Proc. 3rd Berkeley symp. math. statist. and pro. (pp. 197–206). University of California Press], in his seminal paper, came up with the surprising discovery that the sample mean is an inadmissible estimator of the population mean in three or higher dimensions under squared error loss. The past five decades have witnessed multiple extensions and variations of Stein’s results. In this paper we develop Stein-type estimators in a semiparametric framework and prove their coordinatewise asymptotic dominance over the sample mean in terms of Bayes risks.     

Generalized information matrix tests for copulas

We propose a family of goodness-of-fit tests for copulas. The tests use generalizations of the information matrix (IM) equality of White and so relate to the copula test proposed by Huang and Prokhorov. The idea is that eigenspectrum-based statements of the IM equality reduce the degrees of freedom of the test’s asymptotic distribution and lead to better size-power properties, even in high dimensions. The gains are especially pronounced for vine copulas, where additional benefits come from simplifications of score functions and the Hessian. We derive the asymptotic distribution of the generalized tests, accounting for the nonparametric estimation of the marginals and apply a parametric bootstrap procedure, valid when asymptotic critical values are inaccurate. In Monte Carlo simulations, we study the behavior of the new tests, compare them with several Cramer–von Mises type tests and confirm the desired properties of the new tests in high dimensions.     

In Memoriam: Bernard George Greenberg 1919-1985

The bivariate normal density with unit variance and correlation ρ is well known. We show that by integrating out ρ, the result is a function of the maximum norm. The Bayesian interpretation of this result is that if we put a uniform prior over ρ, then the marginal bivariate density depends only on the maximal magnitude of the variables. The square-shaped isodensity contour of this resulting marginal bivariate density can also be regarded as the equally weighted mixture of bivariate normal distributions over all possible correlation coefficients. This density links to the Khintchine mixture method of generating random variables. We use this method to construct the higher dimensional generalizations of this distribution. We further show that for each dimension, there is a unique multivariate density that is a differentiable function of the maximum norm and is marginally normal, and the bivariate density from the integral over ρ is its special case in two dimensions.     

Directional statistics and shape analysis

This paper highlights distributional connections between directional statistics and shape analysis. In particular, we provide a test of uniformity for highly dispersed shapes, using the standard techniques of directional statistics. We exploit the isometric transformation from triangular shapes to a sphere in three dimensions, to provide a rich class of shape distributions. A link between the Fisher distribution and the complex Bingham distribution is re-examined. Some extensions to higher-dimensional shapes are outlined.