Bounds for expectations of concomitants

We present sharp bounds on expected values of concomitants based on a sample of identically distributed random pairs. The dependence between pair components is described by regression functions or modelled by copulas, or generated by sampling without replacement from finite populations.     

Computation of Sharp Bounds on the Expected Value of a Supermodular Function of Risks with Given Marginals

We show that the rearrangement algorithm (RA) introduced in Puccetti and Rüschendorf (2012 Puccetti, G., Rüschendorf, L. (2012). Computation of sharp bounds on the distribution of a function of dependent risks. Journal of Computational and Applied Mathematics 236(7):1833–1840.[Crossref], [Web of Science ®] , [Google Scholar]) to compute distributional bounds can be used also to compute sharp lower and upper bounds on the expected value of a supermodular function of d random variables having fixed marginal distributions. Compared to the analytical methods existing in the literature the algorithm is widely applicable, more easily obtained and gives insight into the dependence structures attaining the bounds.     

Kernel Estimators for Distribution Functions on Dependent Random Fields

Consider observations (representing lifelengths) taken on a random field indexed by lattice points. Estimating the distribution function F(x) = P(X _i  ≤ x) is an important problem in survival analysis. We propose to estimate F(x) by kernel estimators, which take into account the smoothness of the distribution function. Under some general mixing conditions, our estimators are shown to be asymptotically unbiased and consistent. In addition, the proposed estimator is shown to be strongly consistent and sharp rates of convergence are obtained.     

Bounds on the covariate-time transformation for competing-risks survival analysis

A fundamental problem with the latent-time framework in competing risks is the lack of identifiability of the joint distribution. Given observed covariates along with assumptions as to the form of their effect, then identifiability may obtain. However it is difficult to check any assumptions about form since a more general model may lose identifiability. This paper considers a general framework for modelling the effect of covariates, with the single assumption that the copula dependency structure of the latent times is invariant to the covariates. This framework consists of a set of functions: the covariate-time transformations. The main result produces bounds on these functions, which are derived solely from the crude incidence functions. These bounds are a useful model checking tool when considering the covariate-time transformation resulting from any particular set of further assumptions. An example is given where the widely-used assumption of independent competing risks is checked.     

Bounds on the regression coefficients when a covariate is categorized

This paper considers the problem of a linear regression model in which a continuous covariate is categorized. In such a case one can obtain bounds on the OLS-estimate which would be calculated if the complete information were available. These bounds give an indication of the information loss due to grouping and show what can be inferred from the data without additional assumptions.     

Distribution of linear combination of exponential variates

The exact distribution of a linear combination of n independent negative exponential random variables, when the coefficients of the linear combination are distinct and positive quantities, is well-known. This paper extends the above result to the general case, namely when the coefficients are arbitrary real numbers, positive or negative, distinct or coincident.     

Bounds on Bivariate Distribution Functions with Given Margins and Known Values at Several Points

Let H(x, y) be a continuous bivariate distribution function with known marginal distribution functions F(x) and G(y). Suppose the values of H are given at several points, H(x _i , y _i ) = θ _i , i = 1, 2,…, n. We first discuss conditions for the existence of a distribution satisfying these conditions, and present a procedure for checking if such a distribution exists. We then consider finding lower and upper bounds for such distributions. These bounds may be used to establish bounds on the values of Spearman's ρ and Kendall's τ. For n = 2, we present necessary and sufficient conditions for existence of such a distribution function and derive best-possible upper and lower bounds for H(x, y). As shown by a counter-example, these bounds need not be proper distribution functions, and we find conditions for these bounds to be (proper) distribution functions. We also present some results for the general case, where the values of H(x, y) are known at more than two points. In view of the simplification in notation, our results are presented in terms of copulas, but they may easily be expressed in terms of distribution functions.     

Bounds for Moments of the Maximum of Concomitants of Selected Order Statistics With Application

We present sharp bounds for moments of the maximum of concomitants of selected order statistics. The dependence between pair components is modeled by copulas. We use the bounds to compare some insurance premiums.     

Large Sample Theory for the Bounds on the Gini and Related Indices of Inequality Estimated From Grouped Data

The joint asymptotic distribution of the upper and lower bounds for the Gini index derived by Gastwirth for grouped data are obtained. From them a conservative asymptotically distribution-free confidence interval for the population Gini index is presented. The methods also yield similar results for other indices of inequality (e.g., Theil's and Atkinson's).     

Lower confidence bounds for percentiles of weibull and birnbaum-saunders distributions

Approximate lower confidence bounds on percentiles of the Weibull and the Birnbaum-Saunders distributions are investigated. Asymptotic lower confidence bounds based on Bonferroni's inequality and the Fisher information are discussed, and parametric bootstrap methods to provide better bounds are considered. Since the standard percentile bootstrap method typically does not perform well for confidence bounds on quantiles, several other bootstrap procedures are studied via extensive computer simulations. Results of the simulations indicate that the bootstrap methods generally give sharper lower bounds than the Bonferroni bounds but with coverages still near the nominal confidence level. Two illustrative examples are also presented, one for tensile strength of carbon micro-composite specimens and the other for cycles-to-failure data.     

Weighted Multivariate Tests of Independence

We present new tests of marginal independence for ?^d-valued random vectors. Our tests rely upon weighted Cramér–von Mises-type statistics, which are functionals of the empirical copula process based upon a random sample of size n. We establish a decomposition of this process into asymptotically independent components, and describe the tests which follow from these arguments.     

Empirical Likelihood Confidence Intervals for Distribution Functions under Negatively Associated Samples

In this article, we discuss the construction of the confidence intervals for distribution functions under negatively associated samples. It is shown that the blockwise empirical likelihood (EL) ratio statistic for a distribution function is asymptotically χ²-type distributed. The result is used to obtain an EL-based confidence interval for the distribution function.     

Bounds for the Expected Time to Extinction and the Probability of Extinction in the Galton-Watson Process

Upper bounds for the expected time to extinction in the Galton-Watson process are obtained. We also found upper and lower bounds for the probability of extinction of this process. These bounds improve some bounds previously obtained by other authors.     

On the Multivariate Two-Sample Problem Using Strong Approximations of Empirical Copula Processes

In this article, we establish optimal rates for the strong approximation of empirical copula processes in ?² by sequences of Gaussian processes. These results are applied to investigate Cramér–von Mises-type statistics.     

Approximation Multivariate Distribution with Pair Copula Using the Orthonormal Polynomial and Legendre Multiwavelets Basis Functions

We concentrate on constructing higher dimensional distributions using a fast growing graphical model called Vine/ pair-copula model which has been introduced and developed by Joe, Cooke, Bedford, Kurowica, Daneshkhah, and others. They first construct a n-dimensional copula density by stacking together n(n ? 1)/2 bivariate copula density, and they then approximate arbitrarily well these bivariate copulas and the corresponding multivariate distribution using a semi-parametric method. One constructive approach involves the use of minimum information copulas that can be specified to any required degree of precision based on the available data (or possibly based on the experts’ judgments). By using this method, one is able to use a fixed finite dimensional family of copulas to be employed in terms of a vine construction, with the promise of a uniform level of approximation.

The basic idea behind this method is to use a two-dimensional ordinary polynomial series to approximate any log-density of a bivariate copula function by truncating the series at an appropriate point. We make this approximation method more accurate and computationally faster by using the orthonormal polynomial and Legendre multiwavelets (LMW) series as the basis functions. We show the derived approximations are more precise and computationally faster with better properties than the one proposed previous method in the literature. We then apply our method to modeling a dataset of Norwegian financial data that was previously analyzed in the series of articles, and finally compare our results by them. At the end, we present a method to simulate from the approximated models, and validate our approximation using the simulation results to recover the same dependency structure of the original data.     

Multivariate association test for rare variants controlling for cryptic and family relatedness

In genetic studies of complex diseases, multiple measures of related phenotypes are often collected. Jointly analyzing these phenotypes may improve statistical power to detect sets of rare variants affecting multiple traits. In this work, we consider association testing between a set of rare variants and multiple phenotypes in family‐based designs. We use a mixed linear model to express the correlations among the phenotypes and between related individuals. Given the many sources of correlations in this situation, deriving an appropriate test statistic is not straightforward. We derive a vector of score statistics, whose joint distribution is approximated using a copula. This allows us to have closed‐form expressions for the p‐values of several test statistics. A comprehensive simulation study and an application to Genetic Analysis Workshop 18 data highlight the gains associated with joint testing over univariate approaches, especially in the presence of pleiotropy or highly correlated phenotypes. The Canadian Journal of Statistics 47: 90–107; 2019 © 2018 Statistical Society of Canada     

New Families of Copulas Based on Periodic Functions

Abstract

Although there exists a large variety of copula functions, only a few are practically manageable, and often the choice in dependence modeling falls on the Gaussian copula. Furthermore most copulas are exchangeable, thus implying symmetric dependence. We introduce a way to construct copulas based on periodic functions. We study the two-dimensional case based on one dependence parameter and then provide a way to extend the construction to the n-dimensional framework. We can thus construct families of copulas in dimension n and parameterized by n ? 1 parameters, implying possibly asymmetric relations. Such “periodic” copulas can be simulated easily.     

A Note on Continuous Spatial-Temporal Dynamics of Stochastic Processes

We propose a general form to analyze the space-time interdependency of continuous space-time stochastic processes. We present a new space-time approach based on the intensity function of the underlying point process. These formulations can be, to some extent, analytically solved to obtain explicit formulae of interest. We define a general function that controls the space-time interaction and allows for closed forms depending on the particular choice of several mathematical tools playing a role in this interaction function. In particular, we make use of copulas and Laplace transforms to provide interesting examples of the dynamics of the random intensity function and, in turn, of the number of points contained in a given region.     

混合指数几何分布

本文讨论元件的寿命分布服从双参数混合指数分布,元件个数服从几何分布的情形下,定义了混合指数几何分布,并研究了该分布的各种性质,讨论了参数的极大似然估计,并使用EM算法得到参数的近似估计. 最后,本文给出了参数的渐近方差、协方差、置信区间.