Multivariate Chain Ladder

Size: px
Start display at page:

Download "Multivariate Chain Ladder"

Transcription

1 Multivariate Chain Ladder Carsten Pröhl and Klaus D Schmidt Lehrstuhl für Versicherungsmathematik Technische Universität Dresden schmidt@mathtu-dresdende Multivariate Chain Ladder p1/44

2 Multivariate Chain Ladder 1 The Univariate Chain Ladder Method Multivariate Chain Ladder p2/44

3 Multivariate Chain Ladder 1 The Univariate Chain Ladder Method 2 The Univariate Model of Schnaus Multivariate Chain Ladder p2/44

4 Multivariate Chain Ladder 1 The Univariate Chain Ladder Method 2 The Univariate Model of Schnaus 3 The Multivariate Model Multivariate Chain Ladder p2/44

5 Multivariate Chain Ladder 1 The Univariate Chain Ladder Method 2 The Univariate Model of Schnaus 3 The Multivariate Model 4 The Multivariate Chain Ladder Method Multivariate Chain Ladder p2/44

6 Multivariate Chain Ladder 1 The Univariate Chain Ladder Method 2 The Univariate Model of Schnaus 3 The Multivariate Model 4 The Multivariate Chain Ladder Method 5 Conclusion Multivariate Chain Ladder p2/44

7 1 The Univariate Chain Ladder Method Multivariate Chain Ladder p3/44

8 An Example Accident Development Year k Year i Multivariate Chain Ladder p4/44

9 An Example Accident Development Year k Year i Multivariate Chain Ladder p4/44

10 An Example Accident Development Year k Year i Multivariate Chain Ladder p4/44

11 An Example Accident Development Year k Year i Multivariate Chain Ladder p4/44

12 An Example Accident Development Year k Year i Multivariate Chain Ladder p4/44

13 An Example Accident Development Year k Year i Multivariate Chain Ladder p4/44

14 An Example Accident Development Year k Year i Multivariate Chain Ladder p4/44

15 An Example Accident Development Year k Year i , Multivariate Chain Ladder p4/44

16 An Example Accident Development Year k Year i , Multivariate Chain Ladder p4/44

17 An Example Accident Development Year k Year i , Multivariate Chain Ladder p4/44

18 An Example Accident Development Year k Year i ,232 1, Multivariate Chain Ladder p4/44

19 An Example Accident Development Year k Year i ,232 1, Multivariate Chain Ladder p4/44

20 An Example Accident Development Year k Year i ,232 1, Multivariate Chain Ladder p4/44

21 An Example Accident Development Year k Year i Multivariate Chain Ladder p4/44

22 Abstract Run Off Triangle Run off triangle of observable aggregate claims: Accident Development Year Year 0 1 k n i n 1 n 0 S 0,0 S 0,1 S 0,k S 0,n i S 0,n 1 S 0,n 1 S 1,0 S 1,1 S 1,k S 1,n i S 1,n 1 i S i,0 S i,1 S i,k S i,n i n k S n k,0 S n k,1 S n k,k n 1 S n 1,0 S n 1,1 n S n,0 Multivariate Chain Ladder p5/44

23 Chain Ladder Estimation Individual Development Factors: F i,k := S i,k S i,k 1 Multivariate Chain Ladder p6/44

24 Chain Ladder Estimation Individual Development Factors: Chain Ladder Factors: F i,k := S i,k S i,k 1 F CL k := = n k j=0 S j,k n k j=0 S j,k 1 n k j=0 S j,k 1 n k h=0 S h,k 1 F j,k (for k 1) Multivariate Chain Ladder p6/44

25 Chain Ladder Estimation Identity for aggregate claims: S i,k = S i,n i k l=n i+1 F i,l Multivariate Chain Ladder p7/44

26 Chain Ladder Estimation Identity for aggregate claims: S i,k = S i,n i k l=n i+1 F i,l Chain Ladder Estimators: Ŝ CL i,k := S i,n i k k=n i+1 F CL l (for k n i) Multivariate Chain Ladder p7/44

27 Chain Ladder Estimation: Chain Ladder estimation is recursive since Ŝ CL i,n i = S i,n i and Ŝ CL i,k = ŜCL i,k 1 F CL k for k = n i + 1,, n Multivariate Chain Ladder p8/44

28 Modified Chain Ladder Estimation Chain Ladder Factors: F CL k = n k j=0 S j,k 1 n k h=0 S h,k 1 F j,k Multivariate Chain Ladder p9/44

29 Modified Chain Ladder Estimation Chain Ladder Factors: F CL k = n k j=0 S j,k 1 n k h=0 S h,k 1 F j,k Modified Chain Ladder Factors: F k := n k j=0 W j,k F j,k with n k j=0 W j,k = 1 (for k 1) Multivariate Chain Ladder p9/44

30 Modified Chain Ladder Estimation Chain Ladder Estimators: Ŝ CL i,k := S i,n i k k=n i+1 F CL l Multivariate Chain Ladder p10/44

31 Modified Chain Ladder Estimation Chain Ladder Estimators: Ŝ CL i,k := S i,n i k k=n i+1 F CL l Modified Chain Ladder Estimators: Ŝ i,k := S i,n i k k=n i+1 F l = Competition!!! (for k n i) Multivariate Chain Ladder p10/44

32 2 The Univariate Model of Schnaus Multivariate Chain Ladder p11/44

33 Truncated Run Off Square Accident Development Year Year 0 1 k 1 k 0 S 0,0 S 0,1 S 0,k 1 S 0,k 1 S 1,0 S 1,1 S 1,k 1 S 1,k n k S n k,0 S n k,1 S n k,k 1 S n k,k n k + 1 S n k+1,0 S n k+1,1 S n k+1,k 1 S n k+1,k n 1 S n 1,0 S n 1,1 S n 1,k 1 S n 1,k n S n,0 S n 1,1 S n,k 1 S n,k We denote by G k 1 the information provided by development years 0, 1,, k 1 Multivariate Chain Ladder p12/44

34 The Model of Schnaus Assumption: For each k = 1,, n there exist random variables F k and V k such that E(S i,k G k 1 ) = S i,k 1 F k holds for all i = 0, 1,, n Multivariate Chain Ladder p13/44

35 The Model of Schnaus Assumption: For each k = 1,, n there exist random variables F k and V k such that E(S i,k G k 1 ) = S i,k 1 F k var(s i,k G k 1 ) = S i,k 1 V k holds for all i = 0, 1,, n Multivariate Chain Ladder p13/44

36 The Model of Schnaus Assumption: For each k = 1,, n there exist random variables F k and V k such that E(S i,k G k 1 ) = S i,k 1 F k var(s i,k G k 1 ) = S i,k 1 V k cov(s i,k, S j,k G k 1 ) = 0 holds for all i, j = 0, 1,, n such that j i Multivariate Chain Ladder p13/44

37 The Model of Schnaus Assumption: For each k = 1,, n there exist random variables F k and V k such that E(S i,k G k 1 ) = S i,k 1 F k var(s i,k G k 1 ) = S i,k 1 V k cov(s i,k, S j,k G k 1 ) = 0 holds for all i, j = 0, 1,, n such that j i A slightly less general model is the Model of Mack Multivariate Chain Ladder p14/44

38 The Model of Schnaus Equivalent formulation of the assumption: For each k = 1,, n there exist random variables F k and V k such that E(F i,k G k 1 ) = F k var(f i,k G k 1 ) = V k /S i,k 1 cov(f i,k, F j,k G k 1 ) = 0 holds for all i, j = 0, 1,, n such that j i Multivariate Chain Ladder p15/44

39 Estimation of a non observable F i,k (i+k > n): Admissible estimators: F k := n k j=0 W j,k F j,k with random variables W 0,k, W 1,k,, W n k,k G k 1 satisfying n k j=0 W j,k = 1 Multivariate Chain Ladder p16/44

40 Estimation of a non observable F i,k (i+k > n): The Chain Ladder Factor F CL k := n k j=0 S j,k 1 n k h=0 S h,k 1 F j,k is an admissible estimator Multivariate Chain Ladder p17/44

41 Estimation of a non observable F i,k (i+k > n): The Chain Ladder Factor F CL k := n k j=0 S j,k 1 n k h=0 S h,k 1 F j,k is an admissible estimator Every admissible estimator is a Modified Chain Ladder Factor Multivariate Chain Ladder p17/44

42 Estimation of a non observable F i,k (i+k > n): The Chain Ladder Factor F CL k := n k j=0 S j,k 1 n k h=0 S h,k 1 F j,k is an admissible estimator Every admissible estimator is a Modified Chain Ladder Factor Every admissible estimator is unbiased in the sense that ( ) ) E Fk G k 1 = E (F i,k G k 1 Multivariate Chain Ladder p17/44

43 Estimation of a non observable F i,k (i+k > n): Optimization problem: Minimize E( ( Fk F i,k over all admissible estimators ) 2 G k 1 ) Multivariate Chain Ladder p18/44

44 Estimation of a non observable F i,k (i+k > n): Optimization problem: Minimize ) 2 G k 1 ) E( ( Fk F i,k over all admissible estimators Solution: The optimization problem has the unique solution F CL k := n k j=0 S j,k 1 n k h=0 S h,k 1 F j,k Multivariate Chain Ladder p18/44

45 Estimation of S i,n i+1 : Admissible estimators: Ŝ i,n i+1 := S i,n i Fn i+1 with an admissible estimator F n i+1 of F i,n i+1 Multivariate Chain Ladder p19/44

46 Estimation of S i,n i+1 : The Chain Ladder Estimator Ŝ CL i,n i+1 := S i,n i F CL n i+1 is an admissible estimator Multivariate Chain Ladder p20/44

47 Estimation of S i,n i+1 : The Chain Ladder Estimator Ŝ CL i,n i+1 := S i,n i F CL n i+1 is an admissible estimator Every admissible estimator is a Modified Chain Ladder Estimator Multivariate Chain Ladder p20/44

48 Estimation of S i,n i+1 : The Chain Ladder Estimator Ŝ CL i,n i+1 := S i,n i F CL n i+1 is an admissible estimator Every admissible estimator is a Modified Chain Ladder Estimator Every admissible estimator is unbiased in the sense that ) ) E (Ŝi,n i+1 G n i = E (S i,n i+1 G n i Multivariate Chain Ladder p20/44

49 Estimation of S i,n i+1 : Optimization problem: Minimize E( (Ŝi,n i+1 S i,n i+1 over all admissible estimators ) 2 G n i ) Multivariate Chain Ladder p21/44

50 Estimation of S i,n i+1 : Optimization problem: Minimize ) 2 G n i ) E( (Ŝi,n i+1 S i,n i+1 over all admissible estimators Solution: The optimization problem has the unique solution Ŝ CL i,n i+1 := S i,n i F CL n i+1 = Chain Ladder is the Winner!!! Multivariate Chain Ladder p21/44

51 Estimation of S i,n i+1 : Proof: E( (Ŝi,n i+1 S i,n i+1 ) 2 G n i ) = E( ( S i,n i F n i+1 S i,n i F i,n i+1 = S 2 i,n i E( ( Fn i+1 F i,n i+1 ) 2 G n i ) ) 2 G n i ) Multivariate Chain Ladder p22/44

52 Sequential Optimality: As a consequence of the previous result, the chain ladder estimators are optimal in a sequential sense Multivariate Chain Ladder p23/44

53 3 The Multivariate Model Multivariate Chain Ladder p24/44

54 Correlated Lines Of Business In the present section we consider the loss development for m = 2 correlated lines of business Multivariate Chain Ladder p25/44

55 Correlated Lines Of Business In the present section we consider the loss development for m = 2 correlated lines of business All definitions and results can be extended to arbitrary m Multivariate Chain Ladder p25/44

56 Run Off Triangle for Portfolio (1) Accident Development Year Year 0 1 k n i n 1 n 0 S (1) 0,0 S (1) 0,1 S (1) 0,k S (1) 0,n i S (1) 0,n 1 S (1) 0,n 1 S (1) 1,0 S (1) 1,1 S (1) 1,k S (1) 1,n i S (1) 1,n 1 i S (1) i,0 S (1) i,1 S (1) n k S (1) n k,0 S (1) n 1 S (1) n 1,0 S (1) n 1,1 n S (1) n,0 n k,1 S (1) i,n i i,k S (1) n k,k Multivariate Chain Ladder p26/44

57 Run Off Triangle for Portfolio (2) Accident Development Year Year 0 1 k n i n 1 n 0 S (2) 0,0 S (2) 0,1 S (2) 0,k S (2) 0,n i S (2) 0,n 1 S (2) 0,n 1 S (2) 1,0 S (2) 1,1 S (2) 1,k S (2) 1,n i S (2) 1,n 1 i S (2) i,0 S (2) i,1 S (2) n k S (2) n k,0 S (2) n 1 S (2) n 1,0 S (2) n 1,1 n S (2) n,0 n k,1 S (2) i,n i i,k S (2) n k,k Multivariate Chain Ladder p27/44

58 Run Off Triangle for Both Portfolios Accident Development Year Year 0 1 k n i n 1 n 0 S 0,0 S 0,1 S 0,k S 0,n i S 0,n 1 S 0,n 1 S 1,0 S 1,1 S 1,k S 1,n i S 1,n 1 i S i,0 S i,1 S i,k S i,n i n k S n k,0 S n k,1 S n k,k n 1 S n 1,0 S n 1,1 n S n,0 Multivariate Chain Ladder p28/44

59 Notation S i,k := T i,k := F i,k := S(1) i,k S (2) i,k S(1) i,k 0 0 S (2) i,k F (1) i,k F (2) i,k Multivariate Chain Ladder p29/44

60 Truncated Run Off Square Accident Development Year Year 0 1 k 1 k 0 S 0,0 S 0,1 S 0,k 1 S 0,k 1 S 1,0 S 1,1 S 1,k 1 S 1,k n k S n k,0 S n k,1 S n k,k 1 S n k,k n k + 1 S n k+1,0 S n k+1,1 S n k+1,k 1 S n k+1,k n 1 S n 1,0 S n 1,1 S n 1,k 1 S n 1,k n S n,0 S n 1,1 S n,k 1 S n,k We denote by G k 1 the information provided by development years 0, 1,, k 1 Multivariate Chain Ladder p30/44

61 The Multivariate Model Assumption: For each k = 1,, n there exists a random vector F k and a random matrix V k such that E(S i,k G k 1 ) = S i,k 1 F k holds for all i = 0, 1,, n Multivariate Chain Ladder p31/44

62 The Multivariate Model Assumption: For each k = 1,, n there exists a random vector F k and a random matrix V k such that E(S i,k G k 1 ) = T i,k 1 F k var(s i,k G k 1 ) = T 1/2 i,k 1 V kt 1/2 i,k 1 holds for all i = 0, 1,, n Multivariate Chain Ladder p31/44

63 The Multivariate Model Assumption: For each k = 1,, n there exists a random vector F k and a random matrix V k such that E(S i,k G k 1 ) = T i,k 1 F k var(s i,k G k 1 ) = T 1/2 i,k 1 V kt 1/2 i,k 1 cov(s i,k, S j,k G k 1 ) = O holds for all i, j = 0, 1,, n such that j i Multivariate Chain Ladder p31/44

64 The Multivariate Model Assumption: For each k = 1,, n there exists a random vector F k and a random matrix V k such that E(S i,k G k 1 ) = T i,k 1 F k var(s i,k G k 1 ) = T 1/2 i,k 1 V kt 1/2 i,k 1 cov(s i,k, S j,k G k 1 ) = O holds for all i, j = 0, 1,, n such that j i A slightly less general model is the Model of Braun [2004] Multivariate Chain Ladder p32/44

65 The Multivariate Model Equivalent formulation of the assumption: For each k = 1,, n there exists a random vector F k and a random matrix V k such that E(F i,k G k 1 ) = F k var(f i,k G k 1 ) = T 1/2 i,k 1 V kt 1/2 i,k 1 cov(f i,k, F j,k G k 1 ) = O holds for all i, j = 0, 1,, n such that j i Multivariate Chain Ladder p33/44

66 4 The Multivariate Chain Ladder Method Multivariate Chain Ladder p34/44

67 Estimation of a non observable F i,k (i+k > n): Admissible estimators: F k := n k j=0 W j,k F j,k with random matrices W 0,k, W 1,k,, W n k,k G k 1 satisfying n k j=0 W j,k = I Multivariate Chain Ladder p35/44

68 Estimation of a non observable F i,k (i+k > n): Admissible estimators: F k := n k j=0 W j,k F j,k with random matrices W 0,k, W 1,k,, W n k,k G k 1 satisfying n k j=0 W j,k = I Every admissible estimator is unbiased in the sense that ) ) E ( Fk G n i = E (F i,k G n i Multivariate Chain Ladder p35/44

69 Estimation of a non observable F i,k (i+k > n): Optimization problem: Minimize E( ( Fk F i,k over all admissible estimators ) ( Fk F i,k ) G k 1 ) Multivariate Chain Ladder p36/44

70 Estimation of a non observable F i,k (i+k > n): Solution: The optimization problem has the unique solution F CL k := ( n k j=0 n k j=0 T 1/2 j,k 1 V 1 k T1/2 j,k 1) 1 T 1/2 j,k 1 V 1 k T1/2 j,k 1 F j,k The unique solution F CL k is said to be the Multivariate Chain Ladder Factor of development year k Multivariate Chain Ladder p37/44

71 Estimation of S i,n i+1 : Admissible estimators: Ŝ i,n i+1 := T i,n i Fn i+1 with an admissible estimator F n i+1 of F i,n i+1 Multivariate Chain Ladder p38/44

72 Estimation of S i,n i+1 : Admissible estimators: Ŝ i,n i+1 := T i,n i Fn i+1 with an admissible estimator F n i+1 of F i,n i+1 Every admissible estimator is unbiased in the sense that ) ) E (Ŝi,n i+1 G n i = E (S i,n i+1 G n i Multivariate Chain Ladder p38/44

73 Estimation of S i,n i+1 : Optimization problem: Minimize E( (Ŝi,n i+1 S i,n i+1 over all admissible estimators ) 2 G n i ) Multivariate Chain Ladder p39/44

74 Estimation of S i,n i+1 : Optimization problem: Minimize ) 2 G n i ) E( (Ŝi,n i+1 S i,n i+1 over all admissible estimators Solution: The optimization problem has the unique solution Ŝ CL i,n i+1 := T i,n i FCL n i+1 Multivariate Chain Ladder p39/44

75 Estimation of S i,n i+1 : The unique solution ŜCL k is said to be the Multivariate Chain Ladder Estimator of S i,n i+1 Multivariate Chain Ladder p40/44

76 Recursive Chain Ladder Estimation: Define Ŝ CL i,n i := S i,n i and for k = n i + 1,, n Then ŜCL i,k Estimator of S i,k T CL i,k 1 Ŝ CL i,k := := diag(ŝcl i,k 1) T CL i,k 1 F CL k is said to be the Multivariate Chain Ladder Multivariate Chain Ladder p41/44

77 Sequential Optimality: Sequential optimality obtains as in the univariate case Multivariate Chain Ladder p42/44

78 5 Conclusion Multivariate Chain Ladder p43/44

79 Conclusion In the univariate case, the method is there and a stochastic model had to be found for its justification Multivariate Chain Ladder p44/44

80 Conclusion In the univariate case, the method is there and a stochastic model had to be found for its justification In the multivariate case, an extension of the stochastic model of Schnaus leads to an extension of the method Multivariate Chain Ladder p44/44

81 Conclusion In the univariate case, the method is there and a stochastic model had to be found for its justification In the multivariate case, an extension of the stochastic model of Schnaus leads to an extension of the method The coordinates of the Multivariate Chain Ladder Estimators usually differ from the Univariate Chain Ladder Estimators of the Coordinates Multivariate Chain Ladder p44/44

82 Conclusion In the univariate case, the method is there and a stochastic model had to be found for its justification In the multivariate case, an extension of the stochastic model of Schnaus leads to an extension of the method The coordinates of the Multivariate Chain Ladder Estimators usually differ from the Univariate Chain Ladder Estimators of the Coordinates The multivariate approach resolves the problem of additivity Multivariate Chain Ladder p44/44

A Loss Reserving Method for Incomplete Claim Data

A Loss Reserving Method for Incomplete Claim Data A Loss Reserving Method for Incomplete Claim Data René Dahms Bâloise, Aeschengraben 21, CH-4002 Basel renedahms@baloisech June 11, 2008 Abstract A stochastic model of an additive loss reserving method

More information

Various extensions based on Munich Chain Ladder method

Various extensions based on Munich Chain Ladder method Various extensions based on Munich Chain Ladder method Topic 3 Liability Risk - erve models Jedlicka, etr Charles University, Department of Statistics Sokolovska 83 rague 8 Karlin 180 00 Czech Republic.

More information

Various Extensions Based on Munich Chain Ladder Method

Various 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 information

CHAIN LADDER FORECAST EFFICIENCY

CHAIN LADDER FORECAST EFFICIENCY CHAIN LADDER FORECAST EFFICIENCY Greg Taylor Taylor Fry Consulting Actuaries Level 8, 30 Clarence Street Sydney NSW 2000 Australia Professorial Associate, Centre for Actuarial Studies Faculty of Economics

More information

R = µ + Bf Arbitrage Pricing Model, APM

R = µ + Bf Arbitrage Pricing Model, APM 4.2 Arbitrage Pricing Model, APM Empirical evidence indicates that the CAPM beta does not completely explain the cross section of expected asset returns. This suggests that additional factors may be required.

More information

Claims Reserving under Solvency II

Claims 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 information

Prediction Uncertainty in the Bornhuetter-Ferguson Claims Reserving Method: Revisited

Prediction Uncertainty in the Bornhuetter-Ferguson Claims Reserving Method: Revisited Prediction Uncertainty in the Bornhuetter-Ferguson Claims Reserving Method: Revisited Daniel H. Alai 1 Michael Merz 2 Mario V. Wüthrich 3 Department of Mathematics, ETH Zurich, 8092 Zurich, Switzerland

More information

Stochastic Incremental Approach for Modelling the Claims Reserves

Stochastic Incremental Approach for Modelling the Claims Reserves International Mathematical Forum, Vol. 8, 2013, no. 17, 807-828 HIKARI Ltd, www.m-hikari.com Stochastic Incremental Approach for Modelling the Claims Reserves Ilyes Chorfi Department of Mathematics University

More information

ONE-YEAR AND TOTAL RUN-OFF RESERVE RISK ESTIMATORS BASED ON HISTORICAL ULTIMATE ESTIMATES

ONE-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 information

Ross (1976) introduced the Arbitrage Pricing Theory (APT) as an alternative to the CAPM.

Ross (1976) introduced the Arbitrage Pricing Theory (APT) as an alternative to the CAPM. 4.2 Arbitrage Pricing Model, APM Empirical evidence indicates that the CAPM beta does not completely explain the cross section of expected asset returns. This suggests that additional factors may be required.

More information

Individual loss reserving with the Multivariate Skew Normal framework

Individual loss reserving with the Multivariate Skew Normal framework May 22 2013, K. Antonio, KU Leuven and University of Amsterdam 1 / 31 Individual loss reserving with the Multivariate Skew Normal framework Mathieu Pigeon, Katrien Antonio, Michel Denuit ASTIN Colloquium

More information

Dependence. Practitioner Course: Portfolio Optimization. John Dodson. September 10, Dependence. John Dodson. Outline.

Dependence. 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 information

Generalized Mack Chain-Ladder Model of Reserving with Robust Estimation

Generalized Mack Chain-Ladder Model of Reserving with Robust Estimation Generalized Mac Chain-Ladder Model of Reserving with Robust Estimation Przemyslaw Sloma Abstract n the present paper we consider the problem of stochastic claims reserving in the framewor of Development

More information

Math 140A - Fall Final Exam

Math 140A - Fall Final Exam Math 140A - Fall 2014 - Final Exam Problem 1. Let {a n } n 1 be an increasing sequence of real numbers. (i) If {a n } has a bounded subsequence, show that {a n } is itself bounded. (ii) If {a n } has a

More information

= (, ) V λ (1) λ λ ( + + ) P = [ ( ), (1)] ( ) ( ) = ( ) ( ) ( 0 ) ( 0 ) = ( 0 ) ( 0 ) 0 ( 0 ) ( ( 0 )) ( ( 0 )) = ( ( 0 )) ( ( 0 )) ( + ( 0 )) ( + ( 0 )) = ( + ( 0 )) ( ( 0 )) P V V V V V P V P V V V

More information

Received: 2/7/07, Revised: 5/25/07, Accepted: 6/25/07, Published: 7/20/07 Abstract

Received: 2/7/07, Revised: 5/25/07, Accepted: 6/25/07, Published: 7/20/07 Abstract INTEGERS: ELECTRONIC JOURNAL OF COMBINATORIAL NUMBER THEORY 7 2007, #A34 AN INVERSE OF THE FAÀ DI BRUNO FORMULA Gottlieb Pirsic 1 Johann Radon Institute of Computational and Applied Mathematics RICAM,

More information

GARY G. VENTER 1. CHAIN LADDER VARIANCE FORMULA

GARY 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 information

The Chain Ladder Reserve Uncertainties Revisited

The Chain Ladder Reserve Uncertainties Revisited The Chain Ladder Reserve Uncertainties Revisited by Alois Gisler Paper to be presented at the ASTIN Colloquium 6 in Lisbon Abstract: Chain ladder () is still one of the most popular and most used reserving

More information

4 Derivations of the Discrete-Time Kalman Filter

4 Derivations of the Discrete-Time Kalman Filter Technion Israel Institute of Technology, Department of Electrical Engineering Estimation and Identification in Dynamical Systems (048825) Lecture Notes, Fall 2009, Prof N Shimkin 4 Derivations of the Discrete-Time

More information

Greg Taylor. Taylor Fry Consulting Actuaries University of Melbourne University of New South Wales. CAS Spring meeting Phoenix Arizona May 2012

Greg Taylor. Taylor Fry Consulting Actuaries University of Melbourne University of New South Wales. CAS Spring meeting Phoenix Arizona May 2012 Chain ladder correlations Greg Taylor Taylor Fry Consulting Actuaries University of Melbourne University of New South Wales CAS Spring meeting Phoenix Arizona 0-3 May 0 Overview The chain ladder produces

More information

Invertible Matrices over Idempotent Semirings

Invertible Matrices over Idempotent Semirings Chamchuri Journal of Mathematics Volume 1(2009) Number 2, 55 61 http://www.math.sc.chula.ac.th/cjm Invertible Matrices over Idempotent Semirings W. Mora, A. Wasanawichit and Y. Kemprasit Received 28 Sep

More information

DYNAMIC AND COMPROMISE FACTOR ANALYSIS

DYNAMIC AND COMPROMISE FACTOR ANALYSIS DYNAMIC AND COMPROMISE FACTOR ANALYSIS Marianna Bolla Budapest University of Technology and Economics marib@math.bme.hu Many parts are joint work with Gy. Michaletzky, Loránd Eötvös University and G. Tusnády,

More information

Kalman Filter. Predict: Update: x k k 1 = F k x k 1 k 1 + B k u k P k k 1 = F k P k 1 k 1 F T k + Q

Kalman Filter. Predict: Update: x k k 1 = F k x k 1 k 1 + B k u k P k k 1 = F k P k 1 k 1 F T k + Q Kalman Filter Kalman Filter Predict: x k k 1 = F k x k 1 k 1 + B k u k P k k 1 = F k P k 1 k 1 F T k + Q Update: K = P k k 1 Hk T (H k P k k 1 Hk T + R) 1 x k k = x k k 1 + K(z k H k x k k 1 ) P k k =(I

More information

Conditional Least Squares and Copulae in Claims Reserving for a Single Line of Business

Conditional Least Squares and Copulae in Claims Reserving for a Single Line of Business 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

More information

ScienceDirect. ScienceDirec

ScienceDirect. ScienceDirec CRX_ / CR X_ / Q Q Q Q : R J : / // / J : / / N K * Jk k G U U U U N k U NC U : R k R J R H k - - - - - H - K: R - H - - V V V R - V V V - - L L H H - - C L H j q C L H j q k k k X k R k k X L k k k -

More information

Factor Model Risk Analysis

Factor Model Risk Analysis Factor Model Risk Analysis Eric Zivot University of Washington BlackRock Alternative Advisors April 29, 2011 Outline Factor Model Specification Risk measures Factor Risk Budgeting Portfolio Risk Budgeting

More information

P(X 0 = j 0,... X nk = j k )

P(X 0 = j 0,... X nk = j k ) Introduction to Probability Example Sheet 3 - Michaelmas 2006 Michael Tehranchi Problem. Let (X n ) n 0 be a homogeneous Markov chain on S with transition matrix P. Given a k N, let Z n = X kn. Prove that

More information

8.1 Sequences. Example: A sequence is a function f(n) whose domain is a subset of the integers. Notation: *Note: n = 0 vs. n = 1.

8.1 Sequences. Example: A sequence is a function f(n) whose domain is a subset of the integers. Notation: *Note: n = 0 vs. n = 1. 8. Sequences Example: A sequence is a function f(n) whose domain is a subset of the integers. Notation: *Note: n = 0 vs. n = Examples: 6. Find a formula for the general term a n of the sequence, assuming

More information

Iowa State University. Instructor: Alex Roitershtein Summer Exam #1. Solutions. x u = 2 x v

Iowa State University. Instructor: Alex Roitershtein Summer Exam #1. Solutions. x u = 2 x v Math 501 Iowa State University Introduction to Real Analysis Department of Mathematics Instructor: Alex Roitershtein Summer 015 Exam #1 Solutions This is a take-home examination. The exam includes 8 questions.

More information

Game Theory and its Applications to Networks - Part I: Strict Competition

Game Theory and its Applications to Networks - Part I: Strict Competition Game Theory and its Applications to Networks - Part I: Strict Competition Corinne Touati Master ENS Lyon, Fall 200 What is Game Theory and what is it for? Definition (Roger Myerson, Game Theory, Analysis

More information

Recent Development in Claims Reserving

Recent Development in Claims Reserving WDS'06 roceedings o Contributed apers, art, 118 123, 2006 SBN 80-86732-84-3 MATFZRESS Recent Development in Claims Reserving Jedlička Charles University, Faculty o Mathematics and hysics, rague, Czech

More information

Supplement to the paper Accurate distributions of Mallows C p and its unbiased modifications with applications to shrinkage estimation

Supplement to the paper Accurate distributions of Mallows C p and its unbiased modifications with applications to shrinkage estimation To aear in Economic Review (Otaru University of Commerce Vol. No..- 017. Sulement to the aer Accurate distributions of Mallows C its unbiased modifications with alications to shrinkage estimation Haruhiko

More information

PANEL DISCUSSION: THE ROLE OF POTENTIAL OUTPUT IN POLICYMAKING

PANEL DISCUSSION: THE ROLE OF POTENTIAL OUTPUT IN POLICYMAKING PANEL DISCUSSION: THE ROLE OF POTENTIAL OUTPUT IN POLICYMAKING James Bullard* Federal Reserve Bank of St. Louis 33rd Annual Economic Policy Conference St. Louis, MO October 17, 2008 Views expressed are

More information

Chain ladder with random effects

Chain 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 information

Statistics for Business and Economics

Statistics for Business and Economics Statistics for Business and Economics Chapter 6 Sampling and Sampling Distributions Ch. 6-1 6.1 Tools of Business Statistics n Descriptive statistics n Collecting, presenting, and describing data n Inferential

More information

Similar to sequence, note that a series converges if and only if its tail converges, that is, r 1 r ( 1 < r < 1), ( 1) k k. r k =

Similar to sequence, note that a series converges if and only if its tail converges, that is, r 1 r ( 1 < r < 1), ( 1) k k. r k = Infinite Series We say an infinite series a k converges to s if its sequence of initial sums converges to s, that is, lim( n a k : n N) = s. Similar to sequence, note that a series converges if and only

More information

Conditional Tail Expectations for Multivariate Phase Type Distributions

Conditional Tail Expectations for Multivariate Phase Type Distributions Conditional Tail Expectations for Multivariate Phase Type Distributions Jun Cai Department of Statistics and Actuarial Science University of Waterloo Waterloo, ON N2L 3G1, Canada Telphone: 1-519-8884567,

More information

Sample of Ph.D. Advisory Exam For MathFinance

Sample of Ph.D. Advisory Exam For MathFinance Sample of Ph.D. Advisory Exam For MathFinance Students who wish to enter the Ph.D. program of Mathematics of Finance are required to take the advisory exam. This exam consists of three major parts. The

More information

MULTIVARIATE DISCRETE PHASE-TYPE DISTRIBUTIONS

MULTIVARIATE DISCRETE PHASE-TYPE DISTRIBUTIONS MULTIVARIATE DISCRETE PHASE-TYPE DISTRIBUTIONS By MATTHEW GOFF A dissertation submitted in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY WASHINGTON STATE UNIVERSITY Department

More information

Alternative Characterization of Ergodicity for Doubly Stochastic Chains

Alternative Characterization of Ergodicity for Doubly Stochastic Chains Alternative Characterization of Ergodicity for Doubly Stochastic Chains Behrouz Touri and Angelia Nedić Abstract In this paper we discuss the ergodicity of stochastic and doubly stochastic chains. We define

More information

A HOMOTOPY CLASS OF SEMI-RECURSIVE CHAIN LADDER MODELS

A HOMOTOPY CLASS OF SEMI-RECURSIVE CHAIN LADDER MODELS A HOMOTOPY CLASS OF SEMI-RECURSIVE CHAIN LADDER MODELS Greg Taylor Taylor Fry Consulting Actuaries Level, 55 Clarence Street Sydney NSW 2000 Australia Professorial Associate Centre for Actuarial Studies

More information

Week 15-16: Combinatorial Design

Week 15-16: Combinatorial Design Week 15-16: Combinatorial Design May 8, 2017 A combinatorial design, or simply a design, is an arrangement of the objects of a set into subsets satisfying certain prescribed properties. The area of combinatorial

More information

An idea how to solve some of the problems. diverges the same must hold for the original series. T 1 p T 1 p + 1 p 1 = 1. dt = lim

An idea how to solve some of the problems. diverges the same must hold for the original series. T 1 p T 1 p + 1 p 1 = 1. dt = lim An idea how to solve some of the problems 5.2-2. (a) Does not converge: By multiplying across we get Hence 2k 2k 2 /2 k 2k2 k 2 /2 k 2 /2 2k 2k 2 /2 k. As the series diverges the same must hold for the

More information

Modelling Dependence with Copulas and Applications to Risk Management. Filip Lindskog, RiskLab, ETH Zürich

Modelling Dependence with Copulas and Applications to Risk Management. Filip Lindskog, RiskLab, ETH Zürich Modelling Dependence with Copulas and Applications to Risk Management Filip Lindskog, RiskLab, ETH Zürich 02-07-2000 Home page: http://www.math.ethz.ch/ lindskog E-mail: lindskog@math.ethz.ch RiskLab:

More information

9. Multivariate Linear Time Series (II). MA6622, Ernesto Mordecki, CityU, HK, 2006.

9. Multivariate Linear Time Series (II). MA6622, Ernesto Mordecki, CityU, HK, 2006. 9. Multivariate Linear Time Series (II). MA6622, Ernesto Mordecki, CityU, HK, 2006. References for this Lecture: Introduction to Time Series and Forecasting. P.J. Brockwell and R. A. Davis, Springer Texts

More information

EVALUATING THE BIVARIATE COMPOUND GENERALIZED POISSON DISTRIBUTION

EVALUATING THE BIVARIATE COMPOUND GENERALIZED POISSON DISTRIBUTION An. Şt. Univ. Ovidius Constanţa Vol. 92), 2001, 181 192 EVALUATING THE BIVARIATE COMPOUND GENERALIZED POISSON DISTRIBUTION Raluca Vernic Dedicated to Professor Mirela Ştefănescu on the occasion of her

More information

Probability Models of Information Exchange on Networks Lecture 1

Probability Models of Information Exchange on Networks Lecture 1 Probability Models of Information Exchange on Networks Lecture 1 Elchanan Mossel UC Berkeley All Rights Reserved Motivating Questions How are collective decisions made by: people / computational agents

More information

Non-Life Insurance: Mathematics and Statistics

Non-Life Insurance: Mathematics and Statistics Exercise sheet 1 Exercise 1.1 Discrete Distribution Suppose the random variable N follows a geometric distribution with parameter p œ (0, 1), i.e. ; (1 p) P[N = k] = k 1 p if k œ N \{0}, 0 else. (a) Show

More information

Uncertainty Per Krusell & D. Krueger Lecture Notes Chapter 6

Uncertainty Per Krusell & D. Krueger Lecture Notes Chapter 6 1 Uncertainty Per Krusell & D. Krueger Lecture Notes Chapter 6 1 A Two-Period Example Suppose the economy lasts only two periods, t =0, 1. The uncertainty arises in the income (wage) of period 1. Not that

More information

MATH6142 Complex Networks Exam 2016

MATH6142 Complex Networks Exam 2016 MATH642 Complex Networks Exam 26 Solution Problem a) The network is directed because the adjacency matrix is not symmetric. b) The network is shown in Figure. (4 marks)[unseen] 4 3 Figure : The directed

More information

arxiv:math/ v1 [math.co] 13 Jul 2005

arxiv:math/ v1 [math.co] 13 Jul 2005 A NEW PROOF OF THE REFINED ALTERNATING SIGN MATRIX THEOREM arxiv:math/0507270v1 [math.co] 13 Jul 2005 Ilse Fischer Fakultät für Mathematik, Universität Wien Nordbergstrasse 15, A-1090 Wien, Austria E-mail:

More information

Need for Sampling in Machine Learning. Sargur Srihari

Need for Sampling in Machine Learning. Sargur Srihari Need for Sampling in Machine Learning Sargur srihari@cedar.buffalo.edu 1 Rationale for Sampling 1. ML methods model data with probability distributions E.g., p(x,y; θ) 2. Models are used to answer queries,

More information

Individual reserving a survey. R in Insurance, Paris 08-June-2017 Alexandre Boumezoued, joint work with Laurent Devineau

Individual reserving a survey. R in Insurance, Paris 08-June-2017 Alexandre Boumezoued, joint work with Laurent Devineau Individual reserving a survey R in Insurance, Paris 08-June-2017 Alexandre Boumezoued, joint work with Laurent Devineau Agenda 1 Context 2 Model 3 Results 4 Practical illustration 2 Context From aggregate

More information

A CLASS OF THE BEST LINEAR UNBIASED ESTIMATORS IN THE GENERAL LINEAR MODEL

A CLASS OF THE BEST LINEAR UNBIASED ESTIMATORS IN THE GENERAL LINEAR MODEL A CLASS OF THE BEST LINEAR UNBIASED ESTIMATORS IN THE GENERAL LINEAR MODEL GABRIELA BEGANU The existence of the best linear unbiased estimators (BLUE) for parametric estimable functions in linear models

More information

Variance reduced Value at Risk Monte-Carlo simulations

Variance reduced Value at Risk Monte-Carlo simulations Variance reduced Value at Risk Monte-Carlo simulations Armin Müller 1 Discussion paper No. 500 November 24, 2016 Diskussionsbeiträge der Fakultät für Wirtschaftswissenschaft der FernUniversität in Hagen

More information

Autonomous Navigation for Flying Robots

Autonomous Navigation for Flying Robots Computer Vision Group Prof. Daniel Cremers Autonomous Navigation for Flying Robots Lecture 6.2: Kalman Filter Jürgen Sturm Technische Universität München Motivation Bayes filter is a useful tool for state

More information

A Few Special Distributions and Their Properties

A Few Special Distributions and Their Properties A Few Special Distributions and Their Properties Econ 690 Purdue University Justin L. Tobias (Purdue) Distributional Catalog 1 / 20 Special Distributions and Their Associated Properties 1 Uniform Distribution

More information

A Short Course in Basic Statistics

A Short Course in Basic Statistics A Short Course in Basic Statistics Ian Schindler November 5, 2017 Creative commons license share and share alike BY: C 1 Descriptive Statistics 1.1 Presenting statistical data Definition 1 A statistical

More information

Pubh 8482: Sequential Analysis

Pubh 8482: Sequential Analysis Pubh 8482: Sequential Analysis Joseph S. Koopmeiners Division of Biostatistics University of Minnesota Week 7 Course Summary To this point, we have discussed group sequential testing focusing on Maintaining

More information

STA 4273H: Statistical Machine Learning

STA 4273H: Statistical Machine Learning STA 4273H: Statistical Machine Learning Russ Salakhutdinov Department of Statistics! rsalakhu@utstat.toronto.edu! http://www.utstat.utoronto.ca/~rsalakhu/ Sidney Smith Hall, Room 6002 Lecture 11 Project

More information

arxiv: v2 [math-ph] 24 Feb 2016

arxiv: v2 [math-ph] 24 Feb 2016 ON THE CLASSIFICATION OF MULTIDIMENSIONALLY CONSISTENT 3D MAPS MATTEO PETRERA AND YURI B. SURIS Institut für Mathemat MA 7-2 Technische Universität Berlin Str. des 17. Juni 136 10623 Berlin Germany arxiv:1509.03129v2

More information

Competitive Equilibrium and the Welfare Theorems

Competitive Equilibrium and the Welfare Theorems Competitive Equilibrium and the Welfare Theorems Craig Burnside Duke University September 2010 Craig Burnside (Duke University) Competitive Equilibrium September 2010 1 / 32 Competitive Equilibrium and

More information

Institut für Mathematik

Institut für Mathematik RHEINISCH-WESTFÄLISCHE TECHNISCHE HOCHSCHULE AACHEN Institut für Mathematik A General Framework for Portfolio Theory. Part II: drawdown risk measures by S. Maier-Paape Q. J. Zhu Report No. 92 2017 October

More information

Probability and Statistics

Probability and Statistics Probability and Statistics 1 Contents some stochastic processes Stationary Stochastic Processes 2 4. Some Stochastic Processes 4.1 Bernoulli process 4.2 Binomial process 4.3 Sine wave process 4.4 Random-telegraph

More information

DEN: Linear algebra numerical view (GEM: Gauss elimination method for reducing a full rank matrix to upper-triangular

DEN: Linear algebra numerical view (GEM: Gauss elimination method for reducing a full rank matrix to upper-triangular form) Given: matrix C = (c i,j ) n,m i,j=1 ODE and num math: Linear algebra (N) [lectures] c phabala 2016 DEN: Linear algebra numerical view (GEM: Gauss elimination method for reducing a full rank matrix

More information

2.6 Applying the Trigonometric Ratios

2.6 Applying the Trigonometric Ratios 2.6 Applying the Trigonometric atios FOCUS Use trigonometric ratios to solve a right triangle. When we solve a triangle, we find the measures of all the angles and the lengths of all the sides. To do this

More information

. Get closed expressions for the following subsequences and decide if they converge. (1) a n+1 = (2) a 2n = (3) a 2n+1 = (4) a n 2 = (5) b n+1 =

. Get closed expressions for the following subsequences and decide if they converge. (1) a n+1 = (2) a 2n = (3) a 2n+1 = (4) a n 2 = (5) b n+1 = Math 316, Intro to Analysis subsequences. Recall one of our arguments about why a n = ( 1) n diverges. Consider the subsequences a n = ( 1) n = +1. It converges to 1. On the other hand, the subsequences

More information

arxiv: v1 [math.na] 25 Oct 2018

arxiv: v1 [math.na] 25 Oct 2018 Multi-scale control variate methods for uncertainty quantification in kinetic equations arxiv:80.0844v [math.na] 25 Oct 208 Giacomo Dimarco and Lorenzo Pareschi October 26, 208 Abstract Kinetic equations

More information

The Choice of Representative Volumes for Random Materials

The Choice of Representative Volumes for Random Materials The Choice of Representative Volumes for Random Materials Julian Fischer IST Austria 13.4.2018 Effective properties of random materials Material with random microstructure Effective large-scale description?

More information

Financial Econometrics

Financial Econometrics Financial Econometrics Estimation and Inference Gerald P. Dwyer Trinity College, Dublin January 2013 Who am I? Visiting Professor and BB&T Scholar at Clemson University Federal Reserve Bank of Atlanta

More information

MULTIVARIATE TIME SERIES ANALYSIS AN ADAPTATION OF BOX-JENKINS METHODOLOGY Joseph N Ladalla University of Illinois at Springfield, Springfield, IL

MULTIVARIATE TIME SERIES ANALYSIS AN ADAPTATION OF BOX-JENKINS METHODOLOGY Joseph N Ladalla University of Illinois at Springfield, Springfield, IL MULTIVARIATE TIME SERIES ANALYSIS AN ADAPTATION OF BOX-JENKINS METHODOLOGY Joseph N Ladalla University of Illinois at Springfield, Springfield, IL KEYWORDS: Multivariate time series, Box-Jenkins ARIMA

More information

Chapter 30 Design and Analysis of

Chapter 30 Design and Analysis of Chapter 30 Design and Analysis of 2 k DOEs Introduction This chapter describes design alternatives and analysis techniques for conducting a DOE. Tables M1 to M5 in Appendix E can be used to create test

More information

International Competition in Mathematics for Universtiy Students in Plovdiv, Bulgaria 1994

International Competition in Mathematics for Universtiy Students in Plovdiv, Bulgaria 1994 International Competition in Mathematics for Universtiy Students in Plovdiv, Bulgaria 1994 1 PROBLEMS AND SOLUTIONS First day July 29, 1994 Problem 1. 13 points a Let A be a n n, n 2, symmetric, invertible

More information

Variance Component Models for Quantitative Traits. Biostatistics 666

Variance Component Models for Quantitative Traits. Biostatistics 666 Variance Component Models for Quantitative Traits Biostatistics 666 Today Analysis of quantitative traits Modeling covariance for pairs of individuals estimating heritability Extending the model beyond

More information

Calendar Year Dependence Modeling in Run-Off Triangles

Calendar 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 information

Bootstrapping the triangles

Bootstrapping 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 information

The Quantum Heisenberg Ferromagnet

The Quantum Heisenberg Ferromagnet The Quantum Heisenberg Ferromagnet Soon after Schrödinger discovered the wave equation of quantum mechanics, Heisenberg and Dirac developed the first successful quantum theory of ferromagnetism W. Heisenberg,

More information

9.1 Orthogonal factor model.

9.1 Orthogonal factor model. 36 Chapter 9 Factor Analysis Factor analysis may be viewed as a refinement of the principal component analysis The objective is, like the PC analysis, to describe the relevant variables in study in terms

More information

Mid Term-1 : Solutions to practice problems

Mid Term-1 : Solutions to practice problems Mid Term- : Solutions to practice problems 0 October, 06. Is the function fz = e z x iy holomorphic at z = 0? Give proper justification. Here we are using the notation z = x + iy. Solution: Method-. Use

More information

Optimal Randomized Algorithms for Integration on Function Spaces with underlying ANOVA decomposition

Optimal Randomized Algorithms for Integration on Function Spaces with underlying ANOVA decomposition Optimal Randomized on Function Spaces with underlying ANOVA decomposition Michael Gnewuch 1 University of Kaiserslautern, Germany October 16, 2013 Based on Joint Work with Jan Baldeaux (UTS Sydney) & Josef

More information

An introduction to multivariate data

An introduction to multivariate data An introduction to multivariate data Angela Montanari 1 The data matrix The starting point of any analysis of multivariate data is a data matrix, i.e. a collection of n observations on a set of p characters

More information

Summer Jump-Start Program for Analysis, 2012 Song-Ying Li. 1 Lecture 7: Equicontinuity and Series of functions

Summer Jump-Start Program for Analysis, 2012 Song-Ying Li. 1 Lecture 7: Equicontinuity and Series of functions Summer Jump-Start Program for Analysis, 0 Song-Ying Li Lecture 7: Equicontinuity and Series of functions. Equicontinuity Definition. Let (X, d) be a metric space, K X and K is a compact subset of X. C(K)

More information

A Central Limit Theorem for the Sum of Generalized Linear and Quadratic Forms. by R. L. Eubank and Suojin Wang Texas A&M University

A Central Limit Theorem for the Sum of Generalized Linear and Quadratic Forms. by R. L. Eubank and Suojin Wang Texas A&M University A Central Limit Theorem for the Sum of Generalized Linear and Quadratic Forms by R. L. Eubank and Suojin Wang Texas A&M University ABSTRACT. A central limit theorem is established for the sum of stochastically

More information

1 Introduction The purpose of this paper is to illustrate the typical behavior of learning algorithms using stochastic approximations (SA). In particu

1 Introduction The purpose of this paper is to illustrate the typical behavior of learning algorithms using stochastic approximations (SA). In particu Strong Points of Weak Convergence: A Study Using RPA Gradient Estimation for Automatic Learning Felisa J. Vazquez-Abad * Department of Computer Science and Operations Research University of Montreal, Montreal,

More information

460 HOLGER DETTE AND WILLIAM J STUDDEN order to examine how a given design behaves in the model g` with respect to the D-optimality criterion one uses

460 HOLGER DETTE AND WILLIAM J STUDDEN order to examine how a given design behaves in the model g` with respect to the D-optimality criterion one uses Statistica Sinica 5(1995), 459-473 OPTIMAL DESIGNS FOR POLYNOMIAL REGRESSION WHEN THE DEGREE IS NOT KNOWN Holger Dette and William J Studden Technische Universitat Dresden and Purdue University Abstract:

More information

Supermodular ordering of Poisson arrays

Supermodular ordering of Poisson arrays Supermodular ordering of Poisson arrays Bünyamin Kızıldemir Nicolas Privault Division of Mathematical Sciences School of Physical and Mathematical Sciences Nanyang Technological University 637371 Singapore

More information

The purpose of this section is to derive the asymptotic distribution of the Pearson chi-square statistic. k (n j np j ) 2. np j.

The purpose of this section is to derive the asymptotic distribution of the Pearson chi-square statistic. k (n j np j ) 2. np j. Chapter 9 Pearson s chi-square test 9. Null hypothesis asymptotics Let X, X 2, be independent from a multinomial(, p) distribution, where p is a k-vector with nonnegative entries that sum to one. That

More information

A Probability Review

A Probability Review A Probability Review Outline: A probability review Shorthand notation: RV stands for random variable EE 527, Detection and Estimation Theory, # 0b 1 A Probability Review Reading: Go over handouts 2 5 in

More information

Variations. ECE 6540, Lecture 10 Maximum Likelihood Estimation

Variations. ECE 6540, Lecture 10 Maximum Likelihood Estimation Variations ECE 6540, Lecture 10 Last Time BLUE (Best Linear Unbiased Estimator) Formulation Advantages Disadvantages 2 The BLUE A simplification Assume the estimator is a linear system For a single parameter

More information

Chapter 3. Point Estimation. 3.1 Introduction

Chapter 3. Point Estimation. 3.1 Introduction Chapter 3 Point Estimation Let (Ω, A, P θ ), P θ P = {P θ θ Θ}be probability space, X 1, X 2,..., X n : (Ω, A) (IR k, B k ) random variables (X, B X ) sample space γ : Θ IR k measurable function, i.e.

More information

Algorithms: Background

Algorithms: Background Algorithms: Background Amotz Bar-Noy CUNY Amotz Bar-Noy (CUNY) Algorithms: Background 1 / 66 What is a Proof? Definition I: The cogency of evidence that compels acceptance by the mind of a truth or a fact.

More information

Generalized Linear Models (GLZ)

Generalized Linear Models (GLZ) Generalized Linear Models (GLZ) Generalized Linear Models (GLZ) are an extension of the linear modeling process that allows models to be fit to data that follow probability distributions other than the

More information

Boolean Inner-Product Spaces and Boolean Matrices

Boolean Inner-Product Spaces and Boolean Matrices Boolean Inner-Product Spaces and Boolean Matrices Stan Gudder Department of Mathematics, University of Denver, Denver CO 80208 Frédéric Latrémolière Department of Mathematics, University of Denver, Denver

More information

Extended Kalman Filter Derivation Example application: frequency tracking

Extended Kalman Filter Derivation Example application: frequency tracking Extended Kalman Filter Derivation Example application: frequency tracking Nonlinear State Space Model Consider the nonlinear state-space model of a random process x n+ = f n (x n )+g n (x n )u n y n =

More information

We are going to discuss what it means for a sequence to converge in three stages: First, we define what it means for a sequence to converge to zero

We are going to discuss what it means for a sequence to converge in three stages: First, we define what it means for a sequence to converge to zero Chapter Limits of Sequences Calculus Student: lim s n = 0 means the s n are getting closer and closer to zero but never gets there. Instructor: ARGHHHHH! Exercise. Think of a better response for the instructor.

More information

Summary: A Random Walks View of Spectral Segmentation, by Marina Meila (University of Washington) and Jianbo Shi (Carnegie Mellon University)

Summary: A Random Walks View of Spectral Segmentation, by Marina Meila (University of Washington) and Jianbo Shi (Carnegie Mellon University) Summary: A Random Walks View of Spectral Segmentation, by Marina Meila (University of Washington) and Jianbo Shi (Carnegie Mellon University) The authors explain how the NCut algorithm for graph bisection

More information

Higher order moments of the estimated tangency portfolio weights

Higher order moments of the estimated tangency portfolio weights WORKING PAPER 0/07 Higher order moments of the estimated tangency portfolio weights Farrukh Javed, Stepan Mazur, Edward Ngailo Statistics ISSN 403-0586 http://www.oru.se/institutioner/handelshogskolan-vid-orebro-universitet/forskning/publikationer/working-papers/

More information

Matrix representation of a linear map

Matrix representation of a linear map Matrix representation of a linear map As before, let e i = (0,..., 0, 1, 0,..., 0) T, with 1 in the i th place and 0 elsewhere, be standard basis vectors. Given linear map f : R n R m we get n column vectors

More information

Markov Chains CK eqns Classes Hitting times Rec./trans. Strong Markov Stat. distr. Reversibility * Markov Chains

Markov Chains CK eqns Classes Hitting times Rec./trans. Strong Markov Stat. distr. Reversibility * Markov Chains Markov Chains A random process X is a family {X t : t T } of random variables indexed by some set T. When T = {0, 1, 2,... } one speaks about a discrete-time process, for T = R or T = [0, ) one has a continuous-time

More information

Computational Methods CMSC/AMSC/MAPL 460. Fourier transform

Computational Methods CMSC/AMSC/MAPL 460. Fourier transform Computational Methods CMSC/AMSC/MAPL 460 Fourier transform Ramani Duraiswami, Dept. of Computer Science Several slides from Prof. Healy s course at UMD Fourier Methods Fourier analysis ( harmonic analysis

More information