Some tests for uniformity of circular distributions powerful against multimodal alternatives

The author introduces new statistics suited for testing uniformity of circular distributions and powerful against multimodal alternatives. One of them has a simple expression in terms of the geometric mean of the sample of chord lengths. The others belong to a family indexed by a continuous parameter. The asymptotic distributions under the null hypothesis are derived. We compare the power of the new tests against Stephens's alternatives with those of Ajne, Watson, and Hermans‐Rasson's tests. Some of the new tests are the most powerful when the alternative has three or four modes. A heuristic justification of this feature is given. An application to the analysis of archaeological data is provided. The Canadian Journal of Statistics 38:80–96; 2010 © 2010 Statistical Society of Canada     

Testing for order among K populations: theory and examples

Testing for stochastic order among K populations is a common and important problem in statistical practice. It arises in the analysis of both planned experiments and observational studies. The authors develop a new nonparametric test for order among K populations that can accommodate any stochastic ordering. The test is based on a maximally selected chi‐bar‐square statistic. The authors find its limiting distribution and use simulations to derive critical values. Three important examples are used to illustrate the applicability of the general method. The authors find that the new tests outperform the existing methods in many practical cases. The Canadian Journal of Statistics 38: 97–115; 2010 © 2009 Statistical Society of Canada     

Variability explained by covariates in linear mixed‐effect models for longitudinal data

Variability explained by covariates or explained variance is a well‐known concept in assessing the importance of covariates for dependent outcomes. In this paper we study R² statistics of explained variance pertinent to longitudinal data under linear mixed‐effect models, where the R² statistics are computed at two different levels to measure, respectively, within‐ and between‐subject variabilities explained by the covariates. By deriving the limits of R² statistics, we find that the interpretation of explained variance for the existing R² statistics is clear only in the case where the covariance matrix of the outcome vector is compound symmetric. Two new R² statistics are proposed to address the effect of time‐dependent covariate means. In the general case where the outcome covariance matrix is not compound symmetric, we introduce the concept of compound symmetry projection and use it to define level‐one and level‐two R² statistics. Numerical results are provided to support the theoretical findings and demonstrate the performance of the R² statistics. The Canadian Journal of Statistics 38: 352–368; 2010 © 2010 Statistical Society of Canada     

Quantile based tests for exponentiality against DMRQ and NBUE alternatives

Quantile-based reliability analysis has received much attention recently. We propose new quantile-based tests for exponentiality against decreasing mean residual quantile function (DMRQ) and new better than used in expectation (NBUE) classes of alternatives. The exact null distribution of the test statistic is derived when the alternative class is DMRQ. The asymptotic properties of both the test statistics are studied. The performance of the proposed tests with other existing tests in the literature is evaluated through simulation study. Finally, we illustrate our test procedure using real data sets.     

A bivariate distribution with Lomax and geometric margins

We develop a stochastic model describing the joint distribution of $(X,N)$, where $N$ has a geometric distribution while $X$ is the sum of $N$ dependent, heavy-tail Pareto components. Models of this form arise in many applications, ranging from hydro-climatology to finance and insurance. We present fundamental properties of this vector, including marginal and conditional distributions, moments, representations, and parameter estimation. We also include an example from finance, illustrating modeling potential of this new bivariate distribution.     

Mildly explosive autoregression with mixing innovations

In this paper, the limit distribution of the least squares estimator for mildly explosive autoregressive models with strong mixing innovations is established, which is shown to be Cauchy as in the iid case. The result is applied to identify the onset and the end of an explosive period of an econometric time series. Simulations and data analysis are also conducted to demonstrate the usefulness of the result.     

Restricted maximum likelihood estimation of joint mean‐covariance models

The class of joint mean‐covariance models uses the modified Cholesky decomposition of the within subject covariance matrix in order to arrive to an unconstrained, statistically meaningful reparameterisation. The new parameterisation of the covariance matrix has two sets of parameters that separately describe the variances and correlations. Thus, with the mean or regression parameters, these models have three sets of distinct parameters. In order to alleviate the problem of inefficient estimation and downward bias in the variance estimates, inherent in the maximum likelihood estimation procedure, the usual REML estimation procedure adjusts for the degrees of freedom lost due to the estimation of the mean parameters. Because of the parameterisation of the joint mean covariance models, it is possible to adapt the usual REML procedure in order to estimate the variance (correlation) parameters by taking into account the degrees of freedom lost by the estimation of both the mean and correlation (variance) parameters. To this end, here we propose adjustments to the estimation procedures based on the modified and adjusted profile likelihoods. The methods are illustrated by an application to a real data set and simulation studies. The Canadian Journal of Statistics 40: 225–242; 2012 © 2012 Statistical Society of Canada     

A generalized Fleming and Harrington's class of tests for interval‐censored data

The class $G^{\rho,\lambda }$ of weighted log‐rank tests proposed by Fleming & Harrington [Fleming & Harrington (1991) Counting Processes and Survival Analysis, Wiley, New York] has been widely used in survival analysis and is nowadays, unquestionably, the established method to compare, nonparametrically, k different survival functions based on right‐censored survival data. This paper extends the $G^{\rho,\lambda }$ class to interval‐censored data. First we introduce a new general class of rank based tests, then we show the analogy to the above proposal of Fleming & Harrington. The asymptotic behaviour of the proposed tests is derived using an observed Fisher information approach and a permutation approach. Aiming to make this family of tests interpretable and useful for practitioners, we explain how to interpret different choices of weights and we apply it to data from a cohort of intravenous drug users at risk for HIV infection. The Canadian Journal of Statistics 40: 501–516; 2012 © 2012 Statistical Society of Canada