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.     

Kendall's tau for serial dependence

The authors show how Kendall's tau can be adapted to test against serial dependence in a univariate time series context. They provide formulas for the mean and variance of circular and noncircular versions of this statistic, and they prove its asymptotic normality under the hypothesis of independence. They present also a Monte Carlo study comparing the power and size of a test based on Kendall's tau with the power and size of competing procedures based on alternative parametric and nonparametric measures of serial dependence. In particular, their simulations indicate that Kendall's tau outperforms Spearman's rho in detecting first‐order autoregressive dependence, despite the fact that these two statistics are asymptotically equivalent under the null hypothesis, as well as under local alternatives.     

A simulation study on the hybrid nature of Tango's index

Since the early 1990s, there has been an increasing interest in statistical methods for detecting global spatial clustering in data sets. Tango's index is one of the most widely used spatial statistics for assessing whether spatially distributed disease rates are independent or clustered. Interestingly, this statistic can be partitioned into the sum of two terms: one term is similar to the usual chi-square statistic, being based on deviation patterns between the observed and expected values, and the other term, similar to Moran's I, is able to detect the proximity of similar values. In this paper, we examine this hybrid nature of Tango's index. The goal is to evaluate the possibility of distinguishing the spatial sources of clustering: lack of fit or spatial autocorrelation. To comply with the aims of the work, a simulation study is performed, by which examples of patterns driving the goodness-of-fit and spatial autocorrelation components of the statistic are provided. As for the latter aspect, it is worth noting that inducing spatial association among count data without adding lack of fit is not an easy task. In this respect, the overlapping sums method is adopted. The main findings of the simulation experiment are illustrated and a comparison with a previous research on this topic is also highlighted.     

On the distribution of linear combinations of the components of a Dirichlet random vector

On making use of a result of Imhof, an integral representation of the distribution function of linear combinations of the components of a Dirichlet random vector is obtained. In fact, the distributions of several statistics such as Moran and Geary's indices, the Cliff‐Ord statistic for spatial correlation, the sample coefficient of determination, F‐ratios and the sample autocorrelation coefficient can be similarly determined. Linear combinations of the components of Dirichlet random vectors also turn out to be a key component in a decomposition of quadratic forms in spherically symmetric random vectors. An application involving the sample spectrum associated with series generated by ARMA processes is discussed.     

Tests of serial independence based on Kendall's process

The authors propose new rank statistics for testing the white noise hypothesis in a time series. These statistics are Cramér‐von Mises and Kolmogorov‐Smirnov functionals of an empirical distribution function whose mean is related to a serial version of Kendall's tau through a linear transform. The authors determine the asymptotic behaviour of the underlying serial process and the large‐sample distribution of the proposed statistics under the null hypothesis of white noise. They also present simulation results showing the power of their tests.     

Using balanced iterative reducing and clustering hierarchies to compute approximate rank statistics on massive datasets

The balanced iterative reducing and clustering hierarchies (BIRCH) algorithm handles massive datasets by reading the data file only once, clustering the data as it is read, and retaining only a few clustering features to summarize the data read so far. Using BIRCH allows to analyse datasets that are too large to fit in the computer main memory. We propose estimates of Spearman's ρ and Kendall's τ that are calculated from a BIRCH output and assess their performance through Monte Carlo studies. The numerical results show that the BIRCH-based estimates can achieve the same efficiency as the usual estimates of ρ and τ while using only a fraction of the memory otherwise required.     

On Dependence Properties of Random Minima and Maxima

Let X _(n) and X ₍₁₎ be the largest and smallest order statistics, respectively, of a random sample of fixed size n. Quite generally, X ₍₁₎ and X _(n) are approximately independent for n sufficiently large. In this article, we study the dependence properties of random extremes in terms of their copula, when the sample size has a left-truncated binomial distribution and show that they tend to be more dependent in this case. We also give closed-form formulas for the measures of association Kendall's τ and Spearman's ρ to measure the amount of dependence between two extremes.     

Tests for noncorrelation of two multivariate ARMA time series

In many situations, we want to verify the existence of a relationship between multivariate time series. In this paper, we generalize the procedure developed by Haugh (1976) for univariate time series in order to test the hypothesis of noncorrelation between two multivariate stationary ARMA series. The test statistics are based on residual cross-correlation matrices. Under the null hypothesis of noncorrelation, we show that an arbitrary vector of residual cross-correlations asymptotically follows the same distribution as the corresponding vector of cross-correlations between the two innovation series. From this result, it follows that the test statistics considered are asymptotically distributed as chi-square random variables. Two test procedures are described. The first one is based on the residual cross-correlation matrix at a particular lag, whilst the second one is based on a portmanteau type statistic that generalizes Haugh's statistic. We also discuss how the procedures for testing noncorrelation can be adapted to determine the directions of causality in the sense of Granger (1969) between the two series. An advantage of the proposed procedures is that their application does not require the estimation of a global model for the two series. The finite-sample properties of the statistics introduced were studied by simulation under the null hypothesis. It led to modified statistics whose upper quantiles are much better approximated by those of the corresponding chi-square distribution. Finally, the procedures developed are applied to two different sets of economic data.     

基于LSTAR模型的非线性协整检验

大量的经济理论和实践都表明,宏观经济时间序列经常会出现非平稳和非线性特征,因而在统计分析时,需要进行非线性协整检验。基于逻辑平滑转换自回归(LSTAR)模型将传统的线性协整表述方法拓展为非线性形式,构造实用的检验程序及合适的统计量,利用软件R进行蒙特卡洛模拟给出非线性协整检验统计量的临界值,并通过实际数据分析购买力平价动态系统的非线性协整关系,说明方法的有效性。     

A Cauchy estimator test for autocorrelation

This article presents a new test for serial correlation in an observed stationary time series. Rather than using the traditional portmanteau tests based on the sample autocorrelation function, we propose a test based on the Cauchy estimator of correlation. A goodness-of-fit statistic for fitted autoregressive moving average models is also derived and the asymptotic distribution of this statistic is quantified. The test can be employed using either this asymptotic distribution or by using Monte-Carlo quantiles. The small sample behaviour is studied via simulation and the Monte-Carlo-based test seems to be more precise. The method is demonstrated on monthly asset returns for Facebook, Incorporated.     

On estimating the variance of a U-statistic

We compare jackknifing and bootstrapping as methods for estimating the variance of a U-statistic. The use of these estimates in calculating asymptotic confidence intervals is discussed, and the results of a numerical study involving Kendall's tau are reported. For the special case of this statistic, the bootstrap is the estimate of choice.     

Goodness-of-fit tests for binomial AR(1) processes

The binomial AR(1) model describes a nonlinear process with a first-order autoregressive (AR(1)) structure and a binomial marginal distribution. To develop goodness-of-fit tests for the binomial AR(1) model, we investigate the observed marginal distribution of the binomial AR(1) process, and we tackle its autocorrelation structure. Motivated by the family of power-divergence statistics for handling discrete multivariate data, we derive the asymptotic distribution of certain categorized power-divergence statistics for the case of a binomial AR(1) process. Then we consider Bartlett's formula, which is widely used in time series analysis to provide estimates of the asymptotic covariance between sample autocorrelations, but which is not applicable when the underlying process is nonlinear. Hence, we derive a novel Bartlett-type formula for the asymptotic distribution of the sample autocorrelations of a binomial AR(1) process, which is then applied to develop tests concerning the autocorrelation structure. Simulation studies are carried out to evaluate the size and power of the proposed tests under diverse alternative process models. Several real examples are used to illustrate our methods and findings.     

Effects of a single outlier on arma identification

Fox (1972), Box and Tiao (1975), and Abraham and Box (1979) have proposed methods for detecting outliers in time series whose ARMA form is known (or identified). We show that the existence of a single aberrant observation, innovation, or intervention causes an ARMA model to be misidentified using unadjusted autocorrelation (acf) and partial autocorrelation estimates. The magnitude, location, type of outlier, and in some cases the ARMA's parameters, affect the identification outcome. We use variance inflation, signal-to-noise ratios, and acf critical values to determine an ARMA model's susceptibility to misidentifi-cation. Numerical and simulation examples suggest how to iteratively use the outlier detection methods in practice.     

Design of acceptance sampling plan using a modified EWMA statistic

In this paper, a modified exponentially weighted moving average (EWMA) statistic is proposed. The approximate distribution of the proposed modified EWMA statistic is derived. A variable acceptance sampling plan is designed using the proposed EWMA statistic. The plan parameters of the proposed sampling plan are determined such that the given producer's risk and consumer's risk are satisfied. The efficiency of the proposed plan based on the new EWMA statistic is compared with the existing EWMA plan in terms of the sample size required. The application of the proposed plan is given with the help of an example.     

Two-sample inference procedures under nonproportional hazards

We introduce a new two-sample inference procedure to assess the relative performance of two groups over time. Our model-free method does not assume proportional hazards, making it suitable for scenarios where nonproportional hazards may exist. Our procedure includes a diagnostic tau plot to identify changes in hazard timing and a formal inference procedure. The tau-based measures we develop are clinically meaningful and provide interpretable estimands to summarize the treatment effect over time. Our proposed statistic is a U-statistic and exhibits a martingale structure, allowing us to construct confidence intervals and perform hypothesis testing. Our approach is robust with respect to the censoring distribution. We also demonstrate how our method can be applied for sensitivity analysis in scenarios with missing tail information due to insufficient follow-up. Without censoring, Kendall's tau estimator we propose reduces to the Wilcoxon-Mann–Whitney statistic. We evaluate our method using simulations to compare its performance with the restricted mean survival time and log-rank statistics. We also apply our approach to data from several published oncology clinical trials where nonproportional hazards may exist.     

Gluing Copulas

We present a new way of constructing n-copulas, by scaling and gluing finitely many n-copulas. Gluing for bivariate copulas produces a copula that coincides with the independence copula on some grid of horizontal and vertical sections. Examples illustrate how gluing can be applied to build complicated copulas from simple ones. Finally, we investigate the analytical as well as statistical properties of the copulas obtained by gluing, in particular, the behavior of Spearman's ρ and Kendall's τ.     

ASYMPTOTIC TESTS FOR GROWTH CURVE MODELS WITH AUTOREGRESSIVE ERRORS

This paper presents the limit distribution (as the number of time points increase) for the score vector of a growth curve model assuming both stationary and explosive autoregressive (A.R.) errors. Limit distributions of the score statistic and the likelihood-ratio statistic for testing composite hypotheses about the regression parameters of several growth curves, when the autocorrelation parameters are treated as nuisance parameters, are presented.     

An approximated distribution of the gini's rank associatio coefficient

An approximate distribution is proposed for the Gini's rank association coefficient g which is, like Kendall's and Spearman's rank correlation coefficient, a statistic to test independence between two random variables. The purposed distribution can be simply transformed into a Student's T distribution; so, hypothesis testing is made much easier.     

An exponentially weighted moving average chart controlling false discovery rate

A new control chart is developed by using the exponentially weighted moving average (EWMA) statistics and a multiple testing procedure for controlling false discovery rate. The multiple testing procedure considers not only the current EWMA statistic, but also a given number of previous statistics at the same time. Numerical simulations are accomplished to evaluate the performance of the proposed control chart in terms of the average run length and the conditional expected delay. The results are compared with those of the existing control charts including the X-bar chart, EWMA, and cumulative sum control charts. Case studies with real data-sets are also presented.     

On inference for Kendall's τ within a longitudinal data setting

Kendall's τ is a non-parametric measure of correlation based on ranks and is used in a wide range of research disciplines. Although methods are available for making inference about Kendall's τ, none has been extended to modeling multiple Kendall's τs arising in longitudinal data analysis. Compounding this problem is the pervasive issue of missing data in such study designs. In this article, we develop a novel approach to provide inference about Kendall's τ within a longitudinal study setting under both complete and missing data. The proposed approach is illustrated with simulated data and applied to an HIV prevention study.