Introduction to vine copulas

Similar documents
Vine-copula Based Models for Farmland Portfolio Management

Will QE Change the dependence between Baht/Dollar Exchange Rates and Price Returns of AOT and MINT?

INTERNATIONAL JOURNAL FOR INNOVATIVE RESEARCH IN MULTIDISCIPLINARY FIELD ISSN Volume - 3, Issue - 2, Feb

DEPARTMENT OF ECONOMICS AND FINANCE COLLEGE OF BUSINESS AND ECONOMICS UNIVERSITY OF CANTERBURY CHRISTCHURCH, NEW ZEALAND

Application of Vine Copulas to Credit Portfolio Risk Modeling

A Vine Copula Approach for Analyzing Financial Risk and Co-movement of the Indonesian, Philippine and Thailand Stock Markets

Dependence estimation and controlled CVaR portfolio optimization of a highly kurtotic Australian mining sample of stocks

Asymmetric Price Transmission: A Copula Approach

arxiv: v1 [cs.ce] 4 Apr 2017

Volatility Models and Their Applications

MEASURING PORTFOLIO RISKS USING CONDITIONAL COPULA-AR-GARCH MODEL

Uncertainty Analysis with UNICORN

Regime Switching Vine Copula Models for Global Equity and Volatility Indices

Risk Measurement of Multivariate Credit Portfolio based on M-Copula Functions*

ON A PROBLEM BY SCHWEIZER AND SKLAR

MODELING DEPENDENCY RELATIONSHIPS WITH COPULAS

An Empirical Analysis of the Dependence Structure of International Equity and Bond Markets Using Regime-switching Copula Model

An empirical investigation of the short-term relationship between interest rate risk and credit risk

Pair Copula Constructions for Insurance Experience Rating

Risk Factors, Copula Dependence and Risk Sensitivity of a Large Portfolio

FLEXIBLE MODELING OF MULTIVARIATE RISKS IN PRICING MARGIN PROTECTION INSURANCE: MODELING PORTFOLIO RISKS WITH MIXTURES OF MIXTURES

Page 2 Vol. 10 Issue 7 (Ver 1.0) August 2010

Operational Risk Modeling

An Introduction to Copulas with Applications

ASSESSING EFFICIENCY OF D-VINE COPULA ARMA-GARCH METHOD IN VALUE AT RISK FORECASTING: EVIDENCE FROM PSE LISTED COMPANIES

Measuring Risk Dependencies in the Solvency II-Framework. Robert Danilo Molinari Tristan Nguyen WHL Graduate School of Business and Economics

Studying Volatility and Dependency of Chinese Outbound Tourism Demand in Singapore, Malaysia, and Thailand: A Vine Copula Approach

3.4 Copula approach for modeling default dependency. Two aspects of modeling the default times of several obligors

Modeling Asymmetry and Tail Dependence among Multiple Variables by Using Partial Regular Vine

Statistical Assessments of Systemic Risk Measures

Financial Risk Management

Maximum likelihood estimation of skew-t copulas with its applications to stock returns

Subject CS1 Actuarial Statistics 1 Core Principles. Syllabus. for the 2019 exams. 1 June 2018

Dependence Structure and Extreme Comovements in International Equity and Bond Markets

Copulas? What copulas? R. Chicheportiche & J.P. Bouchaud, CFM

2. Copula Methods Background

A Multivariate Analysis of Intercompany Loss Triangles

Market Risk Analysis Volume II. Practical Financial Econometrics

JEL Classification: C32, G13, Q11, Q13, Q41, Q48 Keywords: Agricultural commodity; crude oil; dependence; energy policy; sentiment; vine copula

A Joint Credit Scoring Model for Peer-to-Peer Lending and Credit Bureau

Discussion Papers in Economics

Package ESGtoolkit. February 19, 2015

Lindner, Szimayer: A Limit Theorem for Copulas

Dynamic Copula Methods in Finance

DEFAULT PROBABILITY PREDICTION WITH STATIC MERTON-D-VINE COPULA MODEL

ADVANCED OPERATIONAL RISK MODELLING IN BANKS AND INSURANCE COMPANIES

Model selection of vine copulas with applications

PORTFOLIO MODELLING USING THE THEORY OF COPULA IN LATVIAN AND AMERICAN EQUITY MARKET

Stress Testing German Industry Sectors: Results from a Vine Copula Based Quantile Regression

Statistical Models and Methods for Financial Markets

COMISEF WORKING PAPERS SERIES

OPTIMAL PORTFOLIO OF THE GOVERNMENT PENSION INVESTMENT FUND BASED ON THE SYSTEMIC RISK EVALUATED BY A NEW ASYMMETRIC COPULA

Social Networks, Asset Allocation and Portfolio Diversification

Operational risk Dependencies and the Determination of Risk Capital

A flexible and tractable class of one-factor copulas

Extreme Return-Volume Dependence in East-Asian. Stock Markets: A Copula Approach

PORTFOLIO OPTIMIZATION AND SHARPE RATIO BASED ON COPULA APPROACH

Modeling Co-movements and Tail Dependency in the International Stock Market via Copulae

Modeling Crop prices through a Burr distribution and. Analysis of Correlation between Crop Prices and Yields. using a Copula method

I Preliminary Material 1

Approximating a multifactor di usion on a tree.

Somali Ghosh Department of Agricultural Economics Texas A&M University 2124 TAMU College Station, TX

Algebraic Problems in Graphical Modeling

Can Time-Varying Copulas Generate Profit in Portfolio Optimization?

Can expert knowledge compensate for data scarcity in crop insurance pricing? 1

Modeling of Price. Ximing Wu Texas A&M University

Comparative Analyses of Expected Shortfall and Value-at-Risk under Market Stress

Adaptive Interest Rate Modelling

Innovations in Dependence Modelling for Financial Applications

A New Hybrid Estimation Method for the Generalized Pareto Distribution

Dependence Structure between TOURISM and TRANS Sector Indices of the Stock Exchange of Thailand

Simulation of Extreme Events in the Presence of Spatial Dependence

Molecular Phylogenetics

Modelling Dependence between the Equity and. Foreign Exchange Markets Using Copulas

On the Distribution and Its Properties of the Sum of a Normal and a Doubly Truncated Normal

Contents Part I Descriptive Statistics 1 Introduction and Framework Population, Sample, and Observations Variables Quali

Modeling Dependence in the Design of Whole Farm Insurance Contract A Copula-Based Model Approach

Key Words: emerging markets, copulas, tail dependence, Value-at-Risk JEL Classification: C51, C52, C14, G17

Estimation of Travel Time Distribution. and Travel Time Derivatives

2.1 Mathematical Basis: Risk-Neutral Pricing

Multivariate longitudinal data analysis for actuarial applications

Loss Simulation Model Testing and Enhancement

Limit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies

Pricing bivariate option under GARCH processes with time-varying copula

Essays on Statistical Arbitrage. Der Rechts- und Wirtschaftswissenschaftlichen Fakultät/ dem Fachbereich Wirtschaftswissenschafen

Bivariate Birnbaum-Saunders Distribution

Implied Systemic Risk Index (work in progress, still at an early stage)

Market Risk Analysis Volume I

WORLD FINANCIAL RELATIONS: UNDERSTANDING THE CREDIT DERIVATIVE SWAPS (CDS) DEPENDENCE STRUCTURE 1

Pricing Multi-asset Equity Options Driven by a Multidimensional Variance Gamma Process Under Nonlinear Dependence Structures

Economics 2010c: Lecture 4 Precautionary Savings and Liquidity Constraints

Lecture Note 9 of Bus 41914, Spring Multivariate Volatility Models ChicagoBooth

Comparison results for credit risk portfolios

Web-based Supplementary Materials for. A space-time conditional intensity model. for invasive meningococcal disease occurence

A Markov Chain Monte Carlo Approach to Estimate the Risks of Extremely Large Insurance Claims

Copula information criterion for model selection with two-stage maximum likelihood estimation

Analysis of Energy and Agricultural Commodity Markets with the Policy Mandated: A Vine Copula-based ARMA-EGARCH Model

A Copula-GARCH Model of Conditional Dependencies: Estimating Tehran Market Stock. Exchange Value-at-Risk

Handout 4: Deterministic Systems and the Shortest Path Problem

Stochastic Dual Dynamic Programming

Transcription:

Introduction to vine copulas Nicole Krämer & Ulf Schepsmeier Technische Universität München [kraemer, schepsmeier]@ma.tum.de NIPS Workshop, Granada, December 18, 2011 Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 1 / 21

1 Motivation and background 2 Pair-copula construction (PCC) of vine distribution 3 Model selection and estimation 4 Applications and extensions 5 Summary and Outlook Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 2 / 21

0.1 0.2 0.01 0.01 0.2 0.1 0.05 0.2 0.15 0.01 0.05 0.1 Motivation Copulas model marginal and common dependencies separately. There is a wide range of parametric copula families: Gauss Frank Clayton 3 1 0 1 2 3 0.05 0.15 3 1 0 1 2 3 0.15 3 1 0 1 2 3 3 2 1 0 1 2 3 3 2 1 0 1 2 3 3 2 1 0 1 2 3 But: Standard multivariate copulas can become inflexible in high dimensions. do not allow for different dependency structures between pairs of variables. Vine copulas for higher-dimensional data. Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 3 / 21

Overview Vines Vine pair-copulas Bivariate copulas are building blocks for higher-dimensional distributions. The dependency structure is determined by the bivariate copulas and a nested set of trees. Vine approach is more flexible, as we can select bivariate copulas from a wide range of (parametric) families. Model estimation 1 graph theory to determine the dependency structure of the data 2 statistical inference (maximum-likelihood, Bayesian approach...) to fit bivariate copulas. Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 4 / 21

Background - Bivariate Copulas Bivariate Copula A bivariate copula function C : [0, 1] 2 R is a distribution on [0, 1] 2 with uniform marginals. Let F be a bivariate distribution with marginal distributions F 1, F 2. Sklar s Theorem (1959) There exists a two dimensional copula C(, ), such that (x 1, x 2 ) R 2 : F (x 1, x 2 ) = C(F 1 (x 1 ), F 2 (x 2 )). If F 1 and F 2 are continuous, the copula C is unique. Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 5 / 21

Copula densities Copula density (2-dimensional) This implies joint density conditional density c 12 (u 1, u 2 ) = 2 C 12 (u 1, u 2 ) u 1 u 2 f (x 1, x 2 ) = c 12 (F 1 (x 1 ), F 2 (x 2 )) f 1 (x 1 ) f 2 (x 2 ) f (x 2 x 1 ) = c 12 (F 1 (x 1 ), F 2 (x 2 )) f 2 (x 2 ) Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 6 / 21

Important: pair-copula constructions We can represent a density f (x 1,..., x d ) as a product of pair copula densities and marginal densities! Example: d = 3 dimensions. One possible decomposition of f (x 1, x 2, x 3 ) is: f (x 1, x 2, x 3 ) = f 3 12 (x 3 x 1, x 2 )f 2 1 (x 2 x 1 )f 1 (x 1 ) f 2 1 (x 2 x 1 ) = c 12 (F 1 (x 1 ), F 2 (x 2 ))f 2 (x 2 ) f 3 12 (x 3 x 1, x 2 ) = c 13 2 (F 1 2 (x 1 x 2 ), F 3 2 (x 3 x 2 ))f 3 2 (x 3 x 2 ) f 3 2 (x 3 x 2 ) = c 23 (F 2 (x 2 ), F 3 (x 3 ))f 3 (x 3 ) f (x 1, x 2, x 3 ) = f 3 (x 3 )f 2 (x 2 )f 1 (x 1 ) (marginals) c 12 (F 1 (x 1 ), F 2 (x 2 )) c 23 (F 2 (x 2 ), F 3 (x 3 )) (unconditional pairs) c 13 2 (F 1 2 (x 1 x 2 ), F 3 2 (x 3 x 2 ))(conditional pair) Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 7 / 21

Important: pair-copula constructions We can represent a density f (x 1,..., x d ) as a product of pair copula densities and marginal densities! Example: d = 3 dimensions. One possible decomposition of f (x 1, x 2, x 3 ) is: f (x 1, x 2, x 3 ) = f 3 12 (x 3 x 1, x 2 )f 2 1 (x 2 x 1 )f 1 (x 1 ) f 2 1 (x 2 x 1 ) = c 12 (F 1 (x 1 ), F 2 (x 2 ))f 2 (x 2 ) f 3 12 (x 3 x 1, x 2 ) = c 13 2 (F 1 2 (x 1 x 2 ), F 3 2 (x 3 x 2 ))f 3 2 (x 3 x 2 ) f 3 2 (x 3 x 2 ) = c 23 (F 2 (x 2 ), F 3 (x 3 ))f 3 (x 3 ) f (x 1, x 2, x 3 ) = f 3 (x 3 )f 2 (x 2 )f 1 (x 1 ) (marginals) c 12 (F 1 (x 1 ), F 2 (x 2 )) c 23 (F 2 (x 2 ), F 3 (x 3 )) (unconditional pairs) c 13 2 (F 1 2 (x 1 x 2 ), F 3 2 (x 3 x 2 ))(conditional pair) Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 7 / 21

Important: pair-copula constructions We can represent a density f (x 1,..., x d ) as a product of pair copula densities and marginal densities! Example: d = 3 dimensions. One possible decomposition of f (x 1, x 2, x 3 ) is: f (x 1, x 2, x 3 ) = f 3 12 (x 3 x 1, x 2 )f 2 1 (x 2 x 1 )f 1 (x 1 ) f 2 1 (x 2 x 1 ) = c 12 (F 1 (x 1 ), F 2 (x 2 ))f 2 (x 2 ) f 3 12 (x 3 x 1, x 2 ) = c 13 2 (F 1 2 (x 1 x 2 ), F 3 2 (x 3 x 2 ))f 3 2 (x 3 x 2 ) f 3 2 (x 3 x 2 ) = c 23 (F 2 (x 2 ), F 3 (x 3 ))f 3 (x 3 ) f (x 1, x 2, x 3 ) = f 3 (x 3 )f 2 (x 2 )f 1 (x 1 ) (marginals) c 12 (F 1 (x 1 ), F 2 (x 2 )) c 23 (F 2 (x 2 ), F 3 (x 3 )) (unconditional pairs) c 13 2 (F 1 2 (x 1 x 2 ), F 3 2 (x 3 x 2 ))(conditional pair) Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 7 / 21

Important: pair-copula constructions We can represent a density f (x 1,..., x d ) as a product of pair copula densities and marginal densities! Example: d = 3 dimensions. One possible decomposition of f (x 1, x 2, x 3 ) is: f (x 1, x 2, x 3 ) = f 3 12 (x 3 x 1, x 2 )f 2 1 (x 2 x 1 )f 1 (x 1 ) f 2 1 (x 2 x 1 ) = c 12 (F 1 (x 1 ), F 2 (x 2 ))f 2 (x 2 ) f 3 12 (x 3 x 1, x 2 ) = c 13 2 (F 1 2 (x 1 x 2 ), F 3 2 (x 3 x 2 ))f 3 2 (x 3 x 2 ) f 3 2 (x 3 x 2 ) = c 23 (F 2 (x 2 ), F 3 (x 3 ))f 3 (x 3 ) f (x 1, x 2, x 3 ) = f 3 (x 3 )f 2 (x 2 )f 1 (x 1 ) (marginals) c 12 (F 1 (x 1 ), F 2 (x 2 )) c 23 (F 2 (x 2 ), F 3 (x 3 )) (unconditional pairs) c 13 2 (F 1 2 (x 1 x 2 ), F 3 2 (x 3 x 2 ))(conditional pair) Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 7 / 21

Pair-copula construction (PCC) in d dimensions Joe (1996), Bedford and Cooke (2001),Aas et al. (2009), Czado (2010) f (x 1,..., x d ) = d 1 d j c i,(i+j) (i+1),,(i+j 1) j=1 i=1 } {{ } pair copula densities d f k (x k ) k=1 }{{} marginal densities with c i,j i1,,i k := c i,j i1,,i k (F (x i x i1,, x ik ), (F (x j x i1,, x ik )) for i, j, i 1,, i k with i < j and i 1 < < i k. Remarks: The decomposition is not unique. Bedford and Cooke (2001) introduced a graphical structure called regular vine structure to help organize them. Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 8 / 21

Important: regular vine structure Example: d = 3 dimensions f (x 1, x 2, x 3 ) = f 3 (x 3 )f 2 (x 2 )f 1 (x 1 ) (marginals) c 12 (F 1 (x 1 ), F 2 (x 2 )) c 23 (F 2 (x 2 ), F 3 (x 3 )) (unconditional pairs) c 13 2 (F 1 2 (x 1 x 2 ), F 3 2 (x 3 x 2 ))(conditional pair) Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 9 / 21

Important: regular vine structure Example: d = 3 dimensions f (x 1, x 2, x 3 ) = f 3 (x 3 )f 2 (x 2 )f 1 (x 1 ) (marginals) c 12 (F 1 (x 1 ), F 2 (x 2 )) c 23 (F 2 (x 2 ), F 3 (x 3 )) (unconditional pairs) c 13 2 (F 1 2 (x 1 x 2 ), F 3 2 (x 3 x 2 ))(conditional pair) 1 2 3 Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 9 / 21

Important: regular vine structure Example: d = 3 dimensions f (x 1, x 2, x 3 ) = f 3 (x 3 )f 2 (x 2 )f 1 (x 1 ) (marginals) c 12 (F 1 (x 1 ), F 2 (x 2 )) c 23 (F 2 (x 2 ), F 3 (x 3 )) (unconditional pairs) c 13 2 (F 1 2 (x 1 x 2 ), F 3 2 (x 3 x 2 ))(conditional pair) 12 23 1 2 3 Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 9 / 21

Important: regular vine structure Example: d = 3 dimensions f (x 1, x 2, x 3 ) = f 3 (x 3 )f 2 (x 2 )f 1 (x 1 ) (marginals) c 12 (F 1 (x 1 ), F 2 (x 2 )) c 23 (F 2 (x 2 ), F 3 (x 3 )) (unconditional pairs) c 13 2 (F 1 2 (x 1 x 2 ), F 3 2 (x 3 x 2 ))(conditional pair) 12 23 1 2 3 12 23 Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 9 / 21

Important: regular vine structure Example: d = 3 dimensions f (x 1, x 2, x 3 ) = f 3 (x 3 )f 2 (x 2 )f 1 (x 1 ) (marginals) c 12 (F 1 (x 1 ), F 2 (x 2 )) c 23 (F 2 (x 2 ), F 3 (x 3 )) (unconditional pairs) c 13 2 (F 1 2 (x 1 x 2 ), F 3 2 (x 3 x 2 ))(conditional pair) 12 23 1 2 3 13 2 12 23 Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 9 / 21

R-vine structure (d = 5) formal definition 1 3 4 2 5 T1 Pair-copulas: 1 c 12, c 13, c 34, c 34, c 15 2 proximity condition If two nodes are joined by an edge in tree j + 1, the corresponding edges in tree j share a node. 3 c 23 1, c 14 3, c 35 1 4 c 24 13, c 45 13 5 c 25 134 Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 10 / 21

R-vine structure (d = 5) formal definition 1 3 4 1,3 3,4 1,2 1,5 2 5 T1 Pair-copulas: 1 c 12, c 13, c 34, c 34, c 15 2 proximity condition If two nodes are joined by an edge in tree j + 1, the corresponding edges in tree j share a node. 3 c 23 1, c 14 3, c 35 1 4 c 24 13, c 45 13 5 c 25 134 Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 10 / 21

R-vine structure (d = 5) formal definition 1 3 4 1,3 3,4 1,2 1,5 2 5 1,2 1,3 3,4 1,5 T1 T2 Pair-copulas: 1 c 12, c 13, c 34, c 34, c 15 2 proximity condition If two nodes are joined by an edge in tree j + 1, the corresponding edges in tree j share a node. 3 c 23 1, c 14 3, c 35 1 4 c 24 13, c 45 13 5 c 25 134 Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 10 / 21

R-vine structure (d = 5) formal definition 1 3 4 1,3 3,4 1,2 1,5 2 5 1,2 1,3 3,4 1,5 T1 T2 Pair-copulas: 1 c 12, c 13, c 34, c 34, c 15 2 proximity condition If two nodes are joined by an edge in tree j + 1, the corresponding edges in tree j share a node. 3 c 23 1, c 14 3, c 35 1 4 c 24 13, c 45 13 5 c 25 134 Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 10 / 21

R-vine structure (d = 5) formal definition 1 3 4 1,3 3,4 1,2 T1 1,5 2 5 2,3 1 1,4 3 1,2 1,3 3,4 3,5 1 T2 1,5 2,3 1 1,4 3 3,5 1 T3 Pair-copulas: 1 c 12, c 13, c 34, c 34, c 15 2 proximity condition If two nodes are joined by an edge in tree j + 1, the corresponding edges in tree j share a node. 3 c 23 1, c 14 3, c 35 1 4 c 24 13, c 45 13 5 c 25 134 Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 10 / 21

R-vine structure (d = 5) formal definition 1 3 4 1,3 3,4 1,2 1,5 2 5 2,3 1 1,4 3 1,2 1,3 3,4 3,5 1 T2 1,5 2,4 1,3 4,5 1,3 2,3 1 1,4 3 3,5 1 T1 T3 Pair-copulas: 1 c 12, c 13, c 34, c 34, c 15 2 proximity condition If two nodes are joined by an edge in tree j + 1, the corresponding edges in tree j share a node. 3 c 23 1, c 14 3, c 35 1 2,5 1,3,4 2,4 1,3 4,5 1,3 T4 4 c 24 13, c 45 13 5 c 25 134 Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 10 / 21

C-anonical vines Each tree has a unique node that is connected to all other nodes. f 1234 = f 1 f 2 f 3 f }{{} 4 c 12 c 13 c }{{ 14 } nodes in T 1 edges in T 1 nodes in T 2 c 23 1 c 24 1 }{{} edges in T 2 nodes in T 3 c 34 12 }{{} edge in T 3 12 2 3 13 1 4 14 tree 1 23 1 13 12 14 24 1 tree 2 23 1 24 1 34 12 tree 3 Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 11 / 21

vines Each tree is a path. D-rawable f 1234 = f 1 f 2 f 3 f }{{} 4 c 12 c 23 c }{{ 34 } nodes in T 1 edges in T 1 nodes in T 2 c 13 2 c 24 3 }{{} edges in T 2 nodes in T 3 c 14 23 }{{} edge in T 3 1 2 3 4 12 23 34 tree 1 12 23 34 13 2 24 3 tree 2 13 2 24 3 14 23 tree 3 Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 12 / 21

Preliminary summary: pair-copula decomposition So far Given a d-dimensional density, we can decompose it into products of marginal densities and bivariate copula densities. represent this decomposition with nested set of trees that fulfill a proximity condition. Question Given data, how can we estimate a pair-copula decomposition? Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 13 / 21

Preliminary summary: pair-copula decomposition So far Given a d-dimensional density, we can decompose it into products of marginal densities and bivariate copula densities. represent this decomposition with nested set of trees that fulfill a proximity condition. Question Given data, how can we estimate a pair-copula decomposition? Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 13 / 21

Model selection and parameter estimation Model = structure (trees) + copula families + copula parameters Use our software package CDVine! (Brechmann and Schepsmeier (2011)) Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 14 / 21

Model selection and parameter estimation Model = structure (trees) + copula families + copula parameters 1 2 3 4 Data Use our software package CDVine! (Brechmann and Schepsmeier (2011)) Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 14 / 21

Model selection and parameter estimation Model = structure (trees) + copula families + copula parameters 1 2 3 4 Use our software package CDVine! (Brechmann and Schepsmeier (2011)) Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 14 / 21

Model selection and parameter estimation Model = structure (trees) + copula families + copula parameters Normal 1 2 Clayton Gumbel 3 4 Use our software package CDVine! (Brechmann and Schepsmeier (2011)) Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 14 / 21

Model selection and parameter estimation Model = structure (trees) + copula families + copula parameters Normal,ρ = 0.5 1 2 Clayton,θ = 2.5 Gumbel,θ = 1.7 3 4 Use our software package CDVine! (Brechmann and Schepsmeier (2011)) Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 14 / 21

Model selection and parameter estimation Model = structure (trees) + copula families + copula parameters Normal,ρ = 0.5 1 2 12 13 Clayton,θ = 2.5 Gumbel,θ = 1.7 3 4 14 Use our software package CDVine! (Brechmann and Schepsmeier (2011)) Pseudo observations Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 14 / 21

Model selection and parameter estimation Model = structure (trees) + copula families + copula parameters Problem: Huge number of possible vines d(d 1) 2 pair-copulas structure selection copula selection parameter estimation Use our software package CDVine! (Brechmann and Schepsmeier (2011)) Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 14 / 21

Structure selection Possible edge weights Kendall s τ Spearman s ρ p-values of Goodness-of-Fit tests distances Model selection is done tree by tree via optimal C-vines structure selection (Czado et al. (2011)) Traveling Salesman Problem for D-vines Maximum Spanning Tree for R-vines (Dissmann et al. (2011)) Bayesian approaches (Reversible Jump MCMC) Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 15 / 21

Structure selection Possible edge weights Kendall s τ Spearman s ρ p-values of Goodness-of-Fit tests distances Model selection is done tree by tree via optimal C-vines structure selection (Czado et al. (2011)) Traveling Salesman Problem for D-vines Maximum Spanning Tree for R-vines (Dissmann et al. (2011)) Bayesian approaches (Reversible Jump MCMC) Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 15 / 21

Copula selection Copula selection can be done via Goodness-of-fit tests Independence test AIC/BIC graphical tools like contour plots, λ-function,... Possible copula families Elliptical copulas (Gauss, t-) one-parametric Archimedean copulas (Clayton, Gumbel, Frank, Joe,...) two-parametric Archimedean copulas (BB1, BB7,...) rotated versions of the Archimedean for neg. dependencies... Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 16 / 21

Parameter estimation Estimation approaches: Maximum likelihood estimation Sequential estimation: Parameters are estimated sequentially starting from the top tree. Parameter estimates can be used to define pseudo observations for the next tree Parameter estimation via θ = f (τ) or bivariate MLE Sequential estimates can be used as starting values for maximum likelihood Bayesian estimation Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 17 / 21

Parameter estimation Estimation approaches: Maximum likelihood estimation Sequential estimation: Parameters are estimated sequentially starting from the top tree. Parameter estimates can be used to define pseudo observations for the next tree Parameter estimation via θ = f (τ) or bivariate MLE Sequential estimates can be used as starting values for maximum likelihood Bayesian estimation Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 17 / 21

Parameter estimation Estimation approaches: Maximum likelihood estimation Sequential estimation: Parameters are estimated sequentially starting from the top tree. Parameter estimates can be used to define pseudo observations for the next tree Parameter estimation via θ = f (τ) or bivariate MLE Sequential estimates can be used as starting values for maximum likelihood Bayesian estimation Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 17 / 21

Applications Dimensionality of applications Gaussian vines in arbitrary dimensions (Kurowicka and Cooke 2006) First non Gaussian D-vine models using joint maximum likelihood in 4 dimensions Bayesian D-vines with credible intervals in 7 and 12 dimensions Joint maximum likelihood now feasible in 50 dimensions for R-vines Sequential estimation of R-vines in 100 dimensions Sequential estimation for d 100 dimensions with truncation (i.e. higher order trees only contain independent copulas) Heinen and Valdesogo (2009) sequentially fit a C-vine autoregressive model in 100 dimensions Application areas: finance health insurance images genetics... Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 18 / 21

Extensions (Projects of our research group) Special vine models: vine copulas with time varying parameters regime switching vine models non parametric vine pair copulas Non Gaussian directed acyclic graphical (DAG) models based on PCC s discrete vine copulas truncated and simplified R-vines spatial vines copula discriminant analysis Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 19 / 21

Summary and outlook PCC s such as C-, D- and R-vines allow for very flexible class of multivariate distributions Efficient parameter estimation methods are available for dimensions up to 50 Model selection of vine tree structures and pair copula types for regular vines still needs further work Efficient distance measures between vine distributions would be useful Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 20 / 21

Aas, K., C. Czado, A. Frigessi, and H. Bakken (2009). Pair-copula constructions of multiple dependence. Insurance, Mathematics and Economics 44, 182 198. Bedford, T. and R. M. Cooke (2001). Probability density decomposition for conditionally dependent random variables modeled by vines. Annals of Mathematics and Artificial Intelligence 32, 245 268. Brechmann, E. C. and U. Schepsmeier (2011). Dependence modeling with C- and D-vine copulas: The R-package CDVine. Submitted for publication. Czado, C. (2010). Pair-copula constructions of multivariate copulas. In F. Durante, W. Härdle, P. Jaworki, and T. Rychlik (Eds.), Workshop on Copula Theory and its Applications. Springer, Dortrech. Czado, C., U. Schepsmeier, and A. Min (2011). Maximum likelihood estimation of mixed c-vine pair copula with application to exchange rates. to appear in Statistical Modeling. Joe, H. (1996). Families of m-variate distributions with given margins and m(m-1)/2 bivariate dependence parameters. In L. Rüschendorf and B. Schweizer and M. D. Taylor (Ed.), Distributions with Fixed Marginals and Related Topics. Kurowicka, D. and R. Cooke (2006). Uncertainty analysis with high dimensional dependence modelling. Chichester: Wiley. Reading material, software and current projects: http://www-m4.ma.tum.de/en/research/vine-copula-models Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 21 / 21

Regular vine distribution An d-dimensional regular vine is a sequence of d-1 trees 1 tree 1 d nodes: X 1,..., X d d 1 edges: pair-copula densities between nodes X 1,..., X d 2 tree j d + 1 j nodes: edges of tree j 1 d j edges: conditional pair-copula densities Proximity condition: If two nodes in tree j + 1 are joined by an edge, the corresponding edges in tree j share a node. back to talk Krämer & Schepsmeier (TUM) Introduction to vine copulas December 18, 2011 21 / 21