Generating random variates from a bicompositional Dirichlet distribution

A composition is a vector of positive components summing to a constant. The sample space of a composition is the simplex, and the sample space of two compositions, a bicomposition, is a Cartesian product of two simplices. We present a way of generating random variates from a bicompositional Dirichlet distribution defined on the Cartesian product of two simplices using the rejection method. We derive a general solution for finding a dominating density function and a rejection constant and also compare this solution to using a uniform dominating density function. Finally, some examples of generated bicompositional random variates, with varying number of components, are presented.     

Testing symmetry based on empirical likelihood

In this paper, we propose a general kth correlation coefficient between the density function and distribution function of a continuous variable as a measure of symmetry and asymmetry. We first propose a root-n moment-based estimator of the kth correlation coefficient and present its asymptotic results. Next, we consider statistical inference of the kth correlation coefficient by using the empirical likelihood (EL) method. The EL statistic is shown to be asymptotically a standard chi-squared distribution. Last, we propose a residual-based estimator of the kth correlation coefficient for a parametric regression model to test whether the density function of the true model error is symmetric or not. We present the asymptotic results of the residual-based kth correlation coefficient estimator and also construct its EL-based confidence intervals. Simulation studies are conducted to examine the performance of the proposed estimators, and we also use our proposed estimators to analyze the air quality dataset.     

Correlation analysis with additive distortion measurement errors

This paper studies the estimation of correlation coefficient between unobserved variables of interest. These unobservable variables are distorted in a additive fashion by an observed confounding variable. Two estimators, a direct-plug-in estimator and a residual-based estimator, are proposed. Their asymptotical results are obtained, and the residual-based estimator is shown asymptotically efficient. Moreover, we suggest an asymptotic normal approximation and an empirical likelihood-based statistic to construct the confidence interval. The empirical likelihood statistic is shown to be asymptotically chi-squared. Simulation studies are conducted to examine the performance of the proposed estimators. These methods are applied to analyse the Boston housing price data for an illustration.     

Two new confidence intervals for the coefficient of variation in a normal distribution

In this article we introduce an approximately unbiased estimator for the population coefficient of variation, τ, in a normal distribution. The accuracy of this estimator is examined by several criteria. Using this estimator and its variance, two approximate confidence intervals for τ are introduced. The performance of the new confidence intervals is compared to those obtained by current methods.     

Stationary bootstrapping for realized covariations of high frequency financial data

This paper studies the stationary bootstrap applicability for realized covariations of high frequency asynchronous financial data. The stationary bootstrap method, which is characterized by a block-bootstrap with random block length, is applied to estimate the integrated covariations. The bootstrap realized covariance, bootstrap realized regression coefficient and bootstrap realized correlation coefficient are proposed, and the validity of the stationary bootstrapping for them is established both for large sample and for finite sample. Consistencies of bootstrap distributions are established, which provide us valid stationary bootstrap confidence intervals. The bootstrap confidence intervals do not require a consistent estimator of a nuisance parameter arising from nonsynchronous unequally spaced sampling while those based on a normal asymptotic theory require a consistent estimator. A Monte-Carlo comparison reveals that the proposed stationary bootstrap confidence intervals have better coverage probabilities than those based on normal approximation.     

Estimation for two-way analysis of variance with correlated errors

We consider some estimation and distribution problems encountered in a two way analysis of variance model with only one observation per cell, errors correlated in one level, and the variances are not necessarily equal. The independence criteria for the row and interaction mean sum of squares and distribution of the maximum likelihood estimator of the correlation coefficient are given.     

Estimation of Several Intraclass Correlation Coefficients

An intraclass correlation coefficient observed in several populations is estimated. The basis is a variance-stabilizing transformation. It is shown that the intraclass correlation coefficient from any elliptical distribution should be transformed in the same way. Four estimators are compared. An estimator where the components in a vector consisting of the transformed intraclass correlation coefficients are estimated separately, an estimator based on a weighted average of these components, a pretest estimator where the equality of the components is tested and then the outcome of the test is used in the estimation procedure, and a James-Stein estimator which shrinks toward the mean.     

CONDITIONAL INFERENCE ABOUT AN EQUICORRELATION COEFFICIENT

SenGupta (1987) proposed a locally most powerful test which is globally (one sided) unbiased, and an estimator of p, the equicorrelation coefficient of a standard symmetric multivariate normal (SSMN) distribution. Here we use the idea in Williams (1984) to illustrate the construction and use of ancillary statistics to make inference about p. The test and confidence intervals based on this construction are conditionally optimal.     

Estimating Intraclass Correlation Coefficient and Identifying Influential Observations Under One-Way Random Effects Model

In this paper we consider the estimation of intraclass correlation coefficient and identification of influential observations under one-way random effects model. We introduce an approach to correct negative estimation values induced by the method of moments estimator, and provide an interval estimation for intraclass correlation coefficient. We present the diagnostic tools to identify influential observations through the uncorrected estimate of intraclass correlation coefficient. A simulation study is conducted to investigate the performance of our procedure for identifying influential observations. We also apply the method on a real data of repeated blood pressure measurements.     

Asymptotic distribution of regression correlation coefficient for Poisson regression model

This study treats an asymptotic distribution for measures of predictive power for generalized linear models (GLMs). We focus on the regression correlation coefficient (RCC) that is one of the measures of predictive power. The RCC, proposed by Zheng and Agresti is a population value and a generalization of the population value for the coefficient of determination. Therefore, the RCC is easy to interpret and familiar. Recently, Takahashi and Kurosawa provided an explicit form of the RCC and proposed a new RCC estimator for a Poisson regression model. They also showed the validity of the new estimator compared with other estimators. This study discusses the new statistical properties of the RCC for the Poisson regression model. Furthermore, we show an asymptotic normality of the RCC estimator.     

An Empirical Likelihood Estimate of the Finite Population Correlation Coefficient

In this article, the problem of the estimation of finite population correlation coefficient is considered using the empirical likelihood method. A new estimator that makes the use of both the known mean and variance of an auxiliary variable is proposed. The percent relative bias and percent relative efficiency of the proposed new estimator with respect to the usual estimator of the correlation coefficient is investigated through extensive simulation study for values of the correlation coefficient from ?0.90 to +0.90. The proposed estimator is found to perform better than the simple correlation coefficient from both the bias and relative efficiency points of views, for the population, considered in the investigation. At the end, the proposed estimator has been extended to complex survey designs. Supplementary materials for this article are available online.     

Estimation of the Intraclass Correlation Coefficient for the Analysis of Covariance Model

For the unbalanced analysis of covariance model with one covariate, a simple formula is given for the intraclass correlation coefficient estimator that results from Henderson's Method 3 estimation of variance components. Example calculations and the corresponding interpretations are given for a study of the correlation of iron content among brothers. The example illustrates the manner in which the estimator depends on the pattern of correlation between the covariate and the variable under investigation.     

Canonical Correlation Analysis Using Small Number of Samples

Canonical correlation assesses the relationship between two groups of variables. Although it has been a useful tool in a wide variety of research areas, it is not well known that weaker canonical correlations require larger sample sizes to be correctly inferred. In this article, we investigate small sample bias in canonical correlation analysis and apply the jackknife bias correction to the estimation of canonical correlations. We use bootstrap samples to obtain a better confidence interval for the jackknife canonical correlation estimator.     

Estimation of regression vectors in linear mixed models with Dirichlet process random effects

The Dirichlet process has been used extensively in Bayesian non parametric modeling, and has proven to be very useful. In particular, mixed models with Dirichlet process random effects have been used in modeling many types of data and can often outperform their normal random effect counterparts. Here we examine the linear mixed model with Dirichlet process random effects from a classical view, and derive the best linear unbiased estimator (BLUE) of the fixed effects. We are also able to calculate the resulting covariance matrix and find that the covariance is directly related to the precision parameter of the Dirichlet process, giving a new interpretation of this parameter. We also characterize the relationship between the BLUE and the ordinary least-squares (OLS) estimator and show how confidence intervals can be approximated.     

Bayes methods in the ecological fallacy cntext:estimation of individual correlation from aggregate data

The ecological fallacy is related to Simpson's paradox (1951) where relationships among group means may be counterintuitive and substantially different from relationships within groups, where the groups are usually geographic entities such as census tracts. We consider the problem of estimating the correlation between two jointly normal random variables where only ecological data (group means) are available. Two empirical Bayes estimators and one fully Bayesian estimator are derived and compared with the usual ecological estimator, which is simply the Pearson correlation coefficient of the group sample means. We simulate the bias and mean squared error performance of these estimators, and also give an example employing a dataset where the individual level data are available for model checking. The results indicate superiority of the empirical Bayes estimators in a variety of practical situations where, though we lack individual level data, other relevant prior information is available.     

On correlation analysis of many‐to‐many observations: an alternative to Pearson's correlation coefficient and its application to an ecotoxicological study

Correlation studies are an important hypothesis‐generating and testing tool, and have a wide range of applications in many scientific fields. In ecological studies in particular, multiple environmental variables are often measured in an attempt to determine relationships between chemical, physical and biological factors. For example, one may wish to know whether and how soil properties correlate with plant physiology. Although correlation coefficients are widely used, their properties and limitations are often imperfectly understood. This is especially the case when one is interested in correlations between, say, trace element content in sediments and in marine organisms, where no one‐to‐one correspondence exists. We show that evaluating Pearson's correlation coefficient for either site‐specific means or composite samples results in biased estimates, and we propose an alternative estimator. We use simulation studies to demonstrate that our estimator generally has a much smaller bias and mean squared error. We further illustrate its use in a case study of the correlation between trace element content in sediments and in mussels in Lyttelton Harbour, New Zealand.     

Estimation of the error distribution in a varying coefficient regression model

This paper deals with the estimation of the error distribution function in a varying coefficient regression model. We propose two estimators and study their asymptotic properties by obtaining uniform stochastic expansions. The first estimator is a residual-based empirical distribution function. We study this estimator when the varying coefficients are estimated by under-smoothed local quadratic smoothers. Our second estimator which exploits the fact that the error distribution has mean zero is a weighted residual-based empirical distribution whose weights are chosen to achieve the mean zero property using empirical likelihood methods. The second estimator improves on the first estimator. Bootstrap confidence bands based on the two estimators are also discussed.     

A method of moments estimator of tail dependence in meta-elliptical models

A meta-elliptical model is a distribution function whose copula is that of an elliptical distribution. The tail dependence function in such a bivariate model has a parametric representation with two parameters: a tail parameter and a correlation parameter. The correlation parameter can be estimated by robust methods based on the whole sample. Using the estimated correlation parameter as plug-in estimator, we then estimate the tail parameter applying a modification of the method of moments approach proposed in the paper by Einmahl et al. (2008). We show that such an estimator is consistent and asymptotically normal. Further, we derive the joint limit distribution of the estimators of the two parameters. We illustrate the small sample behavior of the estimator of the tail parameter by a simulation study and on real data, and we compare its performance to that of the competitive estimators.     

Confidence interval estimation for a linear contrast in intraclass correlation coefficients under unequal family sizes for several populations

Confidence intervals [based on F-distribution and (Z) standard normal distribution] for a linear contrast in intraclass correlation coefficients under unequal family sizes for several populations based on several independent multinormal samples have been proposed. It has been found that the confidence interval based on F-distribution consistently and reliably produced better results in terms of shorter average length of the interval than the confidence interval based on standard normal distribution for various combinations of intraclass correlation coefficient values. The coverage probability of the interval based on F-distribution is competitive with the coverage probability of the interval based on standard normal distribution. The interval based on F-distribution can be used for both small sample and large sample situations. An example with real life data has been presented.     

Nonparametric inference for the joint distribution of recurrent marked variables and recurrent survival time

Time between recurrent medical events may be correlated with the cost incurred at each event. As a result, it may be of interest to describe the relationship between recurrent events and recurrent medical costs by estimating a joint distribution. In this paper, we propose a nonparametric estimator for the joint distribution of recurrent events and recurrent medical costs in right-censored data. We also derive the asymptotic variance of our estimator, a test for equality of recurrent marker distributions, and present simulation studies to demonstrate the performance of our point and variance estimators. Our estimator is shown to perform well for a wide range of levels of correlation, demonstrating that our estimators can be employed in a variety of situations when the correlation structure may be unknown in advance. We apply our methods to hospitalization events and their corresponding costs in the second Multicenter Automatic Defibrillator Implantation Trial (MADIT-II), which was a randomized clinical trial studying the effect of implantable cardioverter-defibrillators in preventing ventricular arrhythmia.