Multivariate Kendall's tau for change‐point detection in copulas

Statistical procedures for the detection of a change in the dependence structure of a series of multivariate observations are studied in this work. The test statistics that are proposed are $L_1$ , $L_2$ , and $L_{\infty }$ distances computed from vectors of differences of Kendall's tau; two multivariate extensions of Kendall's measure of association are used. Since the distributions of these statistics under the null hypothesis of no change depend on the unknown underlying copula of the vectors, a procedure based on the multiplier central limit theorem is used for the computation of p‐values; the method is shown to be valid both asymptotically and for moderate sample sizes. Alternative versions of the tests that take into account possible breakpoints in the marginal distributions are also investigated. Monte Carlo simulations show that the tests are powerful under many scenarios of change‐point. In addition, two estimators of the time of change are proposed and their efficiency is carefully studied. The methodologies are illustrated on simulated series from the Canadian Regional Climate Model. The Canadian Journal of Statistics 41: 65–82; 2013 © 2012 Statistical Society of Canada     

A Question of Ethics

This article identifies three consulting roles—helper, leader, and colleague; recommends that statistical consultants ask more questions about basic mechanisms; suggests that statistical consultants regard their primary responsibility as providing guidance about the scientific method itself; discusses the need for continuing education to help improve consulting skills; and describes Wisconsin's Master's degree examination, which was designed to help students become effective practicing statisticians.     

Non-Normal Bivariate Distributions with Normal Marginals

The purpose of this note is to indicate that Fieller's Theorem can be expressed in the matrix formulation of the general linear model. The practical consequence is that one general computer program which can estimate the parameters and test the validity of a pertinent model, can also compute confidence limits for the ratios of any linear combinations of the parameters.     

Pmc-optimal nonparametric quantile estimator

According to Pitman's Measure of Closeness, if T₁and T₂are two estimators of a real parameter $[d], then T₁is better than T₂if Po[d]{T₁-o[d] < T₂-0[d]} > 1/2 for all 0[d]. It may however happen that while T₁is better than T₂and T₂is better than T₃, T₃is better than T₁. Given q ? (0,1) and a sample X₁, X₂, ..., X_nfrom an unknown F ? F, an estimator T* = T*(X₁,X₂...X_n)of the q-th quantile of the distribution F is constructed such that P_F{F(T*)-q <[d] F(T)-q} >[d] 1/2 for all F?F and for all T€T, where F is a nonparametric family of distributions and T is a class of estimators. It is shown that T* =X_j:n'for a suitably chosen jth order statistic.     

Robust methods of estimation of regression coefficients 1

Recently, cumulative residual entropy (CRE) has been found to be a new measure of information that parallels Shannon's entropy (see Rao et al. [Cumulative residual entropy: A new measure of information, IEEE Trans. Inform. Theory. 50(6) (2004), pp. 1220–1228] and Asadi and Zohrevand [On the dynamic cumulative residual entropy, J. Stat. Plann. Inference 137 (2007), pp. 1931–1941]). Motivated by this finding, in this paper, we introduce a generalized measure of it, namely cumulative residual Renyi's entropy, and study its properties. We also examine it in relation to some applied problems such as weighted and equilibrium models. Finally, we extend this measure into the bivariate set-up and prove certain characterizing relationships to identify different bivariate lifetime models.     

Importance of the prior mass for agreement between frequentist and Bayesian approaches in the two-sided test

We study distributional properties of generalized order statistics (gos) related by a random shift or scaling scheme in the continuous and discrete case, respectively. In the continuous case, we obtain new characterizations of distributions relating non-neighbouring gos extending some results given in the literature for the neighbouring cases. On the other hand, in the discrete case, we investigate the existence and uniqueness of a discrete parent distribution supported on the integers whose gos are related by a random translation.     

Optimal orthogonal block designs for four-mixture components in two blocks based on F-squares for Becker's models and K-model

Orthogonal block designs in mixture experiments have been extensively studied by various authors. Aggarwal et al. [M.L. Aggarwal, P. Singh, V. Sarin, and B. Husain, Mixture designs in orthogonal blocks using F-squares, METRON – Int. J. Statist. LXVII(2) (2009), pp. 105–128] considered the case of components assuming the same volume fractions and obtained mixture designs in orthogonal blocks using F-squares. In this paper, we have used the class of designs presented by Aggarwal et al. and have obtained D-, A- and E-optimal orthogonal block designs for four components in two blocks for Becker's mixture models and K-model, respectively. Orthogonality conditions for the considered models are also given.     

Adaptive Test for Periodicity in Self-Exciting Threshold Autoregressive Models

This article is devoted to the study of the periodicity testing problem in a self-exciting threshold autoregressive (SETAR) model. The local asymptotic normality (LAN) property is shown via the adapted sufficient conditions due to Swensen (1985 Swensen , A. R. ( 1985 ). The asymptotic distribution of the likelihood ratio for autoregressive time series with a regression trend . Journal of Multivariate Analysis 16 : 54 – 70 .[Crossref], [Web of Science ®] , [Google Scholar]). Moreover, the LAN of the central sequence is established. First, we consider the case where the innovation density is specified and we obtain a parametric local asymptotic test. Second, we construct an adaptive test in the case where this density is unspecified but symmetric. The performances of these established tests are shown via simulation studies.     

Testing Equality of Mean Vectors in Two Sample Problem with Missing Data

Finding the influence of traffic accident on the road is helpful to analyze the characteristics of traffic flow, and take reasonable and effective control measures. Here, the detrended fluctuation analysis method is applied to investigate the complexity of time series in mixed traffic flow with a blockage induced by an accident. As a parameter to depict the long-term evolutionary behavior of the time series in traffic flow, the scaling exponent is analyzed. According to the scaling exponent, it is shown that the traffic flow time series can display long-range correlation characteristics, short-range correlation characteristics, and non-power-law relation in the long-range correlation characteristics, which is strongly dependent on the entering probability of vehicle, the ratio of slow vehicle and the blockage duration time.     

Dirichlet Processes in Nonlinear Mixed Effects Models

In this article, we use two efficient approaches to deal with the difficulty in computing the intractable integrals when implementing Gibbs sampling in the nonlinear mixed effects model (NLMM) based on Dirichlet processes (DP). In the first approach, we compute the Laplace's approximation to the integral for its high accuracy, low cost, and ease of implementation. The second approach uses the no-gaps algorithm of MacEachern and Müller (1998 MacEachern , S. , Müller , P. ( 1998 ). Estimating mixtures of Dirichlet process models . Journal of Computational and Graphical Statistics 7 : 223 – 238 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) to perform Gibbs sampling without evaluating the difficult integral. We apply both approaches to real problems and simulations. Results show that both approaches perform well in density estimation and prediction and are superior to the parametric analysis in that they can detect important model features, such as skewness, long tails, and multimodality, whereas the parametric analysis cannot.