Ordering results on extremes of scaled random variables with dependence and proportional hazards

This study deals with random variables equipped with Archimedean copulas and following scale proportional hazards (SPHs) or revered hazards models. We build the usual stochastic order both between minimums of two SPHs samples with Archimedean survival copulas and between maximums from two scale proportional reversed hazards (PRHs) samples with Archimedean copulas. The hazard rate order between minimums of independent SPHs samples and the reversed hazard rate order between maximums of independent scale PRHs samples are both derived. Also we have a discussion on the dispersive order between minimums from samples with a common Archimedean survival copula. The present results either generalize or improve some related ones in the recent literature.     

Archimedean copulas with applications to $${{\mathrm{}}}$$ estimation

Assuming absolute continuity of marginals, we give the distribution for sums of dependent random variables from some class of Archimedean copulas and the marginal distribution functions of all order statistics. We use conditional independence structure of random variables from this class of Archimedean copulas and Laplace transform. Additionally, we present an application of our results to \({{\mathrm{VaR}}}\) estimation for sums of data from Archimedean copulas.     

Ordering properties of sample minimum from Kumaraswamy-G random variables

In this paper we compare the minimums of two independent and heterogeneous samples each following Kumaraswamy (Kw)-G distribution with the same and the different parent distribution functions. The comparisons are carried out with respect to usual stochastic ordering and hazard rate ordering with majorized shape parameters of the distributions. The likelihood ratio ordering between the minimum order statistics is established for heterogeneous multiple-outlier Kw-G random variables with the same parent distribution function.     

Comparing Survival Times for Treatments with Those for a Control Under Proportional Hazards

Inferences for survival curves based on right censored data are studied for situations in which it is believed that the treatments have survival times at least as large as the control or at least as small as the control. Testing homogeneity with the appropriate order restricted alternative and testing the order restriction as the null hypothesis are considered. Under a proportional hazards model, the ordering on the survival curves corresponds to an ordering on the regression coefficients. Approximate likelihood methods, which are obtained by applying order restricted procedures to the estimates of the regression coefficients, and ordered analogues to the log rank test, which are based on the score statistics, are considered. Mau's (1988) test, which does not require proportional hazards, is extended to this ordering on the survival curves. Using Monte Carlo techniques, the type I error rates are found to be close to the nominal level and the powers of these tests are compared. Other order restrictions on the survival curves are discussed briefly.     

Some ordering properties of series and parallel systems with dependent component lifetimes

Abstract

In this paper, we consider series systems and parallel systems with the dependence between the component lifetimes modelled by an Archimedean copulas. We obtain sufficient and necessary conditions of relative ageing orders between series (parallel) systems with different component numbers, which partially generalize some main results of Misra and Francis. When the component lifetimes follow the scale model, we also characterize the ordering properties between the series systems and (n–1)-out-of-n systems (parallel systems and 2-out-of-n systems) by mixture distribution.     

Do stock returns have an Archimedean copula?

The flexible class of Archimedean copulas plays an important role in multivariate statistics. While there is a large number of goodness-of-fit tests for copulas and parametric families of copulas, the question if a given data set belongs to an arbitrary Archimedean copula or not has not yet received much attention in the literature. This paper suggests a new, straightforward method to test whether a copula is an Archimedean copula without the need to specify its parametric family. We conduct Monte Carlo simulations to assess the power of the test. The approach is applied to (bivariate) joint distributions of stock asset returns. We find that, in general, stock returns may have Archimedean copulas.     

Hazard rate ordering of the second-order statistics from multiple-outlier PHR samples

In this paper, we compare the hazard rate functions of the second-order statistics arising from two sets of independent multiple-outlier proportional hazard rates (PHR) samples. It is proved that the submajorization order between the sample size vectors together with the supermajorization order between the hazard rate vectors imply the hazard rate ordering between the corresponding second-order statistics from multiple-outlier PHR random variables. The results established here provide theoretical guidance both for the winner's price for the bid in the second-price reverse auction in auction theory and fail-safe system design in reliability. Some numerical examples are also provided for illustration.     

Extended tables for moments of gamma distribution order statistics

Apart from having intrinsic mathematical interest, order statistics are also useful in the solution of many applied sampling and analysis problems. For a general review of the properties and uses of order statistics, see David (1981). This paper provides tabulations of means and variances of certain order statistics from the gamma distribution, for parameter values not previously available. The work was motivated by a particular quota sampling problem, for which existing tables are not adequate. The solution to this sampling problem actually requires the moments of the highest order statistic within a given set; however the calculation algorithm used involves a recurrence relation, which causes all the lower order statistics to be calculated first. Therefore we took the opportunity to develop more extensive tables for the gamma order statistic moments in general. Our tables provide values for the order statistic moments which were not available in previous tables, notably those for higher values of m, the gamma distribution shape parameter. However we have also retained the corresponding statistics for lower values of m, first to allow for checking accuracy of the computtions agtainst previous tables, and second to provide an integrated presentation of our new results with the previously known values in a consistent format     

Testing order restricted hypotheses with proportional hazards

Inferences for survival curves based on right censored continuous or grouped data are studied. Testing homogeneity with an ordered restricted alternative and testing the order restriction as the null hypothesis are considered. Under a proportional hazards model, the ordering on the survival curves corresponds to an ordering on the regression coefficients. Approximate likelihood methods are obtained by applying order restricted procedures to the estimates of the regression coefficients. Ordered analogues to the log rank test which are based on the score statistics are considered also. Chi-bar-squared distributions, which have been studied extensively, are shown to provide reasonable approximations to the null distributions of these tests statistics. Using Monte Carlo techniques, the powers of these two types of tests are compared with those that are available in the literature.     

On structure,family and parameter estimation of hierarchical Archimedean copulas

Research on structure determination and parameter estimation of hierarchical Archimedean copulas (HACs) has so far mostly focused on the case in which all appearing Archimedean copulas belong to the same Archimedean family. The present work addresses this issue and proposes a new approach for estimating HACs that involve different Archimedean families. It is based on employing goodness-of-fit test statistics directly into HAC estimation. The approach is summarized in a simple algorithm, its theoretical justification is given and its applicability is illustrated by several experiments, which include estimation of HACs involving up to five different Archimedean families.     

On a new positive dependence concept based on the conditional mean inactivity time order

In this paper, we introduce a new positive dependence concept between two non negative random variables which is related to a conditional version of the mean inactivity time order. A number of properties and relationship between the new notion and the concept of positive likelihood ratio dependence (PLRD) is discussed. Some results in terms of proposed notions for the Archimedean family of copulas are provided.     

Estimating Archimedean Copulas in High Dimensions

Abstract. This article presents a novel estimation procedure for high‐dimensional Archimedean copulas. In contrast to maximum likelihood estimation, the method presented here does not require derivatives of the Archimedean generator. This is computationally advantageous for high‐dimensional Archimedean copulas in which higher‐order derivatives are needed but are often difficult to obtain. Our procedure is based on a parameter‐dependent transformation of the underlying random variables to a one‐dimensional distribution where a minimum‐distance method is applied. We show strong consistency of the resulting minimum‐distance estimators to the case of known margins as well as to the case of unknown margins when pseudo‐observations are used. Moreover, we conduct a simulation comparing the performance of the proposed estimation procedure with the well‐known maximum likelihood approach according to bias and standard deviation.     

Asymptotic properties of numbers of observations in random regions determined by central order statistics

In this paper, we extend the concept of near order statistic observation by considering observations that fall into a random region determined by a given order statistic and a Borel set. We study asymptotic properties of numbers of such observations as the sample size tends to infinity and the order statistic is a central one. We show that then proportions of these numbers converge in probability to some population probabilities. We also prove that these numbers can be centered and normalized to yield normal limit law. First, we derive results for one order statistic; next we give extensions to the multivariate case of two or more order statistics.     

Limit theorems for intermediate and central order statistics under nonlinear normalization

In this paper the work of Pancheva (1984) for extreme order statistics under nonlinear normalization is extended to order statistics with variable ranks. Two new results are proved. The first is that under nonlinear normalization, the nondegenerate type (family of types) of the distribution functions with two finite growth points is a possible weak limit of any central order statistic with regular rank sequence. The second result is that the possible nondegenerate weak limits of any central order statistic with regular rank under the traditionally linear normalization and under the power normalization are the same. Finally, the class of all possible weak limits for lower and upper intermediate order statistics is derived under power normalization from the corresponding weak limits of extremes under power normalization.     

Applications and asymptotic power of marginal-free tests of stochastic vectorial independence

Fully nonparametric tests for the independence between random vectors are studied in this paper. The test statistics are functionals of an empirical process defined as the difference between the joint empirical copula and the product of the empirical copulas associated to the vectors that are suspected to be independent. The validity of a weighted bootstrap procedure is established, which allows for a quick computation of p-values. A special attention is given to the asymptotic behavior of the tests under contiguous sequences of distributions. Finally, a characteristic of the copulas in the Archimedean class in terms of independence of vectors is exploited in order to propose a new goodness-of-fit procedure.     

Sampling nested Archimedean copulas

We give algorithms for sampling from non-exchangeable Archimedean copulas created by the nesting of Archimedean copula generators, where in the most general algorithm the generators may be nested to an arbitrary depth. These algorithms are based on mixture representations of these copulas using Laplace transforms. While in principle the approach applies to all nested Archimedean copulas, in practice the approach is restricted to certain cases where we are able to sample distributions with given Laplace transforms. Precise instructions are given for the case when all generators are taken from the Gumbel parametric family or the Clayton family; the Gumbel case in particular proves very easy to simulate.     

Sequential order statistics from dependent random variables

Abstract

This paper provides an extension for “sequential order statistics” (SOS) introduced by Kamps. It is called “developed sequential order statistics” (DSOS) and is useful for describing lifetimes of engineering systems when component lifetimes are dependent. Explicit expressions for the joint density function, the marginal distributions and the means of DSOS are derived. Under the well known “conditional proportional hazard rate” (CPHR) model and the Gumbel families of copulas for dependency among component lifetimes, some findings are reported. For example, it is proved that the joint density functions of DSOS and SOS have the same structure. Various illustrative examples are also given.     

LTD and RTI dependence orderings

The authors show how the approach of Capéra à & Genest (The Canadian Journal of Statistics, 1990) can be used to order bivariate distributions with arbitrary marginals by their degree of dependence in the LTD (left‐tail decreasing) or RTI (right‐tail increasing) sense. Some properties of these new orderings are given, along with applications to Archimedean copulas, order statistics and compound random variables.     

Ordering extremes of interdependent random variables

For random variables with Archimedean copula or survival copula, we develop the reversed hazard rate order and the hazard rate order on sample extremes in the context of proportional reversed hazard models and proportional hazard models, respectively. The likelihood ratio order on sample maximum is also investigated for the proportional reversed hazard model. Several numerical examples are presented for illustrations as well.     

A nonparametric test to compare survival distributions with covariate adjustment

Summary.  The analysis of covariance is a technique that is used to improve the power of a k -sample test by adjusting for concomitant variables. If the end point is the time of survival, and some observations are right censored, the score statistic from the Cox proportional hazards model is the method that is most commonly used to test the equality of conditional hazard functions. In many situations, however, the proportional hazards model assumptions are not satisfied. Specifically, the relative risk function is not time invariant or represented as a log-linear function of the covariates. We propose an asymptotically valid k -sample test statistic to compare conditional hazard functions which does not require the assumption of proportional hazards, a parametric specification of the relative risk function or randomization of group assignment. Simulation results indicate that the performance of this statistic is satisfactory. The methodology is demonstrated on a data set in prostate cancer.