A HIERARCHICAL KALMAN FILTER

Size: px
Start display at page:

Download "A HIERARCHICAL KALMAN FILTER"

Transcription

1 A HIERARCHICAL KALMAN FILER Greg aylor aylor Fry Consulng Acuares Level 8, 3 Clarence Sree Sydney NSW Ausrala Professoral Assocae, Cenre for Acuaral Sudes Faculy of Economcs and Commerce Unversy of Melbourne Parkvlle VIC 35 Ausrala Phone: Fax: greg.aylor@aylorfry.com.au March 9

2 A Herarchcal Kalman Fler able of Conens Summary.... Inroducon.... Credbly esmaon Sandard (mul-dmensonal) credbly Kalman fler Herarchcal credbly A herarchcal Kalman fler Numercal example Concluson... 4 References...6

3 A Herarchcal Kalman Fler Summary Sund s herarchcal credbly model s generalsed o a dynamc form, e a form n whch parameers are assumed o evolve over me. hs s done by supermposng a Kalman fler on Sund s model. In he process s shown ha boh Sund s model and he Kalman fler may be derved as drec consequences of Hachemeser s credbly regresson model. A numercal example n Secon 7 llusraes he applcaon of he herarchcal Kalman fler o a herarchy of occupaonal groups. he example shows how parameer esmaes produced by he fler rack he rue values of parameers beer han he esmaes from a sac model when hose parameers evolve over me. Keywords: credbly, Hachemeser model, herarchcal credbly, Kalman fler.

4 A Herarchcal Kalman Fler. Inroducon Herarchcal credbly was nroduced by aylor (979). he parameers were scalars, as were he observaons on hem. hs framework was generalsed by Sund (979, 98) o nclude vecor parameers and observaons. In Sund s framework a lnear regresson model s placed a each level of he herarchy. he model s sac n he sense ha all parameers relae o jus a sngle epoch or me nerval. If, n fac, daa were colleced over mulple me perods, hey would be pooled on he assumpon ha he parameers remaned consan over me. Such poolng of daa renders he model sluggsh n s response o evolvng parameers. One srucure ha recognses parameer evoluon explcly s he Kalman fler (Kalman, 96; Jazwnsk, 97). he purpose of he presen paper s o adap Sund s herarchcal model o a dynamc form by supermposng a Kalman fler on. I s ancpaed ha he parameer esmaes yelded by he dynamc model wll rack he evolvng parameers beer han her sac counerpars. he Kalman fler s derved n Secon 4 and Sudn s resuls n Secon 5. s shown n he process ha boh can be vewed as consequences of Hachemeser s (975) mul-dmensonal credbly model. he herarchcal form of he Kalman fler s derved n Secon 6. Is use s llusraed n Secon 7, where s appled o he orgnal conex of herarchcal credbly (aylor, 979). In hs case he herarchy s a ree of occupaonal groups n a workers compensaon seng.. Credbly esmaon Le m be an unknown vecor random varable, and le x be an observable vecor random varable. A lnear credbly esmaor (or jus credbly esmaor) of m (based on x) s he lnear funcon L(x) of x whch mnmses E L( x) m. he lnear credbly esmaor L(x) s unbased for m. 3. Sandard (mul-dmensonal) credbly Le θ be an unknown random parameer and x an observable random vecor. Assume he followng: (3.) x= Yb( θ ) +ξ for some vecor-valued funcon b (.) Y, and cenred random vecor ξ ; (3.) θ and ξ are ndependen., deermnsc marx

5 A Herarchcal Kalman Fler 3 Le β= E b( θ) (3.) Λ= θ (3.) Var b( ) [ ] Var [ ] Φ= Var x θ = ξ (3.3) hs s he Hachemeser (975) framework, and he lnear credbly esmaor of b( θ ) s ( Z) Zbˆ b = β+ (3.4) where ( Z =ΛY Φ Y +ΛY Φ Y ) (3.5) and ˆb s he classcal regresson esmaor ( ) ˆ = Φ Φ x (3.6) b Y Y Y Here he upper denoes marx ransposon. Exenson. Relax he condon ha ξ be cenred. hen resuls (3.4) and (3.5) connue o hold bu wh (3.6) replaced as follows: ( ) [ ] ˆ b= Y Φ Y Y Φ x E ξ (3.6a) Proof. Wre u = x E[ ξ ] hen, by Assumpon (3.), u = Yb( θ ) +η where s he cenred random vecor η=ξ E[ ξ ] hen resuls (3.4) (3.6) hold wh x replaced by u.

6 A Herarchcal Kalman Fler 4 4. Kalman fler 4. Sandard Kalman fler Le b, =, 3, ec be an unknown vecor parameer ha evolves randomly over me, and le x be an observable random vecor, where b and x sasfy he followng relaons. (4.) + + b = A b +ε +, =,, ec where A + s a deermnsc square marx and ε + a cenred random vecor, and he sysem s naed by a deermnsc vecor b. (4.) x = Y b + ξ, =,, ec where Y s a deermnsc marx and (4.) Each of he ses { ε } and { } ξ a cenred random vecor. ξ s muually sochascally ndependen, and furher he wo ses are sochascally ndependen of each oher. Le s b be he (unbased) credbly esmaor of b Specfcally, on accoun of Assumpon (4.), defne s based on {,,... } x x x. + = A + b (4.) b I s sraghforward o show ha on { x,... x }. s he credbly esmaor of b + based b + Le Φ = Var x b = Var ξ (4.) Ψ = Var ε (4.3) E b b (4.4) Λ = he credbly esmaor Λ. s he lnear funcon of x,, x - ha mnmses b By Assumpon (4.) and (4.), ( ) ( A) E Ab Ab Λ = +ε = AVar b + Ψ (4.5)

7 A Herarchcal Kalman Fler 5 where use has been made of Assumpon (4.) and he fac ha depend on ε. b does no he mahemacs of fndng he credbly esmaor Secon 3 bu wh he followng replacemens: b s he same as n x x Y Y b θ b ( ) ξ ξ β b b b he parallel of Assumpon (3.) here s ndependence of b and guaraneed by Assumpon (4.). ξ, whch s I hen follows ha Λ Λ [from (3.) and (4.4)] Φ Φ [from (3.3) and (4.)] Subsuon of hese replacemens no (3.4) and (3.5) yelds ( ) b = K b + K b ˆ (4.6) where ( ) ( ) = + = + (4.7) K L L L L L ( ) ( ) =Λ Y Φ Y (4.8) ( ) ( ) ( ) ( ) bˆ = Y Φ Y Y Φ x (4.9) Snce ˆ b depends on jus x and does no, (4.6) gves b ( ) ( ) ˆ ( ) Var b K Var b = K + K Var b K (4.) = ( K ) Λ ( K ) + K ( Y ) ( Φ ) Y ( K ) o smplfy hs, emporarly suppress he upper. Subsue (4.7) and (4.8) no (4.), nong ha K = ( + L) -, o oban

8 A Herarchcal Kalman Fler 6 ( ) { }( ) Var b = + L Λ+ L Y Φ Y L + L ( L) = + Λ Wh he upper resored ( K ) (4.) Var b = Λ Equaons (4.), (4.5) (4.9) and (4.) collecvely show how o proceed from b o b and, n dong so, consue a sngle eraon of he Kalman fler (Kalman, 96; Jazwnsk, 97). 4. Exenson + Relax he requremen n Assumpon (4.) ha ε be cenred. hen he resuls of Secon 4. connue o hold wh he sngle change ha (4.) s replaced by he followng: = + ε (4.a) b A b E Proof. Defne he sequence a, =,, ec by he recurson a = A a + E ε a = (4.) hen defne c b a = (4.3) By (4.3) and Assumpon (4.) c A c = +η (4.4) where + η s he cenred random vecor η =ε E ε (4.5) By (4.3) and Assumpon (4.), x = Y c + ζ (4.6) wh ζ = Ya+ξ (4.7)

9 A Herarchcal Kalman Fler 7 Noe ha (4.4) and (4.6) are of he same form as Assumpons (4.) and (4.) wh he followng replacemens: b c ε η ξ ζ and where ζ s no cenred. he dervaon of he Kalman fler for hs sysem s by exacly he same argumen as gven above n relaon o Assumpons (4.) and (4.) excep ha appeal s now made o he exenson n Secon 3 n order o accommodae he non-cenredness of ζ. Noe ha a s deermnsc, and so he varance quanes Φ, Ψ, Λ are all lef unchanged by he above replacemens. he only change requred o he Kalman fler s he replacemen of (4.) by (4.a). 5. Herarchcal credbly Consder he herarchcal credbly model of Sund (979, 98). Here he unknown parameer and observable random vecor x of Secon 3 are replaced by unknown random parameers θ and observable random vecors x, =,,,n, wh denong he level whn a mul-level herarchy. he model s defned by he followng assumpons. (5.) For each =,, n, x θ,..., θ = Yb( θ,..., θ ) +ξ( θ,..., θ ) b., deermnsc marx, and for some vecor-valued funcon ( ) random vecor ξ (,..., θ θ ) sasfyng E ξ ( θ,..., θ ) θ,..., θ =. (5.) For each =,, n, b θ..., θ = W b θ..., θ +η θ..., θ for some deermnsc ( ) ( ) ( ) W and random vecor (..., ) η ( θ,..., θ ) θ,..., θ = marx E. For he case =, b θ =β+η θ ) ( ) ( η θ θ sasfyng wh β some known parameer and ( ) (5.) ξ ( θ θ ) and ( ),..., E η θ =. η,..., θ θ do no depend on θ j, j >. Y

10 A Herarchcal Kalman Fler 8 (5.v) (a) ξ (,..., θ θ ) and ξ (,..., j θ θ j ) are condonally ndependen for j >, gven θ,..., θ η,..., θ θ and η (,..., j θ θ j ) are condonally ndependen for (b) ( ) (c) j >, gven θ,..., θ ξ (,..., θ θ ) and η (,..., j θ θ j ) are ndependen for all,j. Remark 5. I follows from hese assumpons ha (5.v) x and θ j are ndependen for j > (5.v) For each =,..., n, ( x,..., x ) and (,..., x x n ) ndependen gven θ,..., θ. + are condonally he assumpons (5.) (5.v) are n fac sronger han hose of Sund, who assumed (5.v) and (5.v) n place of (5.) and (5.v). Le ( ) Λ = EVar b θ,..., θ θ,..., θ, =,,..., n (5.) [ ] Φ = EVar x θ,..., θ, =,,..., n (5.) Now defne x = x,..., x, =,,..., n ( n ) (5.3) and defne he vecors ξ, η smlarly. When hese vecors are wren as condoned by θ, each componen wll be undersood o be condoned by all relevan nformaon, e n x, x wll be condoned by θ,..., θ ; x + by θ,..., θ+, ec. Equvalenly, each componen s condoned by θ,..., θn. By Assumpons (5.) and (5.), parcularly he lneary of he relaons here, x θ= Y b θ,..., θ + Nη θ+ξ θ (5.4) ( ) for some marces Y and N. o calculae hese marces, noe ha x θ= Y b θ,..., θ + N η θ+ξ θ ( + ) ( + ) + ( + ) + ( + ) ( + ) = Y( ) W b( θ,..., θ ) + Y( ), N ( ) ( +) η θ+ξ θ (5.5) by Assumpon (5.).

11 A Herarchcal Kalman Fler 9 Combnng Assumpon (5.) wh (5.5) and equang he resul o (5.4) yelds Y Y = Y( ) W +, =,,..., n (5.6) N = Y( ) N + +, =,,..., n (5.7) where he zero sub-marces have he same numbers of rows as x. he recurson (5.6) and (5.7) s naed wh Y = Y N = ( ) n, n n (5.8) obaned from (5.4) wh = n. Denoe Φ = EVar,..., x θ θ Π = EVar x,..., θ θ (5.9) (5.) Now [ θ θ ] EVar x,..., Φ = E x( ) θ,..., θ + (5.) Φ = Π ( + ), =,,..., n where he ndependence assumpon (5.v) has been used n he frs equaly. he recurson (5.) s naed wh Φ ( n) =Φ n (5.) Π he marces may be evaluaed by reference o he frs equaly n (5.5) wh + replaced by, whence {,...,,..., } E E x Y b ( ) Π = θ θ θ θ θ,..., θ θ where s convenen for he me beng o show explcly he varaes wh respec o whch expecaons are aken.

12 A Herarchcal Kalman Fler he las relaon may be expanded as follows: {{,...,,..., } Π = E ( ),..., E θ θ x E Y b θ θ θ θ θ { }} + Y E b θ,..., θ θ,..., θ E b θ,..., θ θ,..., θ θ,..., θ ( ) ( ) =Φ + Y Λ Y by (5.) and (5.9). Le (5.3) b denoe he credbly esmaor of b (,..., θ θ ) based on x. hs means ha b s he lnear funcon of x ha mnmses { ( θ,..., θ ) θ,..., θ} E E θ θ b b... Sund (98) shows ha, by Assumpon (5.), W b s he credbly esmaor of E b ( θ,..., θ ) θ,..., θ, and ha, by Assumpon (5.v), b depends on jus he x par of x. he mahemacs of fndng hs esmaor s once agan he same as n Secon 3 bu wh he followng replacemens: x x Y Y b θ b,..., θ θ ξ ξ β W b b b ( ) ( ) he expecaon n he mnmsaon creron s also replaced. he objecve funcon E L( x) m s now replaced by E E b b { (,..., ) θ θ θ,..., θ } θ... θ θ All expecaons on whch he credbly esmaor depends are replaced correspondngly. As a resul

13 A Herarchcal Kalman Fler ( ) Λ Eθ,..., θ Var b θ θ,..., θ θ,..., θ =Λ Φ Φ { } bˆ = Y Φ Y Y Φ x (5.4) Wh hese replacemens, (3.4) and (3.5) become b = Z W b + Z bˆ (5.5) ( ) and Z Y Y =Λ Φ +Λ Y Φ Y (5.6) 6. A herarchcal Kalman fler Secon 5 defned he herarchy of parameers b( ) observaons { :,..., } vecor over me. {,..., :,..., } θ θ = n and x = n. Secon 4 defned he evoluon of a parameer he presen secon combnes he wo conceps no an evoluonary herarchy b : =,..., n; =,,... and observaon vecors of parameer vecors { } { x :,..., n;,,... } = =. he herarchy s subjec o he followng assumpons. (6.) For each =,, and each =,, n, x b,..., b = Y b + ξ for some parameer vecor, deermnsc marx, and cenred random vecor. ξ b (6.) For each =,, and each =,, n, defne c = b W b for some deermnsc marx W. For he case =, defne c = b W β where β = E b and W s some deermnsc marx. he vecors c and β evolve accordng o c = A c +ε, =,..., n; =,, β = A, β +ε =,,... Y for deermnsc marces A + + and cenred random vecors ε. s (6.) he random vecors ξ, ε are muually sochascally ndependen. j

14 A Herarchcal Kalman Fler Remark 6.. he random vecors b are no longer expressed as funcons of absrac laen parameers θ as n Secon 5. Insead, he explc relaons beween hese vecors are spulaed n Assumpon (6.). Remark 6.. By Assumpon (6.), he evoluon of he b s accordng o b = W b + A b A W b +ε, =,..., n; =,,... (6.) For he case =, hs relaon connues o hold f β, β +. b, b + are replaced by he followng wll use a noaon parallel o ha of Secon 5 bu augmened by he superscrp represenng me. For example, x wll denoe he vecor x, defned by (5.4), as occurs n me perod. he followng wll also use he s superscrp noaon from Secon 4 wh he s same meanng as here. For example, b wll denoe he lnearsed esmaor of based on daa x,..., x. s b he objecve s o fnd esmaors b for =,,, n; =, 3,. hs s done by applcaon of he Kalman fler o he vecors x. Begn by nong ha he evoluon of he b may be expressed n he form + = + j j j j j= j= b M b N ε (6.) for marces M + j, N + j defned below, and where b n nerpreed as β. hs relaon may be proved by nducon. I s rue n he case =, by Assumpon (6.). Here M = A, N = (6.3) For he case >, use (6.) o subsue for b + n (6.). hs yelds b W M b N = j j + j j j= j= + A b A W b +ε ε (6.4) Equaly of (6.4) wh (6.) yelds he recurson

15 A Herarchcal Kalman Fler M = W M, j =,,..., j j M + = A M = W M A W N = N j = W N j, j =,,..., (6.5) he relaon (6.) for =,,, n may be summarsed n he form b = M b + N ε (6.6) + ( ) ( ) ( ) ( ) ( ) where he marx M + ( ) akes he block form M M M ( ) M = n + n + n + M M Mn (6.7) and N + s smlarly consruced from he blocks N +. ( ) j Y [] Y Now defne (dsnc from ) as he block dagonal marx wh blocks Y,..., Y. As n Secon 5, where he vecors n x are wren whou as condoned by b, each componen wll be undersood o be condoned by all relevan nformaon, e n x, x wll be condoned by,..., b b ; x + by,..., b b +, ec. Equvalenly, each componen s condoned by b,..., b n. By Assumpon (6.), x b = Y b + (6.8) ξ [] Noe ha (6.6) and (6.8) wh = are of he same form as Assumpons (4.) and (4.) and so he Kalman fler, conssng of (4.), (4.6) (4.9) and (4.), may be appled wh he followng replacemens: b b ε ( ) A M x Y x ξ ξ Y ( ) N [ ] ( ) ( ) ( ) ( ) ε

16 A Herarchcal Kalman Fler 4 { } { } Assumpon (6.) ensures ha he ses N( ) ε ( ) and ξ ( ) are muually sochascally ndependen, as requred by Assumpon (4.). I wll be necessary for hs fler o map an esmaor of o an esmaor of, raher han. hs s because here are no drec observaons Y on b ( ) ˆ and so no esmaor b for use n (6.7). Defne N : ( ) b b b o be he marx obaned by deleng he frs row of M ( : ) M ( ) smlarly. hen (6.6) gves ( ). Defne b = M b + N ε ( ) ( :) ( ) ( :) ( ) (6.9) In correspondence wh (4.) (4.4), defne Φ = Var x b = Var ξ (6.) Ψ = ε (6.) Var Λ = E b b (6.) Noe ha he muual ndependence of he ξ mples ha dagonal srucure. Also defne Φ has block Γ = Var b and noe ha, n correspondence wh (4.5), (6.3) { } Λ = E M b + N ε M b ( ) ( :) ( ) ( :) ( ) ( :) ( ) (6.4) = M Γ M N + Ψ N ( : ) ( ) ( : ) ( : ) ( ) ( : ) (6.5) where, n (6.4), use has been made of (6.9) and (6.6) has been ancpaed. Use has also been made of Assumpon (6.). hen he Kalman fler, wh equaons now wren n operaonal order, s as follows: b M b (6.6) = : ( ) ( )

17 A Herarchcal Kalman Fler 5 ( ) ( ) [] ( []) ( ) ˆ b = Y [] Φ Y Y Φ x (6.7) Λ accordng o (6.5) L [] ( Y ) ( ) [] =Λ Φ Y (6.8) ( ) K L L = + (6.9) ( ) ˆ b = K b + K b (6.) ( K ) Γ = Λ (6.) I remans o calculae b and hence b ( ), and hen Γ ( ). Begn by nong, from (6.), ha ( ) b = W b + A b W b + ε (6.) Now he marx W may no be square. I wll be assumed here o have a leas as many rows as columns and o be of full rank, hs beng he more common case n pracce (eg see example n Secon 7). As shown n he appendx, s always possble o ransform W as follows: R W = U (6.3) where R s defned n he appendx, and U akes he block dagonal form U u = u (6.4) wh each u j a column vecor wh all componens uny. By (6.) and (6.3), ( ) R b = U b + R A b W b + R ε Pre-mulply boh sdes of hs relaon by he block dagonal marx

18 A Herarchcal Kalman Fler 6 Ω = ω ω (6.5) where each ω j s a row vecor and ω u = (6.6) j j hs gves ( ) bo =ΩR b A b W b ε I follows from he unbasedness of ε ha an unbased esmae of b s ( =Ω ) b R b A b W b (6.7) Ω he marx has been consraned by (6.5) and (6.6) bu s oherwse sll free. I may be chosen so ha b s mnmum varance as well as unbased. hs may be done componen of b by componen. By (6.5), each of hese componens akes he form b =ω α (6.8) where he subscrp on ω has been suppressed and α denoes R [ ] n (6.7) and, by (6.6), ω s subjec o he consran ω u = (6.9) By (6.8), Vω (6.3) Var b = ω where V denoes Var [ α ]. I s sraghforward o show ha he choce of ω mnmsng (6.3) subjec o (6.9) s ω= ( uv u) V u (6.3) hs may be evaluaed f V s known. Now V s a submarx of he covarance marx

19 A Herarchcal Kalman Fler 7 ( ) + ( ) ( ) ( ) ( ) R Var b R R AVar b A R + R A W Var b W A R + ( ) ( ) ( ) ( ) ( ) R A Cov b b W A R, RCovb, b A R ( ) ( ) ( ) RCovb, b W A R (6.3) { } where Cov [b,c] denoes E b E[ b] c E[ c] and <.> denoes he symmersaon of he argumen, e, ACov b c B = ACov b, c B + BCov c, b A. [ ] [ ] [ ] Of he sx Var and Cov erms on he rgh, he frs s obanable from Γ n (6.), and he nex hree from Γ ( ). he las wo may be obaned by calculang Cov b, b ( ). By (6.6) and (6.), ( ) ˆ ( ) ( ) ( ) b( ) b = K M b + K : (6.33) ˆ Noe ha depends on daa relang o only me whereas depends on b b ( ) daa relang o only pror perods. he wo are herefore sochascally ndependen, and so (6.33) yelds ( ) ( ) ( ) Cov b, b ( ) = K M Γ (6.34) : Remark 6.3. he ransformaon (6.3) of W no he form (6.4), whle generally applcable, has been seleced here because s convenen for he example n Secon 7. In fac, oher ransformaons could have been chosen, eg (6.3) wh m x n marx U akng he block form U = An alernave choce of U would lead of course o a dfferen R n (6.3) and a dfferen Ω sasfyng (6.6). he remander of he reasonng (6.) (6.34) would be unchanged.

20 A Herarchcal Kalman Fler 8 7. Numercal example 7. Daa he followng example consders workers compensaon clam coss whn he herarchcal occupaonal srucure llusraed n Fgure 7.. Fgure 7. Occupaonal srucure Level oal pool of rsks Level Occupaonal Group A Occupaonal Group B Level Occupaonal sub-group A Occupaonal sub-group A Occupaonal sub-group A3 Occupaonal sub-group B Occupaonal sub-group B Occupaonal sub-group B3 Clams experence has been smulaed over 6 years for each of he occupaonal sub-groups. Each observaon s a realsaon of a gamma, dsrbuon wh parameers as se ou n able 7.. CoV denoes coeffcen of varaon n he able. able 7. Clam cos parameers Segmen Exposure Year Mean CoV Mean CoV Mean CoV Mean CoV Mean CoV Mean CoV $M % % % % % % % % % % % % Pool A B A A A B B Fgure 7. plos he mean clam cos per un exposure for he dfferen occupaonal groups. Noe ha he underlyng clam cos of B undergoes a sep ncrease n Year 4, whle ha of A ramps up over Years 3 and 4.

21 A Herarchcal Kalman Fler 9 Fgure 7. Mean clam cos per un exposure 4.5 Mean clam cos per un exposure (%) Year A A A3 B B he smulaed clams coss are se ou n able 7.. able 7. Clams daa Occupaonal Clam cos per un exposure observed n year group % % % % % % A B A A A B B Sac herarchcal credbly rang For comparson, he occupaonal groups are frs raed by means of Sund s (98) herarchcal credbly, as descrbed n Secon 5. hs s a sac credbly sysem n he sense ha he underlyng clam cos parameers are assumed consan over me. When hey are n fac varable, as n able 7., hs credbly sysem canno be expeced o rack he changes well.

22 A Herarchcal Kalman Fler Snce he sysem s sac, s naural applcaon a he end of years s o he accumulaed clams experence over hose years. Observaons on hs bass are se ou n able 7.3. able 7.3 Accumulang clams coss per un exposure Occupaonal Average clam cos per un exposure observed over years group o o 3 o 4 o 5 o 6 % % % % % % Pool A B A A A B B he paramerc srucure assumed for he occupaonal herarchy s as follows. Le b denoe he vecor of mean clam coss per un exposure a Level of he herarchy. hese vecors wll be of dmensons, and 5 for =, and respecvely. Le x denoe he 5-vecor of observaons a Level. he herarchcal srucure, n he noaon of Secon 5, s as follows: W = W = Y = 5x 5 un marx Y s null as here are no observaons a Level. A leas no drec observaons, hose shown a Levels and n able 7.3 beng merely weghed averages of he observaons a Level. he assumed parameers are as se ou n able 7.4.

23 A Herarchcal Kalman Fler able 7.4 Sac herarchcal credbly parameers Parameer Value β.% (.%) x un marx (x) Λ Λ Φ (.5%) x un marx (5x5) (.5%) x wegh marx (5x5) he wegh marx referred o n able 7.4 s he dagonal marx whose (k,k) elemen s equal o e / ek where ek s he exposure of he k-h occupaonal group a Level. ha s, he varances n Φ are nversely proporonal o he occupaonal group exposures n able 7.. he credbly sysem (5.5) and (5.6) s appled o he observaons, yeldng he resuls se ou n able 7.5. able 7.5 Sac herarchcal credbly rang Occupaonal Sac herarchcal credbly rang group % % % % % % A B A A A B B As seen n Fgure 7., occupaonal groups A and B are subjec o dramac change n rsk over he 6-year perod. Fgure 7.3 plos he credbly esmaes of mean clam cos per un exposure agans he rue parameers. As expeced, he sac naure of he model produces poor rackng of changng rsk parameers.

24 A Herarchcal Kalman Fler Fgure 7.3 Sac herarchcal credbly esmaes rackng parameers Mean clam cos per un exposure (%) Year A parameer A credbly esmae B parameer B credbly esmae 7.3 Dynamc herarchcal credbly rang In hs sub-secon he herarchcal Kalman fler developed n Secon 6 s appled o he observaons se ou n able 7.. he herarchcal srucure s as n Secon 7., wh here akng he same form for each as W n ha earler sub-secon. Each A s a un marx, meanng ha he parameers evolve accordng o ndependen random walks. Furher, Y akes he same form as Y n Secon 7. and R from (6.3) s he 5 5 un marx. W he assumed dsperson parameers are se ou n able 7.6. able 7.6 Dsperson parameers for herarchcal Kalman flers Parameer Value Φ As for Φ n Secon 7. Ψ (.%) Ψ Ψ (.5%) x un marx (x) (.5%) x un marx (5x5) he herarchcal fler (6.6) (6.) s appled o he observaons. I needs o be adaped slghly o he presen herarchy n order o apply o suffx raher han. he resul (6.7) requres smlar adapaon o produce b and b raher han jus he laer.

25 A Herarchcal Kalman Fler 3 he resuls produced by he fler are se ou n able 7.7, whch also reproduces, for comparson, rue parameers, observaons, and sac herarchcal credbly esmaes. able 7.7 Herarchcal Kalman fler credbly rang Occupaonal Observaons and credbly rangs a end of year group % % % % % % Observaons A B A A A B B Sac herarchcal credbly rang A B A A A B B Herarchcal Kalman fler rang A B A A A B B rue parameer A B A A A B B

26 A Herarchcal Kalman Fler 4 Fgure 7.4 provdes a graphcal comparson of rue parameers wh sac and dynamc credbly esmaes for he wo occupaonal groups, A and B, whose parameers have changed over me. he dynamc esmaes (doed lnes) are seen o rack he parameers (sold lnes) far beer han he sac esmaes (dashed lnes). Fgure 7.4 Sac and dynamc herarchcal credbly esmaes rackng parameers Mean clam cos per un exposure (%) Year A parameer A credbly esmae A Kalman esmae B parameer B credbly esmae B Kalman esmae 8. Concluson When he parameers of a herarchy evolve hrough me, he convenonal (sac) herarchcal credbly esmaes of Sund (979, 98) (Secon 5) are lkely o rack poorly. he dynamc equvalen of Sund s model explcly recognses parameer evoluon. hs s done by supermposng a Kalman fler on he herarchcal framework (Secon 6). hs s lkely o lead o much beer rackng of he evolvng parameers by he credbly esmaes. hs s llusraed numercally n Secon 7. Appendx. Proof of (6.3). Suppose W s of dmenson m x n wh m n. Snce W s of full rank, s possble o denfy n lnearly ndependen rows. A row permuaon, represened by m x n marx P, wll poson hese n he frs n rows, e

27 A Herarchcal Kalman Fler 5 PW X = Y where X s he n x n marx conssng of he lnearly ndependen rows, and s herefore non-sngular. hen APW = Y where X A = A sequence of elemenary row operaons, represened by he marx M, may be appled so ha Y s replaced by a marx, each of whose rows consss of a sngle un elemen and he res zero. Fnally, a permuaon of rows, represened by m x m marx Q gves QMAPW of he form (6.4). hen R = QMAP.

28 A Herarchcal Kalman Fler 6 References Hachemeser, C.A. (975). Credbly for regresson models wh applcaon o rend. Appears n Kahn (975, 9-63). Jazwnsk, A.H. (97). Sochasc processes and flerng heory. Academc Press, New York, NY. Kalman, R.E. (96). A new approach o lnear flerng and predcon sysems. Journal of Basc Engneerng, 8, Sund, B. (979). A herarchcal credbly regresson model. Scandnavan Acuaral Journal, 7-4. Sund, B. (98). A mul-level herarchcal credbly regresson model. Scandnavan Acuaral Journal, 5-3. aylor, G.C. (979). Credbly analyss of a general herarchcal model. Scandnavan Acuaral Journal, -.

John Geweke a and Gianni Amisano b a Departments of Economics and Statistics, University of Iowa, USA b European Central Bank, Frankfurt, Germany

John Geweke a and Gianni Amisano b a Departments of Economics and Statistics, University of Iowa, USA b European Central Bank, Frankfurt, Germany Herarchcal Markov Normal Mxure models wh Applcaons o Fnancal Asse Reurns Appendx: Proofs of Theorems and Condonal Poseror Dsrbuons John Geweke a and Gann Amsano b a Deparmens of Economcs and Sascs, Unversy

More information

Panel Data Regression Models

Panel Data Regression Models Panel Daa Regresson Models Wha s Panel Daa? () Mulple dmensoned Dmensons, e.g., cross-secon and me node-o-node (c) Pongsa Pornchawseskul, Faculy of Economcs, Chulalongkorn Unversy (c) Pongsa Pornchawseskul,

More information

Should Exact Index Numbers have Standard Errors? Theory and Application to Asian Growth

Should Exact Index Numbers have Standard Errors? Theory and Application to Asian Growth Should Exac Index umbers have Sandard Errors? Theory and Applcaon o Asan Growh Rober C. Feensra Marshall B. Rensdorf ovember 003 Proof of Proposon APPEDIX () Frs, we wll derve he convenonal Sao-Vara prce

More information

V.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS

V.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS R&RATA # Vol.) 8, March FURTHER AALYSIS OF COFIDECE ITERVALS FOR LARGE CLIET/SERVER COMPUTER ETWORKS Vyacheslav Abramov School of Mahemacal Scences, Monash Unversy, Buldng 8, Level 4, Clayon Campus, Wellngon

More information

( t) Outline of program: BGC1: Survival and event history analysis Oslo, March-May Recapitulation. The additive regression model

( t) Outline of program: BGC1: Survival and event history analysis Oslo, March-May Recapitulation. The additive regression model BGC1: Survval and even hsory analyss Oslo, March-May 212 Monday May 7h and Tuesday May 8h The addve regresson model Ørnulf Borgan Deparmen of Mahemacs Unversy of Oslo Oulne of program: Recapulaon Counng

More information

. The geometric multiplicity is dim[ker( λi. number of linearly independent eigenvectors associated with this eigenvalue.

. The geometric multiplicity is dim[ker( λi. number of linearly independent eigenvectors associated with this eigenvalue. Lnear Algebra Lecure # Noes We connue wh he dscusson of egenvalues, egenvecors, and dagonalzably of marces We wan o know, n parcular wha condons wll assure ha a marx can be dagonalzed and wha he obsrucons

More information

Solution in semi infinite diffusion couples (error function analysis)

Solution in semi infinite diffusion couples (error function analysis) Soluon n sem nfne dffuson couples (error funcon analyss) Le us consder now he sem nfne dffuson couple of wo blocks wh concenraon of and I means ha, n a A- bnary sysem, s bondng beween wo blocks made of

More information

. The geometric multiplicity is dim[ker( λi. A )], i.e. the number of linearly independent eigenvectors associated with this eigenvalue.

. The geometric multiplicity is dim[ker( λi. A )], i.e. the number of linearly independent eigenvectors associated with this eigenvalue. Mah E-b Lecure #0 Noes We connue wh he dscusson of egenvalues, egenvecors, and dagonalzably of marces We wan o know, n parcular wha condons wll assure ha a marx can be dagonalzed and wha he obsrucons are

More information

( ) () we define the interaction representation by the unitary transformation () = ()

( ) () we define the interaction representation by the unitary transformation () = () Hgher Order Perurbaon Theory Mchael Fowler 3/7/6 The neracon Represenaon Recall ha n he frs par of hs course sequence, we dscussed he chrödnger and Hesenberg represenaons of quanum mechancs here n he chrödnger

More information

Department of Economics University of Toronto

Department of Economics University of Toronto Deparmen of Economcs Unversy of Torono ECO408F M.A. Economercs Lecure Noes on Heeroskedascy Heeroskedascy o Ths lecure nvolves lookng a modfcaons we need o make o deal wh he regresson model when some of

More information

F-Tests and Analysis of Variance (ANOVA) in the Simple Linear Regression Model. 1. Introduction

F-Tests and Analysis of Variance (ANOVA) in the Simple Linear Regression Model. 1. Introduction ECOOMICS 35* -- OTE 9 ECO 35* -- OTE 9 F-Tess and Analyss of Varance (AOVA n he Smple Lnear Regresson Model Inroducon The smple lnear regresson model s gven by he followng populaon regresson equaon, or

More information

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4 CS434a/54a: Paern Recognon Prof. Olga Veksler Lecure 4 Oulne Normal Random Varable Properes Dscrmnan funcons Why Normal Random Varables? Analycally racable Works well when observaon comes form a corruped

More information

Dynamic Team Decision Theory. EECS 558 Project Shrutivandana Sharma and David Shuman December 10, 2005

Dynamic Team Decision Theory. EECS 558 Project Shrutivandana Sharma and David Shuman December 10, 2005 Dynamc Team Decson Theory EECS 558 Proec Shruvandana Sharma and Davd Shuman December 0, 005 Oulne Inroducon o Team Decson Theory Decomposon of he Dynamc Team Decson Problem Equvalence of Sac and Dynamc

More information

Robustness Experiments with Two Variance Components

Robustness Experiments with Two Variance Components Naonal Insue of Sandards and Technology (NIST) Informaon Technology Laboraory (ITL) Sascal Engneerng Dvson (SED) Robusness Expermens wh Two Varance Componens by Ana Ivelsse Avlés avles@ns.gov Conference

More information

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Lnear Response Theory: The connecon beween QFT and expermens 3.1. Basc conceps and deas Q: ow do we measure he conducvy of a meal? A: we frs nroduce a weak elecrc feld E, and hen measure

More information

Econ107 Applied Econometrics Topic 5: Specification: Choosing Independent Variables (Studenmund, Chapter 6)

Econ107 Applied Econometrics Topic 5: Specification: Choosing Independent Variables (Studenmund, Chapter 6) Econ7 Appled Economercs Topc 5: Specfcaon: Choosng Independen Varables (Sudenmund, Chaper 6 Specfcaon errors ha we wll deal wh: wrong ndependen varable; wrong funconal form. Ths lecure deals wh wrong ndependen

More information

On One Analytic Method of. Constructing Program Controls

On One Analytic Method of. Constructing Program Controls Appled Mahemacal Scences, Vol. 9, 05, no. 8, 409-407 HIKARI Ld, www.m-hkar.com hp://dx.do.org/0.988/ams.05.54349 On One Analyc Mehod of Consrucng Program Conrols A. N. Kvko, S. V. Chsyakov and Yu. E. Balyna

More information

CS286.2 Lecture 14: Quantum de Finetti Theorems II

CS286.2 Lecture 14: Quantum de Finetti Theorems II CS286.2 Lecure 14: Quanum de Fne Theorems II Scrbe: Mara Okounkova 1 Saemen of he heorem Recall he las saemen of he quanum de Fne heorem from he prevous lecure. Theorem 1 Quanum de Fne). Le ρ Dens C 2

More information

In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!") i+1,q - [(!

In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!) i+1,q - [(! ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL The frs hng o es n wo-way ANOVA: Is here neracon? "No neracon" means: The man effecs model would f. Ths n urn means: In he neracon plo (wh A on he horzonal

More information

Standard Error of Technical Cost Incorporating Parameter Uncertainty

Standard Error of Technical Cost Incorporating Parameter Uncertainty Sandard rror of echncal Cos Incorporang Parameer Uncerany Chrsopher Moron Insurance Ausrala Group Presened o he Acuares Insue General Insurance Semnar 3 ovember 0 Sydney hs paper has been prepared for

More information

ABSTRACT KEYWORDS. Bonus-malus systems, frequency component, severity component. 1. INTRODUCTION

ABSTRACT KEYWORDS. Bonus-malus systems, frequency component, severity component. 1. INTRODUCTION EERAIED BU-MAU YTEM ITH A FREQUECY AD A EVERITY CMET A IDIVIDUA BAI I AUTMBIE IURACE* BY RAHIM MAHMUDVAD AD HEI HAAI ABTRACT Frangos and Vronos (2001) proposed an opmal bonus-malus sysems wh a frequency

More information

HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD

HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD Journal of Appled Mahemacs and Compuaonal Mechancs 3, (), 45-5 HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD Sansław Kukla, Urszula Sedlecka Insue of Mahemacs,

More information

CH.3. COMPATIBILITY EQUATIONS. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

CH.3. COMPATIBILITY EQUATIONS. Continuum Mechanics Course (MMC) - ETSECCPB - UPC CH.3. COMPATIBILITY EQUATIONS Connuum Mechancs Course (MMC) - ETSECCPB - UPC Overvew Compably Condons Compably Equaons of a Poenal Vecor Feld Compably Condons for Infnesmal Srans Inegraon of he Infnesmal

More information

Lecture 6: Learning for Control (Generalised Linear Regression)

Lecture 6: Learning for Control (Generalised Linear Regression) Lecure 6: Learnng for Conrol (Generalsed Lnear Regresson) Conens: Lnear Mehods for Regresson Leas Squares, Gauss Markov heorem Recursve Leas Squares Lecure 6: RLSC - Prof. Sehu Vjayakumar Lnear Regresson

More information

Variants of Pegasos. December 11, 2009

Variants of Pegasos. December 11, 2009 Inroducon Varans of Pegasos SooWoong Ryu bshboy@sanford.edu December, 009 Youngsoo Cho yc344@sanford.edu Developng a new SVM algorhm s ongong research opc. Among many exng SVM algorhms, we wll focus on

More information

Relative controllability of nonlinear systems with delays in control

Relative controllability of nonlinear systems with delays in control Relave conrollably o nonlnear sysems wh delays n conrol Jerzy Klamka Insue o Conrol Engneerng, Slesan Techncal Unversy, 44- Glwce, Poland. phone/ax : 48 32 37227, {jklamka}@a.polsl.glwce.pl Keywor: Conrollably.

More information

Epistemic Game Theory: Online Appendix

Epistemic Game Theory: Online Appendix Epsemc Game Theory: Onlne Appendx Edde Dekel Lucano Pomao Marcano Snscalch July 18, 2014 Prelmnares Fx a fne ype srucure T I, S, T, β I and a probably µ S T. Le T µ I, S, T µ, βµ I be a ype srucure ha

More information

GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS. Youngwoo Ahn and Kitae Kim

GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS. Youngwoo Ahn and Kitae Kim Korean J. Mah. 19 (2011), No. 3, pp. 263 272 GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS Youngwoo Ahn and Kae Km Absrac. In he paper [1], an explc correspondence beween ceran

More information

Advanced Machine Learning & Perception

Advanced Machine Learning & Perception Advanced Machne Learnng & Percepon Insrucor: Tony Jebara SVM Feaure & Kernel Selecon SVM Eensons Feaure Selecon (Flerng and Wrappng) SVM Feaure Selecon SVM Kernel Selecon SVM Eensons Classfcaon Feaure/Kernel

More information

Let s treat the problem of the response of a system to an applied external force. Again,

Let s treat the problem of the response of a system to an applied external force. Again, Page 33 QUANTUM LNEAR RESPONSE FUNCTON Le s rea he problem of he response of a sysem o an appled exernal force. Agan, H() H f () A H + V () Exernal agen acng on nernal varable Hamlonan for equlbrum sysem

More information

Online Appendix for. Strategic safety stocks in supply chains with evolving forecasts

Online Appendix for. Strategic safety stocks in supply chains with evolving forecasts Onlne Appendx for Sraegc safey socs n supply chans wh evolvng forecass Tor Schoenmeyr Sephen C. Graves Opsolar, Inc. 332 Hunwood Avenue Hayward, CA 94544 A. P. Sloan School of Managemen Massachuses Insue

More information

Part II CONTINUOUS TIME STOCHASTIC PROCESSES

Part II CONTINUOUS TIME STOCHASTIC PROCESSES Par II CONTINUOUS TIME STOCHASTIC PROCESSES 4 Chaper 4 For an advanced analyss of he properes of he Wener process, see: Revus D and Yor M: Connuous marngales and Brownan Moon Karazas I and Shreve S E:

More information

New M-Estimator Objective Function. in Simultaneous Equations Model. (A Comparative Study)

New M-Estimator Objective Function. in Simultaneous Equations Model. (A Comparative Study) Inernaonal Mahemacal Forum, Vol. 8, 3, no., 7 - HIKARI Ld, www.m-hkar.com hp://dx.do.org/.988/mf.3.3488 New M-Esmaor Objecve Funcon n Smulaneous Equaons Model (A Comparave Sudy) Ahmed H. Youssef Professor

More information

Lecture 18: The Laplace Transform (See Sections and 14.7 in Boas)

Lecture 18: The Laplace Transform (See Sections and 14.7 in Boas) Lecure 8: The Lalace Transform (See Secons 88- and 47 n Boas) Recall ha our bg-cure goal s he analyss of he dfferenal equaon, ax bx cx F, where we emloy varous exansons for he drvng funcon F deendng on

More information

Time-interval analysis of β decay. V. Horvat and J. C. Hardy

Time-interval analysis of β decay. V. Horvat and J. C. Hardy Tme-nerval analyss of β decay V. Horva and J. C. Hardy Work on he even analyss of β decay [1] connued and resuled n he developmen of a novel mehod of bea-decay me-nerval analyss ha produces hghly accurae

More information

Comparison of Differences between Power Means 1

Comparison of Differences between Power Means 1 In. Journal of Mah. Analyss, Vol. 7, 203, no., 5-55 Comparson of Dfferences beween Power Means Chang-An Tan, Guanghua Sh and Fe Zuo College of Mahemacs and Informaon Scence Henan Normal Unversy, 453007,

More information

Graduate Macroeconomics 2 Problem set 5. - Solutions

Graduate Macroeconomics 2 Problem set 5. - Solutions Graduae Macroeconomcs 2 Problem se. - Soluons Queson 1 To answer hs queson we need he frms frs order condons and he equaon ha deermnes he number of frms n equlbrum. The frms frs order condons are: F K

More information

Lecture VI Regression

Lecture VI Regression Lecure VI Regresson (Lnear Mehods for Regresson) Conens: Lnear Mehods for Regresson Leas Squares, Gauss Markov heorem Recursve Leas Squares Lecure VI: MLSC - Dr. Sehu Vjayakumar Lnear Regresson Model M

More information

Math 128b Project. Jude Yuen

Math 128b Project. Jude Yuen Mah 8b Proec Jude Yuen . Inroducon Le { Z } be a sequence of observed ndependen vecor varables. If he elemens of Z have a on normal dsrbuon hen { Z } has a mean vecor Z and a varancecovarance marx z. Geomercally

More information

Density Matrix Description of NMR BCMB/CHEM 8190

Density Matrix Description of NMR BCMB/CHEM 8190 Densy Marx Descrpon of NMR BCMBCHEM 89 Operaors n Marx Noaon Alernae approach o second order specra: ask abou x magnezaon nsead of energes and ranson probables. If we say wh one bass se, properes vary

More information

Online Supplement for Dynamic Multi-Technology. Production-Inventory Problem with Emissions Trading

Online Supplement for Dynamic Multi-Technology. Production-Inventory Problem with Emissions Trading Onlne Supplemen for Dynamc Mul-Technology Producon-Invenory Problem wh Emssons Tradng by We Zhang Zhongsheng Hua Yu Xa and Baofeng Huo Proof of Lemma For any ( qr ) Θ s easy o verfy ha he lnear programmng

More information

Appendix H: Rarefaction and extrapolation of Hill numbers for incidence data

Appendix H: Rarefaction and extrapolation of Hill numbers for incidence data Anne Chao Ncholas J Goell C seh lzabeh L ander K Ma Rober K Colwell and Aaron M llson 03 Rarefacon and erapolaon wh ll numbers: a framewor for samplng and esmaon n speces dversy sudes cology Monographs

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 5 Lecure 9 Hamlonan Equaons of Moon (Chaper 8) Wha We Dd Las Tme Consruced Hamlonan formalsm H ( q, p, ) = q p L( q, q, ) H p = q H q = p H = L Equvalen o Lagrangan formalsm Smpler, bu

More information

CHAPTER 5: MULTIVARIATE METHODS

CHAPTER 5: MULTIVARIATE METHODS CHAPER 5: MULIVARIAE MEHODS Mulvarae Daa 3 Mulple measuremens (sensors) npus/feaures/arbues: -varae N nsances/observaons/eamples Each row s an eample Each column represens a feaure X a b correspons o he

More information

5th International Conference on Advanced Design and Manufacturing Engineering (ICADME 2015)

5th International Conference on Advanced Design and Manufacturing Engineering (ICADME 2015) 5h Inernaonal onference on Advanced Desgn and Manufacurng Engneerng (IADME 5 The Falure Rae Expermenal Sudy of Specal N Machne Tool hunshan He, a, *, La Pan,b and Bng Hu 3,c,,3 ollege of Mechancal and

More information

Scattering at an Interface: Oblique Incidence

Scattering at an Interface: Oblique Incidence Course Insrucor Dr. Raymond C. Rumpf Offce: A 337 Phone: (915) 747 6958 E Mal: rcrumpf@uep.edu EE 4347 Appled Elecromagnecs Topc 3g Scaerng a an Inerface: Oblque Incdence Scaerng These Oblque noes may

More information

Ordinary Differential Equations in Neuroscience with Matlab examples. Aim 1- Gain understanding of how to set up and solve ODE s

Ordinary Differential Equations in Neuroscience with Matlab examples. Aim 1- Gain understanding of how to set up and solve ODE s Ordnary Dfferenal Equaons n Neuroscence wh Malab eamples. Am - Gan undersandng of how o se up and solve ODE s Am Undersand how o se up an solve a smple eample of he Hebb rule n D Our goal a end of class

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 5 Lecure 9 Hamlonan Equaons of Moon (Chaper 8) Wha We Dd Las Tme Consruced Hamlonan formalsm Hqp (,,) = qp Lqq (,,) H p = q H q = p H L = Equvalen o Lagrangan formalsm Smpler, bu wce as

More information

M. Y. Adamu Mathematical Sciences Programme, AbubakarTafawaBalewa University, Bauchi, Nigeria

M. Y. Adamu Mathematical Sciences Programme, AbubakarTafawaBalewa University, Bauchi, Nigeria IOSR Journal of Mahemacs (IOSR-JM e-issn: 78-578, p-issn: 9-765X. Volume 0, Issue 4 Ver. IV (Jul-Aug. 04, PP 40-44 Mulple SolonSoluons for a (+-dmensonalhroa-sasuma shallow waer wave equaon UsngPanlevé-Bӓclund

More information

Density Matrix Description of NMR BCMB/CHEM 8190

Density Matrix Description of NMR BCMB/CHEM 8190 Densy Marx Descrpon of NMR BCMBCHEM 89 Operaors n Marx Noaon If we say wh one bass se, properes vary only because of changes n he coeffcens weghng each bass se funcon x = h< Ix > - hs s how we calculae

More information

Fall 2010 Graduate Course on Dynamic Learning

Fall 2010 Graduate Course on Dynamic Learning Fall 200 Graduae Course on Dynamc Learnng Chaper 4: Parcle Flers Sepember 27, 200 Byoung-Tak Zhang School of Compuer Scence and Engneerng & Cognve Scence and Bran Scence Programs Seoul aonal Unversy hp://b.snu.ac.kr/~bzhang/

More information

January Examinations 2012

January Examinations 2012 Page of 5 EC79 January Examnaons No. of Pages: 5 No. of Quesons: 8 Subjec ECONOMICS (POSTGRADUATE) Tle of Paper EC79 QUANTITATIVE METHODS FOR BUSINESS AND FINANCE Tme Allowed Two Hours ( hours) Insrucons

More information

Approximate Analytic Solution of (2+1) - Dimensional Zakharov-Kuznetsov(Zk) Equations Using Homotopy

Approximate Analytic Solution of (2+1) - Dimensional Zakharov-Kuznetsov(Zk) Equations Using Homotopy Arcle Inernaonal Journal of Modern Mahemacal Scences, 4, (): - Inernaonal Journal of Modern Mahemacal Scences Journal homepage: www.modernscenfcpress.com/journals/jmms.aspx ISSN: 66-86X Florda, USA Approxmae

More information

NATIONAL UNIVERSITY OF SINGAPORE PC5202 ADVANCED STATISTICAL MECHANICS. (Semester II: AY ) Time Allowed: 2 Hours

NATIONAL UNIVERSITY OF SINGAPORE PC5202 ADVANCED STATISTICAL MECHANICS. (Semester II: AY ) Time Allowed: 2 Hours NATONAL UNVERSTY OF SNGAPORE PC5 ADVANCED STATSTCAL MECHANCS (Semeser : AY 1-13) Tme Allowed: Hours NSTRUCTONS TO CANDDATES 1. Ths examnaon paper conans 5 quesons and comprses 4 prned pages.. Answer all

More information

Volatility Interpolation

Volatility Interpolation Volaly Inerpolaon Prelmnary Verson March 00 Jesper Andreasen and Bran Huge Danse Mares, Copenhagen wan.daddy@danseban.com brno@danseban.com Elecronc copy avalable a: hp://ssrn.com/absrac=69497 Inro Local

More information

2.1 Constitutive Theory

2.1 Constitutive Theory Secon.. Consuve Theory.. Consuve Equaons Governng Equaons The equaons governng he behavour of maerals are (n he spaal form) dρ v & ρ + ρdv v = + ρ = Conservaon of Mass (..a) d x σ j dv dvσ + b = ρ v& +

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 5 Lecure 0 Canoncal Transformaons (Chaper 9) Wha We Dd Las Tme Hamlon s Prncple n he Hamlonan formalsm Dervaon was smple δi δ Addonal end-pon consrans pq H( q, p, ) d 0 δ q ( ) δq ( ) δ

More information

THE PREDICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS

THE PREDICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS THE PREICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS INTROUCTION The wo dmensonal paral dfferenal equaons of second order can be used for he smulaon of compeve envronmen n busness The arcle presens he

More information

TSS = SST + SSE An orthogonal partition of the total SS

TSS = SST + SSE An orthogonal partition of the total SS ANOVA: Topc 4. Orhogonal conrass [ST&D p. 183] H 0 : µ 1 = µ =... = µ H 1 : The mean of a leas one reamen group s dfferen To es hs hypohess, a basc ANOVA allocaes he varaon among reamen means (SST) equally

More information

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Swss Federal Insue of Page 1 The Fne Elemen Mehod for he Analyss of Non-Lnear and Dynamc Sysems Prof. Dr. Mchael Havbro Faber Dr. Nebojsa Mojslovc Swss Federal Insue of ETH Zurch, Swzerland Mehod of Fne

More information

P R = P 0. The system is shown on the next figure:

P R = P 0. The system is shown on the next figure: TPG460 Reservor Smulaon 08 page of INTRODUCTION TO RESERVOIR SIMULATION Analycal and numercal soluons of smple one-dmensonal, one-phase flow equaons As an nroducon o reservor smulaon, we wll revew he smples

More information

Chapter 5. The linear fixed-effects estimators: matrix creation

Chapter 5. The linear fixed-effects estimators: matrix creation haper 5 he lnear fed-effecs esmaors: mar creaon In hs chaper hree basc models and he daa marces needed o creae esmaors for hem are defned. he frs s ermed he cross-secon model: alhough ncorporaes some panel

More information

How about the more general "linear" scalar functions of scalars (i.e., a 1st degree polynomial of the following form with a constant term )?

How about the more general linear scalar functions of scalars (i.e., a 1st degree polynomial of the following form with a constant term )? lmcd Lnear ransformaon of a vecor he deas presened here are que general hey go beyond he radonal mar-vecor ype seen n lnear algebra Furhermore, hey do no deal wh bass and are equally vald for any se of

More information

Robust and Accurate Cancer Classification with Gene Expression Profiling

Robust and Accurate Cancer Classification with Gene Expression Profiling Robus and Accurae Cancer Classfcaon wh Gene Expresson Proflng (Compuaonal ysems Bology, 2005) Auhor: Hafeng L, Keshu Zhang, ao Jang Oulne Background LDA (lnear dscrmnan analyss) and small sample sze problem

More information

Li An-Ping. Beijing , P.R.China

Li An-Ping. Beijing , P.R.China A New Type of Cpher: DICING_csb L An-Png Bejng 100085, P.R.Chna apl0001@sna.com Absrac: In hs paper, we wll propose a new ype of cpher named DICING_csb, whch s derved from our prevous sream cpher DICING.

More information

[Link to MIT-Lab 6P.1 goes here.] After completing the lab, fill in the following blanks: Numerical. Simulation s Calculations

[Link to MIT-Lab 6P.1 goes here.] After completing the lab, fill in the following blanks: Numerical. Simulation s Calculations Chaper 6: Ordnary Leas Squares Esmaon Procedure he Properes Chaper 6 Oulne Cln s Assgnmen: Assess he Effec of Sudyng on Quz Scores Revew o Regresson Model o Ordnary Leas Squares () Esmaon Procedure o he

More information

FI 3103 Quantum Physics

FI 3103 Quantum Physics /9/4 FI 33 Quanum Physcs Aleander A. Iskandar Physcs of Magnesm and Phooncs Research Grou Insu Teknolog Bandung Basc Conces n Quanum Physcs Probably and Eecaon Value Hesenberg Uncerany Prncle Wave Funcon

More information

Existence and Uniqueness Results for Random Impulsive Integro-Differential Equation

Existence and Uniqueness Results for Random Impulsive Integro-Differential Equation Global Journal of Pure and Appled Mahemacs. ISSN 973-768 Volume 4, Number 6 (8), pp. 89-87 Research Inda Publcaons hp://www.rpublcaon.com Exsence and Unqueness Resuls for Random Impulsve Inegro-Dfferenal

More information

Notes on the stability of dynamic systems and the use of Eigen Values.

Notes on the stability of dynamic systems and the use of Eigen Values. Noes on he sabl of dnamc ssems and he use of Egen Values. Source: Macro II course noes, Dr. Davd Bessler s Tme Seres course noes, zarads (999) Ineremporal Macroeconomcs chaper 4 & Techncal ppend, and Hamlon

More information

Advanced Macroeconomics II: Exchange economy

Advanced Macroeconomics II: Exchange economy Advanced Macroeconomcs II: Exchange economy Krzyszof Makarsk 1 Smple deermnsc dynamc model. 1.1 Inroducon Inroducon Smple deermnsc dynamc model. Defnons of equlbrum: Arrow-Debreu Sequenal Recursve Equvalence

More information

Optimal environmental charges under imperfect compliance

Optimal environmental charges under imperfect compliance ISSN 1 746-7233, England, UK World Journal of Modellng and Smulaon Vol. 4 (28) No. 2, pp. 131-139 Opmal envronmenal charges under mperfec complance Dajn Lu 1, Ya Wang 2 Tazhou Insue of Scence and Technology,

More information

NPTEL Project. Econometric Modelling. Module23: Granger Causality Test. Lecture35: Granger Causality Test. Vinod Gupta School of Management

NPTEL Project. Econometric Modelling. Module23: Granger Causality Test. Lecture35: Granger Causality Test. Vinod Gupta School of Management P age NPTEL Proec Economerc Modellng Vnod Gua School of Managemen Module23: Granger Causaly Tes Lecure35: Granger Causaly Tes Rudra P. Pradhan Vnod Gua School of Managemen Indan Insue of Technology Kharagur,

More information

Chapter 6: AC Circuits

Chapter 6: AC Circuits Chaper 6: AC Crcus Chaper 6: Oulne Phasors and he AC Seady Sae AC Crcus A sable, lnear crcu operang n he seady sae wh snusodal excaon (.e., snusodal seady sae. Complee response forced response naural response.

More information

Chapter 8 Dynamic Models

Chapter 8 Dynamic Models Chaper 8 Dnamc odels 8. Inroducon 8. Seral correlaon models 8.3 Cross-seconal correlaons and me-seres crosssecon models 8.4 me-varng coeffcens 8.5 Kalman fler approach 8. Inroducon When s mporan o consder

More information

Bayes rule for a classification problem INF Discriminant functions for the normal density. Euclidean distance. Mahalanobis distance

Bayes rule for a classification problem INF Discriminant functions for the normal density. Euclidean distance. Mahalanobis distance INF 43 3.. Repeon Anne Solberg (anne@f.uo.no Bayes rule for a classfcaon problem Suppose we have J, =,...J classes. s he class label for a pxel, and x s he observed feaure vecor. We can use Bayes rule

More information

Comb Filters. Comb Filters

Comb Filters. Comb Filters The smple flers dscussed so far are characered eher by a sngle passband and/or a sngle sopband There are applcaons where flers wh mulple passbands and sopbands are requred Thecomb fler s an example of

More information

Tight results for Next Fit and Worst Fit with resource augmentation

Tight results for Next Fit and Worst Fit with resource augmentation Tgh resuls for Nex F and Wors F wh resource augmenaon Joan Boyar Leah Epsen Asaf Levn Asrac I s well known ha he wo smple algorhms for he classc n packng prolem, NF and WF oh have an approxmaon rao of

More information

2 Aggregate demand in partial equilibrium static framework

2 Aggregate demand in partial equilibrium static framework Unversy of Mnnesoa 8107 Macroeconomc Theory, Sprng 2009, Mn 1 Fabrzo Perr Lecure 1. Aggregaon 1 Inroducon Probably so far n he macro sequence you have deal drecly wh represenave consumers and represenave

More information

UNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 2017 EXAMINATION

UNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 2017 EXAMINATION INTERNATIONAL TRADE T. J. KEHOE UNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 27 EXAMINATION Please answer wo of he hree quesons. You can consul class noes, workng papers, and arcles whle you are workng on he

More information

CHAPTER 10: LINEAR DISCRIMINATION

CHAPTER 10: LINEAR DISCRIMINATION CHAPER : LINEAR DISCRIMINAION Dscrmnan-based Classfcaon 3 In classfcaon h K classes (C,C,, C k ) We defned dscrmnan funcon g j (), j=,,,k hen gven an es eample, e chose (predced) s class label as C f g

More information

e-journal Reliability: Theory& Applications No 2 (Vol.2) Vyacheslav Abramov

e-journal Reliability: Theory& Applications No 2 (Vol.2) Vyacheslav Abramov June 7 e-ournal Relably: Theory& Applcaons No (Vol. CONFIDENCE INTERVALS ASSOCIATED WITH PERFORMANCE ANALYSIS OF SYMMETRIC LARGE CLOSED CLIENT/SERVER COMPUTER NETWORKS Absrac Vyacheslav Abramov School

More information

Chapter 6 DETECTION AND ESTIMATION: Model of digital communication system. Fundamental issues in digital communications are

Chapter 6 DETECTION AND ESTIMATION: Model of digital communication system. Fundamental issues in digital communications are Chaper 6 DEECIO AD EIMAIO: Fundamenal ssues n dgal communcaons are. Deecon and. Esmaon Deecon heory: I deals wh he desgn and evaluaon of decson makng processor ha observes he receved sgnal and guesses

More information

Technical Appendix to The Equivalence of Wage and Price Staggering in Monetary Business Cycle Models

Technical Appendix to The Equivalence of Wage and Price Staggering in Monetary Business Cycle Models Techncal Appendx o The Equvalence of Wage and Prce Saggerng n Moneary Busness Cycle Models Rochelle M. Edge Dvson of Research and Sascs Federal Reserve Board Sepember 24, 2 Absrac Ths appendx deals he

More information

DEEP UNFOLDING FOR MULTICHANNEL SOURCE SEPARATION SUPPLEMENTARY MATERIAL

DEEP UNFOLDING FOR MULTICHANNEL SOURCE SEPARATION SUPPLEMENTARY MATERIAL DEEP UNFOLDING FOR MULTICHANNEL SOURCE SEPARATION SUPPLEMENTARY MATERIAL Sco Wsdom, John Hershey 2, Jonahan Le Roux 2, and Shnj Waanabe 2 Deparmen o Elecrcal Engneerng, Unversy o Washngon, Seale, WA, USA

More information

CHAPTER FOUR REPEATED MEASURES IN TOXICITY TESTING

CHAPTER FOUR REPEATED MEASURES IN TOXICITY TESTING CHAPTER FOUR REPEATED MEASURES IN TOXICITY TESTING 4. Inroducon The repeaed measures sudy s a very commonly used expermenal desgn n oxcy esng because no only allows one o nvesgae he effecs of he oxcans,

More information

, t 1. Transitions - this one was easy, but in general the hardest part is choosing the which variables are state and control variables

, t 1. Transitions - this one was easy, but in general the hardest part is choosing the which variables are state and control variables Opmal Conrol Why Use I - verss calcls of varaons, opmal conrol More generaly More convenen wh consrans (e.g., can p consrans on he dervaves More nsghs no problem (a leas more apparen han hrogh calcls of

More information

Normal Random Variable and its discriminant functions

Normal Random Variable and its discriminant functions Noral Rando Varable and s dscrnan funcons Oulne Noral Rando Varable Properes Dscrnan funcons Why Noral Rando Varables? Analycally racable Works well when observaon coes for a corruped snle prooype 3 The

More information

Chapter 2 Linear dynamic analysis of a structural system

Chapter 2 Linear dynamic analysis of a structural system Chaper Lnear dynamc analyss of a srucural sysem. Dynamc equlbrum he dynamc equlbrum analyss of a srucure s he mos general case ha can be suded as akes no accoun all he forces acng on. When he exernal loads

More information

Survival Analysis and Reliability. A Note on the Mean Residual Life Function of a Parallel System

Survival Analysis and Reliability. A Note on the Mean Residual Life Function of a Parallel System Communcaons n Sascs Theory and Mehods, 34: 475 484, 2005 Copyrgh Taylor & Francs, Inc. ISSN: 0361-0926 prn/1532-415x onlne DOI: 10.1081/STA-200047430 Survval Analyss and Relably A Noe on he Mean Resdual

More information

RELATIONSHIP BETWEEN VOLATILITY AND TRADING VOLUME: THE CASE OF HSI STOCK RETURNS DATA

RELATIONSHIP BETWEEN VOLATILITY AND TRADING VOLUME: THE CASE OF HSI STOCK RETURNS DATA RELATIONSHIP BETWEEN VOLATILITY AND TRADING VOLUME: THE CASE OF HSI STOCK RETURNS DATA Mchaela Chocholaá Unversy of Economcs Braslava, Slovaka Inroducon (1) one of he characersc feaures of sock reurns

More information

J i-1 i. J i i+1. Numerical integration of the diffusion equation (I) Finite difference method. Spatial Discretization. Internal nodes.

J i-1 i. J i i+1. Numerical integration of the diffusion equation (I) Finite difference method. Spatial Discretization. Internal nodes. umercal negraon of he dffuson equaon (I) Fne dfference mehod. Spaal screaon. Inernal nodes. R L V For hermal conducon le s dscree he spaal doman no small fne spans, =,,: Balance of parcles for an nernal

More information

Machine Learning Linear Regression

Machine Learning Linear Regression Machne Learnng Lnear Regresson Lesson 3 Lnear Regresson Bascs of Regresson Leas Squares esmaon Polynomal Regresson Bass funcons Regresson model Regularzed Regresson Sascal Regresson Mamum Lkelhood (ML)

More information

Appendix to Online Clustering with Experts

Appendix to Online Clustering with Experts A Appendx o Onlne Cluserng wh Expers Furher dscusson of expermens. Here we furher dscuss expermenal resuls repored n he paper. Ineresngly, we observe ha OCE (and n parcular Learn- ) racks he bes exper

More information

arxiv: v1 [cs.sy] 2 Sep 2014

arxiv: v1 [cs.sy] 2 Sep 2014 Noname manuscrp No. wll be nsered by he edor Sgnalng for Decenralzed Roung n a Queueng Nework Y Ouyang Demoshens Tenekezs Receved: dae / Acceped: dae arxv:409.0887v [cs.sy] Sep 04 Absrac A dscree-me decenralzed

More information

Response of MDOF systems

Response of MDOF systems Response of MDOF syses Degree of freedo DOF: he nu nuber of ndependen coordnaes requred o deerne copleely he posons of all pars of a syse a any nsan of e. wo DOF syses hree DOF syses he noral ode analyss

More information

Outline. Probabilistic Model Learning. Probabilistic Model Learning. Probabilistic Model for Time-series Data: Hidden Markov Model

Outline. Probabilistic Model Learning. Probabilistic Model Learning. Probabilistic Model for Time-series Data: Hidden Markov Model Probablsc Model for Tme-seres Daa: Hdden Markov Model Hrosh Mamsuka Bonformacs Cener Kyoo Unversy Oulne Three Problems for probablsc models n machne learnng. Compung lkelhood 2. Learnng 3. Parsng (predcon

More information

Filtrage particulaire et suivi multi-pistes Carine Hue Jean-Pierre Le Cadre and Patrick Pérez

Filtrage particulaire et suivi multi-pistes Carine Hue Jean-Pierre Le Cadre and Patrick Pérez Chaînes de Markov cachées e flrage parculare 2-22 anver 2002 Flrage parculare e suv mul-pses Carne Hue Jean-Perre Le Cadre and Parck Pérez Conex Applcaons: Sgnal processng: arge rackng bearngs-onl rackng

More information

Dual Approximate Dynamic Programming for Large Scale Hydro Valleys

Dual Approximate Dynamic Programming for Large Scale Hydro Valleys Dual Approxmae Dynamc Programmng for Large Scale Hydro Valleys Perre Carpener and Jean-Phlppe Chanceler 1 ENSTA ParsTech and ENPC ParsTech CMM Workshop, January 2016 1 Jon work wh J.-C. Alas, suppored

More information

1 Constant Real Rate C 1

1 Constant Real Rate C 1 Consan Real Rae. Real Rae of Inees Suppose you ae equally happy wh uns of he consumpon good oday o 5 uns of he consumpon good n peod s me. C 5 Tha means you ll be pepaed o gve up uns oday n eun fo 5 uns

More information

ESTIMATIONS OF RESIDUAL LIFETIME OF ALTERNATING PROCESS. COMMON APPROACH TO ESTIMATIONS OF RESIDUAL LIFETIME

ESTIMATIONS OF RESIDUAL LIFETIME OF ALTERNATING PROCESS. COMMON APPROACH TO ESTIMATIONS OF RESIDUAL LIFETIME Srucural relably. The heory and pracce Chumakov I.A., Chepurko V.A., Anonov A.V. ESTIMATIONS OF RESIDUAL LIFETIME OF ALTERNATING PROCESS. COMMON APPROACH TO ESTIMATIONS OF RESIDUAL LIFETIME The paper descrbes

More information