Diagnostic tests for frailty
|
|
- Rosaline Cole
- 6 years ago
- Views:
Transcription
1 Diagnosic ess for fraily P. Economou and C. Caroni Deparmen of Mahemaics School of Applied Mahemaics and Physical Sciences Naional Technical Universiy of Ahens 9 Iroon Polyechniou, Zografou Ahens, Greece ( polikon@mah.nua.gr, ccar@mah.nua.gr) Absrac. A common way of allowing heerogeneiy beween individuals in models for lifeime daa is o inroduce an unobservable individual random effec Z. In a proporional hazards framework, he individual s hazard becomes zh b () where h b () is he baseline hazard. The random variable Z is ofen assumed o follow he Gamma or Inverse Gaussian disribuion. We develop here diagnosic ess for hese assumpions. One simple graphical diagnosic is based on he form of he uncondiional survival funcion when h b () is assumed o be Weibull. Anoher plo uses a closure propery of a family of fraily disribuions, which implies ha he fraily among survivors a ime has he same form as he original disribuion of Z, wih he same shape parameer bu differen scale parameer. In his mehod, we esimae he shape parameer a differen imes and examine graphically wheher i is consan. We give simulaion resuls and examples o illusrae hese mehods. Keywords: Lifeime daa, fraily, proporional hazards, Burr disribuion, Generalized Inverse Gaussian disribuion. 1 Inroducion When modelling daa obained from ime-o-even sudies, i is ofen found ha here is heerogeneiy beween individuals, over and above wha can be accouned for by any available covariaes. One common way of allowing for his heerogeneiy is o inroduce an unobservable individual random effec Z, he so-called fraily. This is usually assumed o operae in a proporional hazards framework, so ha i acs muliplicaively on he baseline hazard funcion h b which is common o all individuals. Thus he hazard funcion for an individual wih fraily Z = z is given by h( z) = zh b () If here are also measured covariaes, he model is usually exended o h( z; X) = ze β x() h b () where x(.) is a q-dimensional vecor of possibly ime-dependen covariaes.
2 1202 Economou and Caroni Inroducing heerogeneiy in lifeime daa by means of an unobserved quaniy in his way was iniiaed by [Clayon, 1978], [Vaupel e al., 1979] and [Hougaard, 1984]. The disribuion of he random variable Z is ofen assumed o be Gamma [Vaupel e al., 1979] or Inverse Gaussian [Hougaard, 1984]. As wih any par of he process of saisical modelling, i is desirable o check ha he assumed disribuion is suppored by he daa. The purpose of he presen paper is o develop diagnosic plos for he fraily disribuion, wih he emphasis on hese wo common choices, he Gamma and Inverse Gaussian. We will be assuming ha he baseline hazard funcion h b has been specified correcly and ha he muliplicaive proporional hazards fraily model is he proper one o describe he daa. 2 Diagnosic plos for mixures From he proporional hazards assumpion, i follows ha he condiional survivor funcion for an individual wih fraily z is S( z) = [S b ()] z where S b is he baseline survivor funcion. In paricular, if he baseline model is aken o be Weibull (η, β), hen he survivor funcion condiional on fraily z is S( z) = exp( zs) where s = (/η) β. If Z has disribuion funcion G on (0, ), hen he uncondiional survivor funcion is given by 0 exp( zs)dg(z) If G is aken o be Gamma wih shape and scale parameers boh equal o ν (so ha he mean is one), hen S() = ( 1 + s ν ) ν (1) (he Burr disribuion), while aking G o be Inverse Gaussian wih scale 1 and shape λ yields he survivor funcion ( ( )) 2s S() = exp λ 1 λ + 1 (2) (In boh cases, a consrain has been applied o he parameers of G o make he model idenifiable. Oher choices of consrain are possible bu make no difference in principle.) We now look for appropriae diagnosic plos o enable us o check ha he assumpion of one of hese disribuions is in fac correc. The idea is o
3 Diagnosic ess for fraily 1203 plo some funcion of he non-parameric Kaplan-Meier esimae of survival, Ŝ(), agains some funcion of, o obain he characerisic shape associaed wih he disribuions (1) and (2). Taking logarihms of (1), and supposing ha s/ν is large enough so ha log(1 + s/ν) log(s/ν), we see ha logŝ() agains log should give a sraigh line. [Wolsenholme, 1999] suggess his plo for he Pareo disribuion. This is he special case of he Burr disribuion when he baseline hazard is exponenial, which is a special case of he Weibull (β = 1). However, we observe ha he above approximaion is poor for he early failures. These give he plo a characerisic horizonal secion, whose lengh depends on ν, disappearing as ν becomes small (high degree of heerogeneiy) (Figure 1). When he fraily Ln(S()) v=0.1 Ln(S()) v=0.5 Ln(S()) v= Fig.1. Diagnosic plo for Burr disribuion for various values of he parameer ν simulaions of samples of size 1000, baseline Weibull parameers 1000 (scale) and 2 (shape). Type I censoring a =3000. disribuion is he Inverse Gaussian, aking logarihms of (2) suggess he plo of log( logŝ()) agains log if λ is large. Under hese circumsances, here is only a small degree of heerogeneiy and he disribuion will no be very differen from he baseline Weibull disribuion. This is he sandard diagnosic plo used o check for he Weibull disribuion. I gives a sraigh line wih slope equal o he shape parameer β. On he oher hand, suppose ha λ is large. I is easy o show ha in his case he Weibull(η, β) - Inverse Gaussian(1,λ) mixure ends o he Weibull wih scale parameer η/(2λ) 1/β and shape β/2. Thus he same plo gives a sraigh line wih slope β/2. For inermediae values of λ, he plo should be curved wih slope falling from β o β/2 as ime increases. Examples are shown in Figure 2. 3 Closure propery of he fraily disribuions To develop anoher kind of diagnosic plo, we sar wih a closure propery of Gamma fraily in he proporional hazards model ([Vaupel e al., 1979]).
4 1204 Economou and Caroni Ln( Ln(S())) lamda=0.001 Ln( Ln(S())) lamda=1 Ln( Ln(S())) lamda= Fig. 2. Diagnosic plo for Weibull-Inverse Gaussian mixure for various values of he parameer λ. Deails of simulaions as in Figure 1. Given ha he fraily disribuion among all individuals is Gamma wih scale parameer κ and shape parameer λ, hen he fraily disribuion among he populaion of survivors a ime is again Gamma wih he same shape parameer λ bu differen scale parameer given by κ + H b (), where H b () is he cumulaive baseline hazard funcion. This propery can be generalized o a whole family of disribuions. 3.1 Generalizaion of he Gamma fraily propery The Gamma fraily propery given by [Vaupel e al., 1979] can be generalized firs o he case of available covariaes. More specifically, given ha he fraily disribuion is Gamma(κ, λ), hen he fraily disribuion among survivors a ime, condiional on he value of he covariaes, is Gamma (κ + Hb x (), λ), where Hb x () is given by H x b () = 0 e β x(u) h b (u)du. The closure propery is no a characerizaion of he Gamma disribuion only. I is quie easy o show ha a similar propery holds also for he Inverse Gaussian and he Generalized Inverse Gaussian (GIG). Furhermore, a similar propery holds for a whole class of disribuions ha belong o he exponenial family [Hougaard, 1984]. Le fraily Z be a random variable wih disribuion F(α) on (0, ), where α is he parameer vecor, wih p.d.f. of he form f z (z) = e [z,g(z)][η1(α),η2(α)] ξ(z) Φ(α) which is an exponenial family disribuion wih canonical saisics z and g(z) [Shao, 1998]. For his fraily disribuion he following heorem holds, which exends Hougaard s resul by including covariaes.
5 Diagnosic ess for fraily 1205 Theorem 1 Given he fraily disribuion F(α) wih p.d.f. as above, hen under he proporional hazards fraily model he fraily disribuion among survivors a ime is again F(.). The value of η 1 (α), he elemen of he parameer vecor corresponding o z, changes, bu he componens of η 2 (α) do no. More specifically, he p.d.f. of fraily among survivors a ime is given by f Z T > (z) = e [z,g(z)][η 1 (α),η2(α)] Φ ξ(z) (α) where η 1 (α) = η 1(α) + H x b () and Φ (α) = Φ(α)S T (). The GIG disribuion (and hence he Gamma and Inverse Gaussian disribuions which i includes) belongs o he above class of he exponenial family. Unforunaely, some oher disribuions, like he lognormal, which are also widely used as fraily disribuions, do no belong o his class because hey do no have z as a canonical saisic. This obsacle can be overcome by considering a generalized disribuion, adding one more parameer [Hougaard, 1986] which will be zero iniially. So, he Theorem can be applied o all disribuions F(α) wih p.d.f given by f z (z) = e T(z)[η(α)] Φ(α) ξ(z) where T(z) does no conain z as a componen, since he above disribuion can be seen as a special case of GF(α, β) wih p.d.f. given by f z (z) = e [z,t(z)][β,η(α)] Φ G (α, β) ξ(z) for β = 0. Φ G (α, β) is he inegral over he range of z, R(z), of he numeraor of he previous relaionship. Applying he Theorem o he disribuion GF(α, 0) shows ha he fraily disribuion among he survivors a ime will be again GF bu wih parameer vecor given by (α, H x b ()). 3.2 Plos Given ha he fraily disribuion has been chosen correcly, hen our above Theorem shows ha he vecor η 2 (α) of he iniial parameers does no change when we resric our aenion o he fraily disribuion among hose unis ha have survived unil ime. Le ˆη 2i (α) T> denoe componen i of he maximum likelihood esimae of his vecor among he survivors a ime. If our assumpion of he fraily disribuion is correc, hen ˆη 2i (α) T> for any ime is an asympoically unbiased esimaor of he same quaniy η 2i (α). Therefore, a plo of ˆη 2i (α) T> agains ime should give a sraigh line parallel o he ime axis.
6 1206 Economou and Caroni For he Gamma(κ, λ) and Inverse Gaussian(κ, λ) disribuions of fraily his plo reduces o ˆλ T> versus since η 2 (κ, λ) = λ. For he GIG(λ, δ, γ), he proposed plos are of ˆλ T> and ˆδ T> agains, since η 2 (λ, δ, γ) = (λ 1, δ 2 ) for his disribuion. In all cases, he maximum likelihood esimaes of he model s parameers ha are required for he plos are obained by maximising he logarihm of he usual likelihood funcion for lifeime daa ( i : i = 1, 2...n) L = n { h(i ) δi S( i ) } i=1 where δ i is he censoring indicaor which akes he value 1 if i is an observed lifeime and zero if i represens a righ censored observaion. The expressions for S() are given above in (1) and (2) for he Gamma and Inverse Gaussian fraily disribuions, respecively, and he hazard funcion h() can be obained as minus he derivaive of logs. We firs carry ou his esimaion using all he daa. Then we selec a sequence of convenien ime poins τ j (j = 1, 2...k) and repea he esimaion k imes, using in he ih esimaion only hose daa poins i saisfying i τ j. 3.3 Simulaions To illusrae he mehod, we simulaed a se of 1000 uncensored daa poins from he Burr disribuion (Weibull-Gamma mixure) and produced he above plo based on repeaed esimaes of ν. Then we fied he incorrec Weibull-Inverse Gaussian mixure o he same daa and produced he corresponding plo (Figure 3). Nex we repeaed he exercise wih he roles of he wo disribuions reversed. Thus we generaed a se of daa from he Weibull- Inverse Gaussian mixure and produced he plos for boh he correc model and for he incorrec Burr disribuion (Figure 4). In boh cases, he plos Gamma shape Inv Gaus shape T T Fig. 3. Plo defined in Secion 3.2, fiing (lef) correc Burr disribuion, (righ) incorrec Weibull-Inverse Gaussian mixure o daa generaed from Burr (ν = 1).
7 Diagnosic ess for fraily 1207 Gamma Shape T Inv Gaus Shape T Fig. 4. Plo defined in Secion 3.2, fiing (lef) incorrec Burr disribuion, (righ) correc Weibull-Inverse Gaussian mixure o daa generaed from Weibull-Inverse Gaussian mixure (λ = 0.5). discriminae exremely well beween he wo fraily disribuion; he plo for he correc disribuion is a horizonal sraigh line as prediced, bu he plo for he incorrec disribuion depars clearly from a horizonal line. 4 Example For a real-daa illusraion of our mehods, we used daa on he duraion of a readmill es underaken by 978 successive paiens a a cardiac clinic in Ahens. Figure 5 shows he simple diagnosic plos ha were developed in Secion 2. The plo for he Burr disribuion, on he righ, has he expeced shape of a sraigh line preceded by a horizonal secion. The plo for he Weibull-Inverse Gaussian mixure, on he lef, is curved as expeced, bu he curvaure is greaer han i should be if his is he correc model. Figure 6 shows he diagnosic plos ha were developed in Secion 3. These indicae Ln( Ln(S())) Ln Ln(S()) Ln Fig. 5. Plos defined in Secion 2 for daa on 978 cardiac paiens. Lef: Weibull- Inverse Gaussian mixure; righ: Burr.
8 1208 Economou and Caroni Inv Gaus Shape Gamma Shape Fig. 6. Plos defined in Secion 3 for daa on 978 cardiac paiens. Lef: Weibull- Inverse Gaussian mixure; righ: Burr. clearly ha he assumpion of a Gamma disribuion for fraily is accepable, because he esimaes of is shape parameer a differen imes are scaered abou a horizonal line, bu he Inverse Gaussian assumpion is no. 5 Furher research Alhough graphical diagnosics have proved very useful in saisical modelling, i can also be valuable o have formal saisical ess for he presence of any fraily, hus showing wheher or no i is necessary o use he models considered here. In hese models, one parameer of he disribuion conrols boh he presence and he degree of fraily. For example, when he fraily disribuion is Gamma(ν, ν), he parameer ν conrols he amoun of fraily since V (Z) = 1/ν 0 (ν ). If he basic disribuion is Weibull and he uncondiional disribuion is herefore Burr, he presence of any fraily can be examined by esing he null hypohesis ν = (or 1/ν = 0) by likelihood-based mehods applied o he Burr disribuion. The heory for one of hese mehods, he score es, was given by [Crowder and Kimber, 1997] for mulivariae lifeime daa and he deails for he univariae case which we are ineresed in by [Kimber, 1996]. (In fac, he es also holds for oher Weibull mixures, no jus for Weibull-Gamma = Burr.) Oher likelihood-based ess ha can be applied include a Wald es and a likelihood raio for his parameer. A difficuly ha arises is ha he null value of he parameer being esed falls on he boundary of he parameer space. In such cases, he disribuion of minus wice he log likelihood raio is no given by he usual chi-squared approximaion. Insead, a mixure of chi-squared disribuions usually applies. We inend o complee a sudy of he likelihood raio es for his model and hen carry ou a simulaion sudy o compare he properies of he differen ess in order o recommend he bes one for use in pracice.
9 Diagnosic ess for fraily 1209 References [Clayon, 1978]D.G. Clayon. A model for associaion in bivariae life ables and is applicaion in epidemiological sudies of familial endency in chronic disease incidence. Biomerika, 65: , [Crowder and Kimber, 1997]M. Crowder and A. Kimber. A score es for he mulivariae burr and oher weibull mixure disribuions. Scandinavian Journal of Saisics, 24: , [Hougaard, 1984]P. Hougaard. Life able mehods for heerogeneous populaions: Disribuions describing he heerogeneiy. Biomerika, 71:75 83, [Hougaard, 1986]P. Hougaard. Survival models for heerogeneous populaions derived from sable disribuions. Biomerika, 73: , [Kimber, 1996]A. Kimber. A weibull-based score es for heerogeneiy. Lifeime Daa Analysis, 2:63 71, [Shao, 1998]J. Shao. Mahemaical Saisics. Springer-Verlag, New York, [Vaupel e al., 1979]J.A. Vaupel, K.G. Manon, and E. Sallard. The impac of heerogeneiy in individual fraily on he dynamics of moraliy. Demography, 16: , [Wolsenholme, 1999]L.C. Wolsenholme. Reliabiliy Modelling: A Saisical Approach. Chapman and Hall/CRC, Boca Raon, 1999.
How to Deal with Structural Breaks in Practical Cointegration Analysis
How o Deal wih Srucural Breaks in Pracical Coinegraion Analysis Roselyne Joyeux * School of Economic and Financial Sudies Macquarie Universiy December 00 ABSTRACT In his noe we consider he reamen of srucural
More informationBias in Conditional and Unconditional Fixed Effects Logit Estimation: a Correction * Tom Coupé
Bias in Condiional and Uncondiional Fixed Effecs Logi Esimaion: a Correcion * Tom Coupé Economics Educaion and Research Consorium, Naional Universiy of Kyiv Mohyla Academy Address: Vul Voloska 10, 04070
More informationA new flexible Weibull distribution
Communicaions for Saisical Applicaions and Mehods 2016, Vol. 23, No. 5, 399 409 hp://dx.doi.org/10.5351/csam.2016.23.5.399 Prin ISSN 2287-7843 / Online ISSN 2383-4757 A new flexible Weibull disribuion
More informationVehicle Arrival Models : Headway
Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where
More informationSemi-Competing Risks on A Trivariate Weibull Survival Model
Semi-Compeing Risks on A Trivariae Weibull Survival Model Jenq-Daw Lee Graduae Insiue of Poliical Economy Naional Cheng Kung Universiy Tainan Taiwan 70101 ROC Cheng K. Lee Loss Forecasing Home Loans &
More informationTesting for a Single Factor Model in the Multivariate State Space Framework
esing for a Single Facor Model in he Mulivariae Sae Space Framework Chen C.-Y. M. Chiba and M. Kobayashi Inernaional Graduae School of Social Sciences Yokohama Naional Universiy Japan Faculy of Economics
More informationKriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Kriging Models Predicing Arazine Concenraions in Surface Waer Draining Agriculural Waersheds Paul L. Mosquin, Jeremy Aldworh, Wenlin Chen Supplemenal Maerial Number
More informationR t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t
Exercise 7 C P = α + β R P + u C = αp + βr + v (a) (b) C R = α P R + β + w (c) Assumpions abou he disurbances u, v, w : Classical assumions on he disurbance of one of he equaions, eg. on (b): E(v v s P,
More informationOn Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature
On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check
More informationDEPARTMENT OF STATISTICS
A Tes for Mulivariae ARCH Effecs R. Sco Hacker and Abdulnasser Haemi-J 004: DEPARTMENT OF STATISTICS S-0 07 LUND SWEDEN A Tes for Mulivariae ARCH Effecs R. Sco Hacker Jönköping Inernaional Business School
More information5.2. The Natural Logarithm. Solution
5.2 The Naural Logarihm The number e is an irraional number, similar in naure o π. Is non-erminaing, non-repeaing value is e 2.718 281 828 59. Like π, e also occurs frequenly in naural phenomena. In fac,
More information2. Nonlinear Conservation Law Equations
. Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear
More informationArticle from. Predictive Analytics and Futurism. July 2016 Issue 13
Aricle from Predicive Analyics and Fuurism July 6 Issue An Inroducion o Incremenal Learning By Qiang Wu and Dave Snell Machine learning provides useful ools for predicive analyics The ypical machine learning
More informationEcon107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)
I. Definiions and Problems A. Perfec Mulicollineariy Econ7 Applied Economerics Topic 7: Mulicollineariy (Sudenmund, Chaper 8) Definiion: Perfec mulicollineariy exiss in a following K-variable regression
More informationPENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD
PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD HAN XIAO 1. Penalized Leas Squares Lasso solves he following opimizaion problem, ˆβ lasso = arg max β R p+1 1 N y i β 0 N x ij β j β j (1.1) for some 0.
More information5.1 - Logarithms and Their Properties
Chaper 5 Logarihmic Funcions 5.1 - Logarihms and Their Properies Suppose ha a populaion grows according o he formula P 10, where P is he colony size a ime, in hours. When will he populaion be 2500? We
More informationRobust estimation based on the first- and third-moment restrictions of the power transformation model
h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,
More informationSTATE-SPACE MODELLING. A mass balance across the tank gives:
B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing
More informationSolutions to Odd Number Exercises in Chapter 6
1 Soluions o Odd Number Exercises in 6.1 R y eˆ 1.7151 y 6.3 From eˆ ( T K) ˆ R 1 1 SST SST SST (1 R ) 55.36(1.7911) we have, ˆ 6.414 T K ( ) 6.5 y ye ye y e 1 1 Consider he erms e and xe b b x e y e b
More information20. Applications of the Genetic-Drift Model
0. Applicaions of he Geneic-Drif Model 1) Deermining he probabiliy of forming any paricular combinaion of genoypes in he nex generaion: Example: If he parenal allele frequencies are p 0 = 0.35 and q 0
More informationPart III: Chap. 2.5,2.6 & 12
Survival Analysis Mah 434 Fall 2011 Par III: Chap. 2.5,2.6 & 12 Jimin Ding Mah Dep. www.mah.wusl.edu/ jmding/mah434/index.hml Jimin Ding, Ocober 4, 2011 Survival Analysis, Fall 2011 - p. 1/14 Jimin Ding,
More informationIntroduction to Probability and Statistics Slides 4 Chapter 4
Inroducion o Probabiliy and Saisics Slides 4 Chaper 4 Ammar M. Sarhan, asarhan@mahsa.dal.ca Deparmen of Mahemaics and Saisics, Dalhousie Universiy Fall Semeser 8 Dr. Ammar Sarhan Chaper 4 Coninuous Random
More informationResearch Article Interval Estimation for Extreme Value Parameter with Censored Data
Inernaional Scholarly Research Nework ISRN Applied Mahemaics Volume 2, Aricle ID 687343, 2 pages doi:.542/2/687343 Research Aricle Inerval Esimaion for Exreme Value Parameer wih Censored Daa Eun-Joo Lee,
More information0.1 MAXIMUM LIKELIHOOD ESTIMATION EXPLAINED
0.1 MAXIMUM LIKELIHOOD ESTIMATIO EXPLAIED Maximum likelihood esimaion is a bes-fi saisical mehod for he esimaion of he values of he parameers of a sysem, based on a se of observaions of a random variable
More informationACE 562 Fall Lecture 4: Simple Linear Regression Model: Specification and Estimation. by Professor Scott H. Irwin
ACE 56 Fall 005 Lecure 4: Simple Linear Regression Model: Specificaion and Esimaion by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Simple Regression: Economic and Saisical Model
More informationInventory Control of Perishable Items in a Two-Echelon Supply Chain
Journal of Indusrial Engineering, Universiy of ehran, Special Issue,, PP. 69-77 69 Invenory Conrol of Perishable Iems in a wo-echelon Supply Chain Fariborz Jolai *, Elmira Gheisariha and Farnaz Nojavan
More informationACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H.
ACE 56 Fall 005 Lecure 5: he Simple Linear Regression Model: Sampling Properies of he Leas Squares Esimaors by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Inference in he Simple
More informationExponential Weighted Moving Average (EWMA) Chart Under The Assumption of Moderateness And Its 3 Control Limits
DOI: 0.545/mjis.07.5009 Exponenial Weighed Moving Average (EWMA) Char Under The Assumpion of Moderaeness And Is 3 Conrol Limis KALPESH S TAILOR Assisan Professor, Deparmen of Saisics, M. K. Bhavnagar Universiy,
More information10. State Space Methods
. Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he
More informationComparing Means: t-tests for One Sample & Two Related Samples
Comparing Means: -Tess for One Sample & Two Relaed Samples Using he z-tes: Assumpions -Tess for One Sample & Two Relaed Samples The z-es (of a sample mean agains a populaion mean) is based on he assumpion
More informationChapter 2. Models, Censoring, and Likelihood for Failure-Time Data
Chaper 2 Models, Censoring, and Likelihood for Failure-Time Daa William Q. Meeker and Luis A. Escobar Iowa Sae Universiy and Louisiana Sae Universiy Copyrigh 1998-2008 W. Q. Meeker and L. A. Escobar. Based
More informationGMM - Generalized Method of Moments
GMM - Generalized Mehod of Momens Conens GMM esimaion, shor inroducion 2 GMM inuiion: Maching momens 2 3 General overview of GMM esimaion. 3 3. Weighing marix...........................................
More informationSTRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN
Inernaional Journal of Applied Economerics and Quaniaive Sudies. Vol.1-3(004) STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN 001-004 OBARA, Takashi * Absrac The
More informationDiebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles
Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance
More informationNotes on Kalman Filtering
Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren
More informationState-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter
Sae-Space Models Iniializaion, Esimaion and Smoohing of he Kalman Filer Iniializaion of he Kalman Filer The Kalman filer shows how o updae pas predicors and he corresponding predicion error variances when
More informationChapter 2. First Order Scalar Equations
Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.
More informationMatrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality
Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]
More informationOn Multicomponent System Reliability with Microshocks - Microdamages Type of Components Interaction
On Mulicomponen Sysem Reliabiliy wih Microshocks - Microdamages Type of Componens Ineracion Jerzy K. Filus, and Lidia Z. Filus Absrac Consider a wo componen parallel sysem. The defined new sochasic dependences
More informationA Group Acceptance Sampling Plans Based on Truncated Life Tests for Type-II Generalized Log-Logistic Distribution
ProbSa Forum, Volume 09, July 2016, Pages 88 94 ISSN 0974-3235 ProbSa Forum is an e-journal. For deails please visi www.probsa.org.in A Group Accepance Sampling Plans Based on Truncaed Life Tess for Type-II
More informationThe Gradient Statistic George R. Terrell Statistics Department, Virginia Polytechnic Institute and State University Blacksburg, Virginia 24061
1 he Gradien Saisic George R. errell Saisics Deparmen, Virginia Polyechnic Insiue and Sae Universiy Blacksburg, Virginia 24061 Absrac: We propose an alernaive o he usual large-sample es saisics, he gradien
More informationClass Meeting # 10: Introduction to the Wave Equation
MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion
More informationSome Basic Information about M-S-D Systems
Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,
More informationMatlab and Python programming: how to get started
Malab and Pyhon programming: how o ge sared Equipping readers he skills o wrie programs o explore complex sysems and discover ineresing paerns from big daa is one of he main goals of his book. In his chaper,
More informationThe Simple Linear Regression Model: Reporting the Results and Choosing the Functional Form
Chaper 6 The Simple Linear Regression Model: Reporing he Resuls and Choosing he Funcional Form To complee he analysis of he simple linear regression model, in his chaper we will consider how o measure
More informationUnit Root Time Series. Univariate random walk
Uni Roo ime Series Univariae random walk Consider he regression y y where ~ iid N 0, he leas squares esimae of is: ˆ yy y y yy Now wha if = If y y hen le y 0 =0 so ha y j j If ~ iid N 0, hen y ~ N 0, he
More information4.1 - Logarithms and Their Properties
Chaper 4 Logarihmic Funcions 4.1 - Logarihms and Their Properies Wha is a Logarihm? We define he common logarihm funcion, simply he log funcion, wrien log 10 x log x, as follows: If x is a posiive number,
More informationMathematical Theory and Modeling ISSN (Paper) ISSN (Online) Vol 3, No.3, 2013
Mahemaical Theory and Modeling ISSN -580 (Paper) ISSN 5-05 (Online) Vol, No., 0 www.iise.org The ffec of Inverse Transformaion on he Uni Mean and Consan Variance Assumpions of a Muliplicaive rror Model
More information3.1 More on model selection
3. More on Model selecion 3. Comparing models AIC, BIC, Adjused R squared. 3. Over Fiing problem. 3.3 Sample spliing. 3. More on model selecion crieria Ofen afer model fiing you are lef wih a handful of
More informationVectorautoregressive Model and Cointegration Analysis. Time Series Analysis Dr. Sevtap Kestel 1
Vecorauoregressive Model and Coinegraion Analysis Par V Time Series Analysis Dr. Sevap Kesel 1 Vecorauoregression Vecor auoregression (VAR) is an economeric model used o capure he evoluion and he inerdependencies
More informationRC, RL and RLC circuits
Name Dae Time o Complee h m Parner Course/ Secion / Grade RC, RL and RLC circuis Inroducion In his experimen we will invesigae he behavior of circuis conaining combinaions of resisors, capaciors, and inducors.
More informationSimulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010
Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid
More information4.1 Other Interpretations of Ridge Regression
CHAPTER 4 FURTHER RIDGE THEORY 4. Oher Inerpreaions of Ridge Regression In his secion we will presen hree inerpreaions for he use of ridge regression. The firs one is analogous o Hoerl and Kennard reasoning
More informationFITTING OF A PARTIALLY REPARAMETERIZED GOMPERTZ MODEL TO BROILER DATA
FITTING OF A PARTIALLY REPARAMETERIZED GOMPERTZ MODEL TO BROILER DATA N. Okendro Singh Associae Professor (Ag. Sa.), College of Agriculure, Cenral Agriculural Universiy, Iroisemba 795 004, Imphal, Manipur
More informationModal identification of structures from roving input data by means of maximum likelihood estimation of the state space model
Modal idenificaion of srucures from roving inpu daa by means of maximum likelihood esimaion of he sae space model J. Cara, J. Juan, E. Alarcón Absrac The usual way o perform a forced vibraion es is o fix
More informationFinal Spring 2007
.615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o
More informationRight tail. Survival function
Densiy fi (con.) Lecure 4 The aim of his lecure is o improve our abiliy of densiy fi and knowledge of relaed opics. Main issues relaed o his lecure are: logarihmic plos, survival funcion, HS-fi mixures,
More informationAir Traffic Forecast Empirical Research Based on the MCMC Method
Compuer and Informaion Science; Vol. 5, No. 5; 0 ISSN 93-8989 E-ISSN 93-8997 Published by Canadian Cener of Science and Educaion Air Traffic Forecas Empirical Research Based on he MCMC Mehod Jian-bo Wang,
More informationThe General Linear Test in the Ridge Regression
ommunicaions for Saisical Applicaions Mehods 2014, Vol. 21, No. 4, 297 307 DOI: hp://dx.doi.org/10.5351/sam.2014.21.4.297 Prin ISSN 2287-7843 / Online ISSN 2383-4757 The General Linear Tes in he Ridge
More informationReliability of Technical Systems
eliabiliy of Technical Sysems Main Topics Inroducion, Key erms, framing he problem eliabiliy parameers: Failure ae, Failure Probabiliy, Availabiliy, ec. Some imporan reliabiliy disribuions Componen reliabiliy
More informationCHERNOFF DISTANCE AND AFFINITY FOR TRUNCATED DISTRIBUTIONS *
haper 5 HERNOFF DISTANE AND AFFINITY FOR TRUNATED DISTRIBUTIONS * 5. Inroducion In he case of disribuions ha saisfy he regulariy condiions, he ramer- Rao inequaliy holds and he maximum likelihood esimaor
More informationOutline. lse-logo. Outline. Outline. 1 Wald Test. 2 The Likelihood Ratio Test. 3 Lagrange Multiplier Tests
Ouline Ouline Hypohesis Tes wihin he Maximum Likelihood Framework There are hree main frequenis approaches o inference wihin he Maximum Likelihood framework: he Wald es, he Likelihood Raio es and he Lagrange
More informationACE 564 Spring Lecture 7. Extensions of The Multiple Regression Model: Dummy Independent Variables. by Professor Scott H.
ACE 564 Spring 2006 Lecure 7 Exensions of The Muliple Regression Model: Dumm Independen Variables b Professor Sco H. Irwin Readings: Griffihs, Hill and Judge. "Dumm Variables and Varing Coefficien Models
More informationChapter 4. Location-Scale-Based Parametric Distributions. William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University
Chaper 4 Locaion-Scale-Based Parameric Disribuions William Q. Meeker and Luis A. Escobar Iowa Sae Universiy and Louisiana Sae Universiy Copyrigh 1998-2008 W. Q. Meeker and L. A. Escobar. Based on he auhors
More informationPhysics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution
Physics 7b: Saisical Mechanics Fokker-Planck Equaion The Langevin equaion approach o he evoluion of he velociy disribuion for he Brownian paricle migh leave you uncomforable. A more formal reamen of his
More informationPhysics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle
Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,
More information) were both constant and we brought them from under the integral.
YIELD-PER-RECRUIT (coninued The yield-per-recrui model applies o a cohor, bu we saw in he Age Disribuions lecure ha he properies of a cohor do no apply in general o a collecion of cohors, which is wha
More informationDEPARTMENT OF ECONOMICS
ISSN 0819-6 ISBN 0 730 609 9 THE UNIVERSITY OF MELBOURNE DEPARTMENT OF ECONOMICS RESEARCH PAPER NUMBER 95 NOVEMBER 005 INTERACTIONS IN REGRESSIONS by Joe Hirschberg & Jenny Lye Deparmen of Economics The
More informationA Specification Test for Linear Dynamic Stochastic General Equilibrium Models
Journal of Saisical and Economeric Mehods, vol.1, no.2, 2012, 65-70 ISSN: 2241-0384 (prin), 2241-0376 (online) Scienpress Ld, 2012 A Specificaion Tes for Linear Dynamic Sochasic General Equilibrium Models
More informationKINEMATICS IN ONE DIMENSION
KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec
More informationThe expectation value of the field operator.
The expecaion value of he field operaor. Dan Solomon Universiy of Illinois Chicago, IL dsolom@uic.edu June, 04 Absrac. Much of he mahemaical developmen of quanum field heory has been in suppor of deermining
More informationMath 10B: Mock Mid II. April 13, 2016
Name: Soluions Mah 10B: Mock Mid II April 13, 016 1. ( poins) Sae, wih jusificaion, wheher he following saemens are rue or false. (a) If a 3 3 marix A saisfies A 3 A = 0, hen i canno be inverible. True.
More informationNew Challenges for Longitudinal Data Analysis Joint modelling of Longitudinal and Competing risks data
New Challenges for Longiudinal Daa Analysis Join modelling of Longiudinal and Compeing risks daa Ruwanhi Kolamunnage-Dona Universiy of Liverpool Acknowledgmen Paula Williamson, Pee Philipson, Tony Marson,
More informationEXERCISES FOR SECTION 1.5
1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler
More informationInternational Journal of Scientific & Engineering Research, Volume 4, Issue 10, October ISSN
Inernaional Journal of Scienific & Engineering Research, Volume 4, Issue 10, Ocober-2013 900 FUZZY MEAN RESIDUAL LIFE ORDERING OF FUZZY RANDOM VARIABLES J. EARNEST LAZARUS PIRIYAKUMAR 1, A. YAMUNA 2 1.
More informationMANAGEMENT SCIENCE doi /mnsc ec pp. ec1 ec20
MANAGEMENT SCIENCE doi.287/mnsc.7.82ec pp. ec ec2 e-companion ONLY AVAILABLE IN ELECTRONIC FORM informs 28 INFORMS Elecronic Companion Saffing of Time-Varying Queues o Achieve Time-Sable Performance by
More informationBayesian inference for spatial parametric proportional hazards model using Spatsurv R
Bayesian inference for spaial parameric proporional hazards model using Spasurv R Sri Asui Thamrin, Amran, Andi Kresna Jaya, Sulvirah Rahmi, and Ansariadi Ciaion: AIP Conference Proceedings 1827, 020015
More informationAn introduction to the theory of SDDP algorithm
An inroducion o he heory of SDDP algorihm V. Leclère (ENPC) Augus 1, 2014 V. Leclère Inroducion o SDDP Augus 1, 2014 1 / 21 Inroducion Large scale sochasic problem are hard o solve. Two ways of aacking
More informationTwo Coupled Oscillators / Normal Modes
Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own
More information23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes
Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals
More informationT L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB
Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal
More informationReferences are appeared in the last slide. Last update: (1393/08/19)
SYSEM IDEIFICAIO Ali Karimpour Associae Professor Ferdowsi Universi of Mashhad References are appeared in he las slide. Las updae: 0..204 393/08/9 Lecure 5 lecure 5 Parameer Esimaion Mehods opics o be
More informationLab 10: RC, RL, and RLC Circuits
Lab 10: RC, RL, and RLC Circuis In his experimen, we will invesigae he behavior of circuis conaining combinaions of resisors, capaciors, and inducors. We will sudy he way volages and currens change in
More informationIntegration Over Manifolds with Variable Coordinate Density
Inegraion Over Manifolds wih Variable Coordinae Densiy Absrac Chrisopher A. Lafore clafore@gmail.com In his paper, he inegraion of a funcion over a curved manifold is examined in he case where he curvaure
More informationU( θ, θ), U(θ 1/2, θ + 1/2) and Cauchy (θ) are not exponential families. (The proofs are not easy and require measure theory. See the references.
Lecure 5 Exponenial Families Exponenial families, also called Koopman-Darmois families, include a quie number of well known disribuions. Many nice properies enjoyed by exponenial families allow us o provide
More informationPredator - Prey Model Trajectories and the nonlinear conservation law
Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories
More informationRecursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems
8 Froniers in Signal Processing, Vol. 1, No. 1, July 217 hps://dx.doi.org/1.2266/fsp.217.112 Recursive Leas-Squares Fixed-Inerval Smooher Using Covariance Informaion based on Innovaion Approach in Linear
More informationLicenciatura de ADE y Licenciatura conjunta Derecho y ADE. Hoja de ejercicios 2 PARTE A
Licenciaura de ADE y Licenciaura conjuna Derecho y ADE Hoja de ejercicios PARTE A 1. Consider he following models Δy = 0.8 + ε (1 + 0.8L) Δ 1 y = ε where ε and ε are independen whie noise processes. In
More informationEffects of Coordinate Curvature on Integration
Effecs of Coordinae Curvaure on Inegraion Chrisopher A. Lafore clafore@gmail.com Absrac In his paper, he inegraion of a funcion over a curved manifold is examined in he case where he curvaure of he manifold
More informationChapter 8 The Complete Response of RL and RC Circuits
Chaper 8 The Complee Response of RL and RC Circuis Seoul Naional Universiy Deparmen of Elecrical and Compuer Engineering Wha is Firs Order Circuis? Circuis ha conain only one inducor or only one capacior
More informationResearch Report Statistical Research Unit Department of Economics University of Gothenburg
Research Repor Saisical Research Uni Deparmen of Economics Universiy of Gohenburg Sweden Hoelling s T Mehod in Mulivariae On-Line Surveillance. On he Delay of an Alarm E. Andersson Research Repor 008:3
More informationDuration Models. Class Notes Manuel Arellano November 1989 Revised: December 3, 2008
Duraion Models Class Noes Manuel Arellano November 1989 Revised: December 3, 28 1 Duraion daa Duraion daa appear in a diversiy of siuaions in economics, bu I will ofen refer o unemploymen duraion o provide
More informationBiol. 356 Lab 8. Mortality, Recruitment, and Migration Rates
Biol. 356 Lab 8. Moraliy, Recruimen, and Migraion Raes (modified from Cox, 00, General Ecology Lab Manual, McGraw Hill) Las week we esimaed populaion size hrough several mehods. One assumpion of all hese
More informationESTIMATION OF DYNAMIC PANEL DATA MODELS WHEN REGRESSION COEFFICIENTS AND INDIVIDUAL EFFECTS ARE TIME-VARYING
Inernaional Journal of Social Science and Economic Research Volume:02 Issue:0 ESTIMATION OF DYNAMIC PANEL DATA MODELS WHEN REGRESSION COEFFICIENTS AND INDIVIDUAL EFFECTS ARE TIME-VARYING Chung-ki Min Professor
More informationACE 562 Fall Lecture 8: The Simple Linear Regression Model: R 2, Reporting the Results and Prediction. by Professor Scott H.
ACE 56 Fall 5 Lecure 8: The Simple Linear Regression Model: R, Reporing he Resuls and Predicion by Professor Sco H. Irwin Required Readings: Griffihs, Hill and Judge. "Explaining Variaion in he Dependen
More informationA Bayesian Approach to Spectral Analysis
Chirped Signals A Bayesian Approach o Specral Analysis Chirped signals are oscillaing signals wih ime variable frequencies, usually wih a linear variaion of frequency wih ime. E.g. f() = A cos(ω + α 2
More informationAppendix to Creating Work Breaks From Available Idleness
Appendix o Creaing Work Breaks From Available Idleness Xu Sun and Ward Whi Deparmen of Indusrial Engineering and Operaions Research, Columbia Universiy, New York, NY, 127; {xs2235,ww24}@columbia.edu Sepember
More informationForecasting optimally
I) ile: Forecas Evaluaion II) Conens: Evaluaing forecass, properies of opimal forecass, esing properies of opimal forecass, saisical comparison of forecas accuracy III) Documenaion: - Diebold, Francis
More informationStochastic Model for Cancer Cell Growth through Single Forward Mutation
Journal of Modern Applied Saisical Mehods Volume 16 Issue 1 Aricle 31 5-1-2017 Sochasic Model for Cancer Cell Growh hrough Single Forward Muaion Jayabharahiraj Jayabalan Pondicherry Universiy, jayabharahi8@gmail.com
More information1. Diagnostic (Misspeci cation) Tests: Testing the Assumptions
Business School, Brunel Universiy MSc. EC5501/5509 Modelling Financial Decisions and Markes/Inroducion o Quaniaive Mehods Prof. Menelaos Karanasos (Room SS269, el. 01895265284) Lecure Noes 6 1. Diagnosic
More information