Prediction Error of the Multivariate Additive Loss Reserving Method for Dependent Lines of Business

Size: px
Start display at page:

Download "Prediction Error of the Multivariate Additive Loss Reserving Method for Dependent Lines of Business"

Transcription

1 Predcton Error of the Mutvarate Addtve Loss Reservng Method for Dependent Lnes of Busness by Mchae Merz and Maro V Wüthrch ABSTRACT Often n non-fe nsurance, cams reserves are the argest poston on the abty sde of the baance sheet Therefore, the predcton of adequate cams reserves for a portfoo consstng of severa run-off subportfoos from dependent nes of busness s of great mportance for every non-fe nsurance company In the present paper, we consder the cams reservng probem n a mutvarate context that s, we study a speca case of the mutvarate addtve oss reservng mode proposed by Hess, Schmdt, and Zocher (26) and Schmdt (26a) Ths mode aows for a smutaneous study of the ndvdua run-off subportfoos and enabes the dervaton of an estmator for the condtona mean square error of predcton (MSEP) for the predctor of the utmate cams of the tota portfoo We ustrate the resuts usng the data gven n Braun (24) and compare them to the resuts derved by the mutvarate chan-adder methods of Braun (24) and Merz and Wüthrch (28) KEYWORDS Cams reservng, sovency, uncertanty, dependent nes of busness, mutvarate addtve oss reservng method, mutvarate chan-adder method, process varance, estmaton error, mean square error of predcton VOLUME 3/ISSUE CASUALTY ACTUARIAL SOCIETY 3

2 Varance Advancng the Scence of Rsk Introducton and motvaton Cams reservng Often n non-fe nsurance, cams reserves are the argest poston on the abty sde of the baance sheet Therefore, gven the avaabe nformaton about the past, the predcton of an adequate amount of cam abty assumed by the non-fe nsurance company, as we as the quantfcaton of the uncertantes n these reserves, s a maor task n actuara practce and scence [eg, Tayor (2); Wüthrch and Merz (28); Casuaty Actuara Socety (2); Teuges and Sundt (24); Engand and Verra (22)] 2 Mutvarate cams reservng methods and ther condtona MSEP In the present paper, we consder the cams reservng probem for a portfoo consstng of severa correated run-off subportfoos Ths smutaneous study of severa ndvdua run-off subportfoos s motvated by the foowng consderatons: ² In practce t s qute natura to subdvde a non-fe run-off portfoo nto severa correated subportfoos, such that each subportfoo satsfes certan homogenety propertes (eg, the chan-adder assumptons or the assumptons of the addtve method) ² It addresses the probem of dependence between the run-off portfoos of dfferent nes of busness (eg, between auto abty and genera abty busness) ² The mutvarate approach has the advantage that by observng one run-off subportfoo we can earn about the behavor of the other runoff subportfoos (eg, subportfoos of sma and arge cams) ² It resoves the probem of addtvty (e, the estmators of the utmate cams for the whoe portfoo are obtaned by summaton over the estmators of the utmate cams for the ndvdua run-off subportfoos) However, n the case of correated run-off subportfoos, the cacuaton of the condtona mean square error of predcton (MSEP) for the predctor of the utmate cam sze of the tota portfoo s more sophstcated than the cacuaton of the condtona MSEP for the predctor of the utmate cam sze of a snge run-off subportfoo An aternatve dea to the smutaneous study of severa ndvdua run-off subportfoos s to cacuate the reserves and ther uncertantes ony for the tota aggregated run-off portfoo However, one shoud pay attenton to the fact that f the subportfoos satsfy, for exampe, the assumptons of the chan-adder or the assumptons of the addtve method, the aggregated run-off portfoo does not n genera satsfy these assumptons (Ane 994; Kemmt 24) Therefore, n most cases t s not a promsng souton to study the aggregated portfoo for the cams reservng probem of severa run-off subportfoos Homberg (994) was probaby the frst one to nvestgate the probem of dependence between run-off portfoos of dfferent nes of busness Later Hawe (997) and Quarg and Mack (24) [see aso Merz and Wüthrch (26)] proposed the frst bvarate modes whch express the dependence between the pad and ncurred osses of a snge run-off subportfoo Braun (24) generazed the we-known unvarate chan-adder mode of Mack (993) to the bvarate case by ncorporatng correatons between two run-off subportfoos In ths setup he derved an estmate for the condtona MSEP for the predctor of the utmate cam sze of two correated run-off subportfoos Usng a mutvarate tme-seres mode for the chan-adder method Merz and Wüthrch (27) gave an estmator for the condtona MSEP n the case of N correated run-off subportfoos However, both the Braun (24) approach and the Merz and Wüthrch (27) approach have the dsadvantage that the chan-adder factors are estmated 32 CASUALTY ACTUARIAL SOCIETY VOLUME 3/ISSUE

3 Predcton Error of the Mutvarate Addtve Loss Reservng Method for Dependent Lnes of Busness n a unvarate way Ths means the estmaton of the chan-adder factors s restrcted to the data of the respectve ndvdua run-off subportfoo and therefore does not take nto account the correaton structure between the dfferent runoff subportfoos Pröh and Schmdt (25) and Schmdt (26a) showed that these unvarate estmates of the chan-adder factors are not optma n terms of a cassca optmaty crteron n the case of correated run-off subportfoos and therefore one shoud repace the unvarate estmators wth mutvarate estmators of the chanadder factors refectng the correaton structure However, ther study dd not go beyond best estmators; that s, they dd not derve an estmator for the condtona MSEP for the predctor of the utmate cam sze of the tota portfoo Fnay, usng a mutvarate chan-adder tmeseres mode, Merz and Wüthrch (28) derved an estmate for the condtona MSEP, n whch the chan-adder factors are estmated n a mutvarate way That s, Merz and Wüthrch (28) studed the condtona MSEP for the mutvarate chan-adder estmates proposed by Pröh and Schmdt (25) and Schmdt (26a) 3 Mutvarate addtve oss reservng method The mutvarate addtve oss reservng method proposed by Hess, Schmdt, and Zocher (26) and Schmdt (26a) s based on a mutvarate near mode whch s sutabe for certan portfoos consstng of severa correated run-off subportfoos The addtve oss reservng method has the foowng features: It s a very smpe cams reservng method whch can easy be mpemented n a spreadsheet 2 Unke the chan-adder method, the addtve oss reservng method combnes past observatons n the upper cams deveopment trange wth externa knowedge from experts or wth a pror nformaton (eg, premum, number of contracts, data from smar run-off portfoos, and market statstcs) 3 It s apped to ncrementa data and thus aows for modeng negatve ncrementa cams n contrast to some other modes such as the (overdspersed) Posson mode [cf Wüthrch and Merz (28)] Ths makes the addtve oss reservng method sutabe for the use of ncurred data, whch often exhbts negatve ncrementa vaues n ater deveopment years due to earer overestmaton of case reserves 4 Unke the chan-adder method, the predcton for the utmate cam does not depend competey on the ast observaton on the dagona Ths means an outer on the dagona w not be proected drecty to the utmate cam Therefore, the addtve oss reservng method s more robust to outers n the ast observatons than the chan-adder method Under the assumptons of ther mutvarate addtve oss reservng mode, Hess, Schmdt, and Zocher (26) and Schmdt (26a) derved a formua for the Gauss-Markov predctor for the nonobservabe ncrementa cam szes whch s optma n terms of a cassca optmaty crteron The components of these predctors are dfferent from the predctors of the unvarate addtve oss reservng method f the subportfoos are correated (eg, see Schmdt (26a; 26b) for the unvarate addtve oss reservng method) Ths means that the predctors of the unvarate method are not optma n the case of correated subportfoos However, Hess, Schmdt, and Zocher (26) and Schmdt (26a) dd not derve an estmator of the condtona MSEP for the mutvarate addtve oss reservng method Snce n actuara practce and scence the condtona MSEP s a very popuar measure to quantfy the uncertantes n cams reserves, ths paper ams to f that gap These studes of uncer- VOLUME 3/ISSUE CASUALTY ACTUARIAL SOCIETY 33

4 Varance Advancng the Scence of Rsk Fgure Cams deveopment trange number n tanty are especay cruca n the deveopment of new sovency gudenes where one exacty quantfes the rsk profe of the dfferent nsurance companes More precsey, we formuate a stochastc mode for the mutvarate addtve oss reservng method to derve an estmator for the condtona MSEP usng the Gauss-Markov predctor proposed by Hess, Schmdt, and Zocher (26) and Schmdt (26a) Furthermore, by means of a detaed exampe, ths estmator s then compared to the estmator for the condtona MSEP of the unvarate predctor (e, f we gnore the correaton structure between ndvdua subportfoos) as we as to the estmator for the condtona MSEP of the mutvarate chan-adder methods consdered by Braun (24) and Merz and Wüthrch (28) 2 Notaton and mutvarate framework In the seque we assume that the data for the N run-off subportfoos consst of run-off tranges of observatons of the same sze However, the mutvarate addtve oss reservng method can aso be apped to other shapes of data (eg, run-off trapezods) In these N tranges the ndces n, n N, refer to subportfoos (tranges),, I, refer to accdent years (rows), and, J, refer to deveopment years (coumns) Fgure shows the cams data structure for the N cams deveopment tranges descrbed above The ncrementa cams (e, ncrementa payments, change of reported cam amount, or number of reported cams wth reportng deay ) of run-off trange n for accdent year and deveopment year are denoted by (n), and cumuatve cams (e, cumuatve payments, cams ncurred, or tota number of reported cams) of accdent year up to deveopment year are gven by C (n), = k= (n),k : () We assume that the ast deveopment year s gven by J, thats (n), = for a >J,and the ast accdent year s gven by I Moreover, our assumpton that we consder run-off tranges mpes I = J 34 CASUALTY ACTUARIAL SOCIETY VOLUME 3/ISSUE

5 Predcton Error of the Mutvarate Addtve Loss Reservng Method for Dependent Lnes of Busness Usuay, at tme I, we have observatons D (n) I = f (n), ; + Ig, (2) for a run-off subportfoos n 2f,:::,Ng Ths meansthatattmei (caendar year I) wehavea tota of observatons over a subportfoos D N I = N[ D (n) I, (3) and we need to predct the random varabes n ts compement D N,c I = f (n), ; I, + >I, n Ng: (4) For the dervaton of the condtona MSEP for severa run-off subportfoos, t s convenent to wrte the data of the N subportfoos n vector form Thus, we defne the N-dmensona random vectors of ncrementa and cumuatve payments by, =( (),,:::,(N), ) and (5) C, =(C (),,:::,C(N), ) for 2f,:::,Ig and 2f,:::,Jg Moreover,we defne the N-dmensona coumn vector consstng of ones by =(,:::,) 2 R N (6) and denote by a D(a)= C A (7) a N the N N-dagona matrx of the vector a = (a,:::,a N ) 2 R N 3 Mutvarate addtve oss reservng method The addtve oss reservng method s easy to appy It s based on the study of ndvdua ncrementa oss ratos We defne for 2f,:::,Ig and 2f,:::,Jg the N-dmensona vector of ndvdua ncrementa oss ratos for accdent year and deveopment year by M, =(M (),,:::,M(N), ) =V,, (8) wthavoumemeasure V = V (,) V (,2) V (,N) V (2,) V (2,2) V (2,N) V (N,) V (N,2) V (N,N), C A (9) whch s a determnstc postve defnte symmetrc N N-matrx The component M (n), of M, denotes the ndvdua ncrementa oss rato (reatve to V ) for accdent year and deveopment year of subportfoo n In the unvarate case N =wehave M, =, =V, () where V s an approprate (determnstc) voume measure If, denotes ncrementa payments and V s the tota premum receved for accdent year, thenm, tes how the tota oss rato s pad over tme 3 Mutvarate addtve oss reservng mode The foowng mutvarate addtve oss reservng mode s a speca case of the mutvarate cams reservng mode studed by Hess, Schmdt, and Zocher (26) and Schmdt (26a) MODEL ASSUMPTIONS 3 (MULTIVARIATE DITIVE MODEL) ² Incrementa payments of dfferent accdent years are ndependent ² There exst N N-dmensona determnstc postve defnte symmetrc matrces V,:::,V I and N-dmensona constants ( =,:::,J) m =(m (),:::,m (N) ) and () ¾ =(¾ (),:::,¾(N) ) VOLUME 3/ISSUE CASUALTY ACTUARIAL SOCIETY 35

6 Varance Advancng the Scence of Rsk wth ¾ (n) > foran =,:::,N as we as N- dmensona random varabes ", =(" (),,:::,"(N), ), (2) such that for a 2f,:::,Ig and 2f,:::,Jg we have, =V m +V =2 D(", ) ¾ : (3) Moreover, the random varabes ", are ndependent wth E[", ]= and Cov(",,", ) = E[", ", ]= ½ (,2) ½ (,N) ½ (2,N) C ½ (2,) ½ (N,) ½ (N,2), C A (4) knowedge for subportfoo m nfuences the ncrementa payments for another subportfoo n n accdent year by choosng V (n,m) 6=Inths case we obtan a nondagona matrx V In the unvarate case N =, the addtve mode satsfes, =V = m + V =2 ¾ ",, (5) wth E[, ]=V m and Var(, )=V ¾ 2 : (6) Hence ths mode can aso be nterpreted as a GLM mode wth Gaussan varance functon (e, V(x) = ), voume measure V and dsperson parameter ¾ 2 [cf McCuagh and Neder (989)] Under Mode Assumptons 3 we have Cov(,,, )=V =2 V =2, (7) where =E[D(", ) ¾ ¾ D(",)] =D(¾ ) Cov(",,", ) D(¾ ) (¾ () )2 ¾ () ¾(2) ½(,2) ¾ () ¾(N) ¾ (2) ¾() ½(2,) (¾ (2) )2 ¾ (2) ¾(N) = ½(,N) ½(2,N) ¾ (N) ¾() ½(N,) ¾ (N) ¾(2) ½(N,2) (¾ (N) )2 : (8) C A where ½ (n,m) 2 (,) for n,m 2f,:::,Ng and n 6= m Ceary, n most practca appcatons V s chosen to be dagona so as to represent a voume measure of accdent year, known a pror (eg, premum, number of contracts, expected number of cams, etc), or an estmate from externa knowedge such as experts, smar portfoos, or market statstcs (see Exampe n Secton 6) However, we can aso take nto account that the voume measure or estmate from externa By Mode Assumptons 3 we restrct any assumpton regardng the correaton between the N run-off subportfoos to each of the correspondng deveopment years ( =,:::,J) nthe N run-off tranges Matrx refects the correaton structure between the ncrementa cams of deveopment year n the N dfferent subportfoos Often correatons between dfferent run-off subportfoos are attrbuted to cams nfaton Under ths pont of vew, t may seem more reasonabe to aow for correaton between 36 CASUALTY ACTUARIAL SOCIETY VOLUME 3/ISSUE

7 Predcton Error of the Mutvarate Addtve Loss Reservng Method for Dependent Lnes of Busness the ncrementa cams of the same caender year (dagonas of the cams deveopment tranges) However, ths woud contradct the assumpton of ndependent accdent years whch s common to most cams reservng methods, and n fact aso necessary to deveop reasonabe estmators from a mathematca pont of vew The Mutvarate Addtve Mode 3 s a speca case of the mutvarate cams reservng mode proposed by Hess, Schmdt, and Zocher (26) and Schmdt (26a), n contrast to whch we assume that ncrementa payments, are ndependent (nstead of ony uncorreated) and generated by the tme seres (3) REMARK 32 ² The ncrementa cams, and k, are ndependent for 6= k or 6= ² The N-dmensona expected ncrementa oss ratos (m ) J can be nterpreted as a mutvarate scaed expected reportng/cashfow pattern over the dfferent deveopment years ² In (7) we use the notaton nstead of snce t smpfes the comparabty wth the dervatons and resuts n Merz and Wüthrch (28) ² Snce we assume that V s a postve defnte symmetrc matrx, there s a we-defned postve defnte symmetrc matrx V =2 (caed square root of V )satsfyngv =V =2 V =2 We obtan for the condtona expectaton (best estmate) E[C,J DI N ] of the utmate cam C,J : PROPERTY 33 Under Mode Assumptons 3 we have for a I J + I E[C,J D N I ]=E[C,J C,I ] = C,I +V =I + m : (9) PROOF Usng the ndependence of the ncrementa cams we obtan 2 3 E[C,J DI N ]=C,I + E 4, D I N 5 = C,I + =I + =I + = C,I +V Ths fnshes the proof E[, ] =I + m = E[C,J C,I ]: (2) QED Ths resut motvates an agorthm for estmatng the expected utmate cams gven the observaton DI N IftheN-dmensona expected ncrementa oss ratos (m ) J are known, the expected outstandng cams abtes of accdent year for the N correated run-off tranges based on the nformaton DI N are estmated by E[C,J D N I ] C,I =V =I + m : (2) However, n most practca appcatons we have to estmate the ratos m from the data n the upper eft trange Hess, Schmdt, and Zocher (26) and Schmdt (26a) propose the foowng mutvarate estmates, for =,:::, J ˆm =(ˆm () I = I =,:::, ˆm (N) ) V =2 V=2 (V =2 V=2 A ) M, : (22) The varabe ˆm (n) denotes the estmated ncrementa oss rato for deveopment year and runoff trange n 2f,:::,Ng basedonthenformaton DI N Note that the covarance structure between the ncrementa cams n the dfferent runoff subportfoos s ncorporated nto the estmaton of m through the matrx VOLUME 3/ISSUE CASUALTY ACTUARIAL SOCIETY 37

8 Varance Advancng the Scence of Rsk Hess, Schmdt, and Zocher (26) and Schmdt (26a) showed the foowng property, whch states that the mutvarate ncrementa oss rato estmates (22) are optma estmators of m wth respect to the crteron of mnma expected squared oss PROPERTY 34 Under Mode Assumptons 3, the estmator ˆm s an unbased estmator for m, whch mnmzes the expected squared oss among a N-dmensona near combnatons of the unbased estmators (M, ) I for m, e, E[(m ˆm ) (m ˆm )] 2 I = mn 4@m W, M A, W, 2R N N E = 3 W, M A, 5: = (23) PROOF See proof of Theorem 4 n Schmdt (26a) QED Note, n Property 34 we assume that the covarance matrx sknownhowever,fwe do not have a reabe estmate for ths covarance matrx t s often more approprate n practce to use the unvarate estmators Property 34 motvates the foowng estmator for the condtonay expected utmate cam: ESTIMATOR 35 (Mutvarate addtve estmator) The mutvarate addtve estmator for E[C, D N I ] s for + I gven by d dc, =( C (),,:::, C d(n), ) = Ê[C, D N I ]=C,I +V =I + ˆm : (24) Ths means that n the mutvarate addtve method we predct the normazed cumuatve cams V C, by the sum of the ast observed normazed cumuatve cams V C,I and the weghted estmated ratos ˆm I +,:::, ˆm,gven the nformaton DI N From (24) we obtan for the ncrementa payments, wth + >Ithe predctors d d, =( (),,:::, d(n), ) REMARK 36 =V ˆm : (25) ² In the case = J (note that we assume I = J) we have ˆm J = M,J ² Estmator (22) s a weghted average of the observed ndvdua normazed ncrementa cams M, InthecaseN = (e, ony one run-off subportfoo), the estmators (22) concde wth the unvarate estmated ncrementa oss ratos I ˆm = = V P I k= V k M, (26) wth determnstc weghts V, whch are used n the unvarate addtve oss reservng method, and from Estmator 35 we obtan the unvarate addtve estmator P I dc,j = k= C,I + k, P I =I + k= V V k (27) [see, for exampe, Schmdt (26a; 26b)] ² If we negect the covarance structure between the ncrementa cams n the dfferent run-off subportfoos [e, n (22) we set =I, where I denotes the dentty matrx], we obtan the foowng (unbased) estmator I ˆm () A V M, : (28) = = V = Moreover, f the voumes V are dagona matrces, then the components of (28) are gven by ˆm (n)() I = = V (n,n) P I k= V(n,n) k M (n), : (29) Ths means that n ths case the components of ˆm () are gven by the estmators of the unvarate addtve oss reservng method Itcaneasybeseenthat ˆm does not depend on the matrx f = J or f and V,:::, V I are dagona In ths case the N components 38 CASUALTY ACTUARIAL SOCIETY VOLUME 3/ISSUE

9 Predcton Error of the Mutvarate Addtve Loss Reservng Method for Dependent Lnes of Busness ˆm (),:::, ˆm (N) of (22) concde wth the unvarate estmators (29) for the N run-off subportfoos Ths means that f,:::, J 2 and V,:::,V I are dagona matrces, the foowng estmates concde: ) the estmaton for the whoe portfoo based on the unvarate estmators (26) for every ndvdua run-off subportfoo, 2) the mutvarate predcton based on the estmators (28), and 3) the mutvarate predcton based on the mutvarate estmators (22) However, Property 34 shows n other cases t s more reasonabe to use the mutvarate estmators (22) Moreover, under Mode Assumptons 3 t hods: PROPERTY 37 Under Mode Assumptons 3 we have a) ˆm and ˆm k are ndependent for 6= k; ³ PI ; b) Var( ˆm )= = V=2 V=2 c) gven C,I,theestmatorC d,j s an unbased estmator for E[C,J DI N]=E[C,J C,I ], e, E[ C d,j C,I ]=E[C,J DI N]; d) C d,j s an unbased estmator for E[C,J ], e, E[ C d,j ]=E[C,J ] PROOF a) Foows from the ndependence of the normazed ncrementa cams M, =V, and M k, =V k k, for 6= b) Usng (7) we obtan Var(M, )=V Var(, ) V =V =2 V =2 : (3) Wth the ndependence of the M, ths eads to where c) We have I A V =2 = E[ d C,J C,I ] = C,I +V = C,I +V =I + =I + V=2 E[ ˆm ] A : (32) m = E[C,J D N I ]: (33) d) Foows mmedatey from c) Ths fnshes the proof QED Observe that Property 37 c) shows that the Estmator 35 s an unbased estmator for E[C,J DI N ] Furthermore, ths mmedatey mpes that the estmator for the aggregated utmate cam of one snge accdent year N d C (n) = C d,j (34),J s, gven C,I, an unbased estmator for P N E[C (n),j DN I ] 4 Condtona MSEP In ths secton we consder the uncertanty n the cams reserves predcted by the estmators P N d C (n),j and P I P N d = C (n),j, gven the ob- I Var( ˆm )=A (V =2 = V=2 ) M A, A 2 3 I =A 4 (V =2 V=2 ) Var(M, ) (V =2 V=2 ) 5 A = 2 I =A 4 V =2 = V=2 3 5 A =A, (3) VOLUME 3/ISSUE CASUALTY ACTUARIAL SOCIETY 39

10 Varance Advancng the Scence of Rsk servatons DI N Ths means our goa s to derve an estmate of the condtona MSEP for ndvdua accdent years 2f,:::,Ig whch s defned as à msep P N d C (n) n C(n),J DN,J I 2à N = E 4 d C (n),j N C (n),j 2 DN I 3 5 = E[( d C,J C,J ) ( d C,J C,J ) D N I ] (35) as we as an estmate of the condtona MSEP for aggregated accdent years msep d C (n) A,n C(n),J DN,J I,n = E 4@ d C (n),j C (n),j A2 D N 7 I 5:,n,n (36) 4 Condtona MSEP for snge accdent years We choose 2f,:::,Ig Snce the estmator P N d C (n),j s known at tme t = I (e, t s based on observatons from DI N ), the condtona MSEP (35) can be decouped nto condtona process varance and condtona estmaton error, that s we derve estmates for both the condtona process varance and the condtona estmaton error for N correated run-off tranges 4 Condtona process varance In ths subsecton we derve an estmate for the condtona process varance of a snge accdent year Var(C,J DI N ) We obtan the foowng resut: PROPERTY 4 (Process varance for a snge accdent year) Under Mode Assumptons 3 the condtona process varance for the utmate cam C,J of accdent year 2f,:::,Ig s gven by Var(C,J DI N ) = V =I + A V =2 : (38) PROOF Usng the ndependence of the ncrementa cam payments, we have Var(C,J DI N ) = A = V =2 =I Var(, ) A =I +, =I + for >I J Ths competes the proof A V =2 (39) QED msep P n C(n),J DN I à N d C (n),j = Var(C,J D N I ) {z } condtona process varance + ( C d,j E[C,J DI N ]) ( C d,j E[C,J DI N ]) : (37) {z } condtona estmaton error The condtona process varance orgnates from the stochastc movement of C,J, whereas the condtona estmaton error refects the uncertanty n the estmaton of the condtona expectaton (best estmate) E[C,J DI N ] In the seque If we repace the parameter n (38) by ts estmate (cf Secton 5), we obtan an estmator of the condtona process varance for accdent year Moreover, from (39) we obtan the recursve formua for the condtona process varance of 4 CASUALTY ACTUARIAL SOCIETY VOLUME 3/ISSUE

11 Predcton Error of the Mutvarate Addtve Loss Reservng Method for Dependent Lnes of Busness PROOF Usng Propertes 37 a) b) we obtan E[( C d,j E[C,J DI N ]) ( C d,j E[C,J DI N ]) ] 2 = E 4@ V ( ˆm m ) V ( ˆm m ) A 3 5 = 2 = 6 V 4 =I + =I + =I + =I + Var( ˆm ) A V (43) V =2 = V=2 A V : (44) accdent year Var(C, D N I ) = (Var(C, D N I ) +V =2 V =2 ), (4) for = I +,:::,J wth Var(C,I DI N)= 42 Condtona estmaton error Now we estmate the uncertanty n the estmaton of E[C,J DI N] by the estmator C d,j Ths means we derve an estmator for the second term on the rght-hand sde of (37) We estmate the condtona estmaton error by ts expected vaue E[( C d,j E[C,J DI N ]) ( d C,J E[C,J D N I ]) ] : (4) We obtan the foowng resut: PROPERTY 42 (Estmator of the estmaton error for a snge accdent year) Under Mode Assumptons 3 the estmator (4) of the condtona estmaton error for P N d C (n),j wth 2 f,:::,ig s gven by E[Var( C d,j C,I )] 2 Ã I 3 = V 4 V =2 V=2 5 =I + = V : (42) On the other hand, usng Property 37 c), we have E[( d C,J E[C,J D N I ]) ( d C,J E[C,J D N I ]) ] = E[Var( d C,J C,I )] : (45) Ths fnshes the proof QED Note, we can rewrte (42) n the recursve form E[Var( d C, C,I )] = E[Var( d C, C,I )] + V Ã I = V =2 V=2 V (46) for = I +,:::,J wth Var( d C,I C,I ) = Fnay, repacng the parameters n (38) and (42) by ther estmates (see Secton 5), we obtan the foowng estmator of the condtona MSEP for a snge accdent year: RESULT 43 (Condtona MSEP for a snge accdent year) Under Mode Assumptons 3 we VOLUME 3/ISSUE CASUALTY ACTUARIAL SOCIETY 4

12 Varance Advancng the Scence of Rsk have the estmator for the condtona MSEP of the utmate cam for a snge accdent year 2 fi J +,:::,Ig à dmsep P N d C (n) n C(n),J DN,J I = V =2 + V 4 =I + 2 =I + ˆ V =2 à I 3 V =2 ˆ V=2 5 = V, (47) where the estmated covarance matrx ˆ s gven n (59), beow For N = formua (47) reduces to the estmator of the condtona MSEP for a snge portfoo n the unvarate addtve oss reservng method msep d C,J D I ( C d,j ) = V =I + ˆ¾ 2 + V2 =I + ˆ¾ 2 P I = V, (48) where V s a known one-dmensona voume measure for accdent year [cf Mack (22)] 42 Condtona MSEP for aggregated accdent years In the foowng we consder the condtona MSEP for aggregated accdent years Our goa s to derve an estmate for (36) From Mode Assumptons 3 we know that the utmate cams C,J and C k,j of two accdent years and k wth <k I are ndependent However, snce the estmators C d,j and C d k,j use the same observatons DI N for estmatng the parameters m,dfferent accdent years are no onger ndependent We start wth the consderaton of two accdent years <k à N msep P P C d (n) N n C(n),J + n C(n) k,j DN,J + C d (n) k,j I 2à N = E 4 ( C d(n),j + C d(n) k,j ) 2 N (C (n),j + C(n) k,j ) DN I 3 5 : (49) We obtan for the condtona MSEP of the sum of two accdent years the decomposton nto process varance and condtona estmaton error whch eads to msep P n C(n),J +P n C(n) k,j DN I =msep P n C(n),J DN I +msep P n C(n) k,j DN I à N d C (n),j à N d + C (n),j à N d C (n) k,j +2 ( d C,J E[C,J D N I ]) N d C (n) k,j ( d C k,j E[Ck,J D N I ]) : (5) Ths shows that we have to derve an estmator for the cross product [thrd term on the rght sde of (5)], whch comes from the dependence descrbed above Anaogousy to (4), we estmate ths cross product by ts expected vaue E[( d C,J E[C,J D N I ]) ( d C k,j E[Ck,J D N I ]) ] (5) and obtan the foowng resut: PROPERTY 44 (Estmator of the cross product) Under Mode Assumptons 3 the estmator (5) of the cross product of aggregated accdent years 42 CASUALTY ACTUARIAL SOCIETY VOLUME 3/ISSUE

13 Predcton Error of the Mutvarate Addtve Loss Reservng Method for Dependent Lnes of Busness and k wth <k I s gven by E[( C,J c E[C,J D N I ]) ( Ck,J d E[Ck,J D N I ]) ] 2 à I 3 = V 4 V =2 V=2 5 Vk : =I + = (52) PROOF Anaogousy to the proof of Property 42 we obtan for <k E[( C,J c E[C,J D N I ]) ( Ck,J d E[Ck,J D N I ]) ] " J # = V Var( ˆm ) V k 2 = V 4 =I + à I =I + = V =2 V=2 3 5 Vk : (53) QED Puttng (47) and (52) n (5) eads to the foowng estmator for the condtona MSEP of the utmate cam for aggregated accdent years: RESULT 45 (Condtona MSEP for aggregated accdent years) Under Mode Assumptons 3 we have the estmator for the condtona MSEP of the utmate cam for aggregated accdent years à I dmsep P N d C Pn (n) C(n),J DN,J I = à I = dmsep P N d C (n) n C(n),J DN,J I = V <k I à I =I + = V =2 ˆ V=2 3 5 V k, (54) where the estmated covarance matrx ˆ s gven n (59), beow For N =, formua (54) reduces to the estmator of the condtona MSEP for aggregated accdent years n the unvarate addtve method à I msep d P C,J D I dc,j = = I msep d C,J D I ( C d,j ) = +2 <k I V V k =I + ˆ¾ 2 P I = V (55) wth known one-dmensona voume measure V for accdent year [cf Mack (22)] 5 Parameter estmaton For the estmaton of the cams reserves and the condtona MSEP we need estmates of the N-dmensona parameters m,:::,m J and of the N N-dmensona covarance parameters,:::, J Estmates of the mutvarate ncrementa oss ratos m are gven n (22) However, estmator (22) s ony an mpct estmator for m snce t depends on parameter, whch on the other hand s estmated by means of ˆm Therefore, as n the mutvarate chan-adder method [cf Merz and Wüthrch (28)], we propose an teratve estmaton of these parameters In ths sprt, the true estmaton error s sghty arger because t shoud aso nvove the uncertantes n the estmate of the varance parameters In order to obtan a feasbe MSEP formua we negect ths term of uncertanty Estmaton of m As startng vaues for the teraton we defne ˆm () by (28) for =,:::,J Estmator ˆm () s an unbased optma estmator for m f the N run-off subportfoos are uncorreated However, f the subportfoos are correated, t s st unbased but no onger optma (cf Property 34) From ˆm () we derve an estmate ˆ () of for =,:::,J [see estmator (59) beow] Then ths estmate s used to determne VOLUME 3/ISSUE CASUALTY ACTUARIAL SOCIETY 43

14 Varance Advancng the Scence of Rsk ˆm () va ˆm (k) =(ˆm ()(k) I = I =,:::, ˆm (N)(k) ) V =2 (k) ( ˆ (V =2 (k) ( ˆ ) V =2 A ) V =2 ) M, (56) for =,:::,J Ths agorthm s then terated unt t has suffcenty converged Estmaton of The N N-dmensona covarance parameters are estmated teratvey from the data for =,:::,J A postve semdefnte estmator of the postve defnte matrx s gven by ˆ = I I V =2 (, V ˆm () ) = (, V ˆm () ) V =2 (57) for =,:::,J If the matrces V are a dagona, the dagona eements of the random matrx (57) are unbased estmators of the correspondng dagona eements (¾ () )2,:::,(¾ (N) )2 (58) of Its nondagona eements sghty underestmate the absoute vaue of the correspondng nondagona eements of However,ths ack of unbasedness s not too mportant snce the random matrx (57) has to be nverted anyway and the nverse of an unbased estmator s n genera not unbased [cf Appendx of Merz and Wüthrch (28)] Ths eads to the foowng teraton for the estmator of : ˆ (k) = I I = V =2 (, V ˆm (k ) ) (, V ˆm (k ) ) V =2 (59) for =,:::,J and k If we have enough data (e, we have a runoff trapezod wth I>J),weareabetoestmate teratvey the parameter J by (59) Otherwse, we can use the estmates b' (n,m)(k) of the eements ' (n,m) of for J n teraton k [e,b' (n,m)(k) s an estmate of ' (n,m) = ¾ (n) ¾(m) ½(n,m) n teraton k, cf (8)] to derve estmates b' (n,m)(k) J of the eements of J for a n m N For exampe, ths can be done by extrapoatng the usuay exponentay decreasng seres b' (n,m)(k),:::,b' (n,m)(k) J 2 (6) by one addtona member b' (n,m)(k) J for n m N and k However, one needs to care- ˆ (k) J fuy check that the estmate s postve defnte In hgher dmensona cases ths s often nontrva, and n fact, many choces are not postve defnte, whch cas for addtona adustments Moreover, observe that the N N- ˆ (k) dmensona estmate s snguar when I N + 2, snce n ths case the dmenson of the near space generated by any reazatons of the (I +) N-dmensona random vectors V =2 (, V ˆm (k ) ) wth 2f,:::,I g (6) s at most I + I (I N +2)+= N Furthermore, the reazatons of (6) may be (neary) neary dependent for some <I N + 2 whch mpes that the correspondng ˆ (k) reazaton of the random matrx s -condtoned or even snguar Therefore, n practca appcaton t s mportant to verfy whether (k) the estmates ˆ are we-condtoned or not and to modfy those estmates (eg, by extrapoaton as n the exampe beow) whch are not we-condtoned Many methods have been suggested to mprove the estmaton of the covarance matrx so that the estmate s postve defnte and we-condtoned By producng a we-condtoned covarance es- 44 CASUALTY ACTUARIAL SOCIETY VOLUME 3/ISSUE

15 Predcton Error of the Mutvarate Addtve Loss Reservng Method for Dependent Lnes of Busness Tabe Genera abty run-off trange (ncrementa cams (), ), source Braun (24) Genera abty run-off trange AY/DY ,966 3,86 9,36 95,2 83,74 42,53 37,882 6,649 7,669,6,738 3,572 6,823,893 49,685 3,659 9,592,43 75,442 44,567 29,257 8,822 4, ,73 2, ,94 8,34 49,56 5,825 78,97 4,77 4,76 7,95,97 2,643,3, ,937 88,246 34,35 39,97 74,45 65,4 49,65 2, ,48 2, ,92 79,48 7,2 3,6 79,64 8,364 2,44,324 6, ,78 73,879 7,295 44,76 93,694 72,6 4,545 25,245 7, ,35 93,57 8,77 53,86 2,96 86,753 45,547 23, ,73 27,43 24,558 22,276,88 4,966 59, ,58 245,7 232,223 93,576 65,86 85,2 9 73, ,73 262,92 232,999 86,45 39,82 297,37 372, ,27 54, ,5 54, ,24 576, ,325 Tabe 2 Auto abty run-off trange (ncrementa cams (2), ), source Braun (24) Auto abty run-off trange AY/DY ,423 33,538 65,2 3,358 27, ,889 4, ,8, ,296 52,879 7,438 4,686 22,9 25,35 7,96 4,843 3, ,383, ,325 62,99 6,365 5,432 55,224 7,95 8,234,49 2, ,94 6,732 79,458 46,642 29,384 5,8 3,598 5,527 2, ,8 4 7,333 7,68 92,6 36,227,872 8,76 3,8 3, ,643 7,48 85,734 6,226 8,479 3,556 7,523, ,22 27,22,8 56,83 28,362 29,79,244 2, ,98 2,39 4,397 45,277 34,888 3,93 7, ,95 25,478 98,68 62,846 52,435 22, , ,796 4,2 82,259 59,29 249,58 33,52 42,26 22,23 258, , , ,762 54, ,997 tmate we automatcay get a we-condtoned estmate for the nverse of the covarance estmate Most of these approaches rey on the concept of shrnkage whch s qute smar to the we-known actuara concept of credbty For more detas and other advanced methods on covarance matrx estmaton we refer to Schäfer and Strmmer (25) 6 Exampe: two correated abty run-off subportfoos To ustrate the methodoogy, we consder two correated run-off portfoos A and B (e, N = 2), whch contan data of genera and auto abty busness, respectvey The data are gven n Tabes and 2 n ncrementa form These are the data used n Braun (24) and aso n Merz and Wüthrch (27; 28) The assumpton that there s a postve correaton between these two nes of busness s ustfed by the fact that both run-off portfoos contan abty busness; that s, certan events (eg, body nury cams) may nfuence both run-off portfoos, and we are abe to earn from the observatons from one portfoo about the behavor of the other portfoo We assume that the 2 2-matrces V are dagona and ther dagona eements V (,) and V (2,2) are pror estmates of the utmate cams n the dfferent accdent years n run-off portfoo A and B, respectvey Such pror estmates are usu- VOLUME 3/ISSUE CASUALTY ACTUARIAL SOCIETY 45

16 Varance Advancng the Scence of Rsk Tabe 3 Pror estmates and chan-adder estmates of the utmate cams Run-off portfoo A V (,) CL C c(),j Run-off portfoo B V (2,2) CL C c(2),j 5,3 549,589 43,23 39, , ,74 537, , ,33 68,4 589,45 54, , , ,49 486, ,32 783,593 5,498 58, , ,88 598, , , ,86 68, ,3 7,56,778,98,2 698, ,4 8,24,569,54,92 74,29 684,79 9,397,23,43,49 93, ,543,832,676,735, , ,734 2,56,78 2,65,99,34,29,69,26 2 2,559,345 2,66,56,538,96,474,54 3 2,456,99 2,274,94,487,234,426,6 Tota 7,758,43 7,498,658,86,268,823,48 ay obtaned from budget fgures, pan vaues or from premum cacuaton parameters Tabe 3 shows these a pror estmates as we as the correspondng cassca unvarate chan-adder estmates C d() CL,J and C d(2) CL,J for comparson purposes We see that the pror estmates and the unvarate chan-adder estmates are cose together [for the unvarate chan-adder method see, eg, Mack (993) or Buchwader, Bühmann, Merz, and Wüthrch (26)] Snce I = J = 3 we do not have enough data to derve an estmate of the 2 2-matrx 2 usng estmator (59) Therefore, we use the extrapoaton b' (n,m) 2 =mnf(b' (n,m) ) 2 =b' (n,m),b' (n,m) g (62) to derve estmates of ts eements ' (n,m) 2 = ¾ (n) 2 ¾ (m) 2 ½(n,m) 2 for n,m =,2 (note ½ (,) 2 = ½ (2,2) 2 =) Moreover, snce estmator (59) woud ead to an -condtoned matrx ˆ,wehaveasoestmated the eements of the 2 2-matrx by b' (n,m) =mnf(b' (n,m) ) 2 =b' (n,m) 9,b' (n,m) 9 g: (63) Tabe 4 shows the estmates for the parameters m, ¾ and ½ (,2) after three teratons k =,2,3 We observe fast convergence of the two-dmensona estmates ˆm (k ), b¾ (k) and the one-dmensona estmates ˆ½ (,2)(k) (k =,2,3) n the sense that there are barey any changes n the estmates after three teratons The frst and second component of the estmates ˆm () and b¾ () are the parameter estmates used n the unvarate addtve method apped to the ndvdua run-off portfoos A and B, respectvey Except for deveopment years, 6, and, we observe postve estmates ˆ½ (,2)(k) for the correaton coeffcents The three negatve estmates shoud not be overstated snce they are cose to zero The frst two coumns of Tabe 5 show for each accdent year the cams reserves for run-off subportfoos A and B estmated by the (unvarate) addtve oss reservng method Coumn portfoo (k = ) shows the reserves for the whoe portfoo consstng of the two run-off subportfoos A and B estmated by the mutvarate addtve oss reservng method These vaues are based on the estmates ˆm () and therefore concde wth the sum of the cams reserves for the two ndvdua subportfoos Coumns portfoo (k = 2) and portfoo (k = 3) contan the cams reserves for the whoe portfoo based on the estmates ˆm () and ˆm (2), respectvey These estmates ead to a tota reserve whch s about 6,9 hgher than the one based on ˆm () The coumn denoted by overa cacuaton shows the estmated reserve when frst aggregatng both run-off tranges to one snge run-off trange and then estmatng the cams reserves wth the (unvarate) addtve oss reservng method Snce n ths approach two run-off tranges wth dfferent deveopment patterns are added together (cf components of estmates ˆm (k) n Tabe 4), ths approach s ony reasonabe f the proporton of exposures from each trange does not change sgnfcanty over the dfferent accdent years In our exampe ths approach eads to a 46 CASUALTY ACTUARIAL SOCIETY VOLUME 3/ISSUE

17 Predcton Error of the Mutvarate Addtve Loss Reservng Method for Dependent Lnes of Busness Tabe 4 Estmates ˆm (k ), b¾ (k) and ˆ½ (,2)(k) for the parameters m, ¾ and ½ (,2) n the frst three teratons k =,2,3 A/B ˆm () :5 355 : :26 b¾ () ˆ½ (:2)() : : : ˆm () :5 354 :97 :26 b¾ (2) ˆ½ (:2)(2) : : : ˆm (2) :5 354 :97 :26 b¾ (3) ˆ½ (:2)(3) : : : Tabe 5 Estmated reserves Addtve method Chan-adder method Mutvarate Unvarate Mutvarate Unvarate Unvarate portfoo portfoo portfoo portfoo portfoo portfoo subportfoo A subportfoo B (k =) (k =2) (k = 3) overa cac Braun (24) MW (28) 2, ,26 2,26 2,26 2,262,8,8 2 5, ,76 5,96 5,96 5,442 4,655 4, ,68,93,8,85,85,356,827, , ,6 4,677 4,677 3,82 6,22 6, ,386 3,54 29,54 29,723 29,723 28,266 29,2 29,49 6 4,96 3,243 44,49 44,749 44,753 4,64 45,793 46, ,946,87 9,32 9,88 9,83 84,45 86,4 87, ,95 2,58 64,973 65,79 65,75 53,693 57,65 58, ,823 55, ,448 34,6 34,66 328,7 344,3 346,42 594,362,5 75,53 76,398 76,45 659,59 679,82 68,729,77,55 235,757,33,272,33,647,33,653,246,294,287,458,287,654 2,86, ,4 2,374,947 2,376,6 2,376,7 2,325,74 2,453,38 2,45,6 3 2,225,22,38,295 3,263,56 3,264,85 3,264,826 3,223,75 3,,679 3,92,98 Tota 6,3,53 2,47,68 8,359,83 8,366,62 8,366,9 8,23,852 8,28,874 8,25,35 tota reserve whch s about 235,3 242,3 ess than the one obtaned by separate cacuaton of the cams reserves n run-off subportfoos A and B The ast two coumns show the vaues cacuated by the mutvarate chan-adder reservng methods proposed by Braun (24) (e, chan-adder factors are estmated n a unvarate way) and Merz and Wüthrch (28) (e, chan-adder factors are estmated n a mutvarate way), respectvey We see that the mutvarate addtve oss reservng method eads to a tota reserve whch s about 47,2 5,8 hgher than the ones obtaned by the two mutvarate chan-adder methods Tabe 6 shows for each accdent year the estmates for the condtona process standard de- VOLUME 3/ISSUE CASUALTY ACTUARIAL SOCIETY 47

18 Varance Advancng the Scence of Rsk Tabe 6 Estmated condtona process standard devatons Addtve method Chan-adder method Mutvarate Unvarate Mutvarate Unvarate Unvarate portfoo portfoo portfoo portfoo overa portfoo portfoo subportfoo A subportfoo B (k =) (k = 2) (k = 3) cacuaton Braun (24) MW (28) 33 57% :% % % % %,289 72%,289 72% %,34 5:8%, %, %, %, % 5, % 5, % 3,64 7% 2,48 227% 2, % 2, % 2, % 2,85 275% 7,29 66% 7,29 66% 4 5,38 392% 2, % 6,42 439% 6,42 437% 6,42 437% 6,96 448% 9,8 65% 9,85 599% 5 2,669 48% 4,743 53% 4,78 5% 4, % 4, % 4,656 58% 6,43 554% 6,49 549% 6 4,763 36% 5,43 555% 7,227 39% 7, % 7, % 7,2 49% 9,2 48% 9,45 49% 7 7,89 22% 6, % 2, % 2, % 2, % 2,33 238% 2,9 255% 2,937 25% 8 23,84 66% 7, % 27,2 64% 27,8 64% 27,8 64% 26,64 73% 28,933 84% 28,966 83% 9 3,227 6% 4, % 36,978 9% 36,985 9% 36,985 9% 37,86 5% 39,28 4% 39,322 4% 43,67 72% 2,49 93% 53,848 76% 53,854 76% 53,854 76% 53,978 82% 63,663 94% 63,724 93% 5,294 48% 28,466 2% 67,39 5% 67,44 5% 67,44 5% 69,957 56% 99,98 78%,4 78% 2 64,43 36% 4,2 7% 9,552 39% 9,569 39% 9,569 39% 94,86 4% 99,543 8% 99,68 8% 3 8,24 36% 5,955 5% 7,567 33% 7,58 33% 7,58 33%,223 34% 36,2 2% 36,2 2% Tota 3,444 2% 77,62 38% 74,596 2% 74,624 2% 74,624 2% 79,43 22% 396,73 48% 396,85 48% Tabe 7 Square roots of estmated condtona estmaton errors Addtve method Chan-adder method Mutvarate Unvarate Mutvarate Unvarate Unvarate portfoo portfoo portfoo wthout portfoo overa portfoo portfoo subportfoo A subportfoo B (k =) (k =2) (k =3) corrn ˆm () cacuaton Braun (24) MW (28) 49 63% :2% % % % % %,32 729%,32 729% % 985 3:9%,3 23%,3 22%,3 22%,3 23%,86 99% 4,533 97,4% 4, % 3,74 2%, %,89 67%,89 67%,89 67%,89 67%,898 83% 6,87 55% 6,87 55% 4 2,96 23%, % 3,55 24% 3,55 239% 3,55 239% 3,56 24% 3, % 7,37 434% 7,34 43% 5 6,7 254% 2,65 829% 7,8 264% 7,8 263% 7,8 263% 7,8 264% 7,64 27% 9, % 9, % 6 7,859 92% 2,75 848% 9,87 26% 9,9 23% 9,9 23% 9,92 26% 8,87 22%, %,742 25% 7,49 3% 3, %,887 3%,89 3%,89 3%,892 3%,283 34% 3,99 63% 3,996 6% 8 2,953 9% 4, 9% 4,5 88% 4,53 88% 4,53 88% 4,56 88% 3,734 89% 6,637 6% 6,644 5% 9 6,473 58% 6,934 25% 9,523 58% 9,527 57% 9,527 57% 9,53 58% 9,446 59% 22,767 66% 22,776 66% 24,583 4% 9,52 86% 28,86 4% 28,865 4% 28,865 4% 28,87 4% 27,84 42% 34,3 5% 34,6 5% 3,469 28% 3,6 56% 36,975 28% 36,982 28% 36,982 28% 36,996 28% 36,798 3% 5,43 4% 5,386 4% 2 38,94 22% 2,38 36% 5,834 2% 5,843 2% 5,843 2% 5,956 2% 5,665 22% 99,933 4% 99,857 4% 3 42,287 9% 23,687 23% 54,274 7% 54,282 7% 54,282 7% 54,38 7% 54,98 7% 3,734 42% 3,59 43% Tota 72,74 27% 74,52 36% 27,9 25% 27,57 25% 27,57 25% 27,3 25% 23,99 25% 33,36 38% 33,74 38% vatons and the correspondng estmates for the coeffcents of varaton The frst two coumns of Tabe 6 contan the vaues for the ndvdua subportfoos A and B cacuated by the unvarate addtve oss reservng method Coumns portfoo (k = ) to portfoo (k = 3) show the estmated condtona process standard devatons for the portfoo consstng of the two subportfoos A and B f we use the mutvarate addtve oss reservng method (frst three teratons) In partcuar ths means that the vaues n coumn k = are based on the parameter estmates ˆm () The coumn denoted by overa cacuaton shows the resuts for the overa cacuaton The ast two coumns show the vaues cacuated by the mutvarate chan-adder reservng methods proposed by Braun (24) and Merz and Wüthrch (28), respectvey Tabe 7 shows the square roots of estmated condtona estmaton errors The frst two coumns contan the estmates for the ndvdua subportfoos A and B cacuated by the unvarate method Coumns portfoo (k = ), portfoo (k = 2) and portfoo (k = 3) show the est- 48 CASUALTY ACTUARIAL SOCIETY VOLUME 3/ISSUE

19 Predcton Error of the Mutvarate Addtve Loss Reservng Method for Dependent Lnes of Busness Tabe 8 Estmated predcton standard errors Addtve method Chan-adder method Mutvarate Unvarate Mutvarate Unvarate Unvarate portfoo portfoo portfoo wthout portfoo overa portfoo portfoo subportfoo A subportfoo B (k =) (k =2) (k =3) corrn ˆm () cacuaton Braun (24) MW (28) 2 85% :% 73 33% 73 33% 73 33% 73 33% 77 34%,845 9%,845 9% %,52 2:%, %, %, %, %,675 38% 7,493 6% 7,493 6% 3,96 24% 2, % 3,39 37% 3,39 37% 3,39 37% 3,39 37% 3,425 33% 9,497 83% 9,497 83% 4 6,2 446% 2, % 7,39 5% 7,32 499% 7,32 499% 7,32 5% 7,59 5% 2,66 744% 2,67 737% 5 4, % 5,46 77% 6,77 566% 6,78 562% 6,78 562% 6,78 566% 6, % 8, % 8, % 6 6,724 49% 5,744 77% 9,477 44% 9, % 9, % 9,479 44% 9,63 46% 22,435 49% 22,459 48% 7 2, % 7, % 23,729 26% 23, % 23, % 23,732 26% 23,79 273% 25,996 32% 26,22 298% 8 27,3 89% 8, % 3,75 86% 3,757 86% 3,757 86% 3,753 86% 29,972 95% 33,376 22% 33,47 2% 9 34,424 2% 5, % 4,85 23% 4,823 23% 4,823 23% 4,88 23% 42,562 29% 45,4 32% 45,442 3% 49,589 83% 23,44 2% 6,94 87% 6,2 86% 6,2 86% 6,99 87% 6,723 92% 72,222 6% 72,282 6% 59,66 55% 3,342 33% 76,868 59% 76,883 59% 76,883 59% 76,878 59% 79,45 63% 2,37 87% 2,434 87% % 44,965 79% 4,78 44% 4,737 44% 4,738 44% 4,777 44% 8,7 46% 223,69 9% 223,92 9% 3 9,67 4% 57, 55% 2,484 37% 2,499 37% 2,499 37% 2,532 37% 23,74 38% 342,377 % 342,322 % Tota 26,63 34% 6,947 52% 27,89 32% 27,938 32% 27,939 32% 27,3 32% 27,358 33% 55,56 62% 55,44 62% mated condtona estmaton errors for the portfoo consstng of the two subportfoos A and B f we use the mutvarate addtve oss reservng method The new coumn wthout corr n ˆm () contans the estmated condtona estmaton errors f we do not take nto account correatons wthn the parameter estmates ˆm and use nstead the estmates ˆm () In contrast to the reserve and the condtona process standard devaton, these estmates do not concde wth the vaues n coumn portfoo (k = ) snce the estmator of the estmaton error for a snge accdent year and the cross product term [e, rght-hand sde of (42) and (52)] are now gven by 2 Ã I Ã I V 4 V V =2 V =2 Ã I and 2 V 4 =I + = Ã I = = V 3 5 V (64) =I + Ã I V = Ã I = V =2 V =2 3 V 5 V k, (65) = respectvey We see (as expected) that the estmaton error s arger (27,3 vs 27,57) f we estmate the parameters on the snge tranges However, the dfference n ths exampe s sma, whch woud ustfy workng wth ˆm () The coumn overa cacuaton shows the estmates for the overa cacuaton The ast two coumns show the vaues cacuated by the mutvarate chan-adder reservng methods proposed by Braun (24) and Merz and Wüthrch (28), respectvey Tabe 8 contans the estmated predcton standard errors and coeffcents of varaton for the same set of modes as above Tabe 9 contans the resuts for the estmated predcton standard errors assumng perfect postve correaton, no correaton, and perfect negatve correaton between the correspondng cams reserves of the two run-off subportfoos A and B These vaues are cacuated by msep d C,J DI N = msep d C (),J DN I +2c msep d =2 + msep d C (2),J DN I C (),J DN I msep d =2 C (2),J DN I (66) wth c =, c = and c =, respectvey Except for accdent year 3, we observe that the estmator n the mutvarate addtve oss reservng method eads to estmates of the predcton standard errors whch are between the ones assumng VOLUME 3/ISSUE CASUALTY ACTUARIAL SOCIETY 49

20 Varance Advancng the Scence of Rsk Tabe 9 Estmated predcton standard errors assumng correaton, and, respectvey Portfoo Portfoo Portfoo dmsep =2 C,J D N I dmsep =2 C,J D N I dmsep =2 C,J D N I correaton = correaton = correaton = ,4, ,826 3, ,5 6,89 3,36 5 9,752 5,325 8, ,469 7,683, ,26 22,24 3, ,66 28,565 8,97 9 5,376 37,94 8,472 73,29 54,85 26,49 9,3 67,392 28,38 2 2,25 87,66 3, ,769 7,5 33,57 Tota 323,56 24,576 9,666 no correaton and a correaton equa to one for a accdent years and a accdent years together (cf coumns 3 5 n Tabe 8) Moreover, we see that an assumed correaton of or woud ead to an estmated predcton standard error that s about 29,5 ower and 52,5 hgher, respectvey, than the one takng the estmated correaton between the two subportfoos nto account 7 Concuson In ths paper we consder the cams reservng probem for a portfoo consstng of severa correated run-off subportfoos The smutaneous study of severa ndvdua run-off subportfoos s motvated by severa mportant facts and s especay cruca n the deveopment of new sovency gudenes However, the cacuaton of the condtona MSEP for the predctor of the utmate cam sze for a whoe portfoo of severa correated run-off subportfoos s more sophstcated snce now mutdmensona matrx cacuatons are nvoved and the mode parameters are nterdependent so that generay an teratve parameter estmaton procedure s requred In the present paper we study a speca case of the mutvarate addtve oss reservng mode proposed by Hess, Schmdt, and Zocher (26) and Schmdt (26a) Our derved formuas for the condtona MSEP n the addtve cams reservng method can be used to quantfy the uncertanty n the cams reserves for a snge runoff portfoo (e, N = ) or a whoe portfoo of severa correated run-off subportfoos (e, N> ) and can easy be mpemented n a spreadsheet By means of a detaed exampe, we compare our mutvarate estmator to the resutng estmator for the condtona MSEP f we gnore the correaton structure between ndvdua subportfoos as we as to the estmator for the condtona MSEP of the mutvarate chan-adder methods consdered by Braun (24) and Merz and Wüthrch (28) We obtan that n our exampe the predcton standard errors are substantay smaer n the mutvarate addtve method than n the mutvarate chan-adder cams reservng methods proposed by Braun (24) and Merz and Wüthrch (28) These fndngs may suggest that n the present case the mutvarate addtve method woud provde a better reserve estmate than the mutvarate chan-adder cams reservng method However, t s mportant to note that such a concuson woud be ony admssbe f we tested that the underyng mode assumptons of the addtve method are fufed Ths coud be done, for exampe, by the technques descrbed n Venter (998) Fnay, we want to emphasze that the condtona MSEP does not provde a compete pcture of the uncertanty assocated wth the predctor of the utmate cams of the tota portfoo Ths can ony be provded by the whoe predctve dstrbuton of the cams reserves [cf Engand and Verra (26) and Wüthrch and Merz (28)] Unfortunatey, n most cases one s not abe to cacuate the predctve dstrbuton anaytcay and one s forced to adopt numerca agorthms such as bootstrappng methods and Markov chan Monte Caro methods [cf Wüthrch and Merz (28)] Endowed wth the smuated predctve dstrbuton, one s not ony abe to cacuate estmates for the frst two moments of the cams reserves but one can aso derve predcton ntervas, quantes (eg, vaue at 5 CASUALTY ACTUARIAL SOCIETY VOLUME 3/ISSUE

Combining Chain-Ladder and Additive Loss Reserving Methods for Dependent Lines of Business

Combining Chain-Ladder and Additive Loss Reserving Methods for Dependent Lines of Business Combnng Chan-Laer an tve Loss Reservng Methos for Depenent Lnes of Busness by Mchae Merz an Maro V Wüthrch BSTRCT Often n non-fe nsurance, cam reserves are the argest poston on the abty se of the baance

More information

Development Pattern and Prediction Error for the Stochastic Bornhuetter-Ferguson Claims Reserving Method

Development Pattern and Prediction Error for the Stochastic Bornhuetter-Ferguson Claims Reserving Method Development Pattern and Predcton Error for the Stochastc Bornhuetter-Ferguson Clams Reservng Method Annna Saluz, Alos Gsler, Maro V. Wüthrch ETH Zurch ASTIN Colloquum Madrd, June 2011 Overvew 1 Notaton

More information

Image Classification Using EM And JE algorithms

Image Classification Using EM And JE algorithms Machne earnng project report Fa, 2 Xaojn Sh, jennfer@soe Image Cassfcaton Usng EM And JE agorthms Xaojn Sh Department of Computer Engneerng, Unversty of Caforna, Santa Cruz, CA, 9564 jennfer@soe.ucsc.edu

More information

COXREG. Estimation (1)

COXREG. Estimation (1) COXREG Cox (972) frst suggested the modes n whch factors reated to fetme have a mutpcatve effect on the hazard functon. These modes are caed proportona hazards (PH) modes. Under the proportona hazards

More information

Supplementary Material: Learning Structured Weight Uncertainty in Bayesian Neural Networks

Supplementary Material: Learning Structured Weight Uncertainty in Bayesian Neural Networks Shengyang Sun, Changyou Chen, Lawrence Carn Suppementary Matera: Learnng Structured Weght Uncertanty n Bayesan Neura Networks Shengyang Sun Changyou Chen Lawrence Carn Tsnghua Unversty Duke Unversty Duke

More information

Associative Memories

Associative Memories Assocatve Memores We consder now modes for unsupervsed earnng probems, caed auto-assocaton probems. Assocaton s the task of mappng patterns to patterns. In an assocatve memory the stmuus of an ncompete

More information

Nested case-control and case-cohort studies

Nested case-control and case-cohort studies Outne: Nested case-contro and case-cohort studes Ørnuf Borgan Department of Mathematcs Unversty of Oso NORBIS course Unversty of Oso 4-8 December 217 1 Radaton and breast cancer data Nested case contro

More information

Linear Approximation with Regularization and Moving Least Squares

Linear Approximation with Regularization and Moving Least Squares Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...

More information

MARKOV CHAIN AND HIDDEN MARKOV MODEL

MARKOV CHAIN AND HIDDEN MARKOV MODEL MARKOV CHAIN AND HIDDEN MARKOV MODEL JIAN ZHANG JIANZHAN@STAT.PURDUE.EDU Markov chan and hdden Markov mode are probaby the smpest modes whch can be used to mode sequenta data,.e. data sampes whch are not

More information

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4) I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes

More information

On the Power Function of the Likelihood Ratio Test for MANOVA

On the Power Function of the Likelihood Ratio Test for MANOVA Journa of Mutvarate Anayss 8, 416 41 (00) do:10.1006/jmva.001.036 On the Power Functon of the Lkehood Rato Test for MANOVA Dua Kumar Bhaumk Unversty of South Aabama and Unversty of Inos at Chcago and Sanat

More information

Neural network-based athletics performance prediction optimization model applied research

Neural network-based athletics performance prediction optimization model applied research Avaabe onne www.jocpr.com Journa of Chemca and Pharmaceutca Research, 04, 6(6):8-5 Research Artce ISSN : 0975-784 CODEN(USA) : JCPRC5 Neura networ-based athetcs performance predcton optmzaton mode apped

More information

Cyclic Codes BCH Codes

Cyclic Codes BCH Codes Cycc Codes BCH Codes Gaos Feds GF m A Gaos fed of m eements can be obtaned usng the symbos 0,, á, and the eements beng 0,, á, á, á 3 m,... so that fed F* s cosed under mutpcaton wth m eements. The operator

More information

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 30 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 2 Remedes for multcollnearty Varous technques have

More information

Composite Hypotheses testing

Composite Hypotheses testing Composte ypotheses testng In many hypothess testng problems there are many possble dstrbutons that can occur under each of the hypotheses. The output of the source s a set of parameters (ponts n a parameter

More information

A finite difference method for heat equation in the unbounded domain

A finite difference method for heat equation in the unbounded domain Internatona Conerence on Advanced ectronc Scence and Technoogy (AST 6) A nte derence method or heat equaton n the unbounded doman a Quan Zheng and Xn Zhao Coege o Scence North Chna nversty o Technoogy

More information

Research on Complex Networks Control Based on Fuzzy Integral Sliding Theory

Research on Complex Networks Control Based on Fuzzy Integral Sliding Theory Advanced Scence and Technoogy Letters Vo.83 (ISA 205), pp.60-65 http://dx.do.org/0.4257/ast.205.83.2 Research on Compex etworks Contro Based on Fuzzy Integra Sdng Theory Dongsheng Yang, Bngqng L, 2, He

More information

Chapter 7 Generalized and Weighted Least Squares Estimation. In this method, the deviation between the observed and expected values of

Chapter 7 Generalized and Weighted Least Squares Estimation. In this method, the deviation between the observed and expected values of Chapter 7 Generalzed and Weghted Least Squares Estmaton The usual lnear regresson model assumes that all the random error components are dentcally and ndependently dstrbuted wth constant varance. When

More information

Time-Varying Systems and Computations Lecture 6

Time-Varying Systems and Computations Lecture 6 Tme-Varyng Systems and Computatons Lecture 6 Klaus Depold 14. Januar 2014 The Kalman Flter The Kalman estmaton flter attempts to estmate the actual state of an unknown dscrete dynamcal system, gven nosy

More information

Journal of Multivariate Analysis

Journal of Multivariate Analysis Journa of Mutvarate Anayss 3 (04) 74 96 Contents sts avaabe at ScenceDrect Journa of Mutvarate Anayss journa homepage: www.esever.com/ocate/jmva Hgh-dmensona sparse MANOVA T. Tony Ca a, Yn Xa b, a Department

More information

Note 2. Ling fong Li. 1 Klein Gordon Equation Probablity interpretation Solutions to Klein-Gordon Equation... 2

Note 2. Ling fong Li. 1 Klein Gordon Equation Probablity interpretation Solutions to Klein-Gordon Equation... 2 Note 2 Lng fong L Contents Ken Gordon Equaton. Probabty nterpretaton......................................2 Soutons to Ken-Gordon Equaton............................... 2 2 Drac Equaton 3 2. Probabty nterpretaton.....................................

More information

Fall 2012 Analysis of Experimental Measurements B. Eisenstein/rev. S. Errede

Fall 2012 Analysis of Experimental Measurements B. Eisenstein/rev. S. Errede Fall 0 Analyss of Expermental easurements B. Esensten/rev. S. Errede We now reformulate the lnear Least Squares ethod n more general terms, sutable for (eventually extendng to the non-lnear case, and also

More information

Analysis of Bivariate Excess Losses

Analysis of Bivariate Excess Losses by Janong Ren ABSTRACT The concept of ecess osses s wey use n rensurance an retrospectve nsurance ratng The mathematcs reate to t has been stue etensvey n the property an casuaty actuara terature However,

More information

Chapter 11: Simple Linear Regression and Correlation

Chapter 11: Simple Linear Regression and Correlation Chapter 11: Smple Lnear Regresson and Correlaton 11-1 Emprcal Models 11-2 Smple Lnear Regresson 11-3 Propertes of the Least Squares Estmators 11-4 Hypothess Test n Smple Lnear Regresson 11-4.1 Use of t-tests

More information

QUARTERLY OF APPLIED MATHEMATICS

QUARTERLY OF APPLIED MATHEMATICS QUARTERLY OF APPLIED MATHEMATICS Voume XLI October 983 Number 3 DIAKOPTICS OR TEARING-A MATHEMATICAL APPROACH* By P. W. AITCHISON Unversty of Mantoba Abstract. The method of dakoptcs or tearng was ntroduced

More information

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U) Econ 413 Exam 13 H ANSWERS Settet er nndelt 9 deloppgaver, A,B,C, som alle anbefales å telle lkt for å gøre det ltt lettere å stå. Svar er gtt . Unfortunately, there s a prntng error n the hnt of

More information

Games of Threats. Elon Kohlberg Abraham Neyman. Working Paper

Games of Threats. Elon Kohlberg Abraham Neyman. Working Paper Games of Threats Elon Kohlberg Abraham Neyman Workng Paper 18-023 Games of Threats Elon Kohlberg Harvard Busness School Abraham Neyman The Hebrew Unversty of Jerusalem Workng Paper 18-023 Copyrght 2017

More information

Lower Bounding Procedures for the Single Allocation Hub Location Problem

Lower Bounding Procedures for the Single Allocation Hub Location Problem Lower Boundng Procedures for the Snge Aocaton Hub Locaton Probem Borzou Rostam 1,2 Chrstoph Buchhem 1,4 Fautät für Mathemat, TU Dortmund, Germany J. Faban Meer 1,3 Uwe Causen 1 Insttute of Transport Logstcs,

More information

Chat eld, C. and A.J.Collins, Introduction to multivariate analysis. Chapman & Hall, 1980

Chat eld, C. and A.J.Collins, Introduction to multivariate analysis. Chapman & Hall, 1980 MT07: Multvarate Statstcal Methods Mke Tso: emal mke.tso@manchester.ac.uk Webpage for notes: http://www.maths.manchester.ac.uk/~mkt/new_teachng.htm. Introducton to multvarate data. Books Chat eld, C. and

More information

Chapter 6. Rotations and Tensors

Chapter 6. Rotations and Tensors Vector Spaces n Physcs 8/6/5 Chapter 6. Rotatons and ensors here s a speca knd of near transformaton whch s used to transforms coordnates from one set of axes to another set of axes (wth the same orgn).

More information

However, since P is a symmetric idempotent matrix, of P are either 0 or 1 [Eigen-values

However, since P is a symmetric idempotent matrix, of P are either 0 or 1 [Eigen-values Fall 007 Soluton to Mdterm Examnaton STAT 7 Dr. Goel. [0 ponts] For the general lnear model = X + ε, wth uncorrelated errors havng mean zero and varance σ, suppose that the desgn matrx X s not necessarly

More information

Uncertainty Specification and Propagation for Loss Estimation Using FOSM Methods

Uncertainty Specification and Propagation for Loss Estimation Using FOSM Methods Uncertanty Specfcaton and Propagaton for Loss Estmaton Usng FOSM Methods J.W. Baer and C.A. Corne Dept. of Cv and Envronmenta Engneerng, Stanford Unversty, Stanford, CA 94305-400 Keywords: Sesmc, oss estmaton,

More information

Multispectral Remote Sensing Image Classification Algorithm Based on Rough Set Theory

Multispectral Remote Sensing Image Classification Algorithm Based on Rough Set Theory Proceedngs of the 2009 IEEE Internatona Conference on Systems Man and Cybernetcs San Antono TX USA - October 2009 Mutspectra Remote Sensng Image Cassfcaton Agorthm Based on Rough Set Theory Yng Wang Xaoyun

More information

Computation of Higher Order Moments from Two Multinomial Overdispersion Likelihood Models

Computation of Higher Order Moments from Two Multinomial Overdispersion Likelihood Models Computaton of Hgher Order Moments from Two Multnomal Overdsperson Lkelhood Models BY J. T. NEWCOMER, N. K. NEERCHAL Department of Mathematcs and Statstcs, Unversty of Maryland, Baltmore County, Baltmore,

More information

A DIMENSION-REDUCTION METHOD FOR STOCHASTIC ANALYSIS SECOND-MOMENT ANALYSIS

A DIMENSION-REDUCTION METHOD FOR STOCHASTIC ANALYSIS SECOND-MOMENT ANALYSIS A DIMESIO-REDUCTIO METHOD FOR STOCHASTIC AALYSIS SECOD-MOMET AALYSIS S. Rahman Department of Mechanca Engneerng and Center for Computer-Aded Desgn The Unversty of Iowa Iowa Cty, IA 52245 June 2003 OUTLIE

More information

Globally Optimal Multisensor Distributed Random Parameter Matrices Kalman Filtering Fusion with Applications

Globally Optimal Multisensor Distributed Random Parameter Matrices Kalman Filtering Fusion with Applications Sensors 2008, 8, 8086-8103; DOI: 10.3390/s8128086 OPEN ACCESS sensors ISSN 1424-8220 www.mdp.com/journa/sensors Artce Gobay Optma Mutsensor Dstrbuted Random Parameter Matrces Kaman Fterng Fuson wth Appcatons

More information

1. Inference on Regression Parameters a. Finding Mean, s.d and covariance amongst estimates. 2. Confidence Intervals and Working Hotelling Bands

1. Inference on Regression Parameters a. Finding Mean, s.d and covariance amongst estimates. 2. Confidence Intervals and Working Hotelling Bands Content. Inference on Regresson Parameters a. Fndng Mean, s.d and covarance amongst estmates.. Confdence Intervals and Workng Hotellng Bands 3. Cochran s Theorem 4. General Lnear Testng 5. Measures of

More information

Uncertainty as the Overlap of Alternate Conditional Distributions

Uncertainty as the Overlap of Alternate Conditional Distributions Uncertanty as the Overlap of Alternate Condtonal Dstrbutons Olena Babak and Clayton V. Deutsch Centre for Computatonal Geostatstcs Department of Cvl & Envronmental Engneerng Unversty of Alberta An mportant

More information

Distributed Moving Horizon State Estimation of Nonlinear Systems. Jing Zhang

Distributed Moving Horizon State Estimation of Nonlinear Systems. Jing Zhang Dstrbuted Movng Horzon State Estmaton of Nonnear Systems by Jng Zhang A thess submtted n parta fufment of the requrements for the degree of Master of Scence n Chemca Engneerng Department of Chemca and

More information

ON AUTOMATIC CONTINUITY OF DERIVATIONS FOR BANACH ALGEBRAS WITH INVOLUTION

ON AUTOMATIC CONTINUITY OF DERIVATIONS FOR BANACH ALGEBRAS WITH INVOLUTION European Journa of Mathematcs and Computer Scence Vo. No. 1, 2017 ON AUTOMATC CONTNUTY OF DERVATONS FOR BANACH ALGEBRAS WTH NVOLUTON Mohamed BELAM & Youssef T DL MATC Laboratory Hassan Unversty MORO CCO

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method Appled Mathematcal Scences, Vol. 7, 0, no. 47, 07-0 HIARI Ltd, www.m-hkar.com Comparson of the Populaton Varance Estmators of -Parameter Exponental Dstrbuton Based on Multple Crtera Decson Makng Method

More information

Exponential Type Product Estimator for Finite Population Mean with Information on Auxiliary Attribute

Exponential Type Product Estimator for Finite Population Mean with Information on Auxiliary Attribute Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 193-9466 Vol. 10, Issue 1 (June 015), pp. 106-113 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) Exponental Tpe Product Estmator

More information

Xin Li Department of Information Systems, College of Business, City University of Hong Kong, Hong Kong, CHINA

Xin Li Department of Information Systems, College of Business, City University of Hong Kong, Hong Kong, CHINA RESEARCH ARTICLE MOELING FIXE OS BETTING FOR FUTURE EVENT PREICTION Weyun Chen eartment of Educatona Informaton Technoogy, Facuty of Educaton, East Chna Norma Unversty, Shangha, CHINA {weyun.chen@qq.com}

More information

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

Lecture 3 Stat102, Spring 2007

Lecture 3 Stat102, Spring 2007 Lecture 3 Stat0, Sprng 007 Chapter 3. 3.: Introducton to regresson analyss Lnear regresson as a descrptve technque The least-squares equatons Chapter 3.3 Samplng dstrbuton of b 0, b. Contnued n net lecture

More information

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur Analyss of Varance and Desgn of Experment-I MODULE VII LECTURE - 3 ANALYSIS OF COVARIANCE Dr Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur Any scentfc experment s performed

More information

3. Stress-strain relationships of a composite layer

3. Stress-strain relationships of a composite layer OM PO I O U P U N I V I Y O F W N ompostes ourse 8-9 Unversty of wente ng. &ech... tress-stran reatonshps of a composte ayer - Laurent Warnet & emo Aerman.. tress-stran reatonshps of a composte ayer Introducton

More information

Quantum Runge-Lenz Vector and the Hydrogen Atom, the hidden SO(4) symmetry

Quantum Runge-Lenz Vector and the Hydrogen Atom, the hidden SO(4) symmetry Quantum Runge-Lenz ector and the Hydrogen Atom, the hdden SO(4) symmetry Pasca Szrftgser and Edgardo S. Cheb-Terrab () Laboratore PhLAM, UMR CNRS 85, Unversté Le, F-59655, France () Mapesoft Let's consder

More information

Predictive Analytics : QM901.1x Prof U Dinesh Kumar, IIMB. All Rights Reserved, Indian Institute of Management Bangalore

Predictive Analytics : QM901.1x Prof U Dinesh Kumar, IIMB. All Rights Reserved, Indian Institute of Management Bangalore Sesson Outlne Introducton to classfcaton problems and dscrete choce models. Introducton to Logstcs Regresson. Logstc functon and Logt functon. Maxmum Lkelhood Estmator (MLE) for estmaton of LR parameters.

More information

Numerical integration in more dimensions part 2. Remo Minero

Numerical integration in more dimensions part 2. Remo Minero Numerca ntegraton n more dmensons part Remo Mnero Outne The roe of a mappng functon n mutdmensona ntegraton Gauss approach n more dmensons and quadrature rues Crtca anass of acceptabt of a gven quadrature

More information

j) = 1 (note sigma notation) ii. Continuous random variable (e.g. Normal distribution) 1. density function: f ( x) 0 and f ( x) dx = 1

j) = 1 (note sigma notation) ii. Continuous random variable (e.g. Normal distribution) 1. density function: f ( x) 0 and f ( x) dx = 1 Random varables Measure of central tendences and varablty (means and varances) Jont densty functons and ndependence Measures of assocaton (covarance and correlaton) Interestng result Condtonal dstrbutons

More information

L-Edge Chromatic Number Of A Graph

L-Edge Chromatic Number Of A Graph IJISET - Internatona Journa of Innovatve Scence Engneerng & Technoogy Vo. 3 Issue 3 March 06. ISSN 348 7968 L-Edge Chromatc Number Of A Graph Dr.R.B.Gnana Joth Assocate Professor of Mathematcs V.V.Vannaperuma

More information

A Robust Method for Calculating the Correlation Coefficient

A Robust Method for Calculating the Correlation Coefficient A Robust Method for Calculatng the Correlaton Coeffcent E.B. Nven and C. V. Deutsch Relatonshps between prmary and secondary data are frequently quantfed usng the correlaton coeffcent; however, the tradtonal

More information

MACHINE APPLIED MACHINE LEARNING LEARNING. Gaussian Mixture Regression

MACHINE APPLIED MACHINE LEARNING LEARNING. Gaussian Mixture Regression 11 MACHINE APPLIED MACHINE LEARNING LEARNING MACHINE LEARNING Gaussan Mture Regresson 22 MACHINE APPLIED MACHINE LEARNING LEARNING Bref summary of last week s lecture 33 MACHINE APPLIED MACHINE LEARNING

More information

Lecture 6: Introduction to Linear Regression

Lecture 6: Introduction to Linear Regression Lecture 6: Introducton to Lnear Regresson An Manchakul amancha@jhsph.edu 24 Aprl 27 Lnear regresson: man dea Lnear regresson can be used to study an outcome as a lnear functon of a predctor Example: 6

More information

Chapter 9: Statistical Inference and the Relationship between Two Variables

Chapter 9: Statistical Inference and the Relationship between Two Variables Chapter 9: Statstcal Inference and the Relatonshp between Two Varables Key Words The Regresson Model The Sample Regresson Equaton The Pearson Correlaton Coeffcent Learnng Outcomes After studyng ths chapter,

More information

Econ Statistical Properties of the OLS estimator. Sanjaya DeSilva

Econ Statistical Properties of the OLS estimator. Sanjaya DeSilva Econ 39 - Statstcal Propertes of the OLS estmator Sanjaya DeSlva September, 008 1 Overvew Recall that the true regresson model s Y = β 0 + β 1 X + u (1) Applyng the OLS method to a sample of data, we estmate

More information

CS286r Assign One. Answer Key

CS286r Assign One. Answer Key CS286r Assgn One Answer Key 1 Game theory 1.1 1.1.1 Let off-equlbrum strateges also be that people contnue to play n Nash equlbrum. Devatng from any Nash equlbrum s a weakly domnated strategy. That s,

More information

Numerical Investigation of Power Tunability in Two-Section QD Superluminescent Diodes

Numerical Investigation of Power Tunability in Two-Section QD Superluminescent Diodes Numerca Investgaton of Power Tunabty n Two-Secton QD Superumnescent Dodes Matta Rossett Paoo Bardea Ivo Montrosset POLITECNICO DI TORINO DELEN Summary 1. A smpfed mode for QD Super Lumnescent Dodes (SLD)

More information

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin Proceedngs of the 007 Wnter Smulaton Conference S G Henderson, B Bller, M-H Hseh, J Shortle, J D Tew, and R R Barton, eds LOW BIAS INTEGRATED PATH ESTIMATORS James M Calvn Department of Computer Scence

More information

STAT 3008 Applied Regression Analysis

STAT 3008 Applied Regression Analysis STAT 3008 Appled Regresson Analyss Tutoral : Smple Lnear Regresson LAI Chun He Department of Statstcs, The Chnese Unversty of Hong Kong 1 Model Assumpton To quantfy the relatonshp between two factors,

More information

Lecture Notes on Linear Regression

Lecture Notes on Linear Regression Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume

More information

Appendix for An Efficient Ascending-Bid Auction for Multiple Objects: Comment For Online Publication

Appendix for An Efficient Ascending-Bid Auction for Multiple Objects: Comment For Online Publication Appendx for An Effcent Ascendng-Bd Aucton for Mutpe Objects: Comment For Onne Pubcaton Norak Okamoto The foowng counterexampe shows that sncere bddng by a bdders s not aways an ex post perfect equbrum

More information

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 31 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 6. Rdge regresson The OLSE s the best lnear unbased

More information

Chapter 2 - The Simple Linear Regression Model S =0. e i is a random error. S β2 β. This is a minimization problem. Solution is a calculus exercise.

Chapter 2 - The Simple Linear Regression Model S =0. e i is a random error. S β2 β. This is a minimization problem. Solution is a calculus exercise. Chapter - The Smple Lnear Regresson Model The lnear regresson equaton s: where y + = β + β e for =,..., y and are observable varables e s a random error How can an estmaton rule be constructed for the

More information

A Dissimilarity Measure Based on Singular Value and Its Application in Incremental Discounting

A Dissimilarity Measure Based on Singular Value and Its Application in Incremental Discounting A Dssmarty Measure Based on Snguar Vaue and Its Appcaton n Incrementa Dscountng KE Xaou Department of Automaton, Unversty of Scence and Technoogy of Chna, Hefe, Chna Ema: kxu@ma.ustc.edu.cn MA Lyao Department

More information

e i is a random error

e i is a random error Chapter - The Smple Lnear Regresson Model The lnear regresson equaton s: where + β + β e for,..., and are observable varables e s a random error How can an estmaton rule be constructed for the unknown

More information

β0 + β1xi. You are interested in estimating the unknown parameters β

β0 + β1xi. You are interested in estimating the unknown parameters β Revsed: v3 Ordnar Least Squares (OLS): Smple Lnear Regresson (SLR) Analtcs The SLR Setup Sample Statstcs Ordnar Least Squares (OLS): FOCs and SOCs Back to OLS and Sample Statstcs Predctons (and Resduals)

More information

APPENDIX A Some Linear Algebra

APPENDIX A Some Linear Algebra APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,

More information

2016 Wiley. Study Session 2: Ethical and Professional Standards Application

2016 Wiley. Study Session 2: Ethical and Professional Standards Application 6 Wley Study Sesson : Ethcal and Professonal Standards Applcaton LESSON : CORRECTION ANALYSIS Readng 9: Correlaton and Regresson LOS 9a: Calculate and nterpret a sample covarance and a sample correlaton

More information

β0 + β1xi. You are interested in estimating the unknown parameters β

β0 + β1xi. You are interested in estimating the unknown parameters β Ordnary Least Squares (OLS): Smple Lnear Regresson (SLR) Analytcs The SLR Setup Sample Statstcs Ordnary Least Squares (OLS): FOCs and SOCs Back to OLS and Sample Statstcs Predctons (and Resduals) wth OLS

More information

Research Article H Estimates for Discrete-Time Markovian Jump Linear Systems

Research Article H Estimates for Discrete-Time Markovian Jump Linear Systems Mathematca Probems n Engneerng Voume 213 Artce ID 945342 7 pages http://dxdoorg/11155/213/945342 Research Artce H Estmates for Dscrete-Tme Markovan Jump Lnear Systems Marco H Terra 1 Gdson Jesus 2 and

More information

9. Binary Dependent Variables

9. Binary Dependent Variables 9. Bnar Dependent Varables 9. Homogeneous models Log, prob models Inference Tax preparers 9.2 Random effects models 9.3 Fxed effects models 9.4 Margnal models and GEE Appendx 9A - Lkelhood calculatons

More information

Deriving the Dual. Prof. Bennett Math of Data Science 1/13/06

Deriving the Dual. Prof. Bennett Math of Data Science 1/13/06 Dervng the Dua Prof. Bennett Math of Data Scence /3/06 Outne Ntty Grtty for SVM Revew Rdge Regresson LS-SVM=KRR Dua Dervaton Bas Issue Summary Ntty Grtty Need Dua of w, b, z w 2 2 mn st. ( x w ) = C z

More information

Chapter 12 Analysis of Covariance

Chapter 12 Analysis of Covariance Chapter Analyss of Covarance Any scentfc experment s performed to know somethng that s unknown about a group of treatments and to test certan hypothess about the correspondng treatment effect When varablty

More information

/ n ) are compared. The logic is: if the two

/ n ) are compared. The logic is: if the two STAT C141, Sprng 2005 Lecture 13 Two sample tests One sample tests: examples of goodness of ft tests, where we are testng whether our data supports predctons. Two sample tests: called as tests of ndependence

More information

Approximate merging of a pair of BeÂzier curves

Approximate merging of a pair of BeÂzier curves COMPUTER-AIDED DESIGN Computer-Aded Desgn 33 (1) 15±136 www.esever.com/ocate/cad Approxmate mergng of a par of BeÂzer curves Sh-Mn Hu a,b, *, Rou-Feng Tong c, Tao Ju a,b, Ja-Guang Sun a,b a Natona CAD

More information

The Insurance Risk in the SST and in Solvency II: Modeling and Parameter Estimators

The Insurance Risk in the SST and in Solvency II: Modeling and Parameter Estimators The Insurance Rsk n the SST and n Solvency II: Modelng and Parameter Estmators ASTIN Colloquum 1-4 June 009, Helsnk Alos Gsler Introducton SST and Solvency II common goal: to nstall a rsk based solvency

More information

Estimation: Part 2. Chapter GREG estimation

Estimation: Part 2. Chapter GREG estimation Chapter 9 Estmaton: Part 2 9. GREG estmaton In Chapter 8, we have seen that the regresson estmator s an effcent estmator when there s a lnear relatonshp between y and x. In ths chapter, we generalzed the

More information

Properties of Least Squares

Properties of Least Squares Week 3 3.1 Smple Lnear Regresson Model 3. Propertes of Least Squares Estmators Y Y β 1 + β X + u weekly famly expendtures X weekly famly ncome For a gven level of x, the expected level of food expendtures

More information

Explaining the Stein Paradox

Explaining the Stein Paradox Explanng the Sten Paradox Kwong Hu Yung 1999/06/10 Abstract Ths report offers several ratonale for the Sten paradox. Sectons 1 and defnes the multvarate normal mean estmaton problem and ntroduces Sten

More information

where I = (n x n) diagonal identity matrix with diagonal elements = 1 and off-diagonal elements = 0; and σ 2 e = variance of (Y X).

where I = (n x n) diagonal identity matrix with diagonal elements = 1 and off-diagonal elements = 0; and σ 2 e = variance of (Y X). 11.4.1 Estmaton of Multple Regresson Coeffcents In multple lnear regresson, we essentally solve n equatons for the p unnown parameters. hus n must e equal to or greater than p and n practce n should e

More information

2.3 Nilpotent endomorphisms

2.3 Nilpotent endomorphisms s a block dagonal matrx, wth A Mat dm U (C) In fact, we can assume that B = B 1 B k, wth B an ordered bass of U, and that A = [f U ] B, where f U : U U s the restrcton of f to U 40 23 Nlpotent endomorphsms

More information

The Application of BP Neural Network principal component analysis in the Forecasting the Road Traffic Accident

The Application of BP Neural Network principal component analysis in the Forecasting the Road Traffic Accident ICTCT Extra Workshop, Bejng Proceedngs The Appcaton of BP Neura Network prncpa component anayss n Forecastng Road Traffc Accdent He Mng, GuoXucheng &LuGuangmng Transportaton Coege of Souast Unversty 07

More information

Some results on a cross-section in the tensor bundle

Some results on a cross-section in the tensor bundle Hacettepe Journa of Matematcs and Statstcs Voume 43 3 214, 391 397 Some resuts on a cross-secton n te tensor bunde ydın Gezer and Murat tunbas bstract Te present paper s devoted to some resuts concernng

More information

F statistic = s2 1 s 2 ( F for Fisher )

F statistic = s2 1 s 2 ( F for Fisher ) Stat 4 ANOVA Analyss of Varance /6/04 Comparng Two varances: F dstrbuton Typcal Data Sets One way analyss of varance : example Notaton for one way ANOVA Comparng Two varances: F dstrbuton We saw that the

More information

Negative Binomial Regression

Negative Binomial Regression STATGRAPHICS Rev. 9/16/2013 Negatve Bnomal Regresson Summary... 1 Data Input... 3 Statstcal Model... 3 Analyss Summary... 4 Analyss Optons... 7 Plot of Ftted Model... 8 Observed Versus Predcted... 10 Predctons...

More information

Polite Water-filling for Weighted Sum-rate Maximization in MIMO B-MAC Networks under. Multiple Linear Constraints

Polite Water-filling for Weighted Sum-rate Maximization in MIMO B-MAC Networks under. Multiple Linear Constraints 2011 IEEE Internatona Symposum on Informaton Theory Proceedngs Pote Water-fng for Weghted Sum-rate Maxmzaton n MIMO B-MAC Networks under Mutpe near Constrants An u 1, Youjan u 2, Vncent K. N. au 3, Hage

More information

ECONOMICS 351*-A Mid-Term Exam -- Fall Term 2000 Page 1 of 13 pages. QUEEN'S UNIVERSITY AT KINGSTON Department of Economics

ECONOMICS 351*-A Mid-Term Exam -- Fall Term 2000 Page 1 of 13 pages. QUEEN'S UNIVERSITY AT KINGSTON Department of Economics ECOOMICS 35*-A Md-Term Exam -- Fall Term 000 Page of 3 pages QUEE'S UIVERSITY AT KIGSTO Department of Economcs ECOOMICS 35* - Secton A Introductory Econometrcs Fall Term 000 MID-TERM EAM ASWERS MG Abbott

More information

x = , so that calculated

x = , so that calculated Stat 4, secton Sngle Factor ANOVA notes by Tm Plachowsk n chapter 8 we conducted hypothess tests n whch we compared a sngle sample s mean or proporton to some hypotheszed value Chapter 9 expanded ths to

More information

STAT 511 FINAL EXAM NAME Spring 2001

STAT 511 FINAL EXAM NAME Spring 2001 STAT 5 FINAL EXAM NAME Sprng Instructons: Ths s a closed book exam. No notes or books are allowed. ou may use a calculator but you are not allowed to store notes or formulas n the calculator. Please wrte

More information

Lower bounds for the Crossing Number of the Cartesian Product of a Vertex-transitive Graph with a Cycle

Lower bounds for the Crossing Number of the Cartesian Product of a Vertex-transitive Graph with a Cycle Lower bounds for the Crossng Number of the Cartesan Product of a Vertex-transtve Graph wth a Cyce Junho Won MIT-PRIMES December 4, 013 Abstract. The mnmum number of crossngs for a drawngs of a gven graph

More information

Continuous vs. Discrete Goods

Continuous vs. Discrete Goods CE 651 Transportaton Economcs Charsma Choudhury Lecture 3-4 Analyss of Demand Contnuous vs. Dscrete Goods Contnuous Goods Dscrete Goods x auto 1 Indfference u curves 3 u u 1 x 1 0 1 bus Outlne Data Modelng

More information

More metrics on cartesian products

More metrics on cartesian products More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of

More information

β0 + β1xi and want to estimate the unknown

β0 + β1xi and want to estimate the unknown SLR Models Estmaton Those OLS Estmates Estmators (e ante) v. estmates (e post) The Smple Lnear Regresson (SLR) Condtons -4 An Asde: The Populaton Regresson Functon B and B are Lnear Estmators (condtonal

More information

UNIVERSITY OF TORONTO Faculty of Arts and Science. December 2005 Examinations STA437H1F/STA1005HF. Duration - 3 hours

UNIVERSITY OF TORONTO Faculty of Arts and Science. December 2005 Examinations STA437H1F/STA1005HF. Duration - 3 hours UNIVERSITY OF TORONTO Faculty of Arts and Scence December 005 Examnatons STA47HF/STA005HF Duraton - hours AIDS ALLOWED: (to be suppled by the student) Non-programmable calculator One handwrtten 8.5'' x

More information