Bayes Estimation of Intraclass Correlation Coefficients Under Unequal Family Sizes

Bayesian inference for the intraclass correlation ρ is considered under unequal family sizes. We obtain the posterior distribution of ρ and then compare the performance of the Bayes estimator (posterior mean of ρ) with that of Srivastava's (1984) estimator through simulation. Simulation study shows that the Bayes estimator performs better than the Srivastava's estimator in terms of lower mean square error. We also obtain large sample posteriors of ρ based on the asymptotic posterior distribution and based on the Laplace approximation.     

Combined estimators for the correlation coefficient

Three combined estimators for the bivariate normal correlation parameter are considered. The data consist of k independent sample correlation coefficients and it is assumed that the underlying correlation parameters are all equal to ρ. Based upon the joint density function of the sample correlations a combined estimator of ρ is obtained as an approximation to the maximum likelihood solution. Two linearly combined estimators are also considered. One of them is based on Fisher's z-transformation of the sample correlations and the other on an unbiased estimator of ρ. The comparison of these three estimators indicates that the combined (approximate) MLE has a slightly smaller estimated mean squared error relative to the other two combined methods of estimation, but it does so at the expense of a relatively larger bias.     

A monte carlo investigation of a statistic for a bivariate missing data problem

Testing the equal means hypothesis of a bivariate normal distribution with homoscedastic varlates when the data are incomplete is considered. If the correlational parameter, ρ, is known, the well-known theory of the general linear model is easily employed to construct the likelihood ratio test for the two sided alternative. A statistic, T, for the case of ρ unknown is proposed by direct analogy to the likelihood ratio statistic when ρ is known. The null and nonnull distribution of T is investigated by Monte Carlo techniques. It is concluded that T may be compared to the conventional t distribution for testing the null hypothesis and that this procedure results in a substantial increase in power-efficiency over the procedure based on the paired t test which ignores the incomplete data. A Monte Carlo comparison to two statistics proposed by Lin and Stivers (1974) suggests that the test based on T is more conservative than either of their statistics.     

Simulated annealing for the bounds of Kendall's τ and Spearman's ρ

It has long been known that, for many joint distributions, Kendall's τ and Spearman's ρ have different values, as they measure different aspects of the dependence structure. Although the classical inequalities between Kendall's τ and Spearman's ρ for pairs of random variables are given, the joint distributions which can attain the bounds between Kendall's τ and Spearman's ρ are difficult to find. We use the simulated annealing method to find the bounds for ρ in terms of τ and its corresponding joint distribution which can attain those bounds. Furthermore, using this same method, we find the improved bounds between τ and ρ, which is different from that given by Durbin and Stuart.     

Tests and confidence intervals for the shape parameter of a gamma distribution

Inference based on the Central Limit Theorem has only first order accuracy. We give tests and confidence intervals (CIs) of second orderaccuracy for the shape parameter ρ of a gamma distribution for both the unscaled and scaled cases.

Tests and CIs based on moment and cumulant estimates are considered as well as those based on the maximum likelihood estimate (MLE).

For the unscaled case the MLE is the moment estimate of order zero; the most efficient moment estimate of integral order is the sample mean, having asymptotic relative efficiency (ARE) .61 when ρ= 1.

For the scaled case the most efficient moment estimate is a functionof the mean and variance. Its ARE is .39 when ρ = 1.

Our motivation for constructing these tests of ρ = 1 and CIs forρ is to provide a simple and convenient method for testing whether a distribution is exponential in situations such as rainfall models where such an assumption is commonly made.     

SB-Robust Estimators of the Parameters of the Wrapped Normal Distribution

In this article, we study the SB-robustness of various estimators of the mean direction (μ) and the concentration parameter (ρ) of the wrapped normal distribution. The functional corresponding to the sample mean direction is seen to be not SB-robust as an estimator of μ at the family of wrapped normal distributions with varying ρ, whereas the γ-trimmed mean direction is SB-robust at the same family of distributions for the different dispersion measures considered in this article. We also study the SB-robustness of the moment estimator of ρ and also that for a newly introduced trimmed estimator of ρ.     

Comparing two poisson parameters: what to do when the optimal isn't done

Suppose two Poisson processes with rates γ₁ and γ₂ are observed for fixed times t_l and t₂. This paper considers hypothesis tests and confidence intervals for the parameter ρ = γ₂/γ₁. Uniformly most powerful unbiased tests and uniformly most accurate unbiased confidence intervals exist for ρ, but they require randomization and so are not used in practice. Several alternative procedures have been proposed. In the context of one-sided hypothesis tests these procedures are reviewed and compared on numerical grounds and by use of the conditionality and repeated sampling principles. It is argued that a conditional binomial test which rejects with conditional level closest to but not necessarily less than, the nominal a is the most reasonable. This test is different from the usual conditional binomial test which rejects with conditional level closeset to but less than or equal to the nominal α Numerical results indicate that an approximate procedure based on the Poisson variance stabilizing transformation has properties similar to the preferred conditional binomial test. Values for λ₁ = t₁λ₁ required to achieve a specified power are considered. These results are also discussed in terms of test inversion to obtain confidence intervals.     

Optimum stratification with two variates

For the 2×2 rectilinear stratification of a bivariate normal distribution with proportional and optimum allocation the dependence of the objective function z(x₁;y₁) on the coefficient of correlation ρ and the sampling fraction q=n/N is investigated. With proportional allocation for great values of ρ (but already for q=0) a so-called ρ-effect arises, which results in a saddle-point of z as “optimum” stratification point in the center of gravity of the distribution and two additional minima. With optimum allocation first for smaller values of q also the ρ-effect arises; for grater values of q a so-called q-effect is superposed, which results in a multitude of minima, saddle-points and maxima of z. All these points satisfy the generalized conditions of Dalenius, but for practical use only the global minimum is of interest.     

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.     

On optimal testing for the equality of equicorrelation: An example of loss in power

Sen Gupta (1988) considered a locally most powerful (LMP) test for testing nonzero values of the equicorrelation coefficient of a standard symmetric multivariate normal distribution. This paper constructs analogous tests for the symmetric multivariate normal distribution. It shows that the new test is uniformly most powerful invariant even in the presence of a nuisance parameter, σ². Further applications of LMP invariant tests to several equicorrelated populations have been considered and an extension to panel data modeling has been suggested.     

Infrence for non-negative autoregressive schemes

The object of this paper is the statistical analysis of Several Closely related models arising in water quality analysis. In particular, concern is with the autoregressive scheme X_r = ρX_r?1 + Y_r where 0 < ρ < 1 and Y's are i.i.d, and non-negative. The estimation and testing problem is considered for three parametric models - Gaussian, uniform and exponential - as well as for the nonparametric case where it is assumed that the Y's have a positive continuous distribution.     

On The Role of a Certain Eigenvalue in Estimating the Growth Rate of a Branching Process

In many situations, the data given on a p-type Galton-Watson process Z_n eP N^p will consist of the total generation sizes |Z_n| only. In that case, the maximum likelihood estimator ρ_ML of the growth rate ρ is not observable, and the asymptotic properties of the most obvious estimators of ρ based on the |Z_n|, as studied by Asmussen & Keiding (1978), show a crucial dependence on |ρ₁|,ρ₁ being a certain other eigenvalue of the offspring mean matrix. In fact, if |ρ₁|²≤ρ, then the speed of convergence compares badly with ρ_ML. In the present note, it is pointed out that recent results of Heyde (1981) on so-called Fibonacci branching processes provide further examples of this phenomenon, and an estimator with the same speed of convergence as ρ_ML and based on the |Z_n| alone is exhibited for the case p= 2, ρ₁²≥ρ.     

ON ASYMPTOTIC OPTIMALITY OF SOME TESTS FOR EXPONENTIAL DISTRIBUTIONS

The probability density function (pdf) of a two parameter exponential distribution is given by f(x; p, s?) =s?^-1 exp {-(x - ρ)/s?} for x≥ρ and 0 elsewhere, where 0 < ρ < ∞ and 0 < s?∞. Suppose we have k independent random samples where the ith sample is drawn from the ith population having the pdf f(x; ρ_i, s?_i), 0 < ρ_i < ∞, 0 < s?_i < s?_i < and f(x; ρ, s?) is as given above. Let X_i1 < X_i2 <… < X_iri denote the first r_i order statistics in a random sample of size n_i, drawn from the ith population with pdf f(x; ρ_i, s?_i), i = 1, 2,…, k. In this paper we show that the well known tests of hypotheses about the parameters ρ_i, s?_i, i = 1, 2,…, k based on the above observations are asymptotically optimal in the sense of Bahadur efficiency. Our results are similar to those for normal distributions.     

Comparative Study of Robust Estimators Based on a Sensitivity Coefficient in Principal Component Analysis

This work is devoted to robust principal component analysis (PCA). We give a comparison between some multivariate estimators of location and scatter by computing the influence functions of the sensitivity coefficient ρ corresponding to these estimators, and the mean squared error (MSE) of estimators of ρ. The coefficient ρ measures the closeness between the subspaces spanned by the initial eigenvectors and their corresponding version derived from an infinitesimal perturbation of the data distribution.     

Locally minimax tests for multiple correlations

Let X be a normally distributed p-dimensional column vector with mean μ and positive definite covariance matrix σ. and let X α, α = 1,…, N, be a random sample of size N from this distribution. Partition X as ( X ₁, X ₍₂₎', X '₍₃₎)', where X₁ is one-dimension, X₍₂₎ is p₂- dimensional, and so 1 + p₁ + p₂ = p. Let ρ₁ and ρ be the multiple correlation coefficients of X₁ with X₍₂₎ and with ( X '₍₂₎, X '₍₃₎)', respectively. Write ρ2/2 = ρ² - ρ2/1. We shall cosider the following two problems     

Some asymptotic properties of point estimation with n-dimensional fuzzy data

Some asymptotic properties of point estimation with n-dimensional fuzzy data with respect to a special L ₂-metric ρ are investigated in this article. It is shown that the collection of all n-dimensional fuzzy data endowed with the ρ-metric is a complete and separable space. Some criterions for point estimation in such fuzzy environments are proposed, and the sample mean, variance and covariance with n-dimensional fuzzy data under these criteria are further studied.     

A statistical test for the hypothesis of Gaussian random function

A Gaussian random function is a functional version of the normal distribution. This paper proposes a statistical hypothesis test to test whether or not a random function is a Gaussian random function. A parameter that is equal to 0 under Gaussian random function is considered, and its unbiased estimator is given. The asymptotic distribution of the estimator is studied, which is used for constructing a test statistic and discussing its asymptotic power. The performance of the proposed test is investigated through several numerical simulations. An illustrative example is also presented.     

THE q-BETA-GEOMETRIC DISTRIBUTION AS A MODEL FOR FECUNDABILITY

The number of sterile couples in a retrospective study of the number of cycles to conception is necessarily zero; this is not so for a prospective study. The paper puts forward a modification of Weinberg and Gladen's beta geometric model for cycles to conception that is suitable for both types of investigation. The probability that a couple achieves conception at the xth cycle, but not earlier, is assumed to take the form R_x = (1 ? ρ)/(1 ? m ^x?1 ρ/u), instead of μ/(1 ? θ + θx). The set of parameter restraints (0 < m < 1, 0< ρ < 1, 1 < u) is appropriate for retrospective data, whilst the alternative set of restraints (1 < m, 1 < ρ, 0 < u < 1) is appropriate for prospective data. The decrease in R_x over time can be interpreted not only as a time effect, but also as a heterogeneity effect by replacing Weinberg and Gladen's beta mixture of geometric distributions by a q-beta mixture.