Modelling dependence using skew t copulas: Bayesian inference and applications

We construct a copula from the skew t distribution of Sahu et al. ( 2003 ). This copula can capture asymmetric and extreme dependence between variables, and is one of the few copulas that can do so and still be used in high dimensions effectively. However, it is difficult to estimate the copula model by maximum likelihood when the multivariate dimension is high, or when some or all of the marginal distributions are discrete‐valued, or when the parameters in the marginal distributions and copula are estimated jointly. We therefore propose a Bayesian approach that overcomes all these problems. The computations are undertaken using a Markov chain Monte Carlo simulation method which exploits the conditionally Gaussian representation of the skew t distribution. We employ the approach in two contemporary econometric studies. The first is the modelling of regional spot prices in the Australian electricity market. Here, we observe complex non‐Gaussian margins and nonlinear inter‐regional dependence. Accurate characterization of this dependence is important for the study of market integration and risk management purposes. The second is the modelling of ordinal exposure measures for 15 major websites. Dependence between websites is important when measuring the impact of multi‐site advertising campaigns. In both cases the skew t copula substantially outperforms symmetric elliptical copula alternatives, demonstrating that the skew t copula is a powerful modelling tool when coupled with Bayesian inference. Copyright © 2010 John Wiley & Sons, Ltd.     

Simulating non‐normal distributions with specified L‐moments and L‐correlations

This paper derives a procedure for simulating continuous non‐normal distributions with specified L‐moments and L‐correlations in the context of power method polynomials of order three. It is demonstrated that the proposed procedure has computational advantages over the traditional product‐moment procedure in terms of solving for intermediate correlations. Simulation results also demonstrate that the proposed L‐moment‐based procedure is an attractive alternative to the traditional procedure when distributions with more severe departures from normality are considered. Specifically, estimates of L‐skew and L‐kurtosis are superior to the conventional estimates of skew and kurtosis in terms of both relative bias and relative standard error. Further, the L‐correlation also demonstrated to be less biased and more stable than the Pearson correlation. It is also shown how the proposed L‐moment‐based procedure can be extended to the larger class of power method distributions associated with polynomials of order five.     

Measuring Inflation Expectations Using Interval‐Coded Data*

To quantify qualitative survey data, the Carlson–Parkin method assumes normality, a time‐invariant symmetric indifference interval, and long‐run unbiased expectations. These assumptions are unnecessary for interval‐coded data. In April 2004, the Monthly Consumer Confidence Survey in Japan started to ask households about their price expectations a year ahead in seven categories with partially known boundaries. Thus one can identify up to six parameters including an indifference interval each month. This paper compares normal, skew normal (SN), skew exponential power (SEP), and skew t (St) distributions, and finds that an St distribution fits the data well. The results help us to better understand the dynamics of heterogeneous expectations.     

A Simple Estimator of Two‐Dimensional Copulas,with Applications1

Copulas are distributions with uniform marginals. Non‐parametric copula estimates may violate the uniformity condition in finite samples. We look at whether it is possible to obtain valid piecewise linear copula densities by triangulation. The copula property imposes strict constraints on design points, making an equi‐spaced grid a natural starting point. However, the mixed‐integer nature of the problem makes a pure triangulation approach impractical on fine grids. As an alternative, we study the ways of approximating copula densities with triangular functions which guarantees that the estimator is a valid copula density. The family of resulting estimators can be viewed as a non‐parametric MLE of B‐spline coefficients on possibly non‐equally spaced grids under simple linear constraints. As such, it can be easily solved using standard convex optimization tools and allows for a degree of localization. A simulation study shows an attractive performance of the estimator in small samples and compares it with some of the leading alternatives. We demonstrate empirical relevance of our approach using three applications. In the first application, we investigate how the body mass index of children depends on that of parents. In the second application, we construct a bivariate copula underlying the Gibson paradox from macroeconomics. In the third application, we show the benefit of using our approach in testing the null of independence against the alternative of an arbitrary dependence pattern.     

A Unifying Framework for Marginalised Random‐Intercept Models of Correlated Binary Outcomes

We demonstrate that many current approaches for marginal modelling of correlated binary outcomes produce likelihoods that are equivalent to the copula‐based models herein. These general copula models of underlying latent threshold random variables yield likelihood‐based models for marginal fixed effects estimation and interpretation in the analysis of correlated binary data with exchangeable correlation structures. Moreover, we propose a nomenclature and set of model relationships that substantially elucidates the complex area of marginalised random‐intercept models for binary data. A diverse collection of didactic mathematical and numerical examples are given to illustrate concepts. Copyright © 2013 John Wiley & Sons, Ltd.     

Modelling Regime‐Specific Stock Price Volatility*

Single‐state generalized autoregressive conditional heteroscedasticity (GARCH) models identify only one mechanism governing the response of volatility to market shocks, and the conditional higher moments are constant, unless modelled explicitly. So they neither capture state‐dependent behaviour of volatility nor explain why the equity index skew persists into long‐dated options. Markov switching (MS) GARCH models specify several volatility states with endogenous conditional skewness and kurtosis; of these the simplest to estimate is normal mixture (NM) GARCH, which has constant state probabilities. We introduce a state‐dependent leverage effect to NM‐GARCH and thereby explain the observed characteristics of equity index returns and implied volatility skews, without resorting to time‐varying volatility risk premia. An empirical study on European equity indices identifies two‐state asymmetric NM‐GARCH as the best fit of the 15 models considered. During stable markets volatility behaviour is broadly similar across all indices, but the crash probability and the behaviour of returns and volatility during a crash depends on the index. The volatility mean‐reversion and leverage effects during crash markets are quite different from those in the stable regime.     

Default models based on scale mixtures of Marshall-Olkin copulas: properties and applications

We present a unification of the Archimedean and the Lévy-frailty copula model for portfolio default models. The new default model exhibits a copula known as scale mixture of Marshall-Olkin copulas and an investigation of the dependence structure reveals that desirable properties of both original models are combined. This allows for a wider range of dependence patterns, while the analytical tractability is retained. Furthermore, simultaneous defaults and default clustering are incorporated. In addition, a hierarchical extension is presented which allows for a heterogeneous dependence structure. Finally, the model is applied to the pricing of CDO contracts. For this purpose, an efficient Laplace transform inversion approach is developed. Supporting a separation of marginal default probabilities and dependence structure, the model can be calibrated to CDS contracts in a first step. In a second step, the calibration of several parametric families to CDO contracts demonstrates a good fitting quality, which further emphasizes the suitability of the approach.     

Dynamic stochastic copula models: estimation,inference and applications

We propose a new dynamic copula model in which the parameter characterizing dependence follows an autoregressive process. As this model class includes the Gaussian copula with stochastic correlation process, it can be viewed as a generalization of multivariate stochastic volatility models. Despite the complexity of the model, the decoupling of marginals and dependence parameters facilitates estimation. We propose estimation in two steps, where first the parameters of the marginal distributions are estimated, and then those of the copula. Parameters of the latent processes (volatilities and dependence) are estimated using efficient importance sampling. We discuss goodness‐of‐fit tests and ways to forecast the dependence parameter. For two bivariate stock index series, we show that the proposed model outperforms standard competing models. Copyright © 2010 John Wiley & Sons, Ltd.     

Tests for zero inflation in a bivariate zero‐inflated Poisson model

The score test statistics for testing zero inflation and covariance parameter are proposed in the bivariate zero‐inflated Poisson (BZIP) regression model. The Monte Carlo studies show that the score test and likelihood ratio test for testing zero inflation underestimate the nominal significance level, while the score test for covariance parameter keeps the significance level close to the nominal one. To overcome this nominal level underestimation, we propose a bootstrap method of the score test for the testing problem of zero inflation. An empirical example with covariates is provided to illustrate the results. In addition, score test for zero inflation is also proposed in the BZIP model, which allows a flexible dependence structure using copula.     

A social cognitive interpretation of person‐organization fitting: The maintenance and development of professional technical competency

Information technology professionals must continually align their competencies with new technological innovation and changing organizational technological climates. Existing theory and practice is deficient with respect to dynamic models of person‐organization (P‐O) fit. Social cognitive theory proposes a rich and well‐defined theoretical framework for understanding dynamic, interactive behavioral processes, and so is adopted as a basis for developing a model of P‐O fitting focused on the development of professional technical competencies. P‐O fitting is so named because it refers to organizational fit as an ongoing process of adaptation. The article begins by reviewing existing literature on P‐O fit and social cognitive theory to produce a set of propositions and a theoretical model of P‐O fitting as an adaptive, dynamic phenomenon. Applications of the model for both future research and practice are suggested. © 2007 Wiley Periodicals, Inc.     

Exact Skewness–Kurtosis Tests for Multivariate Normality and Goodness‐of‐Fit in Multivariate Regressions with Application to Asset Pricing Models*

We study the problem of testing the error distribution in a multivariate linear regression (MLR) model. The tests are functions of appropriately standardized multivariate least squares residuals whose distribution is invariant to the unknown cross‐equation error covariance matrix. Empirical multivariate skewness and kurtosis criteria are then compared with a simulation‐based estimate of their expected value under the hypothesized distribution. Special cases considered include testing multivariate normal and stable error distributions. In the Gaussian case, finite‐sample versions of the standard multivariate skewness and kurtosis tests are derived. To do this, we exploit simple, double and multi‐stage Monte Carlo test methods. For non‐Gaussian distribution families involving nuisance parameters, confidence sets are derived for the nuisance parameters and the error distribution. The tests are applied to an asset pricing model with observable risk‐free rates, using monthly returns on New York Stock Exchange (NYSE) portfolios over 5‐year subperiods from 1926 to 1995.     

Trivariate non‐Gaussian copulas to analyze the spatial behavior of rat sightings

The brown rat lives with man in a wide variety of environmental contexts and adversely affects public health by transmission of diseases, bites, and allergies. Understanding behavioral and spatial correlation aspects of pest species can contribute to their effective management and control. Rat sightings can be described by spatial coordinates in a particular region of interest defining a spatial point pattern. In this paper, we investigate the spatial structure of rat sightings in the Latina district of Madrid (Spain) and its relation to a number of distance‐based covariates that relate to the proliferation of rats. Given a number of locations, biologically considered as attractor points, the spatial dependence is modeled by distance‐based covariates and angular orientations through copula functions. We build a particular spatial trivariate distribution using univariate margins coming from the covariate information and provide predictive distributions for such distances and angular orientations.     

Adjusting for Non‐Response in Population‐Based Case‐Control Studies

In this paper we discuss the analysis of data from population‐based case‐control studies when there is appreciable non‐response. We develop a class of estimating equations that are relatively easy to implement. For some important special cases, we also provide efficient semi‐parametric maximum‐likelihood methods. We compare the methods in a simulation study based on data from the Women's Cardiovascular Health Study discussed in Arbogast et al. (Estimating incidence rates from population‐based case‐control studies in the presence of non‐respondents, Biometrical Journal 44, 227–239, 2002).     

Combining Significance of Correlated Statistics with Application to Panel Data*

The inverse normal method, which is used to combine P‐values from a series of statistical tests, requires independence of single test statistics in order to obtain asymptotic normality of the joint test statistic. The paper discusses the modification by Hartung (1999, Biometrical Journal, Vol. 41, pp. 849–855) , which is designed to allow for a certain correlation matrix of the transformed P‐values. First, the modified inverse normal method is shown here to be valid with more general correlation matrices. Secondly, a necessary and sufficient condition for (asymptotic) normality is provided, using the copula approach. Thirdly, applications to panels of cross‐correlated time series, stationary as well as integrated, are considered. The behaviour of the modified inverse normal method is quantified by means of Monte Carlo experiments.     

Ultimate 100‐m world records through extreme‐value theory

We use extreme‐value theory to estimate the ultimate world records for the 100‐m running, for both men and women. For this aim we collected the fastest personal best times set between January 1991 and June 2008. Estimators of the extreme‐value index are based on a certain number of upper order statistics. To optimize this number of order statistics we minimize the asymptotic mean‐squared error of the moment estimator. Using the thus obtained estimate for the extreme‐value index, the right endpoint of the speed distribution is estimated. The corresponding time can be interpreted as the estimated ultimate world record: the best possible time that could be run in the near future. We find 9.51 seconds for the 100‐m men and 10.33 seconds for the women.     

Estimating Mann–Whitney‐type Causal Effects

Mann–Whitney‐type causal effects are generally applicable to outcome variables with a natural ordering, have been recommended for clinical trials because of their clinical relevance and interpretability and are particularly useful in analysing an ordinal composite outcome that combines an original primary outcome with death and possibly treatment discontinuation. In this article, we consider robust and efficient estimation of such causal effects in observational studies and clinical trials. For observational studies, we propose and compare several estimators: regression estimators based on an outcome regression (OR) model or a generalised probabilistic index (GPI) model, an inverse probability weighted estimator based on a propensity score model and two doubly robust (DR), locally efficient estimators. One of the DR estimators involves a propensity score model and an OR model, is consistent and asymptotically normal under the union of the two models and attains the semiparametric information bound when both models are correct. The other DR estimator has the same properties with the OR model replaced by a GPI model. For clinical trials, we extend an existing augmented estimator based on a GPI model and propose a new one based on an OR model. The methods are evaluated and compared in simulation experiments and applied to a clinical trial in cardiology and an observational study in obstetrics.     

A Simple Estimator of Two-Dimensional Copulas,with Applications1

Copulas are distributions with uniform marginals. Non-parametric copula estimates may violate the uniformity condition in finite samples. We look at whether it is possible to obtain valid piecewise linear copula densities by triangulation. The copula property imposes strict constraints on design points, making an equi-spaced grid a natural starting point. However, the mixed-integer nature of the problem makes a pure triangulation approach impractical on fine grids. As an alternative, we study the ways of approximating copula densities with triangular functions which guarantees that the estimator is a valid copula density. The family of resulting estimators can be viewed as a non-parametric MLE of B-spline coefficients on possibly non-equally spaced grids under simple linear constraints. As such, it can be easily solved using standard convex optimization tools and allows for a degree of localization. A simulation study shows an attractive performance of the estimator in small samples and compares it with some of the leading alternatives. We demonstrate empirical relevance of our approach using three applications. In the first application, we investigate how the body mass index of children depends on that of parents. In the second application, we construct a bivariate copula underlying the Gibson paradox from macroeconomics. In the third application, we show the benefit of using our approach in testing the null of independence against the alternative of an arbitrary dependence pattern.     

Analysis of failure time using threshold regression with semi‐parametric varying coefficients

Many new statistical models may enjoy better interpretability and numerical stability than traditional models in survival data analysis. Specifically, the threshold regression (TR) technique based on the inverse Gaussian distribution is a useful alternative to the Cox proportional hazards model to analyse lifetime data. In this article we consider a semi‐parametric modelling approach for TR and contribute implementational and theoretical details for model fitting and statistical inferences. Extensive simulations are carried out to examine the finite sample performance of the parametric and non‐parametric estimates. A real example is analysed to illustrate our methods, along with a careful diagnosis of model assumptions.     

Estimation of multivariate models for time series of possibly different lengths

We consider the problem of estimating parametric multivariate density models when unequal amounts of data are available on each variable. We focus in particular on the case that the unknown parameter vector may be partitioned into elements relating only to a marginal distribution and elements relating to the copula. In such a case we propose using a multi‐stage maximum likelihood estimator (MSMLE) based on all available data rather than the usual one‐stage maximum likelihood estimator (1SMLE) based only on the overlapping data. We provide conditions under which the MSMLE is not less asymptotically efficient than the 1SMLE, and we examine the small sample efficiency of the estimators via simulations. The analysis in this paper is motivated by a model of the joint distribution of daily Japanese yen–US dollar and euro–US dollar exchange rates. We find significant evidence of time variation in the conditional copula of these exchange rates, and evidence of greater dependence during extreme events than under the normal distribution. Copyright © 2006 John Wiley & Sons, Ltd.