Multivariate geometric expectiles

A generalization of expectiles for d-dimensional multivariate distribution functions is introduced. The resulting geometric expectiles are unique solutions to a convex risk minimization problem and are given by d-dimensional vectors. They are well behaved under common data transformations and the corresponding sample version is shown to be a consistent estimator. We exemplify their usage as risk measures in a number of multivariate settings, highlighting the influence of varying margins and dependence structures.     

Multivariate extremes for models with constant conditional correlations

Models with constant conditional correlations are versatile tools for describing the behavior of multivariate time series of financial returns. Mathematically speaking, they are solutions of a special class of stochastic recurrence equations (SRE). The extremal behavior of general solutions of SRE has been studied in detail by Kesten [Kesten, H., 1973. Random difference equations and renewal theory for products of random matrices. Acta Mathematica 131, 207–248] and Perfekt [Perfekt, R., 1997. Extreme value theory for a class of Markov chains with values in ^d. Advances in Applied Probability 29, 138–164]. The central concept to understanding the joint extremal behavior of such multivariate time series is the multivariate regular variation spectral measure. In this paper, we propose an estimator for the spectral measure associated with solutions of SRE and prove its consistency. Our estimator is the tail empirical measure of the multivariate time series. Successful use of the estimator depends on a good choice of k, the number of upper order statistics contributing to the empirical measure. We introduce a new criteria for the choice of k based on a scaling property of the spectral measure. We investigate the performance of our estimation technique on exchange rate time series from HFDF96 data set. The estimated spectral measure is used to calculate probabilities of joint extreme returns and probabilities of large movements in an exchange rate conditional on the occurrence of extreme returns in another exchange rate. We find a high level of dependence between the extreme movements of most of the currencies in the EU. We also investigate the changes in the level of dependence between the extreme returns of pairs of currencies as the sampling frequency decreases. When at least one return is extreme, a strong dependence between the components is present already at the 4-hour level for most of the European currencies.     

Prediction Error of the Multivariate Chain Ladder Reserving Method

Abstract

In this paper we consider the claims reserving problem in a multivariate context: that is, we study the multivariate chain-ladder (CL) method for a portfolio of N correlated runoff triangles based on multivariate age-to-age factors. This method allows for a simultaneous study of individual runoff subportfolios and facilitates the derivation of an estimator for the mean square error of prediction (MSEP) for the CL predictor of the ultimate claim of the total portfolio. However, unlike the already existing approaches we replace the univariate CL predictors with multivariate ones. These multivariate CL predictors reflect the correlation structure between the subportfolios and are optimal in terms of a classical optimality criterion, which leads to an improvement of the estimator for the MSEP. Moreover, all formulas are easy to implement on a spreadsheet because they are in matrix notation. We illustrate the results by means of an example.     

A semi-parametric approach to risk management

Abstract

The benchmark theory of mathematical finance is the Black–Scholes–Merton (BSM) theory, based on Brownian motion as the driving noise process for stock prices. Here the distributions of financial returns of the stocks in a portfolio are multivariate normal. Risk management based on BSM underestimates tails. Hence estimation of tail behaviour is often based on extreme value theory (EVT). Here we discuss a semi-parametric replacement for the multivariate normal involving normal variance–mean mixtures. This allows a more accurate modelling of tails, together with various degrees of tail dependence, while (unlike EVT) the whole return distribution can be modelled. We use a parametric component, incorporating the mean vector μ and covariance matrix Σ, and a non-parametric component, which we can think of as a density on [0,∞), modelling the shape (in particular the tail decay) of the distribution. We work mainly within the family of elliptically contoured distributions, focusing particularly on normal variance mixtures with self-decomposable mixing distributions. We discuss efficient methods to estimate the parametric and non-parametric components of our model and provide an algorithm for simulating from such a model. We fit our model to several financial data series. Finally, we calculate value at risk (VaR) quantities for several portfolios and compare these VaRs to those obtained from simple multivariate normal and parametric mixture models.     

On optimal stratifications for multivariate distributions

Summary

In the present paper we study the problem of optimal stratifications for estimating the mean vector y of a given multivariate distribution F(x) with covariance matrix ζ both in cases of proportionate and of optimal (or generalized Neyman) allocations. It is noted that an “optimal stratification” is meant for one to make the covariance matrix of an unbiased estimator X for μ minimal, in the sense of semi-order defined below, in the symmetric matrix space. We show the existence of an optimal stratification and the necessary conditions for a stratification to be optimal. Besides we prove that an optimal stratification can be represented by a “hyperplane stratification” or a “quadratic hypersurface stratification” according to the proportionate or optimal (or generalized Neyman) allocation, and that the set of all optimal (or admissible) stratifications is a minimal complete class in the analogous sense of decision theory. Further we discuss the optimal stratification when a criterion based on a suitable real-valued function is adopted instead of the semi-order.     

Alternative constructions of Tobin's q: An empirical comparison

Although Tobin's q is an attractive theoretical firm performance measure, its empirical construction is subject to considerable measurement error. In this paper we compare five estimators of q that range from a simple-to-construct estimator based on book-values to a relatively complex estimator based upon the methodology developed by Lindenberg and Ross (1981). We present comparisons of the means, medians and variances of the q estimates, and examine how robust sorting and regression results are to changes in the construction of q. We find that empirical results are sensitive to the method used to estimate Tobin's q. The simple-to-construct estimator produces empirical results that differ significantly from the alternative estimators. Among the other four estimators, one developed by Hall (1990) produces means that are higher and variances that are larger than the three alternative estimators, but does approximate those estimators in most of the empirical comparisons. Those three alternative q ratio estimators, furthermore, produce empirical results that are robust.     

A global consistency result for the two-dimensional Pareto distribution in the presence of misspecified inflation

A global consistency result for the ML estimator of a misspecified two-parameter Pareto distribution is proved. The misspecification is due to the assumption of a wrong inflation rate, which violates the i.i.d. assumption in the model. We also investigate how far away from the true parameters one finds the ML estimator of the misspecified model (asymptotically for a small misspecification r). Finally, for the case where the misspecification depends on the number of observations n, i.e., r=r _n , and where $r_{n}\stackrel{n\to \infty}{\longrightarrow}0A global consistency result for the ML estimator of a misspecified two-parameter Pareto distribution is proved. The misspecification is due to the assumption of a wrong inflation rate, which violates the i.i.d. assumption in the model. We also investigate how far away from the true parameters one finds the ML estimator of the misspecified model (asymptotically for a small misspecification r). Finally, for the case where the misspecification depends on the number of observations n, i.e., r=r _n , and where r_n? ^{n? ￥}0r_{n}\stackrel{n\to \infty}{\longrightarrow}0, we prove a central limit theorem for the ML estimator.     

Large Sample Behavior of the CTE and VaR Estimators under Importance Sampling

Abstract

The α-level value at risk (Var) and the α-level conditional tail expectation (CTE) of a continuous random variable X are defined as its α-level quantile (denoted by q_α ) and its conditional expectation given the event {X > q_α }, respectively. Var is a popular risk measure in the banking sector, for both external and internal reporting purposes, while the CTE has recently become the risk measure of choice for insurance regulation in North America. Estimation of the CTE for company assets and liabilities is becoming an important actuarial exercise, and the size and complexity of these liabilities make inference procedures with good small sample performance very desirable. A common situation is one in which the CTE of the portfolio loss is estimated using simulated values, and in such situations use of variance reduction techniques such as importance sampling have proved to be fruitful. Construction of confidence intervals for the CTE relies on the availability of the asymptotic distribution of the normalized CTE estimator, and although such a result has been available to actuaries, it has so far been supported only by heuristics. The main goal of this paper is to provide an honest theorem establishing the convergence of the normalized CTE estimator under importance sampling to a normal distribution. In the process, we also provide a similar result for the Var estimator under importance sampling, which improves upon an earlier result. Also, through examples we motivate the practical need for such theoretical results and include simulation studies to lend insight into the sample sizes at which these asymptotic results become meaningful.     

Asymmetries in Price‐Setting Behavior: New Microeconometric Evidence from Switzerland

In this paper, we follow the recent empirical literature that has specified reduced‐form models for price setting that are closely tied to (S, s) ‐pricing rules. Our contribution to the literature is twofold. First, we propose an estimator that relaxes distributional assumptions on the unobserved heterogeneity. Second, we use the estimator to examine the prevalence of positive price changes in a low‐inflation environment. Our model estimates suggest that, if inflation falls from 0.9% to zero, the share of positive price changes in all price changes falls from 63.6% to 56.2%.     

Model misspecification,Bayesian versus credibility estimation,and Gibbs posteriors

ABSTRACT

In the context of predicting future claims, a fully Bayesian analysis – one that specifies a statistical model, prior distribution, and updates using Bayes's formula – is often viewed as the gold-standard, while Bühlmann's credibility estimator serves as a simple approximation. But those desirable properties that give the Bayesian solution its elevated status depend critically on the posited model being correctly specified. Here we investigate the asymptotic behavior of Bayesian posterior distributions under a misspecified model, and our conclusion is that misspecification bias generally has damaging effects that can lead to inaccurate inference and prediction. The credibility estimator, on the other hand, is not sensitive at all to model misspecification, giving it an advantage over the Bayesian solution in those practically relevant cases where the model is uncertain. This begs the question: does robustness to model misspecification require that we abandon uncertainty quantification based on a posterior distribution? Our answer to this question is No, and we offer an alternative Gibbs posterior construction. Furthermore, we argue that this Gibbs perspective provides a new characterization of Bühlmann's credibility estimator.     

On optimal portfolio diversification with respect to extreme risks

Extreme losses of portfolios with heavy-tailed components are studied in the framework of multivariate regular variation. Asymptotic distributions of extreme portfolio losses are characterized by a functional γ _ξ =γ _ξ (α,Ψ) of the tail index α, the spectral measure Ψ, and the vector ξ of portfolio weights. Existence, uniqueness, and location of the optimal portfolio are analysed and applied to the minimization of risk measures. It is shown that diversification effects are positive for α>1 and negative for α<1. Strong consistency and asymptotic normality are established for a semiparametric estimator of the mapping ξ ↦ γ _ξ . Strong consistency is also established for the estimated optimal portfolio.     

Multivariate Shrinkage for Optimal Portfolio Weights

Abstract

This paper proposes a multivariate shrinkage estimator for the optimal portfolio weights. The estimated classical Markowitz weights are shrunk to the deterministic target portfolio weights. Assuming log asset returns to be i.i.d. Gaussian, explicit solutions are derived for the optimal shrinkage factors. The properties of the estimated shrinkage weights are investigated both analytically and using Monte Carlo simulations. The empirical study compares the competing portfolio selection approaches. Both simulation and empirical studies show that the proposed shrinkage estimator is robust and provides significant gains to the investor compared to benchmark procedures.     

Analytical studies in stop-loss reinsurance. II

Abstract

I

In an earlier paper [5] we discussed the problem of finding an unbiased estimator of where p (x, 0) is a given frequency density and 0 is a (set of) parameter(s). In general, will not be an unbiased estimator of (1), when Ô is an unbiased estimate of O. In [5] it was shown that is an unbiased estimator of (1), if we define y_i , as the larger of 0 and X _j - c. It was emphasized that the resulting estimate may very well be zero, even when it is unreasonable to assume that the premium for a stop.loss reinsurance. defined by a frequency p (x, 0) of claims x and a critical limit c, should be zero when the critical limit has not been exceeded during the n years considered for the determination of the premium.     

Value at risk linear exponent (VARLINEX) forecasts

Abstract

The VARLINEX (value at risk linear exponent) forecasting procedure is presented in this paper, which explicitly adjusts the forecasts when the loss functions of the forecaster are asymmetric. The theory of order statistics is applied to derive the VARLINEX forecasts and their corresponding confidence intervals, which are distribution-free. An empirical study based on our method is carried out for the S&P 500 returns and compared with the RiskMetrics? and the EVT method. It is found that our method can perform very well in relation to EVT and always performs much better than RiskMetrics?.     

Nonparametric methods with applications to hedonic models

Current real estate statistical valuation involves the estimation of parameters within a posited specification. Suchparametric estimation requires judgment concerning model (1) variables; and (2) functional form. In contrast,nonparametric regression estimation requires attention to (1) but permits greatly reduced attention to (2). Parametric estimators functionally model the parameters and variables affectingE(y¦x) while nonparametric estimators directly modelpdf(y, x) and henceE(y¦x).This article applies the kernel nonparametric regression estimator to two different data sets and specifications. The article shows the nonparametric estimator outperforms the standard parametric estimator (OLS) across variable transformations and across data subsets differing in quality. In addition, the article reviews properties of nonparametric estimators, presents the history of nonparametric estimators in real estate, and discusses a representation of the kernel estimator as a nonparametric grid method.     

A maximum likelihood approach to volatility estimation for a Brownian motion using high,low and close price data

Abstract

Volatility plays an important role in derivatives pricing, asset allocation, and risk management, to name but a few areas. It is therefore crucial to make the utmost use of the scant information typically available in short time windows when estimating the volatility. We propose a volatility estimator using the high and the low information in addition to the close price, all of which are typically available to investors. The proposed estimator is based on a maximum likelihood approach. We present explicit formulae for the likelihood of the drift and volatility parameters when the underlying asset is assumed to follow a Brownian motion with constant drift and volatility. Our approach is to then maximize this likelihood to obtain the estimator of the volatility. While we present the method in the context of a Brownian motion, the general methodology is applicable whenever one can obtain the likelihood of the volatility parameter given the high, low and close information. We present simulations which indicate that our estimator achieves consistently better performance than existing estimators (that use the same information and assumptions) for simulated data. In addition, our simulations using real price data demonstrate that our method produces more stable estimates. We also consider the effects of quantized prices and discretized time.     

Measuring price discovery: The variance ratio, the R, and the weighted price contribution

We analyze the statistical properties of three price discovery measures: The variance ratio, the weighted price contribution (WPC), and the R² of unbiasedness regressions. We find that, if the price process is a driftless martingale, only the WPC is an unbiased estimator for the return variance explained during a time interval. For autocorrelated processes with a drift, only the R² of the unbiasedness regression is consistent, but it is biased for small samples.     

Estimating cost of capital in firm valuations with arithmetic or geometric mean – or better use the Cooper estimator?

We test the extent and determinants of bias effects of the arithmetic as well as the geometric mean estimator and the estimator of Cooper [1996. Arithmetic versus geometric mean estimators: Setting discount rates for capital budgeting. European Financial Management 2 (July): 157–67] regarding discount rate estimation for firm valuation by way of a bootstrap approach for 13 different countries. The Cooper estimator is superior to both the geometric and the (conventional) arithmetic mean estimator. However, a ‘truncated’ version of the arithmetic mean estimator leads generally to better estimation outcomes than the Cooper estimator. This means that, in order to reduce problems of upward-biased firm value estimates, expected cash flows beyond a certain time horizon are completely neglected in terminal value estimation. Such an approach seems particularly reasonable for the valuation of young growth companies as well as for companies from quickly developing countries such as Brazil, China, or Thailand, because the bias in terminal value estimation is increasing in the growth rate of future expected cash flows.     

Sampling distribution of the relative risk aversion estimator: Theory and applications

Brown and Gibbons (1985) developed a theory of relative risk aversion estimation in terms of average market rates of return and the variance of market rates of return. However, the exact sampling distributions of the relative risk aversion estimators have not been derived. The main purpose of this paper is to derive the exact sampling distribution of an appropriate relative risk aversion estimator. First, we have derived theoretically the density of Brown and Gibbons' maximum likelihood estimator. It is shown that the centralt is not appropriate for testing the significance of estimated relative risk aversion distribution. Then we derived the minimum variance unbiased estimator by a linear transformation of the Brown and Gibbons' maximum likelihood estimator. The density function is neither a central nor a noncentralt distribution. The density function of this new distribution has been tabulated. There is an empirical example to illustrate the application of this new sampling distribution.     

A parametric alternative to the Hill estimator for heavy-tailed distributions

Despite its wide use, the Hill estimator and its plot remain to be difficult to use in Extreme Value Theory (EVT) due to substantial sampling variations in extreme sample quantiles. In this paper, we propose a new plot we call the eigenvalue plot which can be seen as a generalization of the Hill plot. The theory behind the plot is based on a heavy-tailed parametric distribution class called the scaled Log phase-type (LogPH) distributions, a generalization of the ordinary LogPH distribution class which was previously used to model insurance claims data. We show that its tail property and moment condition are well aligned with EVT. Based on our findings, we construct the eigenvalue plot from fitting a shifted PH distribution to the excess log data with a minimal phase size. Through various numerical examples we illustrate and compare our method against the Hill plot.