Conditional Least Squares and Copulae in Claims Reserving for a Single Line of Business
|
|
- Victor Watkins
- 6 years ago
- Views:
Transcription
1 Conditional Least Squares and Copulae in Claims Reserving for a Single Line of Business Michal Pešta Charles University in Prague Faculty of Mathematics and Physics Ostap Okhrin Dresden University of Technology Institute of Economics and Transport
2 Introduction 2 30 Overview Motivated by claims reserving in non-life insurance Joint work with Ostap Okhrin (TU Dresden) Triangular data models Published in Insurance: Mathematics and Economics, 56 (May 2014):
3 Data structure 3 30 Triangular data i\j 1 2 n 1 n 1 Y 1,1 Y 1,2 Y 1,n 1 Y 1,n 2 Y 2,1 Y 2,2 Y 2,n Y i,n+1 i n 1 Y n 1,1 Y n 1,2 n Y n,1 n copies of stochastic process The first realization consists of n observations The last one has only one observation
4 Main aims 4 30 Terminology and goals Y i,j... cumulative payments in origin year i after j development periods (accounting year i + j) n... current year corresponds to the most recent accident year and development period Our data history consists of right-angled isosceles triangles Y i,j, where i + j n + 1 Predict Y i,n and R i = Y i,n Y i,n+1 i (claims reserve) Estimate distribution of the reserves
5 CMV model 5 30 Conditional mean and variance (CMV) model CMV model Y i,j = µ(y i,j 1, α, j) + σ(y i,j 1, β, j)ε i,j (α, β) α and β are unknown parameters, which dimensions do not depend on n µ is a continuous function in α σ is a positive and continuous function in β Errors ε i,j (α, β)
6 CMV model 6 30 CMV model s errors Disturbances {ε i,j (α, β)} n+1 i j=1 are independent sample copies of a stationary first-order Markov process for all i All ε i,j (α, β) have the common true invariant distribution G α,β which is absolutely continuous with respect to Lebesgue measure on the real line Filtration F i,j = σ(y k,l : l j, k i + 1 j) denotes the information set generated by that trapezoid E[ε i,j (α, β) F i,j 1 ] = 0 var[ε i,j (α, β) F i,j 1 ] = s(α, β)
7 CMV model 7 30 Properties of the CMV model Unknown true values [α, β ] of parameters [α, β ] set (due to identifiability purposes): s(α, β ) = 1 Model s name come from the fact that E[Y i,j F i,j 1 ] = µ(y i,j 1, α, j) var[y i,j F i,j 1 ] = σ 2 (Y i,j 1, β, j)s(α, β) Conditional mean models: types of ARMA models, vector autoregressions, linear and nonlinear regressions,... Conditional variance models: ARCH and any of its numerous parametric extensions (GARCH, EGARCH, GJR-GARCH, etc.), stochastic volatility models,...
8 CMV model 8 30 Candidates for the mean and variance function From the nature of data: Y i,j C i R + almost surely as j, i (stabilizing property) One may propose, e.g., µ(y i,j 1, α, j) = η(α, j)y i,j 1 σ(y i,j 1, β, j) = ν(β, j) Y i,j 1 η(α, j) should be decreasing in j with limit 1 as j ν(β, j) should be decreasing in j with limit 0 as j
9 CMV model 9 30 Dependence modeling Since the mean and variance trends are removed by the CMV model, the rest of the relationship among claim amounts Y i,j can be additionally captured by modeling dependent errors {ε i,j (α, β)} n+1 i j=1 are independent sample copies of a stationary first-order Markov process for all i generated from (G α,β ( ), C(, ; γ)) C(, ; γ) is the true parametric copula for [ε i,j 1 (α, β), ε i,j (α, β)], which is given and fixed up to unknown parameter γ
10 CMV model Copula-based model It is believed that there exist a kind of information overlap between the claims from consecutive development periods Joint bivariate distribution of [ε i,j 1 (α, β), ε i,j (α, β)] has distribution function H(e 1, e 2 ) = C(G α,β (e 1 ), G α,β (e 2 ); γ) Conditional copula density can be derived as h(e 2 e 1 ) = g α,β (e 2 )c(g α,β (e 1 ), G α,β (e 2 ); γ) where c is the copula density and g α,β is the marginal density corresponding to G α,β Play an important role in making" the dependent errors conditionally independent
11 Nonlinear generalized semiparametric regression Parameter estimation CMV model with copula assume three vector parameters to be estimated Estimation process consists of two stages In the first one, mean and variance parameters α and β are estimated in a distribution-free fashion, since no specific distributional assumptions are proposed nor required for the claims The second stage concerns estimation of the dependence structure, mainly the copula parameter γ, in a likelihood based way
12 Nonlinear generalized semiparametric regression Conditional least squares (CLS) Denote M n (α, β) = 1 n n+1 j 1 n 1 n + 1 j j=2 i=1 [ Yi,j µ(y i,j 1, α, j) ] 2 σ 2 (Y i,j 1, β, j) n n+1 j i=1 V n (α, β) = 1 1 n 1 n + 1 j j=2 { [Yi,j µ(y i,j 1, α, j) ]2 2 σ 2 (Y i,j 1, β, j)}
13 Nonlinear generalized semiparametric regression CLS estimates CLS estimate of the mean parameter α for a fixed value of parameter β Θ 2 is defined as α(β) = arg min α Θ 1 M n (α, β) and CLS estimate of the variance parameter β for a fixed value of parameter α Θ 1 is defined as β(α) = arg min β Θ 2 V n (α, β) Computationally not feasible to find the global minimum of M n and V n with respect to [α, β ] simultaneously
14 Nonlinear generalized semiparametric regression Consistency Under regularity conditions α(β) P n α P (β), β; β(α) n β (α), α Mixingales are to mixing processes as martingale differences are to independent processes
15 Nonlinear generalized semiparametric regression Iterative CLS What is the connection between the true unknown parameter values α and β of the CMV model and true unknown parameter values α (β) and β (α)? [ α(β ) β(α ) ] [ P α n β Iteratively estimate α given the fixed value of β and, consequently, estimate β given the fixed value of α (obtained from previous step) Repeat in turns until almost no change in consecutive estimates of [α, β ] ]
16 Nonlinear generalized semiparametric regression Estimation of dependence structure Estimate the unknown marginal distribution function G α,β of CMV model errors ε i,j (α, β) non-parametrically by the empirical distribution function Ĝ n (e) = 1 n(n 1)/2 + 1 of the fitted residuals n 1 i=1 ε i,j ( α, β) = Y i,j µ(y i,j 1, α, j) σ(y i,j 1, β, j) n+1 i I{ ε i,j ( α, β) e} j=2
17 Nonlinear generalized semiparametric regression Likelihood for copula Full log-likelihood for copula parameter γ L(γ) = + n 2 i=1 n 2 i=1 n+1 i log g α,β (ε i,j (α, β)) j=2 n+1 i log c(g α,β (ε i,j 1 (α, β)), G α,β (ε i,j (α, β)); γ) j=3
18 Nonlinear generalized semiparametric regression Psuedo (quasi) likelihood Ignoring the first term in L(γ) and replacing ε s and G α,β by their estimated counterparts ε s and Ĝ n, parameter γ can be estimated by the so-called canonical maximum likelihood, i.e., maximizing the partial (pseudo) log-likelihood γ = arg max γ L(γ) = n 2 i=1 L(γ) n+1 i log c(ĝ n ( ε i,j 1 ( α, β)), Ĝ n ( ε i,j ( α, β)); γ) j=3
19 Nonlinear generalized semiparametric regression Prediction Predictor for reserve R (n) i R (n) i = Ŷ i,n Y i,n+1 i can be defined as Prediction of unobserved claims may be done in a telescopic way based on the CMV model formulation: start with the diagonal element Y i,n+1 i and predict Y i,j, j > n + 1 i stepwise in each row Ŷ i,j = Y i,j, i + j n + 1 Ŷ i,j = µ(ŷ i,j 1, α, j) + σ(ŷ i,j 1, β, j) ε j, i + j > n + 1
20 Nonlinear generalized semiparametric regression Semiparametric bootstrap Errors ε j are simulated from the fitted residuals Takes advantage of the fact that ε i,j (α, β) = G 1 α,β (X j) for all i (due to the independent rows), where {X j } n j=2 is a stationary first-order Markov process with the copula C(x 1, x 2 ; γ) being the joint distribution of [X j 1, X j ]
21 Nonlinear generalized semiparametric regression Resampling algorithm Generate n 1 independent Un(0, 1) rvs {X j } n j=2 Repeat b = 1,..., B (b) U 2 X 2 (b) ε 2 Ĝn ( (b) U 2 ) (b) U j C (X j (b) U j 1 ; γ), j = 3,..., n (b) ε j Ĝn ( (b) U j ), j = 3,..., n Center bootstrap residuals (b) ε j (b) ε j 1 n 1 n l=2 (b) ε l (b) Ŷ i,n+1 j Y i,n+1 i (b) Ŷ i,j µ( (b) Ŷ i,j 1, α, j) + σ( (b) Ŷ i,j 1, β, j) (b) ε j, j = n + 2 i,..., n (b) R (n) i (b) Ŷ i,n Y i,n+1 i
22 Real data example Development period j Y ij a
23 Real data example Real data Data set from Zehnwirth and Barnett (2000) ( µ(y i,j 1, α, j) = 1 + α 1 α 2 j 1 α 2 exp { α 1 j α }) 2 Y i,j 1 σ(y i,j 1, β, j) = β 1 exp{ β 2 j} Y i,j 1 CLS estimates: α 1 = 2.033, α 2 = 1.106, β 1 = 109.8, β 2 =
24 Real data example Residuals ε^i,j Development period j Still some slight pattern (trend) not captured by mean and variance parametric part
25 Real data example Copula godness-of-fit Kendall τ for the pairs of consecutive residuals {[ ε i,j 1 ( α, β), n 2,n+1 i ε i,j ( α, β)]} i=1,j=3 equals 0.43, which indicates at least mild dependence Three Archimedean copulae (Clayton, Frank, and Gumbel) together with Gaussian and Student t 5 -copula considered S (C) n goodness-of-fit test proposed by Genest et al. (2009) Gumbel copula ( γ = 1.776) was chosen Exhibits strong right tail dependence and relatively weak left tail dependence Transformed residuals (by the residuals marginal edf Ĝ n ; having uniform margins) seem to be strongly correlated at high values but less correlated at low values
26 Real data example Residuals ε^i,j 1 ε^i,j Simulated Residuals ε^j 1 ε^j Transformed Residuals G^ n(ε^i,j 1) G^ n(ε^i,j) Gumbel(γ^) Copula Density Un[0,1] Un[0,1]
27 Real data example Results Benchmark: traditional bootstrapped chain ladder (BCL) Disadvantages, which can be overcome by our approach: Number of parameters depending on the sample size Some parameters estimated by just ratio of two numbers (yielding zero sample variance) Questionable consistency of the estimates Non-realistic assumption of independence of the residuals Our approach: Slightly smaller predictions of reserves But even more important is that the estimates of the reserves distribution are less volatile (interesting finding in terms of CoV: GEE AR(1) 1.95; CLSC 1.99)
28 Reserves BCL+CLSC i=2 BCL+CLSC BCL+CLSC i=3 i=4 0k 25k Reserves 25k 50k 75k Reserves 50k 75k 100k BCL+CLSC i=5 BCL+CLSC i=6 75k 100k 125k 150k Reserves 100k 150k 200k Reserves 150k 200k 250k BCL+CLSC i=8 BCL+CLSC i=7 BCL+CLSC i=9 300k 350k 400k 450k 500k Reserves 600k 675k 750k 825k 900k Reserves 1100k 1250k 1400k 1550k BCL+CLSC i=11 BCL+CLSC i= k 5000k 5250k 5500k 5750k 6000k Total Reserves 1775k 2000k 2225k 2450k 2675k Reserves Reserves BCL+CLSC Reserves Reserves
29 Conclusions Summary Conditional mean and variance (CMV) time series model for triangular data with innovations being a stationary first-order Markov process Framework is demonstrated to be suitable for stochastic claims reserving in general insurance Very flexible modeling approach, relatively smaller number of model parameters not depending on the number of development periods, and time series innovations not considered as independent Increase in precision of the claims reserves prediction Theoretical justification of the proposed approached shown
30 Questions? Thank you!
arxiv: v1 [stat.ap] 19 Jun 2013
Conditional Least Squares and Copulae in Claims Reserving for a Single Line of Business Michal Pešta a,, Ostap Okhrin b arxiv:1306.4529v1 [stat.ap] 19 Jun 2013 a Charles University in Prague, Faculty of
More informationBootstrapping the triangles
Bootstrapping the triangles Prečo, kedy a ako (NE)bootstrapovat v trojuholníkoch? Michal Pešta Charles University in Prague Faculty of Mathematics and Physics Actuarial Seminar, Prague 19 October 01 Overview
More informationForecasting with the age-period-cohort model and the extended chain-ladder model
Forecasting with the age-period-cohort model and the extended chain-ladder model By D. KUANG Department of Statistics, University of Oxford, Oxford OX1 3TG, U.K. di.kuang@some.ox.ac.uk B. Nielsen Nuffield
More informationClaims Reserving under Solvency II
Claims Reserving under Solvency II Mario V. Wüthrich RiskLab, ETH Zurich Swiss Finance Institute Professor joint work with Michael Merz (University of Hamburg) April 21, 2016 European Congress of Actuaries,
More informationONE-YEAR AND TOTAL RUN-OFF RESERVE RISK ESTIMATORS BASED ON HISTORICAL ULTIMATE ESTIMATES
FILIPPO SIEGENTHALER / filippo78@bluewin.ch 1 ONE-YEAR AND TOTAL RUN-OFF RESERVE RISK ESTIMATORS BASED ON HISTORICAL ULTIMATE ESTIMATES ABSTRACT In this contribution we present closed-form formulas in
More informationGARCH Models Estimation and Inference
GARCH Models Estimation and Inference Eduardo Rossi University of Pavia December 013 Rossi GARCH Financial Econometrics - 013 1 / 1 Likelihood function The procedure most often used in estimating θ 0 in
More informationVarious Extensions Based on Munich Chain Ladder Method
Various Extensions Based on Munich Chain Ladder Method etr Jedlička Charles University, Department of Statistics 20th June 2007, 50th Anniversary ASTIN Colloquium etr Jedlička (UK MFF) Various Extensions
More informationConvolution Based Unit Root Processes: a Simulation Approach
International Journal of Statistics and Probability; Vol., No. 6; November 26 ISSN 927-732 E-ISSN 927-74 Published by Canadian Center of Science and Education Convolution Based Unit Root Processes: a Simulation
More informationJuly 31, 2009 / Ben Kedem Symposium
ing The s ing The Department of Statistics North Carolina State University July 31, 2009 / Ben Kedem Symposium Outline ing The s 1 2 s 3 4 5 Ben Kedem ing The s Ben has made many contributions to time
More informationEfficient estimation of a semiparametric dynamic copula model
Efficient estimation of a semiparametric dynamic copula model Christian Hafner Olga Reznikova Institute of Statistics Université catholique de Louvain Louvain-la-Neuve, Blgium 30 January 2009 Young Researchers
More informationGARY G. VENTER 1. CHAIN LADDER VARIANCE FORMULA
DISCUSSION OF THE MEAN SQUARE ERROR OF PREDICTION IN THE CHAIN LADDER RESERVING METHOD BY GARY G. VENTER 1. CHAIN LADDER VARIANCE FORMULA For a dozen or so years the chain ladder models of Murphy and Mack
More informationFinancial Econometrics and Volatility Models Copulas
Financial Econometrics and Volatility Models Copulas Eric Zivot Updated: May 10, 2010 Reading MFTS, chapter 19 FMUND, chapters 6 and 7 Introduction Capturing co-movement between financial asset returns
More informationGeneralized Autoregressive Score Models
Generalized Autoregressive Score Models by: Drew Creal, Siem Jan Koopman, André Lucas To capture the dynamic behavior of univariate and multivariate time series processes, we can allow parameters to be
More information13. Estimation and Extensions in the ARCH model. MA6622, Ernesto Mordecki, CityU, HK, References for this Lecture:
13. Estimation and Extensions in the ARCH model MA6622, Ernesto Mordecki, CityU, HK, 2006. References for this Lecture: Robert F. Engle. GARCH 101: The Use of ARCH/GARCH Models in Applied Econometrics,
More informationMultivariate Non-Normally Distributed Random Variables
Multivariate Non-Normally Distributed Random Variables An Introduction to the Copula Approach Workgroup seminar on climate dynamics Meteorological Institute at the University of Bonn 18 January 2008, Bonn
More informationOn the Systemic Nature of Weather Risk
Martin Odening 1 Ostap Okhrin 2 Wei Xu 1 Department of Agricultural Economics 1 Ladislaus von Bortkiewicz Chair of Statistics C.A.S.E. Center for Applied Statistics and Economics 2 Humboldt Universität
More informationLecture Quantitative Finance Spring Term 2015
on bivariate Lecture Quantitative Finance Spring Term 2015 Prof. Dr. Erich Walter Farkas Lecture 07: April 2, 2015 1 / 54 Outline on bivariate 1 2 bivariate 3 Distribution 4 5 6 7 8 Comments and conclusions
More informationGARCH Models Estimation and Inference. Eduardo Rossi University of Pavia
GARCH Models Estimation and Inference Eduardo Rossi University of Pavia Likelihood function The procedure most often used in estimating θ 0 in ARCH models involves the maximization of a likelihood function
More informationSwitching Regime Estimation
Switching Regime Estimation Series de Tiempo BIrkbeck March 2013 Martin Sola (FE) Markov Switching models 01/13 1 / 52 The economy (the time series) often behaves very different in periods such as booms
More informationIdentification of the age-period-cohort model and the extended chain ladder model
Identification of the age-period-cohort model and the extended chain ladder model By D. KUANG Department of Statistics, University of Oxford, Oxford OX TG, U.K. di.kuang@some.ox.ac.uk B. Nielsen Nuffield
More informationTime Varying Hierarchical Archimedean Copulae (HALOC)
Time Varying Hierarchical Archimedean Copulae () Wolfgang Härdle Ostap Okhrin Yarema Okhrin Ladislaus von Bortkiewicz Chair of Statistics C.A.S.E. Center for Applied Statistics and Economics Humboldt-Universität
More informationAppendix A: The time series behavior of employment growth
Unpublished appendices from The Relationship between Firm Size and Firm Growth in the U.S. Manufacturing Sector Bronwyn H. Hall Journal of Industrial Economics 35 (June 987): 583-606. Appendix A: The time
More informationNonlinear Parameter Estimation for State-Space ARCH Models with Missing Observations
Nonlinear Parameter Estimation for State-Space ARCH Models with Missing Observations SEBASTIÁN OSSANDÓN Pontificia Universidad Católica de Valparaíso Instituto de Matemáticas Blanco Viel 596, Cerro Barón,
More informationExtreme Value Analysis and Spatial Extremes
Extreme Value Analysis and Department of Statistics Purdue University 11/07/2013 Outline Motivation 1 Motivation 2 Extreme Value Theorem and 3 Bayesian Hierarchical Models Copula Models Max-stable Models
More informationDELTA METHOD and RESERVING
XXXVI th ASTIN COLLOQUIUM Zurich, 4 6 September 2005 DELTA METHOD and RESERVING C.PARTRAT, Lyon 1 university (ISFA) N.PEY, AXA Canada J.SCHILLING, GIE AXA Introduction Presentation of methods based on
More informationA copula goodness-of-t approach. conditional probability integral transform. Daniel Berg 1 Henrik Bakken 2
based on the conditional probability integral transform Daniel Berg 1 Henrik Bakken 2 1 Norwegian Computing Center (NR) & University of Oslo (UiO) 2 Norwegian University of Science and Technology (NTNU)
More informationOn the Importance of Dispersion Modeling for Claims Reserving: Application of the Double GLM Theory
On the Importance of Dispersion Modeling for Claims Reserving: Application of the Double GLM Theory Danaïl Davidov under the supervision of Jean-Philippe Boucher Département de mathématiques Université
More informationCorrelation: Copulas and Conditioning
Correlation: Copulas and Conditioning This note reviews two methods of simulating correlated variates: copula methods and conditional distributions, and the relationships between them. Particular emphasis
More informationGARCH Models. Eduardo Rossi University of Pavia. December Rossi GARCH Financial Econometrics / 50
GARCH Models Eduardo Rossi University of Pavia December 013 Rossi GARCH Financial Econometrics - 013 1 / 50 Outline 1 Stylized Facts ARCH model: definition 3 GARCH model 4 EGARCH 5 Asymmetric Models 6
More informationExplicit Bounds for the Distribution Function of the Sum of Dependent Normally Distributed Random Variables
Explicit Bounds for the Distribution Function of the Sum of Dependent Normally Distributed Random Variables Walter Schneider July 26, 20 Abstract In this paper an analytic expression is given for the bounds
More informationRegression Models - Introduction
Regression Models - Introduction In regression models, two types of variables that are studied: A dependent variable, Y, also called response variable. It is modeled as random. An independent variable,
More informationThe Instability of Correlations: Measurement and the Implications for Market Risk
The Instability of Correlations: Measurement and the Implications for Market Risk Prof. Massimo Guidolin 20254 Advanced Quantitative Methods for Asset Pricing and Structuring Winter/Spring 2018 Threshold
More informationEcon 423 Lecture Notes: Additional Topics in Time Series 1
Econ 423 Lecture Notes: Additional Topics in Time Series 1 John C. Chao April 25, 2017 1 These notes are based in large part on Chapter 16 of Stock and Watson (2011). They are for instructional purposes
More informationModelling and Estimation of Stochastic Dependence
Modelling and Estimation of Stochastic Dependence Uwe Schmock Based on joint work with Dr. Barbara Dengler Financial and Actuarial Mathematics and Christian Doppler Laboratory for Portfolio Risk Management
More informationVolatility. Gerald P. Dwyer. February Clemson University
Volatility Gerald P. Dwyer Clemson University February 2016 Outline 1 Volatility Characteristics of Time Series Heteroskedasticity Simpler Estimation Strategies Exponentially Weighted Moving Average Use
More informationThe Use of Copulas to Model Conditional Expectation for Multivariate Data
The Use of Copulas to Model Conditional Expectation for Multivariate Data Kääri, Meelis; Selart, Anne; Kääri, Ene University of Tartu, Institute of Mathematical Statistics J. Liivi 2 50409 Tartu, Estonia
More informationSongklanakarin Journal of Science and Technology SJST R1 Sukparungsee
Songklanakarin Journal of Science and Technology SJST-0-0.R Sukparungsee Bivariate copulas on the exponentially weighted moving average control chart Journal: Songklanakarin Journal of Science and Technology
More informationModelling Dropouts by Conditional Distribution, a Copula-Based Approach
The 8th Tartu Conference on MULTIVARIATE STATISTICS, The 6th Conference on MULTIVARIATE DISTRIBUTIONS with Fixed Marginals Modelling Dropouts by Conditional Distribution, a Copula-Based Approach Ene Käärik
More informationECONOMICS 7200 MODERN TIME SERIES ANALYSIS Econometric Theory and Applications
ECONOMICS 7200 MODERN TIME SERIES ANALYSIS Econometric Theory and Applications Yongmiao Hong Department of Economics & Department of Statistical Sciences Cornell University Spring 2019 Time and uncertainty
More informationSemi-parametric estimation of non-stationary Pickands functions
Semi-parametric estimation of non-stationary Pickands functions Linda Mhalla 1 Joint work with: Valérie Chavez-Demoulin 2 and Philippe Naveau 3 1 Geneva School of Economics and Management, University of
More informationDiscussion of Bootstrap prediction intervals for linear, nonlinear, and nonparametric autoregressions, by Li Pan and Dimitris Politis
Discussion of Bootstrap prediction intervals for linear, nonlinear, and nonparametric autoregressions, by Li Pan and Dimitris Politis Sílvia Gonçalves and Benoit Perron Département de sciences économiques,
More informationSemi-parametric predictive inference for bivariate data using copulas
Semi-parametric predictive inference for bivariate data using copulas Tahani Coolen-Maturi a, Frank P.A. Coolen b,, Noryanti Muhammad b a Durham University Business School, Durham University, Durham, DH1
More informationStatistical Inference and Methods
Department of Mathematics Imperial College London d.stephens@imperial.ac.uk http://stats.ma.ic.ac.uk/ das01/ 31st January 2006 Part VI Session 6: Filtering and Time to Event Data Session 6: Filtering and
More informationBayesian Semiparametric GARCH Models
Bayesian Semiparametric GARCH Models Xibin (Bill) Zhang and Maxwell L. King Department of Econometrics and Business Statistics Faculty of Business and Economics xibin.zhang@monash.edu Quantitative Methods
More informationComputer Intensive Methods in Mathematical Statistics
Computer Intensive Methods in Mathematical Statistics Department of mathematics johawes@kth.se Lecture 16 Advanced topics in computational statistics 18 May 2017 Computer Intensive Methods (1) Plan of
More informationBayesian Semiparametric GARCH Models
Bayesian Semiparametric GARCH Models Xibin (Bill) Zhang and Maxwell L. King Department of Econometrics and Business Statistics Faculty of Business and Economics xibin.zhang@monash.edu Quantitative Methods
More informationMarginal Specifications and a Gaussian Copula Estimation
Marginal Specifications and a Gaussian Copula Estimation Kazim Azam Abstract Multivariate analysis involving random variables of different type like count, continuous or mixture of both is frequently required
More informationCalendar Year Dependence Modeling in Run-Off Triangles
Calendar Year Dependence Modeling in Run-Off Triangles Mario V. Wüthrich RiskLab, ETH Zurich May 21-24, 2013 ASTIN Colloquium The Hague www.math.ethz.ch/ wueth Claims reserves at time I = 2010 accident
More informationModel Specification Testing in Nonparametric and Semiparametric Time Series Econometrics. Jiti Gao
Model Specification Testing in Nonparametric and Semiparametric Time Series Econometrics Jiti Gao Department of Statistics School of Mathematics and Statistics The University of Western Australia Crawley
More informationSymmetric btw positive & negative prior returns. where c is referred to as risk premium, which is expected to be positive.
Advantages of GARCH model Simplicity Generates volatility clustering Heavy tails (high kurtosis) Weaknesses of GARCH model Symmetric btw positive & negative prior returns Restrictive Provides no explanation
More informationGARCH Models Estimation and Inference
Università di Pavia GARCH Models Estimation and Inference Eduardo Rossi Likelihood function The procedure most often used in estimating θ 0 in ARCH models involves the maximization of a likelihood function
More informationAutoregressive and Moving-Average Models
Chapter 3 Autoregressive and Moving-Average Models 3.1 Introduction Let y be a random variable. We consider the elements of an observed time series {y 0,y 1,y2,...,y t } as being realizations of this randoms
More informationWeb-based Supplementary Material for A Two-Part Joint. Model for the Analysis of Survival and Longitudinal Binary. Data with excess Zeros
Web-based Supplementary Material for A Two-Part Joint Model for the Analysis of Survival and Longitudinal Binary Data with excess Zeros Dimitris Rizopoulos, 1 Geert Verbeke, 1 Emmanuel Lesaffre 1 and Yves
More informationGaussian kernel GARCH models
Gaussian kernel GARCH models Xibin (Bill) Zhang and Maxwell L. King Department of Econometrics and Business Statistics Faculty of Business and Economics 7 June 2013 Motivation A regression model is often
More informationReserving for multiple excess layers
Reserving for multiple excess layers Ben Zehnwirth and Glen Barnett Abstract Patterns and changing trends among several excess-type layers on the same business tend to be closely related. The changes in
More information1. The OLS Estimator. 1.1 Population model and notation
1. The OLS Estimator OLS stands for Ordinary Least Squares. There are 6 assumptions ordinarily made, and the method of fitting a line through data is by least-squares. OLS is a common estimation methodology
More informationGENERALIZED LINEAR MODELING APPROACH TO STOCHASTIC WEATHER GENERATORS
GENERALIZED LINEAR MODELING APPROACH TO STOCHASTIC WEATHER GENERATORS Rick Katz Institute for Study of Society and Environment National Center for Atmospheric Research Boulder, CO USA Joint work with Eva
More informationMarkov Switching Regular Vine Copulas
Int. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin (Session CPS057) p.5304 Markov Switching Regular Vine Copulas Stöber, Jakob and Czado, Claudia Lehrstuhl für Mathematische Statistik,
More informationCalibration Estimation of Semiparametric Copula Models with Data Missing at Random
Calibration Estimation of Semiparametric Copula Models with Data Missing at Random Shigeyuki Hamori 1 Kaiji Motegi 1 Zheng Zhang 2 1 Kobe University 2 Renmin University of China Econometrics Workshop UNC
More informationBivariate Rainfall and Runoff Analysis Using Entropy and Copula Theories
Entropy 2012, 14, 1784-1812; doi:10.3390/e14091784 Article OPEN ACCESS entropy ISSN 1099-4300 www.mdpi.com/journal/entropy Bivariate Rainfall and Runoff Analysis Using Entropy and Copula Theories Lan Zhang
More information1 Class Organization. 2 Introduction
Time Series Analysis, Lecture 1, 2018 1 1 Class Organization Course Description Prerequisite Homework and Grading Readings and Lecture Notes Course Website: http://www.nanlifinance.org/teaching.html wechat
More informationGARCH processes probabilistic properties (Part 1)
GARCH processes probabilistic properties (Part 1) Alexander Lindner Centre of Mathematical Sciences Technical University of Munich D 85747 Garching Germany lindner@ma.tum.de http://www-m1.ma.tum.de/m4/pers/lindner/
More informationAnalytics Software. Beyond deterministic chain ladder reserves. Neil Covington Director of Solutions Management GI
Analytics Software Beyond deterministic chain ladder reserves Neil Covington Director of Solutions Management GI Objectives 2 Contents 01 Background 02 Formulaic Stochastic Reserving Methods 03 Bootstrapping
More informationarxiv: v1 [stat.me] 9 Feb 2012
Modeling high dimensional time-varying dependence using D-vine SCAR models Carlos Almeida a, Claudia Czado b, Hans Manner c, arxiv:1202.2008v1 [stat.me] 9 Feb 2012 a Georges Lemaitre Centre for Earth and
More informationVine copula specifications for stationary multivariate Markov chains
Vine copula specifications for stationary multivariate Markov chains Brendan K. Beare and Juwon Seo Department of Economics, University of California, San Diego January 14, 2014 Abstract Vine copulae provide
More informationRobustness of a semiparametric estimator of a copula
Robustness of a semiparametric estimator of a copula Gunky Kim a, Mervyn J. Silvapulle b and Paramsothy Silvapulle c a Department of Econometrics and Business Statistics, Monash University, c Caulfield
More informationChain ladder with random effects
Greg Fry Consulting Actuaries Sydney Australia Astin Colloquium, The Hague 21-24 May 2013 Overview Fixed effects models Families of chain ladder models Maximum likelihood estimators Random effects models
More informationA measure of radial asymmetry for bivariate copulas based on Sobolev norm
A measure of radial asymmetry for bivariate copulas based on Sobolev norm Ahmad Alikhani-Vafa Ali Dolati Abstract The modified Sobolev norm is used to construct an index for measuring the degree of radial
More informationForecasting time series with multivariate copulas
Depend. Model. 2015; 3:59 82 Research Article Open Access Clarence Simard* and Bruno Rémillard Forecasting time series with multivariate copulas DOI 10.1515/demo-2015-0005 5 Received august 17, 2014; accepted
More informationECON 3150/4150, Spring term Lecture 6
ECON 3150/4150, Spring term 2013. Lecture 6 Review of theoretical statistics for econometric modelling (II) Ragnar Nymoen University of Oslo 31 January 2013 1 / 25 References to Lecture 3 and 6 Lecture
More informationConsistency of the maximum likelihood estimator for general hidden Markov models
Consistency of the maximum likelihood estimator for general hidden Markov models Jimmy Olsson Centre for Mathematical Sciences Lund University Nordstat 2012 Umeå, Sweden Collaborators Hidden Markov models
More informationSimulating Exchangeable Multivariate Archimedean Copulas and its Applications. Authors: Florence Wu Emiliano A. Valdez Michael Sherris
Simulating Exchangeable Multivariate Archimedean Copulas and its Applications Authors: Florence Wu Emiliano A. Valdez Michael Sherris Literatures Frees and Valdez (1999) Understanding Relationships Using
More informationUsing copulas to model time dependence in stochastic frontier models
Using copulas to model time dependence in stochastic frontier models Christine Amsler Michigan State University Artem Prokhorov Concordia University November 2008 Peter Schmidt Michigan State University
More informationMultivariate Random Variable
Multivariate Random Variable Author: Author: Andrés Hincapié and Linyi Cao This Version: August 7, 2016 Multivariate Random Variable 3 Now we consider models with more than one r.v. These are called multivariate
More informationSystemic Weather Risk and Crop Insurance: The Case of China
and Crop Insurance: The Case of China Ostap Okhrin 1 Martin Odening 2 Wei Xu 3 Ladislaus von Bortkiewicz Chair of Statistics C.A.S.E. Center for Applied Statistics and Economics 1 Department of Agricultural
More information2.5 Forecasting and Impulse Response Functions
2.5 Forecasting and Impulse Response Functions Principles of forecasting Forecast based on conditional expectations Suppose we are interested in forecasting the value of y t+1 based on a set of variables
More informationSTAT 100C: Linear models
STAT 100C: Linear models Arash A. Amini June 9, 2018 1 / 56 Table of Contents Multiple linear regression Linear model setup Estimation of β Geometric interpretation Estimation of σ 2 Hat matrix Gram matrix
More informationRecent Advances in the Field of Trade Theory and Policy Analysis Using Micro-Level Data
Recent Advances in the Field of Trade Theory and Policy Analysis Using Micro-Level Data July 2012 Bangkok, Thailand Cosimo Beverelli (World Trade Organization) 1 Content a) Classical regression model b)
More informationDependence. Practitioner Course: Portfolio Optimization. John Dodson. September 10, Dependence. John Dodson. Outline.
Practitioner Course: Portfolio Optimization September 10, 2008 Before we define dependence, it is useful to define Random variables X and Y are independent iff For all x, y. In particular, F (X,Y ) (x,
More informationReducing Model Risk With Goodness-of-fit Victory Idowu London School of Economics
Reducing Model Risk With Goodness-of-fit Victory Idowu London School of Economics Agenda I. An overview of Copula Theory II. Copulas and Model Risk III. Goodness-of-fit methods for copulas IV. Presentation
More informationGenerated Covariates in Nonparametric Estimation: A Short Review.
Generated Covariates in Nonparametric Estimation: A Short Review. Enno Mammen, Christoph Rothe, and Melanie Schienle Abstract In many applications, covariates are not observed but have to be estimated
More informationMultivariate GARCH models.
Multivariate GARCH models. Financial market volatility moves together over time across assets and markets. Recognizing this commonality through a multivariate modeling framework leads to obvious gains
More informationSynchronous bootstrapping of loss reserves
Synchronous bootstrapping of loss reserves Greg Taylor Taylor Fry Consulting Actuaries University of Melbourne University of New South Wales Gráinne McGuire Taylor Fry Consulting Actuaries ASTIN Colloquium
More informationFORECASTING IN FINANCIAL MARKET USING MARKOV R MARKOV REGIME SWITCHING AND PRINCIPAL COMPONENT ANALYSIS.
FORECASTING IN FINANCIAL MARKET USING MARKOV REGIME SWITCHING AND PRINCIPAL COMPONENT ANALYSIS. September 13, 2012 1 2 3 4 5 Heteroskedasticity Model Multicollinearily 6 In Thai Outline การพยากรณ ในตลาดทางการเง
More informationImputation Algorithm Using Copulas
Metodološki zvezki, Vol. 3, No. 1, 2006, 109-120 Imputation Algorithm Using Copulas Ene Käärik 1 Abstract In this paper the author demonstrates how the copulas approach can be used to find algorithms for
More informationDefine y t+h t as the forecast of y t+h based on I t known parameters. The forecast error is. Forecasting
Forecasting Let {y t } be a covariance stationary are ergodic process, eg an ARMA(p, q) process with Wold representation y t = X μ + ψ j ε t j, ε t ~WN(0,σ 2 ) j=0 = μ + ε t + ψ 1 ε t 1 + ψ 2 ε t 2 + Let
More informationAnalysing geoadditive regression data: a mixed model approach
Analysing geoadditive regression data: a mixed model approach Institut für Statistik, Ludwig-Maximilians-Universität München Joint work with Ludwig Fahrmeir & Stefan Lang 25.11.2005 Spatio-temporal regression
More informationThe Realized Hierarchical Archimedean Copula in Risk Modelling
econometrics Article The Realized Hierarchical Archimedean Copula in Risk Modelling Ostap Okhrin 1 and Anastasija Tetereva 2, * 1 Chair of Econometrics and Statistics esp. Transportation, Institute of
More informationTIME SERIES ANALYSIS. Forecasting and Control. Wiley. Fifth Edition GWILYM M. JENKINS GEORGE E. P. BOX GREGORY C. REINSEL GRETA M.
TIME SERIES ANALYSIS Forecasting and Control Fifth Edition GEORGE E. P. BOX GWILYM M. JENKINS GREGORY C. REINSEL GRETA M. LJUNG Wiley CONTENTS PREFACE TO THE FIFTH EDITION PREFACE TO THE FOURTH EDITION
More informationGoodness-of-Fit Tests for Time Series Models: A Score-Marked Empirical Process Approach
Goodness-of-Fit Tests for Time Series Models: A Score-Marked Empirical Process Approach By Shiqing Ling Department of Mathematics Hong Kong University of Science and Technology Let {y t : t = 0, ±1, ±2,
More informationTime Series Models and Inference. James L. Powell Department of Economics University of California, Berkeley
Time Series Models and Inference James L. Powell Department of Economics University of California, Berkeley Overview In contrast to the classical linear regression model, in which the components of the
More informationStatistics 910, #15 1. Kalman Filter
Statistics 910, #15 1 Overview 1. Summary of Kalman filter 2. Derivations 3. ARMA likelihoods 4. Recursions for the variance Kalman Filter Summary of Kalman filter Simplifications To make the derivations
More informationApplied Econometrics (QEM)
Applied Econometrics (QEM) The Simple Linear Regression Model based on Prinicples of Econometrics Jakub Mućk Department of Quantitative Economics Jakub Mućk Applied Econometrics (QEM) Meeting #2 The Simple
More informationCopulas, a novel approach to model spatial and spatio-temporal dependence
Copulas, a novel approach to model spatial and spatio-temporal dependence Benedikt Gräler 1, Hannes Kazianka 2, Giovana Mira de Espindola 3 1 Institute for Geoinformatics, University of Münster, Germany
More informationUsing all observations when forecasting under structural breaks
Using all observations when forecasting under structural breaks Stanislav Anatolyev New Economic School Victor Kitov Moscow State University December 2007 Abstract We extend the idea of the trade-off window
More informationHow to select a good vine
Universitetet i Oslo ingrihaf@math.uio.no International FocuStat Workshop on Focused Information Criteria and Related Themes, May 9-11, 2016 Copulae Regular vines Model selection and reduction Limitations
More informationClass 1: Stationary Time Series Analysis
Class 1: Stationary Time Series Analysis Macroeconometrics - Fall 2009 Jacek Suda, BdF and PSE February 28, 2011 Outline Outline: 1 Covariance-Stationary Processes 2 Wold Decomposition Theorem 3 ARMA Models
More informationTime Series Models for Measuring Market Risk
Time Series Models for Measuring Market Risk José Miguel Hernández Lobato Universidad Autónoma de Madrid, Computer Science Department June 28, 2007 1/ 32 Outline 1 Introduction 2 Competitive and collaborative
More informationGeneralized Linear Models. Kurt Hornik
Generalized Linear Models Kurt Hornik Motivation Assuming normality, the linear model y = Xβ + e has y = β + ε, ε N(0, σ 2 ) such that y N(μ, σ 2 ), E(y ) = μ = β. Various generalizations, including general
More informationMath 180B Problem Set 3
Math 180B Problem Set 3 Problem 1. (Exercise 3.1.2) Solution. By the definition of conditional probabilities we have Pr{X 2 = 1, X 3 = 1 X 1 = 0} = Pr{X 3 = 1 X 2 = 1, X 1 = 0} Pr{X 2 = 1 X 1 = 0} = P
More information