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.     

The Flexible Fourier Form and Local Generalised Least Squares De‐trended Unit Root Tests*

In two recent papers Enders and Lee (2009) and Becker, Enders and Lee (2006) provide Lagrange multiplier and ordinary least squares de‐trended unit root tests, and stationarity tests, respectively, which incorporate a Fourier approximation element in the deterministic component. Such an approach can prove useful in providing robustness against a variety of breaks in the deterministic trend function of unknown form and number. In this article, we generalize the unit root testing procedure based on local generalized least squares (GLS) de‐trending proposed by Elliott, Rothenberg and Stock (1996) to allow for a Fourier approximation to the unknown deterministic component in the same way. We show that the resulting unit root tests possess good finite sample size and power properties and the test statistics have stable non‐standard distributions, despite the curious result that their limiting null distributions exhibit asymptotic rank deficiency.     

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.     

Are output growth‐rate distributions fat‐tailed? some evidence from OECD countries

This work explores some distributional properties of aggregate output growth‐rate time series. We show that, in the majority of OECD countries, output growth‐rate distributions are well approximated by symmetric exponential power densities with tails much fatter than those of a Gaussian (but with finite moments of any order). Fat tails robustly emerge in output growth rates independently of: (i) the way we measure aggregate output; (ii) the family of densities employed in the estimation; (iii) the length of time lags used to compute growth rates. We also show that fat tails still characterize output growth‐rate distributions even after one washes away outliers, autocorrelation and heteroscedasticity. Copyright © 2008 John Wiley & Sons, Ltd.     

Semi‐nonparametric competing risks analysis of recidivism

In this paper we specify a semi‐nonparametric competing risks (SNP‐CR) model of recidivism, for misdemeanors and felonies. The model is a bivariate mixed proportional hazard model with Weibull baseline hazards and common unobserved heterogeneity. The distribution of the latter is modeled semi‐nonparametrically, using orthonormal Legendre polynomials on the unit interval, and integrated out to make the two durations dependent, conditional on the covariates. The SNP‐CR model involved corresponds to a Logit model for felony arrest; hence the validity of the SNP‐CR model can be tested by testing the validity of the implied Logit model. The latter will be done by using the integrated conditional moment (ICM) test. In the first instance we have estimated and tested two versions of the SNP‐CR model, without and with fixed state effects. However, the ICM test rejects these models. Therefore, we have estimated and tested the model for each state separately. These state models are not rejected by the ICM test. Indeed, the estimation results vary substantially per state. Copyright © 2007 John Wiley & Sons, Ltd.     

Fitting Regional Income Distributions in the European Union*

Cross‐sectional distribution of per capita (log) GDP across the European Union regions from 1977 to 1996 is analysed. Kernel density estimates reveal a multimodal structure of the distribution during the 1970s and early 1980s, and a tendency towards unimodality since the mid‐1980s. The distribution is further analysed by a mixture of normal densities. A two well‐separated component mixture fits the distributions in the 1970s and early 1980s. These two clusters tend to converge, supporting the idea of a process of catching up. In the mid‐1990s, a small group of very rich regions is generated by a separated component.     

On Quantile‐based Asymmetric Family of Distributions: Properties and Inference

In this paper, we provide a detailed study of a general family of asymmetric densities. In the general framework, we establish expressions for important characteristics of the distributions and discuss estimation of the parameters via method‐of‐moments as well as maximum likelihood estimation. Asymptotic normality results for the estimators are provided. The results under the general framework are then applied to some specific examples of asymmetric densities. The use of the asymmetric densities is illustrated in a real‐data analysis.     

Partial Horizontal Inequity Orderings: A Non‐parametric Approach*

We provide a partial ordering view of horizontal inequity (HI), based on the Lorenz criterion, associated with different post‐tax income distributions and a (bistochastic) non‐parametric estimated benchmark distribution. As a consequence, several measures consistent with the Lorenz criterion can be rationalized. In addition, we establish the so‐called HI transfer principle, which imposes a normative minimum requirement that any HI measure must satisfy. Our proposed HI ordering is consistent with this principle. Moreover, we adopt a cardinal view to decompose the total effect of a tax system into a welfare gain caused by HI‐free income redistribution and a welfare loss caused by HI, without any additive decomposable restriction on the indices. Hence, more robust tests can be applied. Other decompositions in the literature are seen as particular cases.     

Testing Stationarity in Small‐ and Medium‐Sized Samples when Disturbances are Serially Correlated*

In this article, we study the size distortions of the KPSS test for stationarity when serial correlation is present and samples are small‐ and medium‐sized. It is argued that two distinct sources of the size distortions can be identified. The first source is the finite‐sample distribution of the long‐run variance estimator used in the KPSS test, while the second source of the size distortions is the serial correlation not captured by the long‐run variance estimator because of a too narrow choice of truncation lag parameter. When the relative importance of the two sources is studied, it is found that the size of the KPSS test can be reasonably well controlled if the finite‐sample distribution of the KPSS test statistic, conditional on the time‐series dimension and the truncation lag parameter, is used. Hence, finite‐sample critical values, which can be applied to reduce the size distortions of the KPSS test, are supplied. When the power of the test is studied, it is found that the price paid for the increased size control is a lower raw power against a non‐stationary alternative hypothesis.     

A Simple Panel‐CADF Test for Unit Roots*

In this paper, we propose a simple extension to the panel case of the covariate‐augmented Dickey–Fuller (CADF) test for unit roots developed in Hansen (1995) . The panel test we propose is based on a P values combination approach that takes into account cross‐section dependence. We show that the test has good size properties and gives power gains with respect to other popular panel approaches. An empirical application is carried out for illustration purposes on international data to test the purchasing power parity (PPP) hypothesis.     

Value‐at‐risk for long and short trading positions

In this paper we model Value‐at‐Risk (VaR) for daily asset returns using a collection of parametric univariate and multivariate models of the ARCH class based on the skewed Student distribution. We show that models that rely on a symmetric density distribution for the error term underperform with respect to skewed density models when the left and right tails of the distribution of returns must be modelled. Thus, VaR for traders having both long and short positions is not adequately modelled using usual normal or Student distributions. We suggest using an APARCH model based on the skewed Student distribution (combined with a time‐varying correlation in the multivariate case) to fully take into account the fat left and right tails of the returns distribution. This allows for an adequate modelling of large returns defined on long and short trading positions. The performances of the univariate models are assessed on daily data for three international stock indexes and three US stocks of the Dow Jones index. In a second application, we consider a portfolio of three US stocks and model its long and short VaR using a multivariate skewed Student density. Copyright © 2003 John Wiley & Sons, Ltd.     

APPLYING BETA‐TYPE SIZE DISTRIBUTIONS TO HEALTHCARE COST REGRESSIONS

This paper extends the literature on modelling healthcare cost data by applying the generalised beta of the second kind (GB2) distribution to English hospital inpatient cost data. A quasi‐experimental design, estimating models on a sub‐population of the data and evaluating performance on another sub‐population, is used to compare this distribution with its nested and limiting cases. While for these data the beta of the second kind (B2) distribution and generalised gamma (GG) distribution outperform the GB2, our results illustrate that the GB2 can be used as a device for choosing among competing parametric distributions for healthcare cost data. Copyright © 2013 John Wiley & Sons, Ltd.     

CAUTIONARY NOTES ON THE RANK‐SIZE DISTRIBUTION

In the analysis of center populations at both the national and regional levels considerable use has been made of the rank‐size distribution. Within the extensive literature on this distribution, however, there are a number areas that are unclear, neglected, or misleading. Consideration is initially given to the form of the rank‐size distribution and to the manner in which it is related to other distributions, notably the Pareto and the lognormal. Attention is also focused on the application of the rank‐size distribution to data on center size and on the difficulties that arise in this connection, including the problem of establishing rank‐size regularity. The discussion concludes by taking issue with the contention that the rank‐size distribution is of relevance in the design and formulation of urban policy.     

Bayesian Estimation of Willingness‐to‐pay Where Respondents Mis‐report Their Preferences*

We introduce a modified conditional logit model that takes account of uncertainty associated with mis‐reporting in revealed preference experiments estimating willingness‐to‐pay (WTP). Like Hausman et al. [Journal of Econometrics (1988) Vol. 87, pp. 239–269], our model captures the extent and direction of uncertainty by respondents. Using a Bayesian methodology, we apply our model to a choice modelling (CM) data set examining UK consumer preferences for non‐pesticide food. We compare the results of our model with the Hausman model. WTP estimates are produced for different groups of consumers and we find that modified estimates of WTP, that take account of mis‐reporting, are substantially revised downwards. We find a significant proportion of respondents mis‐reporting in favour of the non‐pesticide option. Finally, with this data set, Bayes factors suggest that our model is preferred to the Hausman model.     

Density Nowcasts and Model Combination: Nowcasting Euro‐Area GDP Growth over the 2008–09 Recession*

Combined density nowcasts for quarterly Euro‐area GDP growth are produced based on the real‐time performance of component models. Components are distinguished by their use of ‘hard’ and ‘soft’, aggregate and disaggregate, indicators. We consider the accuracy of the density nowcasts as within‐quarter indicator data accumulate. We find that the relative utility of ‘soft’ indicators surged during the recession. But as this instability was hard to detect in real‐time it helps, when producing density nowcasts unknowing any within‐quarter ‘hard’ data, to weight the different indicators equally. On receipt of ‘hard’ data for the second month in the quarter better calibrated densities are obtained by giving a higher weight in the combination to ‘hard’ indicators.     

Regression Models with Data‐based Indicator Variables*

Ordinary least squares estimation of an impulse‐indicator coefficient is inconsistent, but its variance can be consistently estimated. Although the ratio of the inconsistent estimator to its standard error has a t‐distribution, that test is inconsistent: one solution is to form an index of indicators. We provide Monte Carlo evidence that including a plethora of indicators need not distort model selection, permitting the use of many dummies in a general‐to‐specific framework. Although White's (1980) heteroskedasticity test is incorrectly sized in that context, we suggest an easy alteration. Finally, a possible modification to impulse ‘intercept corrections’ is considered.     

Multi‐state Stochastic Processes: A Statistical Flowgraph Perspective

Two‐state models (working/failed or alive/dead) are widely used in reliability and survival analysis. In contrast, multi‐state stochastic processes provide a richer framework for modeling and analyzing the progression of a process from an initial to a terminal state, allowing incorporation of more details of the process mechanism. We review multi‐state models, focusing on time‐homogeneous semi‐Markov processes (SMPs), and then describe the statistical flowgraph framework, which comprises analysis methods and algorithms for computing quantities of interest such as the distribution of first passage times to a terminal state. These algorithms algebraically combine integral transforms of the waiting time distributions in each state and invert them to get the required results. The estimated transforms may be based on parametric distributions or on empirical distributions of sample transition data, which may be censored. The methods are illustrated with several applications.     

Dynamics of state price densities

State price densities (SPDs) are an important element in applied quantitative finance. In a Black–Scholes world they are lognormal distributions, but in practice volatility changes and the distribution deviates from log-normality. In order to study the degree of this deviation, we estimate SPDs using EUREX option data on the DAX index via a nonparametric estimator of the second derivative of the (European) call pricing function. The estimator is constrained so as to satisfy no-arbitrage constraints and corrects for the intraday covariance structure in option prices. In contrast to existing methods, we do not use any parametric or smoothness assumptions.     

General‐to‐Specific Model Selection Procedures for Structural Vector Autoregressions*

Structural vector autoregressive (SVAR) models have emerged as a dominant research strategy in empirical macroeconomics, but suffer from the large number of parameters employed and the resulting estimation uncertainty associated with their impulse responses. In this paper, we propose general‐to‐specific (Gets) model selection procedures to overcome these limitations. It is shown that single‐equation procedures are generally efficient for the reduction of recursive SVAR models. The small‐sample properties of the proposed reduction procedure (as implemented using PcGets) are evaluated in a realistic Monte Carlo experiment. The impulse responses generated by the selected SVAR are found to be more precise and accurate than those of the unrestricted VAR. The proposed reduction strategy is then applied to the US monetary system considered by Christiano, Eichenbaum and Evans (Review of Economics and Statistics, Vol. 78, pp. 16–34, 1996) . The results are consistent with the Monte Carlo and question the validity of the impulse responses generated by the full system.