Luhm, Pruscha: Semi-parametric Inference for Regression Models Based on Marked Point Processes

Size: px
Start display at page:

Download "Luhm, Pruscha: Semi-parametric Inference for Regression Models Based on Marked Point Processes"

Transcription

1 Luhm, Pruscha: Sem-paramerc Inference for Regresson Models Based on Marked Pon Processes Sonderforschungsberech 386, Paper 78 (1997) Onlne uner: hp://epub.ub.un-muenchen.de/ Projekparner

2 Sem-paramerc Inference for Regresson Models Based on Marked Pon Processes by A. Luhm and H. Pruscha Unversy of Munch SUMMARY. We sudy marked pon processes (MPP's) wh an arbrary mark space. Frs we develop some sascally relevan opcs n he heory of MPP's admng an nensy kernel (dz), namely marngale resuls, cenral lm heorems for boh he number n of objecs under observaon and he me endng o nny, he decomposon no a local characersc ( (dz)) and a lkelhood approach. Then we presen sem-paramerc sascal nference n a class of Aalen (1975)-ype mulplcave regresson models for MPP's as n!1, usng paral lkelhood mehods. Furhermore, consderng he case!1,we sudy purely paramerc M-esmaors. KYWORDS: Marked pon process, nensy kernel, (locally square negrable) marngale, local characersc, paral lkelhood, M-esmaor. 1 Inroducon and Basc Denons The monography by Andersen e al. (1993) presens a knd of canoncal approach o he sascal analyss of pon process models. I deals wh mulvarae pon processes where each random even carres nformaon on he occurrence me and he ype of even, he laer beng from a ne se of alernaves. The heorecal fundamen o mulvarae pon processes was lad - among ohers - by Jacod (1975), Bremaud (1981) and Dellachere & Meyer (1982). There are applcaons, however, where an uncounable se (e.g., he se of real numbers) of alernaves - now called marks - s more approprae, see Scheke (1994a,b), Murphy (1995) and Pruscha (1997). A mahemacal foundaon of marked pon processes (MPP's) s gven by Las and Brand (1995), bu hs work does no conan all ools necessary for sascal analyss. The rs goal of he presen paper s o ll hs gap. We presen () resuls on MPPnegrals, () lkelhood funcons of an MPP observaon, () cenral lm laws for wo deren suaons denoed as I and II below. These ools are hen used for he asympoc sascal nference n he case of wo deren knds of daa schemes. Scheme I conans n realsaons over a xed me nerval [,T], wh n endng o nny. Here he sem-paramerc analyss of a wde class of Aalen (1975)-ype models can be presened. Scheme II has one sngle realsaon over a longer me nerval [,T], wh T endng o nny. Here we presen a purely paramerc analyss only, he work on he nonparamerc par of he problem s sll n progress. The proofs are skeched only, for complee versons we refer o a forhcomng paper by Luhm. Mahema. Ins. der Ludwg-Maxmlans-Unversa Munchen, Theresensr.39, D-8333 Munchen 1

3 Unless menoned oherwse, we suppose all random elemens and hus all processes o be dened on some probably space ( F ) beng equpped wh a complee, rgh connuous lraon F. Followng Bremaud (1981), we x an arbrary measurable space ( ) and dene an MPP o be a double sequence ( k k ) k2in such ha () ( k ) k2in s a pon process (wh =) () ( k ) k2in s a sequence of random elemens n : The double sequence ( k k ) k2in shall be dened wh s assocaed counng measure N(d dz) whch s dened by N(] ] A) :=N (A) = 1X k=1 1( k )1( k 2 A) A 2: The lraon F N =(F N ) ha consss of he sgma algebras F N := (N s (A) : s A 2) s called nernal hsory of he MPP. Now le P(F ) denoe he F -predcable -algebra on ] 1[ and ~ P(F ):=P(F ) he F -predcable -algebra. ach mappng H : [ 1[! IR such ha () () Hj fg s F B-measurable, Hj ] 1[ s ~P(F ) B-measurable, we call F-predcable -ndexed process or shorly F-P. Fnally, we dene an nensy kernel (dz) o be a ranson measure from ( [ 1[ FB + )no ( ) such ha 8 A 2 he pon process (N (A)) adms he F -predcable nensy ( (A)).We pu 8 A 2 (A) := s (A) ds = s (dz) ds seng especally := () and := (). Furhermore, we suppose ha Denng < 1 ; a:s: 8 : (dz) := (dz) we oban a probably measure on he mark space gven he hsory unl he occurrence me (cf. Jacod (1975), Bremaud (1981) and Las & Brand (1995)). The par ( (dz)) s hen called local characersc. Throughou he paper, we denoe by jaj and jaj he ukldan norm of a vecor a and a marx A, respecvely a 2 s he marx wh ( j)-enry a a j. For a counable se I, we wre jij for he number of elemens of I, whle I d sands for he d-dmensonal un marx, d 2 IN. A 2

4 2 Resuls on Marked Pon Processes 2.1 Marngale Theory As a connuaon of Bremaud (1981), we consder d-dmensonal sochasc processes (d 2 IN) whch are generaed by henegraon of a d-dmensonal F-P H( z) w.r.. he measure M(d dz) :=N(d dz) ; (dz)d.e., we deal wh erms of he form M = Obvously, all hese processes sasfy M =. Theorem 1. The followng mplcaons hold: (a) (1) ) M s a local marngale, (b) (2) ) M s a marngale, H(s z) M(ds dz) : (c) (3) ) M s a unformly negrable marngale, where I I jh (s z)j s (dz) ds < 1 ; a:s: 8 1 d 8 (1) 1 jh (s z)j s (dz) ds jh (s z)j s (dz) ds < d 8 (2) < d: (3) Proof: As o (a),(b), see Bremaud (1981, VIII Corollary 4). Par (c) can be shown by applyng (b), Bremaud (1981, VIII Theorem 3) and Kopp (1984, Theorem 3.3.8). 2 Observe ha all processes we deal wh are of bounded varaon on ne nervals. For square negrable marngales we oban a resul smlar o Theorem 1, connung Bremaud (1972) and Boel e al. (1975). Theorem 2. (a) M s a square negrable marngale () 1 I H 2 (s z) s (dz) ds < d: (4) (b) M s a marngale whch s locally square negrable (= I H 2 (s z) s(dz) ds < d 8 : (5) (c) M s a locally square negrable marngale () H 2 (s z) s(dz) ds < 1 ; a:s: 8 1 d 8 : (6) = R R Proof: Par (a) can be obaned usng he oponal covaraon process [M] H 2 (s z) N(ds dz),, whch s derved accordng o Dellachere & Meyer (1982, VII Theorem 36), whle (c) follows from (a) generalzng he lnes of Lpser & Shryayev 3

5 (1978, Theorem 18.8). Par (b) s a drec consequence of (c) and Theorem 1(b) and was rs formulaed by Scheke (1994a). 2 Noce ha a locally square negrable marngale need no be a marngale as such. To close hs secon, we gve an explc formula for he characersc hm of M. Theorem 3. Le (6) hold. Then he characersc of M s gven by hm = H 2 (s z) s (dz) ds : (7) Proof: Apply Theorem 1 (a) o he oponal varaon process [M] and use he unqueness of he compensaor Cenral Lm Theorems Frs we formulae a Rebolledo (198)-ype cenral lm heorem for he number n of objecs under observaon endng o nny. Theorem 4. For each n 2 IN, le( ) be ameasurable space andn (ddz) be an MPP on ( F F ) wh he nensy kernel (dz). Le furher H be a d-dmensonal F -P for each n 2 IN, fulllng (s z)) 2 s (dz) ds < 1 ; a:s: 8 1 d 8 : (8) (H Fnally, le he followng wo condons hold 8 ">, (n!1) jh (s z)j 2 1(jH (s z)j ") s (dz) ds (H (s z)) 2 s (dz) ds ;! (9) ;! G (1) where G =(g j ) 1 jd s 8 > aposve dene d d-marx, connuous n, and g = 8 1 d. Then we have he followng convergence n dsrbuon n he space D d [ 1[: M D ;! M (n!1) := R R H (s z) M (ds dz),, andm are locally square negrable where M marngales. Furhermore, M s Gaussan wh characersc hm = G. Proof: Apply Meser (1991, Theorem 2.6) n combnaon wh Lndvall (1973, Theorem 3 ), usng Theorems 1,2,3. 2 In he nex heorem, we consder lms of he form!1nsead of n!1. We wll use a famly of non-sngular d d-marces (; ), fulllng (T) () ;! (elemen-wse) as!1 () 9 a famly of d d-marces (C ) beng non-sngular for each > and connuous n such ha for each xed s we have ; ; ;1 s! C s as!1: 4

6 Theorem 5. Le N(d dz) be an MPP wh nensy kernel (dz). Le H be a d-dmensonal F-P sasfyng (6) and (; ) be a famly of non-sngular d d-marces fulllng (T). Le furher 8 "> and for!1he followng wo condons hold j; H(s z)j 2 1(j; H(s z)j ") s (dz) ds (; H(s z)) 2 s (dz) ds where G s a posve dene d d-marx. Then ; M D ;! N d ( G) (!1): ;! (11) ;! G (12) Proof: Generalze he lnes of Pruscha (1984, Theorem 2.4.9), where Rebolledo (198) was used Lkelhood Followng Jacod (1975), Lpser & Shryayev (1978), Bremaud (1981), Pruscha (1984) and Las & Brand (1995), we presen a process L, whch can be nerpreed as he Radon- Nkodym dervave of an MPP w.r.. anoher MPP, especally a marked Posson process, and whch opens he door o he lkelhood approach. Le wo probably measures on ( F) be gven sasfyng. If N(d dz) adms an (F )-nensy kernel (dz), hen here exss a unque (up o - ndsngushably) F-P h = h( z), called Jacod-process, such ha N(d dz) adms he (F )-nensy kernel (dz) =h( z) (dz) z 2 (13) (see Jacod (1975, Theorem 4.1) and Las & Brand (1995, Theorem 1.2.1)). Consderng he decomposons ( (dz)) and ( (dz)) as local characerscs,.e., (dz) = (dz) (dz) = (dz) we ge an F-P g = g( z) sasfyng (dz) =g( z) (dz) by denng g( z) := h( z) z2 : (14) Now le be he probablyonf admng a local characersc (1 (dz)) ndependen of, sasfyng (dz) (dz) 8. Then N(d dz) s a marked sandard Posson process (see Bremaud (1981), VIII xercse 3), and denng ~ := ~ j F for any probably measure ~ and any soppng me, we formulae Theorem 6. Le h denoe he Jacod process nroduced n(13) and L () be he densy of w.r... If s an F -soppng me wh <1 -and-a.s., hen wh d d = L 5

7 where he process (L ) s dened by L := L Y 1kN h( k k ) 1 A exp (1 ; h(s z)) (dz) ds : Proof: The Theorem can be derved from Jacod (1975, Proposon 4.3), Bremaud (1981, VIII Theorem 1) and Las & Brand (1995, Theorems and 1.2.6), followng he lnes of Lpser & Shryayev (1978, Theorem 19.9) and Pruscha (1984, Theorem 2.3.1). 2 Le he nensy kernel depend on some - no necessarly ne dmensonal - unknown parameer of he form ( dz) =h( z ) (dz). Then L L() can be wren as L () = exp log h(s z ) N(ds dz) ; ()+R 2 (15) where R := log L () +, does no depend on 2. 3 Sem-paramerc Mulplcave Models 3.1 A Class of Mulplcave Models As a generalzaon of Andersen e al. (1993, sec.vii.2), we consder a wde class of mulplcave regresson models. Le I be a counable se and N (d dz) for each 2 I be an MPP wh nensy kernel (dz) of he ype ( dz) = ( 1 ) ( 2 dz) ( 1 ) = () r( 1 X ) Y 2 I z 2 where (wh c d 1 d 2 2 IN, d := d 1 + d 2, B = B 1 B 2, B j open IR d j, j =1 2) Y 2 I s an observable, bounded, (ofen f 1g-valued) nonnegave F-P, : IR +! ] 1[ s R an unknown baselne hazard funcon sasfyng (s) ds < 1 8, r : IR d 1 IR c! ] 1[ s a known regresson funcon, =( 1 2 ) T 2 B s an unknown d-dmensonal parameer, X 2 I s an observable, c-dmensonal F -predcable process of covaraes. By N we denoe he superposon P 2I N. Assume ha here exss a probably measure (dz) on such ha ( 2 dz) (dz) 8 2 I : Then (dz) nduces a probably measure on ( F) such ha he MPP N(d dz) on ( F F ) wh he (F )-local characersc (1 (dz)) s a marked sandard Posson process. Consequenly, here exss 8 2 I an F-P g as n (14) sasfyng 6

8 ( 2 dz)=g ( z 2 ) (dz) 2 I z 2 leadng o an nensy kernel of he form ( dz) = () r( 1 X ) Y g ( z 2 ) (dz) 2 I z 2 : Thus, he Jacod process s gven by h ( z )=() r( 1 X ) Y g ( z 2 ) 2 I z 2 : Smlarly o Andersen e al. (1993, p.482), we dene S ( 1 )= X 2I r( 1 X ) Y and by vrue of formula (15), wh =( ), we ge log L ( ) = ; [log (s) + log r( 1 X s ) + log g (s z 2 )] N (ds dz) ; (16) S s ( 1 ) (s) ds + R where R neher depends on nor on. In case of a mulvarae pon process (N ), 2 I (.e. when jj = 1), wh c = d 1 = d and r( X ):=exp( T X ) 2 I we oban he classc Cox' regresson model. 3.2 Paral Log-Lkelhood Guded by he Nelson-Aalen esmaor (cf. Andersen e al. (1993, p.482)), we subsue 1 log () by log and () S d by 1 ( 1 ) S dn ( 1 ) n formula (16), obanng he paral loglkelhood l () = X 2I ; [log r( 1 X s )+logg (s z 2 )] N (ds dz) ; log S s ( 1 ) dn s + R whch only depends on he unknown parameer 2 B. We wll suppose ha all processes are connuously derenable as ofen as needed, and ha he order of summaon, negraon and derenaon may always be changed. For sake of smpler noaon, we dene for 2 I z 2 he followng processes whch weassumeobef-p's. 7

9 r (1) ( 1 X ):= d d 1 r( 1 X ), (d 1 -dmensonal vecor), r (2) ( 1 X ):= d2 d 1 d T 1 r( 1 X ), (d 1 d 1 -marx), g (1) ( z 2 ):= d d 2 g ( z 2 ), g (2) ( z 2 ):= d2 g d 2 d T ( z 2 ), 2 S (1) ( 1 ):= d d 1 S ( 1 ) S (2) ( 1 ):= d2 S d 1 d T ( 1 ) 1 (d 2 -dmensonal vecor), (d 2 d 2 -marx), (d 1 -dmensonal vecor), (d 1 d 1 -marx). Furhermore, S (1) () :=((S (1) ( 1 )) T ::: ) T (d-dmensonal vecor), () := S(2) ( 1 ) S (2)! (d d-marx). Fnally, wenroduce for j =1 2hevecors and marces, respecvely, (j) ( 1 X ):= r(j) ( 1 X ), r( 1 X ) (j) ( z ( z 2 ):= g(j) 2 ) g ( z 2 ). Usng mehods of purely paramerc nference (cf. Pruscha (1984), Andersen e al. (1993, sec. VI.2.2)), we ge he followng hree Lemmas. Lemma 1. Denng U () := d d l (), we oban! U () = X 2I = X 2I (1) ( 1 X s ) (s z 2 ) (1) N(ds dz) ; K (s z ) M ( ds dz) where M ( d dz) =N (d dz) ; ( dz) d and K ( z ):= (1) ( 1 X ) (1) ( z 2 ) R! ; S(1) S (1) () s S s () dn s = () 2 I z 2 : S () R K (s z ) s ( dz) ds = 8. 2 Proof: Observe ha P 2I Lemma 2. Assume ha condon X (A) K 2 j (s z ) s( dz) ds < 1 -a.s. 8 1 j d 2I holds. Then hu( ) = X 2I K 2 (s z ) s ( dz) ds : Proof: Apply Theorems 2 (c) and

10 For he nex Lemma, we pu for 2 I z 2 2 B C z () := D z () := (2) ( 1 X ) ; (1) 2 ( 1 X ) (1) (1) ( 1 X ) 2 (2) ( z 2 ) ; (1) 2 ( z 2 ) ( z 2 )( (1) ( 1 X )) T (1) (1) ( 1 X )( (1) ( z 2 )) T 2 ( z 2 ) Lemma 3. For he d d-marces W () := d d (U ()) T, we ge under (A) X 2! 2 3 W () = C s z () N (ds dz) ; 4 S(2) s () 2I S s () ; S(1) s () 5 dns = S s () where for = w ( ) ; w ( ) w ( ) = X 2I 8 2 < : C s z() ; () S s () ; 4 S(2) s X w ( ) = D s z () s ( dz) ds ; 2I = hu( ) S(1) s () S s ()! 2 39 = 5 1 C A : 1 C A M ( ds dz) (s) (S(1) s ()) 2 S s () ds = Noce ha, n conras o hu( ), neher U() nor W () depend on he funcon. Inroducng he furher abbrevaon R () := X 2I D z () r( 1 X ) Y g ( z 2 ) (dz) we oban w ( ) = hu( ) = ()) 2 R s () ; (S(1) s S s ()! (s) ds : Observe ha n he purely mulvarae case (where jj =1and = 1 ) R () = X 2I (r (1) ( X )) 2 r( X ) Y : In he even more specal Cox' case we have R() S (2) () andc z =. 3.3 Asympocal Inference as n!1 Nowwe subsue he counable se I of he las secons by a sequence (I ) n2in of counable ses, where I mgh label he objecs under observaon (.e., I = f1 ::: ng). Accordngly,we consder sequences ((N (ddz) Y X ), 2 I, z 2, ) n2in 9

11 of he correspondng processes such ha for each 2 I he MPP N (d dz) adms he local characersc ( ( 1 ) ( 2 dz)) wh ( 1 ) = () r( 1 X ) Y ( 2 dz) = g ( z 2 ) (dz) z2 n 2 IN: Whn hs seng we denoe condon (A) by (A n ). Observe ha he funcons r and he unknown parameer do no depend on n 2 IN. For he res of hs secon, we x an arbrary T> o descrbe he end of he observaon nervall [ T]. All lms n hs secon are aken as n ends o nny. Nowwenroduce a sequence of non-sngular d d-normng marces (; n ) n2in sasfyng (N) () ; n! (elemen-wse) () 9 C ; 2 ] 1[ such ha j; n j 2 ji jc ; 8 n 2 IN: Adopng he erms of he prevous secons and equppng hem wh an addonal ndex n 2 IN f necessary, we presen a furher se of condons: There exs mappngs s (1) (2) R (1) : [ T] B! IR IR d IR dd IR dd,respecvely, such ha as n!1 (Bn ) () () () (v) (v) sup 2[ T ] 2B sup 2[ T ] 2B sup 2[ T ] 2B sup 2[ T ] 2B sup 2[ T ] 2B j; n j 2 S S (1 n) () S () ; s () () ; (1) () S (2 n) () ; n S () ; (2) () (S (1 n) ()) 2 S () ; n R (); T n ; R (1) ;! ;! ;! ; T n ; s () ( (1) ()) 2 ;! () ;! (C n ) (D n ) () s (1) (2) and R (1) are bounded n [ T] B () s (1) (2) and R (1) are connuous funcons n 2 B unformly n 2 [ T]: The d d-marx () := s posve dene 8 2 ] T]: R (1) s () ; s s () ( (1) s ()) 2 (s) ds ( n ) There exss a > such ha sup 2I 2[ T ] 1 r( 1 X ; n )g (1) ( 1 X (1 n) ) ( z 2 ) ( z 2 ) > exp ( ;! Y (dz) 1 (1) ( 1 X (1 n) ) ( z 2 ) ;! :! )!

12 In case of ; n = 1 p an I d wh <a n!1, he followng condon s sucen for (B n )()-(v) (B n) () () sup 2[ T ] 2B sup 2[ T ] 2B 1 S (m n) () ; s (m) () a n 1 R () ; R (1) () a n ;! ;! (m = 1 2) where s (m) () s() (m) () m=1 2, demandng addonally ha s () : [ T]! IR s bounded away from. In he purely mulvarae Cox' case, (B n)() conans (B n)() wh R 1 () s (2) (), see Andersen e al. (1993 p.497, condon VII.2.1). A sequence ( ^ ) n2in of d-dmensonal random vecors s called conssen paral lkelhood esmaor or shorly conssen PML, wehaveforn!1 (j ^ ; j < U ( ^ )=)! 1 8 > where U s gven by Lemma1. Proposon 1. Under (A n ),(B n )(),(),(v),(v),(c n ),(D n ),( n ) we have (U n ) ; n U () D ;! N d ( ()) 8 2 [ T]: Proof: One can show ha(b n )(),(), (C n ) and ( n ) mply (9), whle (B n )(v),(v) and (D n ) yeld (1). Condon (A n )allows he applcaon of Theorem 4 whch complees he proof. 2 Now we consder sequences of d-dmensonal random vecors ( n) n2in fulllng (B n) ; ;T n ( n ; ) n 2 IN s -sochascally bounded: Proposon 2. Under (A n ),(B n ),(C n ) we have (W n) ;; n W ; T n ;! () 8 2 [ T] and 8 sequences of d-dmensonal random vecors ( ) n n2in sasfyng (B ): n Proof: Decompose j;; n W ; T n ; ()j smlar o Andersen e al. (1993, p.5-1). 2 Theorem 7. Le (A n ),(B n ),(C n ),(D n ) hold. Then here exss for each 2 [ T]a conssen PML ( ^ ) n2in for fulllng (B ). n Proof: See Pruscha (1996 sec.vi, Saz 1.2). The basc deas are due o Achson & Slvey (1958), Bllngsley (1971) and Fegn (1975). 2 Theorem 8. (A n ),(B n ),(C n ),(D n ) and ( n ) mply (U n), (W n), and under (U n), (W n) we have for any conssen PML ( ^ ) n2in sasfyng (B n) ; ;T n ( ^ ; ) D ;! Nd ( ;1 ()) 8 2 [ T]: Proof: Apply Proposons 1,2 and Pruscha (1996 sec.vi, Saz 1.6). 2 11

13 Now we can esmae as n he classcal pon process heory (see, e.g., Andersen e al. (1993, sec.vii.2)). Frs we approxmae A() := R (s) ds by he Breslow esmaor ^A ( ^ ):= 1(Y S s s > ) ( ^ ) dn s 2 [ T] n2 IN where Y := P 2I Y. Fnally, we oban he kernel esmaor ^ = 1 b N X k=1 ; k b where K s a kernel funcon and b he bandwdh. 1 A 1 S k 2 [ T] n2 IN 4 Paramerc Inference as!1 In hs chaper, we deal wh purely paramerc problems as me ends o nny when here s jus one objec under observaon, develloppng opcs from Pruscha (1984), Hjor (1985) and Andersen e al. (1993). 4.1 M-esmaors We consder an MPP N(d dz) wh nensy kernel ( dz) admng he decomposon ( dz) =h( z ) (dz) z2 where 2 open IR d d2 IN, (dz) s a probably measure on sasfyng ( dz) (dz) 8 2, and h s he accompanyng Jacod process whch we assume o be connuously derenable w.r.. An M-esmaor ^ of 2 s a soluon of he equaons U () =, where U () := K(s z ) M( ds dz) and M( d dz) =N(d dz) ; ( dz) d. For he d-dmensonal process K, we assume jk(s z )j s ( dz) ds < 1 ; a:s: 8 : If K( z ) = d h( z ) d,we refer o he lkelhood case, and^ h( z ) s called maxmum lkelhood esmaor (ML). Now we sae a rs se of basc condons (A ) () K has connuous rs order dervaves w.r.. he processes K K (1) := d d K T and h (1) := d h are F-P's. d () K and h have connuous second order dervaves w.r.. he processes K j K (1)! 1jd 12 and h (2) := d2 h are F-P's. ddt

14 Noe ha he Jacod-process h s by denon an F-P. The followng Lemma can easly be proven usng derenaon rules. Lemma 4. (a) Assumng (A )(), we have for where W () := d d (U ()) T w () := w () := = w () ; w () K (1) (s z ) M( ds dz) K(s z ) (h(1) (s z )) T h(s z ) s ( dz) ds: (b) In he lkelhood case, (A )() mples for wh W () = v () ; v () v () := v () := h (2) (s z ) h(s z ) (h (1) (s z )) 2 h 2 (s z ) M( ds dz) N(ds dz): ven n he lkelhood case, we have w () 6= v () andw() 6= v() however, w() she compensaor of v(). Boh w () andv () have marngale srucure. We furher dene he d d-marx and he d 2 d 2 -marx I () := I () := For hese marces we sae a furher condon: K 2 (s z ) s ( dz) ds (K (1) (s z )) 2 s ( dz) ds : (B ) () Le he dagonal elemens of I () be ne -a.s. 8 : () Le he dagonal elemens of I () be ne -a.s. 8 : Noe ha (B ) mples va Cauchy-Schwarz nequalyhaall he elemens of he marces are ne -a.s. Now Theorems 2(c) and 3 yeld Lemma 5. Le (A )() hold. Then we have (a) (B )() ) U() s a locally square negrable marngale wh characersc hu() = K 2 (s z ) s ( dz) ds : (b) (B )() ) w () s a locally square negrable marngale wh characersc hw () = (K (1) (s z )) 2 s ( dz) ds : 13

15 4.2 Asympocal Inference as!1 As a generalzaon of Pruscha (1984) and Andersen e al. (1993, sec.vi.2), we presen a se of condons under whch he asympoc heory of Pruscha (1996, sec.vi.1) can be appled for M-esmaors. Consderng a famly of non-sngular d d-marces (; ) sasfyng (T), we formulae a furher condon whch has an ergodc [(),()] and a Lndeberg-lke par [()]. All he lms n hs secon are aken as!1. (C ) () () () ; hu() ; T ;! () where () s a posve dene, symmerc d d-marx ; w (); T ;! B() where B() s a posve dene d d-marx L ( ") ;! 8 "> where L ( ") := j; K(s z )j 2 1(j; K(s z )j >") s ( dz) ds: In he lkelhood case, condons (C )() and (C )() are dencal, wh B() =(). Proposon 3. Le condons (A )(), (B )(), (C )(),() be sased. Then we have he followng convergence n dsrbuon (U ) ; U () D ;! N d ( ()) : Proof: (B )(), (C )() and (C )() yeld (6), (11) and (12), respecvely, such ha he applcaon of Theorem 5 complees he proof. 2 To sae asympoc resuls for he process W (), we formulae a furher se of condons (D () kg jh hw hg ) () s -sochascally bounded 8 1 g h j k d where ; =( jk ) 1j kd and w () =(wjk ()) 1j kd : () There exss a neghborhood of and 9 M < 1 such ha lm!1 (j; ~ R ( ~ ); T j <M 8 ~ 2 )=1 where ~ R () =( ~ R jk ()) 1j kd ~ R jk () := dx l=1 () R jkl () W jk l R jkl Proposon 4. (A ), (B )(), (C )() and (D )() mply (W ) ; ; W (); T ;! B() Proof: Apply Lemma 5(b) and use Lenglar's nequaly. 2 Now we consder sequences of d-dmensonal random vecors ( ) fulllng (B ) ; ;T ( ; ) s -sochascally bounded: Proposon 5. Under condons (A ), (B )(), (C )() and (D ) we have as!1 ;; W ( );T (W ) ;! B() 8 sequences of d-dmensonal random vecors ( ) sasfyng (B ): 14

16 Proof: xpand W () n a Taylor seres around he rue parameer and apply Proposon 4. 2 Theorem 9. (A ), (B ), (C ) and (D ) mply (U ), (W), and under (U), (W),he followng holds: There exss a conssen M-esmaor (^ ) for fulllng (B ). If here exss an esmaon funcon l : IR +! IR such ha U() = d l(), hen wh d -probably endng o one, l akes a local maxmum a ^. Furhermore, we have for conssen M-esmaors (^ ) fulllng (B ) ; ;T (^ ; ) D ;! N d ( B ;1 ()()B ;T ()): Proof: Apply Proposons 3,5 and Pruscha (1996 sec.vi, Saz 1.2 and Saz 1.6, usng deas due o Achson & Slvey (1958), Bllngsley (1961) and Fegn (1975)). 2 References [1] Aalen, O.O. (1975): Sascal Inference for a Famly of Counng Processes. PhD hess, Unversy of Calforna, Berkeley. [2] Achson, J. and Slvey, S.D. (1958): Maxmum Lkelhood smaon of Parameers subjec o Resrans. Ann. Mah. Sas. 29, [3] Andersen, P.K. and Borgan,. (1985): Counng Process Models for Lfe Hsory Daa: A Revew. Scand. J. Sascs 12, [4] Andersen, P.K., Borgan,., Gll, R.D. and Kedng, N. (1993): Sascal Models Based on Counng Processes. Sprnger, New York. [5] Bllngsley, P. (1961): Sascal Inference for Markov Processes. The Unversy of Chcago Press, Chcago. [6] Boel, R., Varaya, P. and Wong,. (1975): Marngales on Jump Processes I, II. Sam J. Conrol 13, and [7] Bremaud, P. (1972): A Marngale Approach o Pon Processes. lecronc Res. Lab., Memo no. M-345. Unversy of Calforna, Berkeley. [8] Bremaud, P. (1981): Pon Processes and Queues Marngale Dynamcs. Sprnger, New York. [9] Dellachere, C. and Meyer, P.-A. (1982): Probables and Poenal B. Norh- Holland, Amserdam. 15

17 [1] Fegn, P.D. (1975): Maxmum Lkelhood smaon for Sochasc Processes - a Marngale Approach. Ph. D. Thess, Ausralan Naonal Unversy. [11] Hjor, N. (1985): Dscusson of Andersen & Borgan (1985). Scand. J. Sascs 12, [12] Jacod, J. (1975): Mulvarae Pon Processes: Predcable Projecon, Radon- Nkodym Dervaves, Represenaon of Marngales.. Wahrsch. verw. Gebee 31, [13] Kopp, P.. (1984): Marngales and Sochasc Inegrals. Cambrdge Unversy Press, Cambrdge. [14] Las, G. and Brand, A. (1995): Marked Pon Processes on he Real Lne. Sprnger, New York. [15] Lndvall, T. (1973): Weak Convergence of Probably Measures and Random Funcons n he Funcon Space D[ 1[. Journal of Appled Probably 1, [16] Lpser, R.S. and Shryayev, A.N. (1977): Sascs of Random Processes I General Theory. Sprnger, New York. [17] Lpser, R.S. and Shryayev, A.N. (1978): Sascs of Random Processes II Applcaons. Sprnger, New York. [18] Meser, H. (1991): Funkonale Grenzwersaze fur Marngale. Dplomarbe am mahemaschen Insu der Ludwg-Maxmlans-Unversa n Munchen. [19] Murphy, S.A. (1995): A Cenral Lm Theorem for Local Marngales wh Applcaons o he Analyss of Longudnal Daa. Scand. J. Sascs , Blackwell, Cambrdge (USA). [2] Pruscha, H. (1984): Sascal Inference n Mulvarae Pon Processes wh Paramerc Inensy. Hablaonsschrf, Fachberech Mahemak der Ludwg- Maxmlans-Unversa n Munchen. [21] Pruscha, H. (1996): Angewande Mehoden der Mahemaschen Sask. Teubner, Sugar. [22] Pruscha, H. (1997): Semparamerc smaon n Regresson Models for Pon Processes based on One Realzaon. Proc. 12h In. Workshop Sa. Mod., Schrfenrehe d. Oserr. Sa. Ges. 5, Wen. [23] Rebolledo, R. (198): Cenral Lm Theorems for Local Marngales.. Wahrsch. verw. Gebee 51, [24] Scheke, T.H. (1994a): Paramerc Regresson for Longudnal Daa wh Counng Process Measuremen Tmes. Scand. J. Sascs , Blackwell, Cambrdge (USA). [25] Scheke, T.H. (1994b): Nonparamerc Kernel Regresson when he Regressor follows a Counng Process. Research Repor 94/5, Unv. of Copenhagen. 16

( 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

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

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

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

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

ON THE WEAK LIMITS OF SMOOTH MAPS FOR THE DIRICHLET ENERGY BETWEEN MANIFOLDS

ON THE WEAK LIMITS OF SMOOTH MAPS FOR THE DIRICHLET ENERGY BETWEEN MANIFOLDS ON THE WEA LIMITS OF SMOOTH MAPS FOR THE DIRICHLET ENERGY BETWEEN MANIFOLDS FENGBO HANG Absrac. We denfy all he weak sequenal lms of smooh maps n W (M N). In parcular, hs mples a necessary su cen opologcal

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

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

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

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

SOME NOISELESS CODING THEOREMS OF INACCURACY MEASURE OF ORDER α AND TYPE β

SOME NOISELESS CODING THEOREMS OF INACCURACY MEASURE OF ORDER α AND TYPE β SARAJEVO JOURNAL OF MATHEMATICS Vol.3 (15) (2007), 137 143 SOME NOISELESS CODING THEOREMS OF INACCURACY MEASURE OF ORDER α AND TYPE β M. A. K. BAIG AND RAYEES AHMAD DAR Absrac. In hs paper, we propose

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

Track Properities of Normal Chain

Track Properities of Normal Chain In. J. Conemp. Mah. Scences, Vol. 8, 213, no. 4, 163-171 HIKARI Ld, www.m-har.com rac Propes of Normal Chan L Chen School of Mahemacs and Sascs, Zhengzhou Normal Unversy Zhengzhou Cy, Hennan Provnce, 4544,

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

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

( ) [ ] MAP Decision Rule

( ) [ ] MAP Decision Rule Announcemens Bayes Decson Theory wh Normal Dsrbuons HW0 due oday HW o be assgned soon Proec descrpon posed Bomercs CSE 90 Lecure 4 CSE90, Sprng 04 CSE90, Sprng 04 Key Probables 4 ω class label X feaure

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

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

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

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

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

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

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

( ) () 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

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

. 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

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

. 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

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

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

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

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

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

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

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

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

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

Sampling Procedure of the Sum of two Binary Markov Process Realizations

Sampling Procedure of the Sum of two Binary Markov Process Realizations Samplng Procedure of he Sum of wo Bnary Markov Process Realzaons YURY GORITSKIY Dep. of Mahemacal Modelng of Moscow Power Insue (Techncal Unversy), Moscow, RUSSIA, E-mal: gorsky@yandex.ru VLADIMIR KAZAKOV

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

ESTIMATION METHODS FOR DURATION MODELS. Brian P. McCall. University of Minnesota. and. John J. McCall. University of California, Los Angeles

ESTIMATION METHODS FOR DURATION MODELS. Brian P. McCall. University of Minnesota. and. John J. McCall. University of California, Los Angeles Draf: Commens Welcome ESTIMATION METHODS FOR DURATION MODELS by Bran P. McCall Unversy of Mnnesoa and John J. McCall Unversy of Calforna, Los Angeles and Unversy of Calforna, Sana Barbara JANUARY 25 24

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

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

[ ] 2. [ ]3 + (Δx i + Δx i 1 ) / 2. Δx i-1 Δx i Δx i+1. TPG4160 Reservoir Simulation 2018 Lecture note 3. page 1 of 5

[ ] 2. [ ]3 + (Δx i + Δx i 1 ) / 2. Δx i-1 Δx i Δx i+1. TPG4160 Reservoir Simulation 2018 Lecture note 3. page 1 of 5 TPG460 Reservor Smulaon 08 page of 5 DISCRETIZATIO OF THE FOW EQUATIOS As we already have seen, fne dfference appromaons of he paral dervaves appearng n he flow equaons may be obaned from Taylor seres

More information

Performance Analysis for a Network having Standby Redundant Unit with Waiting in Repair

Performance Analysis for a Network having Standby Redundant Unit with Waiting in Repair TECHNI Inernaonal Journal of Compung Scence Communcaon Technologes VOL.5 NO. July 22 (ISSN 974-3375 erformance nalyss for a Nework havng Sby edundan Un wh ang n epar Jendra Sngh 2 abns orwal 2 Deparmen

More information

A HIERARCHICAL KALMAN FILTER

A HIERARCHICAL KALMAN FILTER 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

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

On elements with index of the form 2 a 3 b in a parametric family of biquadratic elds

On elements with index of the form 2 a 3 b in a parametric family of biquadratic elds On elemens wh ndex of he form a 3 b n a paramerc famly of bquadrac elds Bora JadrevĆ Absrac In hs paper we gve some resuls abou prmve negral elemens p(c p n he famly of bcyclc bquadrac elds L c = Q ) c;

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

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

SELFSIMILAR PROCESSES WITH STATIONARY INCREMENTS IN THE SECOND WIENER CHAOS

SELFSIMILAR PROCESSES WITH STATIONARY INCREMENTS IN THE SECOND WIENER CHAOS POBABILITY AD MATEMATICAL STATISTICS Vol., Fasc., pp. SELFSIMILA POCESSES WIT STATIOAY ICEMETS I TE SECOD WIEE CAOS BY M. M A E J I M A YOKOAMA AD C. A. T U D O LILLE Absrac. We sudy selfsmlar processes

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

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

ANOVA FOR DIFFUSIONS AND ITO PROCESSES. By Per Aslak Mykland and Lan Zhang

ANOVA FOR DIFFUSIONS AND ITO PROCESSES. By Per Aslak Mykland and Lan Zhang Submed o he Annals of Sascs ANOVA FOR DIFFUSIONS AND ITO PROCESSES By Per Aslak Mykland and Lan Zhang The Unversy of Chcago, and Carnege Mellon Unversy and Unversy of Illnos a Chcago Io processes are he

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

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

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

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

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

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

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

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

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

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

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

Cubic Bezier Homotopy Function for Solving Exponential Equations

Cubic Bezier Homotopy Function for Solving Exponential Equations Penerb Journal of Advanced Research n Compung and Applcaons ISSN (onlne: 46-97 Vol. 4, No.. Pages -8, 6 omoopy Funcon for Solvng Eponenal Equaons S. S. Raml *,,. Mohamad Nor,a, N. S. Saharzan,b and M.

More information

Method of upper lower solutions for nonlinear system of fractional differential equations and applications

Method of upper lower solutions for nonlinear system of fractional differential equations and applications Malaya Journal of Maemak, Vol. 6, No. 3, 467-472, 218 hps://do.org/1.26637/mjm63/1 Mehod of upper lower soluons for nonlnear sysem of fraconal dfferenal equaons and applcaons D.B. Dhagude1 *, N.B. Jadhav2

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

Discrete Time Approximation and Monte-Carlo Simulation of Backward Stochastic Differential Equations

Discrete Time Approximation and Monte-Carlo Simulation of Backward Stochastic Differential Equations Dscree Tme Approxmaon and Mone-Carlo Smulaon of Backward Sochasc Dfferenal Equaons Bruno Bouchard Unversé Pars VI, PMA, and CREST Pars, France bouchard@ccrjusseufr Nzar Touz CREST Pars, France ouz@ensaefr

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

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

@FMI c Kyung Moon Sa Co.

@FMI c Kyung Moon Sa Co. Annals of Fuzzy Mahemacs and Informacs Volume 8, No. 2, (Augus 2014), pp. 245 257 ISSN: 2093 9310 (prn verson) ISSN: 2287 6235 (elecronc verson) hp://www.afm.or.kr @FMI c Kyung Moon Sa Co. hp://www.kyungmoon.com

More information

Existence of Time Periodic Solutions for the Ginzburg-Landau Equations. model of superconductivity

Existence of Time Periodic Solutions for the Ginzburg-Landau Equations. model of superconductivity Journal of Mahemacal Analyss and Applcaons 3, 3944 999 Arcle ID jmaa.999.683, avalable onlne a hp:www.dealbrary.com on Exsence of me Perodc Soluons for he Gnzburg-Landau Equaons of Superconducvy Bxang

More information

Analysis And Evaluation of Econometric Time Series Models: Dynamic Transfer Function Approach

Analysis And Evaluation of Econometric Time Series Models: Dynamic Transfer Function Approach 1 Appeared n Proceedng of he 62 h Annual Sesson of he SLAAS (2006) pp 96. Analyss And Evaluaon of Economerc Tme Seres Models: Dynamc Transfer Funcon Approach T.M.J.A.COORAY Deparmen of Mahemacs Unversy

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

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

First-order piecewise-linear dynamic circuits

First-order piecewise-linear dynamic circuits Frs-order pecewse-lnear dynamc crcus. Fndng he soluon We wll sudy rs-order dynamc crcus composed o a nonlnear resse one-por, ermnaed eher by a lnear capacor or a lnear nducor (see Fg.. Nonlnear resse one-por

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

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

TitleA random walk analogue of Levy's th. Studia scientiarum mathematicarum H Citation

TitleA random walk analogue of Levy's th. Studia scientiarum mathematicarum H Citation TleA random walk analogue of Levy's h Auhor(s) Fuja, Takahko Suda scenarum mahemacarum H Caon 3-33 Issue 008-06 Dae Type Journal Arcle Tex Verson auhor URL hp://hdlhandlene/10086/15876 Ths s an auhor's

More information

arxiv: v1 [math.pr] 6 Mar 2019

arxiv: v1 [math.pr] 6 Mar 2019 Local law and Tracy Wdom lm for sparse sochasc block models Jong Yun Hwang J Oon Lee Wooseok Yang arxv:1903.02179v1 mah.pr 6 Mar 2019 March 7, 2019 Absrac We consder he specral properes of sparse sochasc

More information

A NEW TECHNIQUE FOR SOLVING THE 1-D BURGERS EQUATION

A NEW TECHNIQUE FOR SOLVING THE 1-D BURGERS EQUATION S19 A NEW TECHNIQUE FOR SOLVING THE 1-D BURGERS EQUATION by Xaojun YANG a,b, Yugu YANG a*, Carlo CATTANI c, and Mngzheng ZHU b a Sae Key Laboraory for Geomechancs and Deep Underground Engneerng, Chna Unversy

More information

This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and

This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and hs arcle appeared n a journal publshed by Elsever. he aached copy s furnshed o he auhor for nernal non-commercal research and educaon use, ncludng for nsrucon a he auhors nsuon and sharng wh colleagues.

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

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

This document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore.

This document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore. Ths documen s downloaded from DR-NTU, Nanyang Technologcal Unversy Lbrary, Sngapore. Tle A smplfed verb machng algorhm for word paron n vsual speech processng( Acceped verson ) Auhor(s) Foo, Say We; Yong,

More information

On the numerical treatment ofthenonlinear partial differentialequation of fractional order

On the numerical treatment ofthenonlinear partial differentialequation of fractional order IOSR Journal of Mahemacs (IOSR-JM) e-iss: 2278-5728, p-iss: 239-765X. Volume 2, Issue 6 Ver. I (ov. - Dec.26), PP 28-37 www.osrjournals.org On he numercal reamen ofhenonlnear paral dfferenalequaon of fraconal

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

Advanced time-series analysis (University of Lund, Economic History Department)

Advanced time-series analysis (University of Lund, Economic History Department) Advanced me-seres analss (Unvers of Lund, Economc Hsor Dearmen) 3 Jan-3 Februar and 6-3 March Lecure 4 Economerc echnues for saonar seres : Unvarae sochasc models wh Box- Jenns mehodolog, smle forecasng

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

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

Including the ordinary differential of distance with time as velocity makes a system of ordinary differential equations.

Including the ordinary differential of distance with time as velocity makes a system of ordinary differential equations. Soluons o Ordnary Derenal Equaons An ordnary derenal equaon has only one ndependen varable. A sysem o ordnary derenal equaons consss o several derenal equaons each wh he same ndependen varable. An eample

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

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

Testing a new idea to solve the P = NP problem with mathematical induction

Testing a new idea to solve the P = NP problem with mathematical induction Tesng a new dea o solve he P = NP problem wh mahemacal nducon Bacground P and NP are wo classes (ses) of languages n Compuer Scence An open problem s wheher P = NP Ths paper ess a new dea o compare he

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

3. OVERVIEW OF NUMERICAL METHODS

3. OVERVIEW OF NUMERICAL METHODS 3 OVERVIEW OF NUMERICAL METHODS 3 Inroducory remarks Ths chaper summarzes hose numercal echnques whose knowledge s ndspensable for he undersandng of he dfferen dscree elemen mehods: he Newon-Raphson-mehod,

More information

CS 536: Machine Learning. Nonparametric Density Estimation Unsupervised Learning - Clustering

CS 536: Machine Learning. Nonparametric Density Estimation Unsupervised Learning - Clustering CS 536: Machne Learnng Nonparamerc Densy Esmaon Unsupervsed Learnng - Cluserng Fall 2005 Ahmed Elgammal Dep of Compuer Scence Rugers Unversy CS 536 Densy Esmaon - Cluserng - 1 Oulnes Densy esmaon Nonparamerc

More information

The Shapley value for fuzzy games on vague sets

The Shapley value for fuzzy games on vague sets WEA TRANACTIN on INFRMATIN CIENCE APPLICATIN The hapley value or uzzy games on vague ses Fan-Yong Meng* (Correspondng Auhor chool o Managemen Qngdao Technologcal nversy Qngdao 266520 hong Provnce P R Chna

More information

FTCS Solution to the Heat Equation

FTCS Solution to the Heat Equation FTCS Soluon o he Hea Equaon ME 448/548 Noes Gerald Reckenwald Porland Sae Unversy Deparmen of Mechancal Engneerng gerry@pdxedu ME 448/548: FTCS Soluon o he Hea Equaon Overvew Use he forward fne d erence

More information

Childhood Cancer Survivor Study Analysis Concept Proposal

Childhood Cancer Survivor Study Analysis Concept Proposal Chldhood Cancer Survvor Sudy Analyss Concep Proposal 1. Tle: Inverse probably censored weghng (IPCW) o adjus for selecon bas and drop ou n he conex of CCSS analyses 2. Workng group and nvesgaors: Epdemology/Bosascs

More information