Simulating from Exchangeable Archimedean Copulas

This article considers identifying the existence of the local likelihood estimator in binary regression. A simple method for the identification is proposed, which is derived by recognizing the problem as a linear complementarity problem through a support vector machine problem with soft margins.     

A stochastic representation and sampling algorithm for nested Archimedean copulas

A general sampling algorithm for nested Archimedean copulas was recently suggested. It is given in two different forms, a recursive or an explicit one. The explicit form allows for a simpler version of the algorithm which is numerically more stable and faster since less function evaluations are required. The algorithm can also be given in general form, not being restricted to a particular nesting such as fully nested Archimedean copulas. Further, several examples are given.     

Copula representation of bivariate L-moments: a new estimation method for multiparameter two-dimensional copula models

Serfling and Xiao [A contribution to multivariate L-moments, L-comoment matrices. J Multivariate Anal. 2007;98:1765–1781] extended the L-moment theory to the multivariate setting. In the present paper, we focus on the two-dimensional random vectors to establish a link between the bivariate L-moments (BLM) and the underlying bivariate copula functions. This connection provides a new estimate of dependence parameters of bivariate statistical data. Extensive simulation study is carried out to compare estimators based on the BLM, the maximum likelihood, the minimum distance and a rank approximate Z-estimation. The obtained results show that, when the sample size increases, BLM-based estimation performs better as far as the bias and computation time are concerned. Moreover, the root-mean-squared error is quite reasonable and less sensitive in general to outliers than those of the above cited methods. Further, the proposed BLM method is an easy-to-use tool for the estimation of multiparameter copula models. A generalization of the BLM estimation method to the multivariate case is discussed.     

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.     

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.     

A self-consistent estimator of marginal survival functions based on dependent competing risk data and an assumed copula

When the time to death, X, and the time to censoring, Y, are associated some additional information is need to identify the marginal survival functions. A natural function which provides this additional information is the copula of X and Y. Assuming that the copula is known, we use the notion of self consistency to construct an estimator of the marginal survival functions based on dependent competing risk data. Results of a small simulation study are shown to compare this estimator to other estimators of the marginal survival function based on an assumed copula.     

SOME RESULTS ON TRUNCATION DEPENDENCE INVARIANT CLASS OF COPULAS

ABSTRACT

In this paper, m-dimensional distribution functions with truncation invariant dependence structure are studied. Some of the properties of generalized Archimedean class of copulas under this dependence structure are presented including some results on the conditions of compatibility. It has been shown that Archimedean copula generalized as it is described by Jouini and Clemen[1] Jouini, M.N. and Clemen, R.T. 1996. Copula Models for Aggregating Expert Opinions. Operations Research, 44(3): 444–457.  [Google Scholar] which has the truncation invariant dependence structure has to have the form of independence or Cook-Johnson copula. We also consider a multi-parameter class of copulas derived from one-parameter Archimedean copulas. It has been shown that this class has a probabilistic meaning as a connecting copula of the truncated random pair with a right truncation region on the third variable. Multi-parameter copulas generated in this paper stays in the Archimedean class. We provide formulas to compute Kendall's tau and explore the dependence behavior of this multi-parameter class through examples.     

A CHARACTERIZATION OF GEOMETRIC DISTRIBUTIONS THROUGH CONDITIONAL INDEPENDENCE

Let X₁,…, X_n be mutually independent non-negative integer-valued random variables with probability mass functions f_i(x) > 0 for z= 0,1,…. Let E denote the event that {X₁≥X₂≥…≥X_n}. This note shows that, conditional on the event E, X_i-X_i+ 1 and X_i+ 1 are independent for all t = 1,…, k if and only if X_i (i= 1,…, k) are geometric random variables, where 1 ≤k≤n-1. The k geometric distributions can have different parameters θ_i, i= 1,…, k.     

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.     

Characterization theorems for some discrete distributions based on conditional structure

The joint distribution of (X,Y) is determined if the conditional expectation E {g(X)|Y = y} is given and the conditional distribution of Y|(X = x) is a conditional power series distribution, where g(·) is a function satisfying some minor conditions.     

A characterization based on conditional expectations

We obtain the necessary and sufficient conditions so that any real function (x) is the conditional expectation E(h(X)/X≥x) of a random variable X with continuous distribution function, where h is a given real, continuous and strictly monotonic function.     

Measures of biomarker dependence using a copula-based multivariate epsilon–skew–normal family of distributions

In this note we develop a new multivariate copula model based on epsilon–skew–normal marginal densities for the purpose of examining biomarker dependency structures. We illustrate the flexibility and utility of this model via a variety of graphical tools and a data analysis example pertaining to salivary biomarker. The multivariate normal model is a sub-model of the multivariate epsilon–skew–normal distribution.     

THE REDUCTION OF RISKS FOR SET-VALUED ESTIMATORS

An extension of the Theorem of Rao-Blackwell will be given for the case of set-valued estimators, its images being nonempty, compact-convex subsets of a certain Banach space. Moreover sufficiency will be characterized in connection with set-valued estimators by related point-valued functions.     

Non-random shuffling strategies in blackjack

Numerous playing and betting strategies for the game of twenty-one have been computed assuming the deck or decks are randomly shuffled. In practice, dealers do not spend the time necessary (it takes too long) to completely randomly shuffle the decks used. Hence, there is information not only from the current round of play, but potentially from the previous round of play. We present a model for a non-random shuffle and assert ways in which this information can be used. Rules are derived using a normal approximation which updates the current strategies utilizing information from a non-random shuffle.     

TWO-MOMENT SCALING FORMULAE FOR AGGREGATING EXAMINATION MARKS

Three different methods are presented for deriving linear transformations of raw examination marks in order to construct ostensibly fair’ average or aggregate marks for students taking possibly different subsets of subjects from the curriculum. All methods are based on the same linear model, which in the case of a balanced dataset (that is, if all students have marks in all subjects) is closely related to a simple principal component representation.     

Risk analysis in the brazilian stock market: copula-APARCH modeling for value-at-risk

Risk management of stock portfolios is a fundamental problem for the financial analysis since it indicates the potential losses of an investment at any given time. The objective of this study is to use bivariate static conditional copulas to quantify the dependence structure and to estimate the risk measure Value-at-Risk (VaR). There were selected stocks that have been performing outstandingly on the Brazilian Stock Exchange to compose pairs trading portfolios (B3, Gerdau, Magazine Luiza, and Petrobras). Due to the flexibility that this methodology offers in the construction of multivariate distributions and risk aggregation in finance, we used the copula-APARCH approach with the Normal, T-student, and Joe-Clayton copula functions. In most scenarios, the results showed a pattern of dependence at the extremes. Moreover, the copula form seems not to be relevant for VaR estimation, since in most portfolios the appropriate copulas lead to significant VaR estimates. It has found that the best models fitted provided conservative risk measures, estimates at 5% and 1%, in a scenario more aggressive.     

Optimal Weight Functions for Marginal Proportional Hazards Analysis of Clustered Failure Time Data

The choice of weights in estimating equations for multivariate survival data is considered. Specifically, we consider families of weight functions which are constant on fixed time intervals, including the special case of time-constant weights. For a fixed set of time intervals, the optimal weights are identified as the solution to a system of linear equations. The optimal weights are computed for several scenarios. It is found that for the scenarios examined, the gains in efficiency using the optimal weights are quite small relative to simpler approaches except under extreme dependence, and that a simple estimator of an exchangeable approximation to the weights also performs well.     

On the Characterization of Copulas by Differential Equations

We study the semigroup action induced by univariate conditioning of copulas. Based on this, we give a new characterization of bivariate copulas in terms of flows generated by solutions of ordinary differential equations with not necessary continuous right side. Several applications, related to concordance ordering of copulas, illustrate the usefulness of this result.     

A NONPARAMETRIC TEST OF CHANGING CONDITIONAL VARIANCES IN AUTOREGRESSIVE TIME SERIES

A nonparametric test for detecting changing conditional variances in stationary AR(p) time series is proposed in this paper. For AR(1) models, the test statistic is a Kolmogorov-Smirnov type statistic and the asymptotic theory is developed under both the null and the alternative hypotheses. For AR(p) models (p ≥ 2), an approximate test procedure is proposed. The empirical upper percentage points for our test are tabulated for both p = 1 and p = 2 cases and a bootstrap procedure is suggested for the p ≥ 3 case. Monte Carlo simulations demonstrate that the test has very good powers for finite samples under both normal and non-normal errors.     

AN EXACT CONDITIONAL TEST FOR COVARIANCE SELECTION MODELS

An exact conditional test is developed for testing the absence of an edge in a graphical covariance selection model and is shown to be equivalent to a test based on the partial correlation coefficient. An example is given.