A Bivariate Generalization of Gamma Distribution

In this article, a bivariate generalisation of the gamma distribution is proposed by using an unsymmetrical bivariate characteristic function; an extension to the non central case also receives attention. The probability density functions of the product and ratio of the correlated components of this distribution are also derived. The benefits of introducing this generalized bivariate gamma distribution and the distributions of the product and the ratio of its components will be demonstrated by graphical representations of their density functions. An example of this generalized bivariate gamma distribution to rainfall data for two specific districts in the North West province is also given to illustrate the greater versatility of the new distribution.     

Depth Measures for Multivariate Functional Data

In this article, we address the problem of mining and analyzing multivariate functional data. That is, data where each observation is a set of possibly correlated functions. Complex data of this kind is more and more common in many research fields, particularly in the biomedical context. In this work, we propose and apply a new concept of depth measure for multivariate functional data. With this new depth measure it is possible to generalize robust statistics, such as the median, to the multivariate functional framework, which in turn allows the application of outlier detection, boxplots construction, and nonparametric tests also in this more general framework. We present an application to Electrocardiographic (ECG) signals.     

Time-varying clustering of multivariate longitudinal observations

Abstract

We propose a statistical method for clustering multivariate longitudinal data into homogeneous groups. This method relies on a time-varying extension of the classical K-means algorithm, where a multivariate vector autoregressive model is additionally assumed for modeling the evolution of clusters' centroids over time. Model inference is based on a least-squares method and on a coordinate descent algorithm. To illustrate our work, we consider a longitudinal dataset on human development. Three variables are modeled, namely life expectancy, education and gross domestic product.     

Tukey-Type Distributions in the Context of Financial Data

Abstract

In one-parameter (θ) families, we were not aware of explicit hypothesis testing scenarios where maximal invariant statistics failed to distinguish the models. We start with a concrete example (Sec. 2.2) to highlight such a hypothesis testing problem involving markedly different models. In this problem, because of the absence of a nontrivial uniformly most powerful invariant (UMPI) test, we briefly suggest two approaches to test the hypothesis. The first resolution (Sec. 3.1) is frequentist in nature. It utilizes a weight function on the parameter space and compares “average” distributions obtained under the null and alternative models in the sense of Wald (1947 Wald , A. ( 1947 ). Sequential Analysis . New York : Wiley . [Google Scholar] 1950 Wald , A. ( 1950 ). Statistical Decision Functions . New York : Wiley . [Google Scholar]). In contrast, a fully Bayesian resolution (Sec. 3.2) is also included. The note ends with a series of other interesting examples involving one-parameter families where maximal invariant statistics fail to distinguish the hypothesized models. The examples include easy-to-construct families of probability models involving only a single location or scale parameter θ.     

Hypothesis Testing Problems in an Unbalanced Longitudinal Ophthalmology Study

ABSTRACT

The analysis of clustered data in a longitudinal ophthalmology study is complicated by correlations between repeatedly measured visual outcomes of paired eyes in a participant and missing observations due to the loss of follow-up. In the present article we consider hypothesis testing problems in an ophthalmology study, where eligible eyes are randomized to two treatments (when two eyes of a participant are eligible, the paired eyes are assigned to different treatments), and vision function outcomes are repeatedly measured over time. A large sample-based nonparametric test statistic and a nonparametric Bootstrap test analog are proposed for testing an interaction effect of two factors and testing an effect of a eye-specific factor within a level of the other person-specific factor on visual function outcomes. Both test statistics allow for missing observations, correlations between repeatedly measured outcomes on individual eyes, and correlations between repeatedly measured outcomes on both eyes of each participant. A simulation study shows that these proposed test statistics maintain nominal significance levels approximately and comparable powers to each other, as well as higher powers than the naive test statistic ignoring correlations between repeated bilateral measurements of both eyes in the same person. For illustration, we apply the proposed test statistics to the changes of visual field defect score in the Advanced Glaucoma Intervention Study.     

Goodness-of-Fit Tests for Multivariate Distributions

We propose three new statistics, Z _p , C _p , and R _p for testing a p-variate (p ≥ 2) normal distribution and compare them with the prominent test statistics. We show that C _p is overall most powerful and is effective against skew, long-tailed as well as short-tailed symmetric alternatives. We show that Z _p and R _p are most powerful against skew and long-tailed alternatives, respectively. The Z _p and R _p statistics can also be used for testing an assumed p-variate nonnormal distribution.     

Estimation of Survivorship Function Based on Doubly Censored Data with Application to HIV/AIDS

Left censoring concept has been defined in different ways in statistical applications. Turnbull (1974 Turnbull , B. W. ( 1974 ). Nonparametric estimation of a survivorship function with doubly censored data . J. Amer. Statist. Assoc. 69 : 169 – 173 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) defines it in a particular way. Whereas in recent literature, especially in epidemiological studies, it has been defined in another way. This difference between the two approaches is the main reason that despite simplicity, Turnbull method cannot be applicable in all cases of doubly censored data. In this article we present a modified Turnbull method for analysis of doubly censored data adequate with recent definition. Comparison has been done with other statistical methods, including imputation estimator, full likelihood-based and conditional likelihood-based approach using Iranian HIV data.     

Strong Gaussian Approximations of Product-Limit and Quantile Processes for Strong Mixing and Censored Data

In this article, we consider the product-limit quantile estimator of an unknown quantile function under a censored dependent model. This is a parallel problem to the estimation of the unknown distribution function by the product-limit estimator under the same model. Simultaneous strong Gaussian approximations of the product-limit process and product-limit quantile process are constructed with rate O[(log n)^?λ] for some λ > 0. The strong Gaussian approximation of the product-limit process is then applied to derive the laws of the iterated logarithm for product-limit process.     

Estimation Bias in Complete-Case Analysis in Crossover Studies with Missing Data

Crossover designs are used often in clinical trials. It is not uncommon that subjects discontinue before completing all treatment periods in a crossover study. Despite availability of statistical methodologies utilizing all available data and software for obtaining valid inferences under the assumption of missing at random (MAR), naïve approaches, such as the complete case (CC) analysis, which is only valid with a strong assumption of missing completely at random are still widely used in practice. In this article, we obtain the analytical form of the estimation bias of treatment effects with CC for linear mixed models. We use simulation studies to examine the inflation of Type I error and efficiency loss in the inferences with CC under MAR. Invalidity and inefficiency of two other commonly used approaches for defining analyzed data in the presence of missing data, including data from at least two periods in three period crossover and available cases for a specific comparison of interest, are also demonstrated through simulation studies.