Odds ratio for a single 2 × 2 table with correlated binomials for two margins

While analyzing 2 × 2 contingency tables, the log odds ratio for measuring the strength of association is often approximated by a normal distribution with some variance. We show that the expression of that variance needs to be modified in the presence of correlation between two binomial distributions of the contingency table. In the present paper, we derive a correlation-adjusted variance of the limiting normal distribution of log odds ratio. We also propose a correlation adjusted test based on the standard odds ratio for analyzing matched-pair studies and any other study settings that induce correlated binary outcomes. We demonstrate that our proposed test outperforms the classical McNemar’s test. Simulation studies show the gains in power are especially manifest when sample size is small and strong correlation is present. Two examples of real data sets are used to demonstrate that the proposed method may lead to conclusions significantly different from those reached using McNemar’s test.     

Data-driven smooth tests of the proportional hazards assumption

A new test of the proportional hazards assumption in the Cox model is proposed. The idea is based on Neyman’s smooth tests. The Cox model with proportional hazards (i.e. time-constant covariate effects) is embedded in a model with a smoothly time-varying covariate effect that is expressed as a combination of some basis functions (e.g., Legendre polynomials, cosines). Then the smooth test is the score test for significance of these artificial covariates. Furthermore, we apply a modification of Schwarz’s selection rule to choosing the dimension of the smooth model (the number of the basis functions). The score test is then used in the selected model. In a simulation study, we compare the proposed tests with standard tests based on the score process.     

A test of concordance based on Gini’s mean difference

A new rank correlation index, which can be used to measure the extent of concordance or discordance between two rankings, is proposed. This index is based on Gini’s mean difference computed on the totals ranks corresponding to each unit and it turns out to be a special case of a more general measure of the agreement of m rankings. The proposed index can be used in a test for the independence of two criteria used to rank the units of a sample, against their concordance/discordance. It can then be regarded as a competitor of other classical methods, such as Kendall’s tau. The exact distribution of the proposed test-statistic under the null hypothesis of independence is studied and its expectation and variance are determined; moreover, the asymptotic distribution of the test-statistic is derived. Finally, the implementation of the proposed test and its performance are discussed. Both the authors contributed equally to this work; however, the actual writing of the paper was as follows: Sects. 2 and 3 are due to C. G. Borroni, Sects. 1 and 4 are due to M. Zenga.     

Statistical regularity of firm size distribution: the Pareto IV and truncated Yule for Italian SCI manufacturing

In this paper we model the firm size distribution (FSD) of Italian manufacturing firms of SCI, the GDP survey of ISTAT, by a continuous and a discrete distribution: the Pareto IV distribution on total assets and the Yule distribution on Number of Employees. The Pareto IV distribution is characterized by four parameters and shows a better fit than both the Lognormal and Pareto I, which are the distributions more frequently applied to model firm size. The Pareto IV is inconsistent with Gibrat’s Law according to which the different segments of an Industry are characterized by proportionate growth and the distribution of size is Lognormal. A truncation of the Yule distribution has been necessary because the dataset is characterized by firms with at least 20 employees. The truncated Yule distribution shows a good fit for medium–large firms (firms with more than 50 employees). The partition of the dataset in innovative and non-innovative firms – both of which are well described by the Pareto IV – reveals a beneficial effect of scale on innovation. Finally, the good fit of both distributions holds not only for the composite industry, but for the single sectors too. The present work is part of a more general research project: “Industry evolution: innovation, profitability and firm’s growth”, conducted within the Department of Economic and Social Sciences of the Università Cattolica del Sacro Cuore (UCSC), Piacenza, coordinated by Professor Maurizio Baussola in cooperation with ISTAT (Italian Statistical Office, regional office for Lombardy). Part of this research was done when Lisa Crosato was a visiting research student at the LSE Statistics Department, during her Ph.D program in “Quantitative Models for Policy Analysis” at the UCSC of Piacenza.     

Bayesian analysis of paired survival data using a bivariate exponential distribution

We consider a Bayesian analysis method of paired survival data using a bivariate exponential model proposed by Moran (1967, Biometrika 54:385–394). Important features of Moran’s model include that the marginal distributions are exponential and the range of the correlation coefficient is between 0 and 1. These contrast with the popular exponential model with gamma frailty. Despite these nice properties, statistical analysis with Moran’s model has been hampered by lack of a closed form likelihood function. In this paper, we introduce a latent variable to circumvent the difficulty in the Bayesian computation. We also consider a model checking procedure using the predictive Bayesian P-value.     

A quadratic approximation for Jackknife estimators of the variance of sample mean functions

Jackknife estimators of the variance of estimators which are functions of the sample mean are considered. A quadratic approximation of them is proposed and compared with a linear approximation by Monte Carlo experiments carried out by statistical software Minitab.     

Bivariate semi α-Laplace distribution and processes

The Laplace distribution is considered as a better choice for modeling whenever data exhibit high kurtosis and heavier tails than Gaussian tails. Even though this is the case, not much work has been done on bivariate Laplace distribution. In this work, we introduce and study a new class of bivariate distributions called bivariate semi α-Laplace distribution, containing bivariate Laplace distributions. Three characterizations of bivariate semi α-Laplace distribution are obtained. Relation with bivariate semi stable distribution is established. An autoregressive model with bivariate semi α-Laplace marginal distributions is developed.     

Inequalities on partial correlations in Gaussian graphical models containing star shapes

ABSTRACT

This short paper proves inequalities that restrict the magnitudes of the partial correlations in star-shaped structures in Gaussian graphical models. These inequalities have to be satisfied by distributions that are used for generating simulated data to test structure-learning algorithms, but methods that have been used to create such distributions do not always ensure that they are. The inequalities are also noteworthy because stars are common and meaningful in real-world networks.     

A bivariate geometric distribution allowing for positive or negative correlation

In this paper, we propose a new bivariate geometric model, derived by linking two univariate geometric distributions through a specific copula function, allowing for positive and negative correlations. Some properties of this joint distribution are presented and discussed, with particular reference to attainable correlations, conditional distributions, reliability concepts, and parameter estimation. A Monte Carlo simulation study empirically evaluates and compares the performance of the proposed estimators in terms of bias and standard error. Finally, in order to demonstrate its usefulness, the model is applied to a real data set.     

An Improved Boxplot for Univariate Data

The boxplot is an effective data-visualization tool useful in diverse applications and disciplines. Although more sophisticated graphical methods exist, the boxplot remains relevant due to its simplicity, interpretability, and usefulness, even in the age of big data. This article highlights the origins and developments of the boxplot that is now widely viewed as an industry standard as well as its inherent limitations when dealing with data from skewed distributions, particularly when detecting outliers. The proposed Ratio-Skewed boxplot is shown to be practical and suitable for outlier labeling across several parametric distributions.     

CAMCR: Computer-Assisted Mixture model analysis for Capture–Recapture count data

Population size estimation with discrete or nonparametric mixture models is considered, and reliable ways of construction of the nonparametric mixture model estimator are reviewed and set into perspective. Construction of the maximum likelihood estimator of the mixing distribution is done for any number of components up to the global nonparametric maximum likelihood bound using the EM algorithm. In addition, the estimators of Chao and Zelterman are considered with some generalisations of Zelterman’s estimator. All computations are done with CAMCR, a special software developed for population size estimation with mixture models. Several examples and data sets are discussed and the estimators illustrated. Problems using the mixture model-based estimators are highlighted.     

On families of beta- and generalized gamma-generated distributions and associated inference

A general family of univariate distributions generated by beta random variables, proposed by Jones, has been discussed recently in the literature. This family of distributions possesses great flexibility while fitting symmetric as well as skewed models with varying tail weights. In a similar vein, we define here a family of univariate distributions generated by Stacy’s generalized gamma variables. For these two families of univariate distributions, we discuss maximum entropy characterizations under suitable constraints. Based on these characterizations, an expected ratio of quantile densities is proposed for the discrimination of members of these two broad families of distributions. Several special cases of these results are then highlighted. An alternative to the usual method of moments is also proposed for the estimation of the parameters, and the form of these estimators is particularly amenable to these two families of distributions.     

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.     

Comment on testing for spurious and cointegrated regressions: a wavelet approach

In a recent paper, Leong and Huang [6] proposed a wavelet-correlation-based approach to test for cointegration between two time series. However, correlation and cointegration are two different concepts even when wavelet analysis is used. It is known that statistics based on non-stationary integrated variables have non-standard asymptotic distributions. However, wavelet analysis offsets the integrating order of non-stationary series so that traditional asymptotics on stationary variables suffices to ascertain the statistical properties of wavelet-based statistics. Based on this, this note shows that wavelet correlations cannot be used as a test of cointegration.     

Computing the Point-biserial Correlation under Any Underlying Continuous Distribution

The connection between the point-biserial and biserial correlations is well-established when the underlying distribution is bivariate normal. For many other bivariate distributions, the formula that links these two quantities is not straightforward to derive or does not have a closed form. We propose a simple technique that enables researchers to compute one of these correlations when the other is specified. For this, we take advantage of the constancy of their ratio, which can be easily approximated for any distribution. We illustrate the proposed method using several examples and discuss its extension to the ordinal case. We believe that this approach is potentially useful in stochastic simulation..     

A Comparison of Some Methods for the Selection of a Common Eigenvector Model for the Covariance Matrices of Two Groups

Two new nonparametric common principal component model selection procedures based on bootstrap distributions of the vector correlations of all combinations of the eigenvectors from two groups are proposed. The performance of these methods is compared in a simulation study to the two parametric methods previously suggested by Flury in 1988, as well as modified versions of two nonparametric methods proposed by Klingenberg in 1996 and then by Klingenberg and McIntyre in 1998. The proposed bootstrap vector correlation distribution (BVD) method is shown to outperform all of the existing methods in most of the simulated situations considered.     

The effects on the distributions of sample canonical correlations under nonnormality

In this paper, we study the effects of nonnormality on the distributions of sample canonical correlations when the population canonical correlations are simple. In order to achieve the purpose, we derive asymptotic expansion formulas for the distributions of a function of the canonical correlations as well as the individual canonical correlations under nonnormal populations. We particularly discuss the distribution of sample canonical correlations under the class of elliptical population. These expansions are given by using a perturbation method. Simulation results are also given.     

Small sample comparison of six bivariate survival curve estimators 1

We study the performance of six proposed bivariate survival curve estimators on simulated right censored data. The performance of the estimators is compared for data generated by three bivariate models with exponential marginal distributions. The estimators are compared in their ability to estimate correlations and survival functions probabilities. Simulated data results are presented so that the proposed estimators in this relatively new area of analysis can be explicitly compared to the known distribution of the data and the parameters of the underlying model. The results show clear differences in the performance of the estimators.     

Permutation Tests for Linear Models

Several approximate permutation tests have been proposed for tests of partial regression coefficients in a linear model based on sample partial correlations. This paper begins with an explanation and notation for an exact test. It then compares the distributions of the test statistics under the various permutation methods proposed, and shows that the partial correlations under permutation are asymptotically jointly normal with means 0 and variances 1. The method of Freedman & Lane (1983) is found to have asymptotic correlation 1 with the exact test, and the other methods are found to have smaller correlations with this test. Under local alternatives the critical values of all the approximate permutation tests converge to the same constant, so they all have the same asymptotic power. Simulations demonstrate these theoretical results.