Nonparametric estimation of the cross ratio function

The cross ratio function (CRF) is a commonly used tool to describe local dependence between two correlated variables. Being a ratio of conditional hazards, the CRF can be rewritten in terms of (first and second derivatives of) the survival copula of these variables. Bernstein estimators for (the derivatives of) this survival copula are used to define a nonparametric estimator of the cross ratio, and asymptotic normality thereof is established. We consider simulations to study the finite sample performance of our estimator for copulas with different types of local dependency. A real dataset is used to investigate the dependence between food expenditure and net income. The estimated CRF reveals that families with a low net income relative to the mean net income will spend less money to buy food compared to families with larger net incomes. This dependence, however, disappears when the net income is large compared to the mean income.

     

Approximation of bivariate copulas by patched bivariate Fréchet copulas

Bivariate Fréchet (BF) copulas characterize dependence as a mixture of three simple structures: comonotonicity, independence and countermonotonicity. They are easily interpretable but have limitations when used as approximations to general dependence structures. To improve the approximation property of the BF copulas and keep the advantage of easy interpretation, we develop a new copula approximation scheme by using BF copulas locally and patching the local pieces together. Error bounds and a probabilistic interpretation of this approximation scheme are developed. The new approximation scheme is compared with several existing copula approximations, including shuffle of min, checkmin, checkerboard and Bernstein approximations and exhibits better performance, especially in characterizing the local dependence. The utility of the new approximation scheme in insurance and finance is illustrated in the computation of the rainbow option prices and stop-loss premiums.     

Nonparametric estimation of an extreme-value copula in arbitrary dimensions

Inference on an extreme-value copula usually proceeds via its Pickands dependence function, which is a convex function on the unit simplex satisfying certain inequality constraints. In the setting of an i.i.d. random sample from a multivariate distribution with known margins and an unknown extreme-value copula, an extension of the Capéraà-Fougères-Genest estimator was introduced by D. Zhang, M. T. Wells and L. Peng [Nonparametric estimation of the dependence function for a multivariate extreme-value distribution, Journal of Multivariate Analysis 99 (4) (2008) 577-588]. The joint asymptotic distribution of the estimator as a random function on the simplex was not provided. Moreover, implementation of the estimator requires the choice of a number of weight functions on the simplex, the issue of their optimal selection being left unresolved.A new, simplified representation of the CFG-estimator combined with standard empirical process theory provides the means to uncover its asymptotic distribution in the space of continuous, real-valued functions on the simplex. Moreover, the ordinary least-squares estimator of the intercept in a certain linear regression model provides an adaptive version of the CFG-estimator whose asymptotic behavior is the same as if the variance-minimizing weight functions were used. As illustrated in a simulation study, the gain in efficiency can be quite sizable.     

Extreme value copula estimation based on block maxima of a multivariate stationary time series

The core of the classical block maxima method consists of fitting an extreme value distribution to a sample of maxima over blocks extracted from an underlying series. In asymptotic theory, it is usually postulated that the block maxima are an independent random sample of an extreme value distribution. In practice however, block sizes are finite, so that the extreme value postulate will only hold approximately. A more accurate asymptotic framework is that of a triangular array of block maxima, the block size depending on the size of the underlying sample in such a way that both the block size and the number of blocks within that sample tend to infinity. The copula of the vector of componentwise maxima in a block is assumed to converge to a limit, which, under mild conditions, is then necessarily an extreme value copula. Under this setting and for absolutely regular stationary sequences, the empirical copula of the sample of vectors of block maxima is shown to be a consistent and asymptotically normal estimator for the limiting extreme value copula. Moreover, the empirical copula serves as a basis for rank-based, nonparametric estimation of the Pickands dependence function of the extreme value copula. The results are illustrated by theoretical examples and a Monte Carlo simulation study.     

Dependence modeling in non-life insurance using the Bernstein copula

This paper illustrates the modeling of dependence structures of non-life insurance risks using the Bernstein copula. We conduct a goodness-of-fit analysis and compare the Bernstein copula with other widely used copulas. Then, we illustrate the use of the Bernstein copula in a value-at-risk and tail-value-at-risk simulation study. For both analyses we utilize German claims data on storm, flood, and water damage insurance for calibration. Our results highlight the advantages of the Bernstein copula, including its flexibility in mapping inhomogeneous dependence structures and its easy use in a simulation context due to its representation as mixture of independent Beta densities. Practitioners and regulators working toward appropriate modeling of dependences in a risk management and solvency context can benefit from our results.     

Empirical likelihood based confidence intervals for copulas

Copula as an effective way of modeling dependence has become more or less a standard tool in risk management, and a wide range of applications of copula models appear in the literature of economics, econometrics, insurance, finance, etc. How to estimate and test a copula plays an important role in practice, and both parametric and nonparametric methods have been studied in the literature. In this paper, we focus on interval estimation and propose an empirical likelihood based confidence interval for a copula. A simulation study and a real data analysis are conducted to compare the finite sample behavior of the proposed empirical likelihood method with the bootstrap method based on either the empirical copula estimator or the kernel smoothing copula estimator.     

Robust estimation of the Pickands dependence function under random right censoring

We consider robust nonparametric estimation of the Pickands dependence function under random right censoring. The estimator is obtained by applying the minimum density power divergence criterion to properly transformed bivariate observations. The asymptotic properties are investigated by making use of results for Kaplan–Meier integrals. We investigate the finite sample properties of the proposed estimator with a simulation experiment and illustrate its practical applicability on a dataset of insurance indemnity losses.     

Nonparametric rank-based tests of bivariate extreme-value dependence

A new class of tests of extreme-value dependence for bivariate copulas is proposed. It is based on the process comparing the empirical copula with a natural nonparametric rank-based estimator of the unknown copula under extreme-value dependence. A multiplier technique is used to compute approximate p-values for several candidate test statistics. Extensive Monte Carlo experiments were carried out to compare the resulting procedures with the tests of extreme-value dependence recently studied in Ben Ghorbal et al. (2009) [1] and Kojadinovic and Yan (2010) [19]. The finite-sample performance study of the tests is complemented by local power calculations.     

Nonparametric estimation of the dependence function for a multivariate extreme value distribution

Understanding and modeling dependence structures for multivariate extreme values are of interest in a number of application areas. One of the well-known approaches is to investigate the Pickands dependence function. In the bivariate setting, there exist several estimators for estimating the Pickands dependence function which assume known marginal distributions [J. Pickands, Multivariate extreme value distributions, Bull. Internat. Statist. Inst., 49 (1981) 859-878; P. Deheuvels, On the limiting behavior of the Pickands estimator for bivariate extreme-value distributions, Statist. Probab. Lett. 12 (1991) 429-439; P. Hall, N. Tajvidi, Distribution and dependence-function estimation for bivariate extreme-value distributions, Bernoulli 6 (2000) 835-844; P. Capéraà, A.-L. Fougères, C. Genest, A nonparametric estimation procedure for bivariate extreme value copulas, Biometrika 84 (1997) 567-577]. In this paper, we generalize the bivariate results to p-variate multivariate extreme value distributions with p?2. We demonstrate that the proposed estimators are consistent and asymptotically normal as well as have excellent small sample behavior.     

An estimator of the stable tail dependence function based on the empirical beta copula

The replacement of indicator functions by integrated beta kernels in the definition of the empirical tail dependence function is shown to produce a smoothed version of the latter estimator with the same asymptotic distribution but superior finite-sample performance. The link of the new estimator with the empirical beta copula enables a simple but effective resampling scheme.     

Bivariate extreme-value copulas with discrete Pickands dependence measure

A two-parametric family of bivariate extreme-value copulas (EVCs), which corresponds to precisely the bivariate EVCs whose Pickands dependence measure is discrete with at most two atoms, is introduced and analyzed. It is shown how bivariate EVCs with arbitrary discrete Pickands dependence measure can be represented as the geometric mean of such basis copulas. General bivariate EVCs can thus be represented as the limit of this construction when the number of involved basis copulas tends to infinity. Besides the theoretical value of such a representation, it is shown how several properties of the represented copula can be deduced from properties of the involved basis copulas. An algorithm for the computation of the representation is given.     

����Copula������Ӧ��-ǿ��ģ�Ϳɿ��ȵķDz�������

In this paper, a nonparametric method for reliability of the stress-strength model is proposed when the dependent stress variable and strength variable are subject to right censoring. The dependence between variables is measured by the common Farlie-Gumbel-Morgenstern copula function and Clayton copula function. Using the empirical process theory, consistency and asymptotic normality of the proposed estimator is established in this paper. The results of numerical simulation show that the proposed method performs well in the case of finite sample. The method proposed in this paper has a wide application prospect in practice.     

Bernstein Polynomial Estimation in the Partially Linear Model under Monotonicity Constraints

In this paper, a Bernstein-polynomial-based likelihood method is proposed for the partially linear model under monotonicity constraints. Monotone Bernstein polynomials are employed to approximate the monotone nonparametric function in the model. The estimator of the regression parameter is shown to be asymptotically normal and efficient, and the rate of convergence of the estimator of the nonparametric component is established, which could be the optimal under the smooth assumptions. A simulation study and a real data analysis are conducted to evaluate the finite sample performance of the proposed method.     

Comparison of three semiparametric methods for estimating dependence parameters in copula models

Three semiparametric methods for estimating dependence parameters in copula models are compared, namely maximum pseudo-likelihood estimation and the two method-of-moment approaches based on the inversion of Spearman’s rho and Kendall’s tau. For each of these three asymptotically normal estimators, an estimator of their asymptotic (co)variance is stated in three different situations, namely the bivariate one-parameter case, the multivariate one-parameter case and the multivariate multiparameter case. An extensive Monte Carlo study is carried out to compare the finite-sample performance of the three estimators under consideration in these three situations. In the one-parameter case, it involves up to six bivariate and four-variate copula families, and up to five levels of dependence. In the multiparameter case, attention is restricted to trivariate and four-variate normal and t copulas. The maximum pseudo-likelihood estimator appears as the best choice in terms of mean square error in all situations except for small and weakly dependent samples. It is followed by the method-of-moment estimator based on Kendall’s tau, which overall appears to be significantly better than its analogue based on Spearman’s rho. The simulation results are complemented by asymptotic relative efficiency calculations. The numerical computation of Spearman’s rho, Kendall’s tau and their derivatives in the case of copula families for which explicit expressions are not available is also investigated.     

Copula parameter estimation by maximum-likelihood and minimum-distance estimators: a simulation study

The purpose of this paper is to present a comprehensive Monte Carlo simulation study on the performance of minimum-distance (MD) and maximum-likelihood (ML) estimators for bivariate parametric copulas. In particular, I consider Cramér-von-Mises-, Kolmogorov-Smirnov- and L ₁-variants of the CvM-statistic based on the empirical copula process, Kendall’s dependence function and Rosenblatt’s probability integral transform. The results presented in this paper show that regardless of the parametric form of the copula, the sample size or the location of the parameter, maximum-likelihood yields smaller estimation biases at less computational effort than any of the MD-estimators. The MD-estimators based on copula goodness-of-fit metrics, on the other hand, suffer from large biases especially when used for estimating the parameters of archimedean copulas. Moreover, the results show that the bias and efficiency of the minimum-distance estimators are strongly influenced by the location of the parameter. Conversely, the results for the maximum-likelihood estimator are relatively stable over the parameter interval of the respective parametric copula.     

Lévy-frailty copulas

A parametric family of n-dimensional extreme-value copulas of Marshall–Olkin type is introduced. Members of this class arise as survival copulas in Lévy-frailty models. The underlying probabilistic construction introduces dependence to initially independent exponential random variables by means of first-passage times of a Lévy subordinator. Jumps of the subordinator correspond to a singular component of the copula. Additionally, a characterization of completely monotone sequences via the introduced family of copulas is derived. An alternative characterization is given by Hausdorff’s moment problem in terms of random variables with compact support. The resulting correspondence between random variables, Lévy subordinators, and copulas is studied and illustrated with several examples. Finally, it is used to provide a general methodology for sampling the copula in many cases. The new class is shown to share some properties with Archimedean copulas regarding construction and analytical form. Finally, the parametric form allows us to compute different measures of dependence and the Pickands representation.     

Estimation of nonstrict Archimedean copulas and its application to quantum networks

Bivariate nonstrict Archimedean copulas form a subclass of Archimedean copulas and are able to model the dependence structure of random variables that do not take on low quantiles simultaneously; i.e. their domain includes a set, the so‐called zero set, with positive Lebesgue measure but zero probability mass. Standard methods to fit a parametric Archimedean copula, e.g. classical maximum likelihood estimation, are either getting computationally more involved or even fail when dealing with this subclass. We propose an alternative method for estimating the parameter of a nonstrict Archimedean copula that is based on the zero set and the functional form of its boundary curve. This estimator is fast to compute and can be applied to absolutely continuous copulas but also allows singular components. In a simulation study, we compare its performance to that of the standard estimators. Finally, the estimator is applied when modeling the dependence structure of quantities describing the quality of transmission in a quantum network, and it is shown how this model can be used effectively to detect potential intruders in this network. Copyright © 2014 John Wiley & Sons, Ltd.     

Applying copula models to individual claim loss reserving methods

The estimation of loss reserves for incurred but not reported (IBNR) claims presents an important task for insurance companies to predict their liabilities. Recently, individual claim loss models have attracted a great deal of interest in the actuarial literature, which overcome some shortcomings of aggregated claim loss models. The dependence of the event times with the delays is a crucial issue for estimating the claim loss reserving. In this article, we propose to use semi-competing risks copula and semi-survival copula models to fit the dependence structure of the event times with delays in the individual claim loss model. A nonstandard two-step procedure is applied to our setting in which the associate parameter and one margin are estimated based on an ad hoc estimator of the other margin. The asymptotic properties of the estimators are established as well. A simulation study is carried out to evaluate the performance of the proposed methods.     

Constructing fixed rank optimal estimators with method of best recurrent approximations

We propose a new approach which generalizes and improves principal component analysis (PCA) and its recent advances. The approach is based on the following underlying ideas. PCA can be reformulated as a technique which provides the best linear estimator of the fixed rank for random vectors. By the proposed method, the vector estimate is presented in a special quadratic form aimed to improve the error of estimation compared with customary linear estimates. The vector is first pre-estimated from the special iterative procedure such that each iterative loop consists of a solution of the unconstrained nonlinear best approximation problem. Then, the final vector estimate is obtained from a solution of the constrained best approximation problem with the quadratic approximant. We show that the combination of these techniques allows us to provide a new nonlinear estimator with a significantly better performance compared with that of PCA and its known modifications.