Statistical procedures searching for suspicious measurements and their sequences in measurement series. Ing. Marek Brabec, PhD CHMU, SZU, UI Praha
|
|
- Ezra Underwood
- 5 years ago
- Views:
Transcription
1 Saisical procedures searching for suspicious measuremens and heir sequences in measuremen series Ing. Marek Brabec, PhD CHMU, SZU, UI Praha 8h AQ EIONET meeing Oslo, November 2003
2 Goals: Propose procedures for checking suspicious spos in ime series of air qualiy measuremens. Look: - no only for ouliers (large values) - bu also for more complicaed scenarios (e.g. aypical sequences, unusual dynamics) Formalize ad hoc procedures uilized by air polluion expers for manual checking Suppress subjeciviy in human exper appraisal - make rules for finding suspicious values explici, objecive and open o conrol Techniques applicable in a sofware for rouine use - help expers o find problemaic spos - wihou going hrough masses of normal daa - proposed mehods find only suspicious obs. - ulimae decision is reserved for human expers Use mehods ha can be employed in on-line mode - possible applicaions in real ime (conrol sofware of he measuremen devices)
3 Our approach is: based on saisical modelling measuremen series we use several flexible ime series models of DLM (Dynamic Linear Model) ype unforunaely, we do no have pure raining daa available wih reliable a priori classificaion problem / no- problem for all daa poins (e.g. gold sandard classificaion by expers) if such caegorizaion were available, model for problem free series would be more precise and deecion of problemaic daa would be more sharp from he heoreical poin of view, we have o derive a classifier wihou having raining classificaion consequenly, he fied model describes ypical properies of an essenially problem free series
4 Subsequenly: we look for measuremens and sequences ha are: - in some sense non-conforming wih he model developed for ypical daa, - i.e. hose ha are suspicious we check many pre-specified scenarios (i.e. pre-specified suspicious siuaions o look for) lis of suspicious scenarios - was prepared according o exper opinion abou wha are he mos frequen and mos imporan problems occurring in measuremen daabases - can be expanded based on furher exper advice o localize dubious observaions, or heir paches, we apply saisical qualiy conrol mehods (o residuals afer he all righ model)
5 Series of air qualiy measuremens have raher complex saisical properies - call for more complicaed ime series mehods when subsanial feaures are no prop. modelled, subsequen ess/checks are invalid: - inflaion of false alarms (compared o nominal) - increasing # of undeeced real problems Subsanial saisical feaures of observed series: auocorrelaion - measuremens closer in ime are more alike - more auocorrelaion for frequen measuremens (half hourly daa) - auoregressive models show close o uni roos (Dickey-Fuller problem) long erm rends - conneced o seasonal changes, long erm evoluion of he polluion level and oher facors - measuremens are ypically nonsaionary - shape of he rend is ofen complicaed (canno be described by simple parameric f.) periodiciy (day/nigh, weekdays/weekend changes, ohers) polluan concenraions are no observed direcly blurred by a measuremen error from: - device imprecision - local influences in space and ime
6 General saisical model of measuremen series: Measuremen equaion: T Y = F θ + ε Sae equaion(s): θ Gθ + η = 1 Where: Y value observed a ime. θ sae vecor (m x 1) G ime-invarian ransiion marix (m x m). T F (possibly) ime-varying vecor (1x m) 2 ε ~ N 0; σ ε observaion error, ( ) σ 2 is unknown, will be esimaed η sae error (srucural disurbance), η ~ N ( 0; ) W is unknown W To complee he model: ε and η are assumed independen (muually and across differen s) o ge sared, a priori iniial sae disribuion is θ ~ N 0 C assumed, 0 ( ; 0 ) W is specified by discouning, 1 δ = δ W GC 1 for a specific discoun facor 0 < δ < 1 G T
7 Imporan model feaures: T rue air qualiy sae ( F θ ) is no dir. observable, - wha is observed is only a noisy version - blurred by measuremen error ( ε ) decomposiion of he oal variabiliy ino: 2 - measuremen error variance ( σ ) - srucural variabiliy ( W ) 2 σ relaive amoun of srucural variabiliy w.r.. - amouns o signal-o-noise-raio - when S/N is larger, measuremens reflec rue sae of he naure closer and vice versa rue concenraion dynamics is linear (DLM) - vecor auoregressive model (VAR1) - simple Markov sysem modelled observaion Y can consis from - eiher he original measuremens - or heir ransformaion, like Y = log( Z a) effecively, shifed LN disribuion is assumed for original measuremens: Z T ( F θ ). exp( ) a = exp ε + - useful e.g. for modeling skewed daa (like O 3 )
8 Some concree model specificaions: he general DLM: - is flexible and modular in naure - easily adaped o paricular polluan s properies - unified saisical and programming framework modulariy G11 G12 T T T - block marix specificaion of G =, F = ( F 1, F 2 ) - each block can address a specific modeled feaure G flexibiliy T - locally consan model: G = 1, F = T - local linear rend: G =, F = ( 1 0) 0 i.e. linear rend, parameers randomly changing Y = µ + ε µ = µ β = β β + η η 1 1 (adapable alernaive o global linear regression) - wihin day periodiciy for half-hourly daa harmonics can be chosen o represen paricular polluan (srong d. periodiciy for O 3, vs. PM10) randomly changing rigonomeric coefficiens e.g. for k-h harmonic: a k 2πk sin + b 48 k 2πk cos 48 a b k k = a = b k, 1 k, 1 + η + η - relaionship beween (spaially close) saions T G = 1, F = Y 21 k1 k 2 G 22
9 Applying model o daa: esimae srucural parameers 2 - like σ ; δ, or elemens of W - hese are needed for subsequen seps - choices: maximum likelihood fiing Bayesian esimaion T esimae rue concenraions µ = F θ, or θ several asks: - when observaions Y1, hy are available filraion - when only Y1, l Y 1 are available one sep ahead (e.g. 30 min) predicion - when Y1,f Y, fyt for some T > are available smoohing smoohing produces revises esimaes when hisory of he process beyond ime of ineres is available we use only filraion, (one sep ahead) predicion - hese can be compued on-line (hisory up o he curren ime is available) - hese are mos useful for inspecion of suspicious measuremens in rouine regimen
10 Compuaionally, filraion, predicion can be done via Kalman algorihm swiching beween predicion/updaing seps relaively simple, linear calculaions advanageous marix formulaion dimension of marices (and inversions involved) is small (compared o sandard LS approach) due o he recursive naure of he Kalman filraion, - sorage requiremens are small - implemenaion easy and fas (yearlong, 30-min) can easily handle missing daa (jus skipping updaing sep) no only esimaes ( a, f sae and observaion predicions) - bu also heir variances can be easily obained ( R, Q sae predicion error covariance and observaion predicion error variance) - when normal disribuion of Y, θ is assumed, confidence inervals can be compued easily e.g. f ± q0.95( n 1) Q
11 Kalman filering Iniialisaion: m 0 = 0 (m x 1) vecor C 0 = I (m x m) marix S 0 = n ( Yi Y ) i= 1 n 1 2 scalar Predicion ( parameers of an a priori disribuion before a new observaion comes): Variance of one sep ahead predicion error ( Y predicion error, fromy1,, Y 1 ): T R = GC 1 G + W (m x m) marix see below for C calculaion 1 δ T W = GC 1G (due o he discoun simplificaion), hence δ 1 T R = GC 1G δ Predicion of he sae vecor one sep ahead: a = Gm 1 (m x 1) vecor predicion of he observaion one sep ahead: T f = F a scalar esimae of is variance is: T Q F R F + S scalar = 1 Updaing he parameers afer a new observaion comes (when we know obs. Y, 1,Y ): The residual afer predicing one sep ahead: e = Y f scalar auxiliary vecor: A = R F / Q (m x 1) vecor degrees of freedom for he esimae of he observaion variance: n = n 1 +1 scalar esimae of he observaion variance: 2 S 1 e S = + S 1 1 scalar n Q esimae of he sae vecor: m = a + A e esimae of is covariance marix: S T C = ( R A A Q ) S 1 This concludes one enire cycle of compuaions: calculaion reurns o R + 1 again (for he subsequen ime +1).
12 To check: how he model fis daa overall wheher paricular daum conform o he model look a: residuals afer he model or sandardized residuals e U = y y f = Q for various saisical reasons, we work wih sepahead residuals (orhogonaliy, independence) Idea: Any subsanial sysemaic feaures of he daa: - should be capured by he model - and hence should be removed from he residuals herefore, residuals should: - behave (approximaely) as a simple whie noise sequence (-disribued, under normaliy) - be ideal for checking suspicious deparures form he all righ model subsanial saisical heory available for he checks saisical qualiy conrol mehods are useful some ess amenable o graphical implemenaion, - - suiable for scanning large amouns of daa by human expers f
13 Ouline of he sraegy for suspicious runs checks: fi a model ha describes daa behavior under all righ siuaion reasonably well inspec residuals and sandardized residuals - o check he proposed model qualiy as he normal siuaion descripor (overall residual behavior should no show subsanial sysemaic feaures) - o localize suspicious observaions (showing unusually large or small residuals) - sd. res. are comparable o heor. criical values - residual checks can deec broad range of problems and deparures, bu wih smaller power han a focused es inspec model coefficiens (esimaed saes) o address more specific problems (changing/loosing periodiciy, loosing relaionship o geographically close saion, presence of a srong unidirecional rend, lack of any rend, ec.) o deec deparures from a ypical paern, use mehods analogous o saisical qualiy conrol - conrol chars wih rules, EWMA, CUSUM - applied o eiher residuals or coefficiens
14 Example of real daa (30-min. ozone, saion 777) Original daa: log(o3+0.1) cas (pulhodina) Sandardized residuals (from one-sep-ahead predicion): sandardizovane reziduum predikce o krok dopredu cas
15 ACF from residuals: ACF original daa ACF (wihou model correcing sysemaic feaures): lag ACF lag
16 Hisogram of sandardized residuals: sandardizovane reziduum z predikce o krok dopredu In addiion o residual checks, look a sae vecor elemens (behavior of model coefficiens): Variabiliy conneced wih differen harmonics (periodic model, ozone, 1011): cislo harmonicke slozky
17 Checking residuals for unusual behavior: 1. Tesing for gross, isolaed deparures from ypical behavior (large U or e ): Shewhar conrol scheme: - If U > ku, give an alarm signal (prediced value is oo small, compared o ypical rajecory) - If U < k L, give an alarm signal (prediced value is oo large, compared o ypical rajecory) - If wu < U ku or wl > U k L, give a warning signal - If wl U wu, do no do anyhing suiable choice of consans k L < wl < wu < ku - ofen k L = ku = 3, w L = wu = 2 (normaliy) - oher choices for heavy-ailed disribuions - larger absolue values of k L, ku, wl, wu number of false alarms, undeeced problems for consans seleced o saisfy prescribed false alarm rae, his is a sandard saisical es can be easily implemened in graphical form
18 Shewhar conrol char: Limis: warning, w L = wu = 2 decision, k L = ku = 3 sandardizovane reziduum cas
19 2. Tesing for slow, bu consisen deparures from ypical behavior (slow rends in U or e ): EWMA(Exponenially Weighed Moving Average) Defined recursively: E = 0 = 1,2, 0 E = ( 1 λ ) E 1 + λu for 0 < λ < 1 can be viewed as a degree of innovaion λ means subsanial E upgrade by U each ime λ leads o E smooher han U Decision scheme (wo-sided): - If E > k, give alarm signal (large residual; predicions are oo low) - If E < k, give alarm signal (small residual) - If k < E < k, do no do anyhing Choice of k : - ARL (Average Run Lengh) beween wo alarms - should be for all OK siuaion - should be for a specified problemaic scenario - analogous o ype I and ype II error raes, ( α, β ) - heoreical compuaions, ables, Crowder (1987) - Markov chain compuaions for simple siuaions - Mone Carlo simulaions for complicaed schemes CUSUM - ypically yields similar resuls, Wes, Harrison (1997)
20 EWMA: Procedure parameers: λ = 0. 5, k = 5 Decisions: +1 for large posiive residual, curren measuremen is much smaller han ypical -1 for large negaive residual Theoreical ARL s: for all OK siuaion 7.31 for sysemaic under/over-esimaion by 2 se s (0.0062% and respecively) Observed alarm rae: 0.80% wih posiive residues, 0.83% wih negaive
21 Checks can be applied no only o residuals, bu also o sae componens (esimaed model coefficiens) o address specific problems: sandardized i-h sae vecor componen, U i mi mi = C (local) lin. regr. coefficien, saion 771 on 777: ii sandardizovany regresni koeficien cas warning and decision limis as before unusually large pos./neg. coefficiens show unusual posiive/negaive relaionship beween wo geographically close saion zeros and oher values can be esed (deparure from long erm average relaionship)
22 3. Tesing for specific scenarios: Shewhar char wih addiional rules Example: implemened in CHMI for rouine checks wih doubly exponenial disribuion, wih quaniles: ( ) Q p log 2 p = for p 0. 5 = 1 ( 2( p) ) log for p 0. 5 use a simplified (random walk) model: Y = Y 1 + η work wih differences: e = Y Y 1 and heir sandardized forms: U = e σ
23 define following evens: P1 U ( Q ; ) P3 U ( Q 0.999;Q ] P2 U ( Q 0.99;Q0.999 ] Z U ( Q 0.499;Q0.501] N2 U ( Q 0.001;Q0.01] N3 U ( Q ;Q0.001] N1 ( ) U ;Q furher, define 8 conrolled siuaions A1 P1 A2 N1 A3 P3-1 and N3 A4 N3-1 and P3 A5 Z -3 and P2-2 and N2-1 and Z A6 Z -3 and N2-2 and P2-1 and Z A7 P2-1 and N2 and U U 1 k 1 A8 N2-1 and P2 and k inerpreaion in erms of original series ( Y ): A1 jump A2 jump A3 jump, followed by jump A4 jump, followed by jump A5 almos consan, jump, jump, almos cons. A6 almos consan, jump, jump, almos cons. A7 jump, jump wih seep decline A8 jump, jump wih seep rise U U 1
24 Example of an applicaion on real daa: Saion 771 (Prague, Náměsí Republiky) yearlong (2002), ozone, 30-minue measuremens Theoreical all OK ARL=11880 (abou % of checked daa should be alarms of any kind) The procedure found 6 alarms in daa (0.0356%) Observed frequencies of various alarm ypes: A1 A2 A3 A4 A5 A6 A7 A Exremely simple implemenaion in a sofware for rouine applicaions Furher rules can be added, according o air polluion specialiss experise Underlying saisical model is raher simple, can be improved: - srucured DLM, alernaive o simple differencing - more deailed model leads o improvemen in false alarm and undeeced error raes ( α, β )
25 Examples of siuaions labeled as suspicious by he procedure among real daa: Horizonal axis: ime since he alarm signalled, verical axis: O 3 O cas od poplachu O cas od poplachu
26 O cas od poplachu Beware of he imporan feaure: by design, procedure produces some false alarms in order no o miss problem spos hence, i is imporan ha an exper will check suggesed suspicious spos manually before declaring hem o be in error checking procedures produce lis of suspicious values/sequences, hey are no inended as auomaic ools for deleion of measured values he lis of suspicious values is raher shor, so ha he exper can concenrae on heir inspecion, wihou going hrough masses of normal daa
27 Special procedures for deecing specific problems: example of device jamming conrol Problem descripion provided by A-P expers: occasionally, measuring device freezes on a single value jamming occurs for various reasons, e.g. daa pre-processing, ec., he final daa are no compleely consan, bu oscillae somehow around he jammed value hence, he jamming deecion is no as rivial as checking for zero differences Formally, we consider wo models: Model 1: normal measuremens - This is a more complicaed DLM (e.g. local linear rend) Model 2: jamming - This is a simpler alernaive (locally consan model) Alhough we know abou possibiliy of: normal o jammed change (M1->M2) and he reverse change (M2->M1) For a paricular ime : one isn sure which model holds
28 Afer inroducing formal mixure m. (dynamic linear mixure model) having apriori probabiliy π ha M1 holds (and 1 π for jammed M2) one can compue a poseriori probabiliy p ha normal model (M1) holds, given he observed daa up o ime (using Bayesian approach, see e.g. Wes-Harrison (1997)) he mos imporan inspecion ool in his conex is p and is changes hroughou he ime whenever p drops below a specified hreshold P, jamming is suspeced Example: Saion 1011(Úsí nad Labem, Kočkov) wih parameers: π = 0. 95, P = % of he 2002 daa swiched he alarm on p() cas
29 Real daa example of spos where alarm was se off Horizonal axis: ime since he alarm signalled, verical axis: log(o ) log(o3+0.1) cas sandardizovane reziduum cas
30 Behavior of he procedure for simulaed daa 2 M1 is locally linear, Var η = diag( 1, 0.01 ) M2 is locally consan σ = 0. 1 for boh of hem simulaed 200 M1 poins and 200 M2 (jammed) y cas Mixure model was applied wih: = 5 π, P = , δ = 0. 5 p cas
31 Conclusions: Inroduced saisical checks for problemaic spos The ess are broader han jus checking for ouliers include ess for: - unusually large/small values - unusual rajecories - specific feaures ( spaial cor.) - paricular problems (online jamming monior) Model srucure and subsequen esing is: - explici - open o criicism - open o furher improvemen (AP exp. suggesion) - modular - flexible Procedures are based on unified framework: - saisically - from programming poin of view The goal is o provide a ool for air polluion expers when checking he daa - provide lis of suspicious siuaions - o be checked manually - replaces he need o go hrough los of normal daa (usable for years-long 30-min daa) - removes parially subjeciviy in exper judgemen Presened mehods, procedures and ools - can be programmed easily - ino a rouinely usable sofware
R 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 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 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 informationMethodology. -ratios are biased and that the appropriate critical values have to be increased by an amount. that depends on the sample size.
Mehodology. Uni Roo Tess A ime series is inegraed when i has a mean revering propery and a finie variance. I is only emporarily ou of equilibrium and is called saionary in I(0). However a ime series ha
More informationOBJECTIVES OF TIME SERIES ANALYSIS
OBJECTIVES OF TIME SERIES ANALYSIS Undersanding he dynamic or imedependen srucure of he observaions of a single series (univariae analysis) Forecasing of fuure observaions Asceraining he leading, lagging
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 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 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 informationNon-parametric techniques. Instance Based Learning. NN Decision Boundaries. Nearest Neighbor Algorithm. Distance metric important
on-parameric echniques Insance Based Learning AKA: neares neighbor mehods, non-parameric, lazy, memorybased, or case-based learning Copyrigh 2005 by David Helmbold 1 Do no fi a model (as do LTU, decision
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 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 informationNature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time.
Supplemenary Figure 1 Spike-coun auocorrelaions in ime. Normalized auocorrelaion marices are shown for each area in a daase. The marix shows he mean correlaion of he spike coun in each ime bin wih he spike
More informationLecture 3: Exponential Smoothing
NATCOR: Forecasing & Predicive Analyics Lecure 3: Exponenial Smoohing John Boylan Lancaser Cenre for Forecasing Deparmen of Managemen Science Mehods and Models Forecasing Mehod A (numerical) procedure
More informationNon-parametric techniques. Instance Based Learning. NN Decision Boundaries. Nearest Neighbor Algorithm. Distance metric important
on-parameric echniques Insance Based Learning AKA: neares neighbor mehods, non-parameric, lazy, memorybased, or case-based learning Copyrigh 2005 by David Helmbold 1 Do no fi a model (as do LDA, logisic
More informationGINI MEAN DIFFERENCE AND EWMA CHARTS. Muhammad Riaz, Department of Statistics, Quaid-e-Azam University Islamabad,
GINI MEAN DIFFEENCE AND EWMA CHATS Muhammad iaz, Deparmen of Saisics, Quaid-e-Azam Universiy Islamabad, Pakisan. E-Mail: riaz76qau@yahoo.com Saddam Akbar Abbasi, Deparmen of Saisics, Quaid-e-Azam Universiy
More informationSmoothing. Backward smoother: At any give T, replace the observation yt by a combination of observations at & before T
Smoohing Consan process Separae signal & noise Smooh he daa: Backward smooher: A an give, replace he observaion b a combinaion of observaions a & before Simple smooher : replace he curren observaion wih
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 informationEstimation of Poses with Particle Filters
Esimaion of Poses wih Paricle Filers Dr.-Ing. Bernd Ludwig Chair for Arificial Inelligence Deparmen of Compuer Science Friedrich-Alexander-Universiä Erlangen-Nürnberg 12/05/2008 Dr.-Ing. Bernd Ludwig (FAU
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 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 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 informationWednesday, November 7 Handout: Heteroskedasticity
Amhers College Deparmen of Economics Economics 360 Fall 202 Wednesday, November 7 Handou: Heeroskedasiciy Preview Review o Regression Model o Sandard Ordinary Leas Squares (OLS) Premises o Esimaion Procedures
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 informationChapter 5. Heterocedastic Models. Introduction to time series (2008) 1
Chaper 5 Heerocedasic Models Inroducion o ime series (2008) 1 Chaper 5. Conens. 5.1. The ARCH model. 5.2. The GARCH model. 5.3. The exponenial GARCH model. 5.4. The CHARMA model. 5.5. Random coefficien
More informationTime series Decomposition method
Time series Decomposiion mehod A ime series is described using a mulifacor model such as = f (rend, cyclical, seasonal, error) = f (T, C, S, e) Long- Iner-mediaed Seasonal Irregular erm erm effec, effec,
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 informationEnsamble methods: Bagging and Boosting
Lecure 21 Ensamble mehods: Bagging and Boosing Milos Hauskrech milos@cs.pi.edu 5329 Senno Square Ensemble mehods Mixure of expers Muliple base models (classifiers, regressors), each covers a differen par
More informationModeling Economic Time Series with Stochastic Linear Difference Equations
A. Thiemer, SLDG.mcd, 6..6 FH-Kiel Universiy of Applied Sciences Prof. Dr. Andreas Thiemer e-mail: andreas.hiemer@fh-kiel.de Modeling Economic Time Series wih Sochasic Linear Difference Equaions Summary:
More informationZürich. ETH Master Course: L Autonomous Mobile Robots Localization II
Roland Siegwar Margaria Chli Paul Furgale Marco Huer Marin Rufli Davide Scaramuzza ETH Maser Course: 151-0854-00L Auonomous Mobile Robos Localizaion II ACT and SEE For all do, (predicion updae / ACT),
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 informationInnova Junior College H2 Mathematics JC2 Preliminary Examinations Paper 2 Solutions 0 (*)
Soluion 3 x 4x3 x 3 x 0 4x3 x 4x3 x 4x3 x 4x3 x x 3x 3 4x3 x Innova Junior College H Mahemaics JC Preliminary Examinaions Paper Soluions 3x 3 4x 3x 0 4x 3 4x 3 0 (*) 0 0 + + + - 3 3 4 3 3 3 3 Hence x or
More informationHow 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 informationModule 2 F c i k c s la l w a s o s f dif di fusi s o i n
Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms
More informationTime series model fitting via Kalman smoothing and EM estimation in TimeModels.jl
Time series model fiing via Kalman smoohing and EM esimaion in TimeModels.jl Gord Sephen Las updaed: January 206 Conens Inroducion 2. Moivaion and Acknowledgemens....................... 2.2 Noaion......................................
More informationDistribution of Estimates
Disribuion of Esimaes From Economerics (40) Linear Regression Model Assume (y,x ) is iid and E(x e )0 Esimaion Consisency y α + βx + he esimaes approach he rue values as he sample size increases Esimaion
More informationEnsamble methods: Boosting
Lecure 21 Ensamble mehods: Boosing Milos Hauskrech milos@cs.pi.edu 5329 Senno Square Schedule Final exam: April 18: 1:00-2:15pm, in-class Term projecs April 23 & April 25: a 1:00-2:30pm in CS seminar room
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 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 informationChapter 15. Time Series: Descriptive Analyses, Models, and Forecasting
Chaper 15 Time Series: Descripive Analyses, Models, and Forecasing Descripive Analysis: Index Numbers Index Number a number ha measures he change in a variable over ime relaive o he value of he variable
More informationMonitoring and data filtering II. Dynamic Linear Models
Ouline Monioring and daa filering II. Dynamic Linear Models (Wes and Harrison, chaper 2 Updaing equaions: Kalman Filer Discoun facor as an aid o choose W Incorporae exernal informaion: Inervenion General
More informationStationary Time Series
3-Jul-3 Time Series Analysis Assoc. Prof. Dr. Sevap Kesel July 03 Saionary Time Series Sricly saionary process: If he oin dis. of is he same as he oin dis. of ( X,... X n) ( X h,... X nh) Weakly Saionary
More informationTypes of Exponential Smoothing Methods. Simple Exponential Smoothing. Simple Exponential Smoothing
M Business Forecasing Mehods Exponenial moohing Mehods ecurer : Dr Iris Yeung Room No : P79 Tel No : 788 8 Types of Exponenial moohing Mehods imple Exponenial moohing Double Exponenial moohing Brown s
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 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 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 informationAugmented Reality II - Kalman Filters - Gudrun Klinker May 25, 2004
Augmened Realiy II Kalman Filers Gudrun Klinker May 25, 2004 Ouline Moivaion Discree Kalman Filer Modeled Process Compuing Model Parameers Algorihm Exended Kalman Filer Kalman Filer for Sensor Fusion Lieraure
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 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 informationCointegration and Implications for Forecasting
Coinegraion and Implicaions for Forecasing Two examples (A) Y Y 1 1 1 2 (B) Y 0.3 0.9 1 1 2 Example B: Coinegraion Y and coinegraed wih coinegraing vecor [1, 0.9] because Y 0.9 0.3 is a saionary process
More informationLesson 2, page 1. Outline of lesson 2
Lesson 2, page Ouline of lesson 2 Inroduce he Auocorrelaion Coefficien Undersand and define saionariy Discuss ransformaion Discuss rend and rend removal C:\Kyrre\sudier\drgrad\Kurs\series\lecure 02 03022.doc,
More informationLecture 5. Time series: ECM. Bernardina Algieri Department Economics, Statistics and Finance
Lecure 5 Time series: ECM Bernardina Algieri Deparmen Economics, Saisics and Finance Conens Time Series Modelling Coinegraion Error Correcion Model Two Seps, Engle-Granger procedure Error Correcion Model
More informationAn recursive analytical technique to estimate time dependent physical parameters in the presence of noise processes
WHAT IS A KALMAN FILTER An recursive analyical echnique o esimae ime dependen physical parameers in he presence of noise processes Example of a ime and frequency applicaion: Offse beween wo clocks PREDICTORS,
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 informationDistribution of Least Squares
Disribuion of Leas Squares In classic regression, if he errors are iid normal, and independen of he regressors, hen he leas squares esimaes have an exac normal disribuion, no jus asympoic his is no rue
More informationSequential Importance Resampling (SIR) Particle Filter
Paricle Filers++ Pieer Abbeel UC Berkeley EECS Many slides adaped from Thrun, Burgard and Fox, Probabilisic Roboics 1. Algorihm paricle_filer( S -1, u, z ): 2. Sequenial Imporance Resampling (SIR) Paricle
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 informationThe electromagnetic interference in case of onboard navy ships computers - a new approach
The elecromagneic inerference in case of onboard navy ships compuers - a new approach Prof. dr. ing. Alexandru SOTIR Naval Academy Mircea cel Bărân, Fulgerului Sree, Consanţa, soiralexandru@yahoo.com Absrac.
More informationGeorey E. Hinton. University oftoronto. Technical Report CRG-TR February 22, Abstract
Parameer Esimaion for Linear Dynamical Sysems Zoubin Ghahramani Georey E. Hinon Deparmen of Compuer Science Universiy oftorono 6 King's College Road Torono, Canada M5S A4 Email: zoubin@cs.orono.edu Technical
More informationA Robust Exponentially Weighted Moving Average Control Chart for the Process Mean
Journal of Modern Applied Saisical Mehods Volume 5 Issue Aricle --005 A Robus Exponenially Weighed Moving Average Conrol Char for he Process Mean Michael B. C. Khoo Universii Sains, Malaysia, mkbc@usm.my
More informationModeling and Forecasting Volatility Autoregressive Conditional Heteroskedasticity Models. Economic Forecasting Anthony Tay Slide 1
Modeling and Forecasing Volailiy Auoregressive Condiional Heeroskedasiciy Models Anhony Tay Slide 1 smpl @all line(m) sii dl_sii S TII D L _ S TII 4,000. 3,000.1.0,000 -.1 1,000 -. 0 86 88 90 9 94 96 98
More informationDynamic Econometric Models: Y t = + 0 X t + 1 X t X t k X t-k + e t. A. Autoregressive Model:
Dynamic Economeric Models: A. Auoregressive Model: Y = + 0 X 1 Y -1 + 2 Y -2 + k Y -k + e (Wih lagged dependen variable(s) on he RHS) B. Disribued-lag Model: Y = + 0 X + 1 X -1 + 2 X -2 + + k X -k + e
More informationLecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.
Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in
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 informationAn Adaptive Generalized Likelihood Ratio Control Chart for Detecting an Unknown Mean Pattern
An adapive GLR conrol char... 1/32 An Adapive Generalized Likelihood Raio Conrol Char for Deecing an Unknown Mean Paern GIOVANNA CAPIZZI and GUIDO MASAROTTO Deparmen of Saisical Sciences Universiy of Padua
More information14 Autoregressive Moving Average Models
14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class
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 informationFinancial Econometrics Jeffrey R. Russell Midterm Winter 2009 SOLUTIONS
Name SOLUTIONS Financial Economerics Jeffrey R. Russell Miderm Winer 009 SOLUTIONS You have 80 minues o complee he exam. Use can use a calculaor and noes. Try o fi all your work in he space provided. If
More informationFailure of the work-hamiltonian connection for free energy calculations. Abstract
Failure of he work-hamilonian connecion for free energy calculaions Jose M. G. Vilar 1 and J. Miguel Rubi 1 Compuaional Biology Program, Memorial Sloan-Keering Cancer Cener, 175 York Avenue, New York,
More informationTwo Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017
Two Popular Bayesian Esimaors: Paricle and Kalman Filers McGill COMP 765 Sep 14 h, 2017 1 1 1, dx x Bel x u x P x z P Recall: Bayes Filers,,,,,,, 1 1 1 1 u z u x P u z u x z P Bayes z = observaion u =
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 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 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 informationCHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK
175 CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK 10.1 INTRODUCTION Amongs he research work performed, he bes resuls of experimenal work are validaed wih Arificial Neural Nework. From he
More informationL07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms
L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS NA568 Mobile Roboics: Mehods & Algorihms Today s Topic Quick review on (Linear) Kalman Filer Kalman Filering for Non-Linear Sysems Exended Kalman Filer (EKF)
More information2.160 System Identification, Estimation, and Learning. Lecture Notes No. 8. March 6, 2006
2.160 Sysem Idenificaion, Esimaion, and Learning Lecure Noes No. 8 March 6, 2006 4.9 Eended Kalman Filer In many pracical problems, he process dynamics are nonlinear. w Process Dynamics v y u Model (Linearized)
More information3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon
3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of
More informationSpeaker Adaptation Techniques For Continuous Speech Using Medium and Small Adaptation Data Sets. Constantinos Boulis
Speaker Adapaion Techniques For Coninuous Speech Using Medium and Small Adapaion Daa Ses Consaninos Boulis Ouline of he Presenaion Inroducion o he speaker adapaion problem Maximum Likelihood Sochasic Transformaions
More informationChapter 16. Regression with Time Series Data
Chaper 16 Regression wih Time Series Daa The analysis of ime series daa is of vial ineres o many groups, such as macroeconomiss sudying he behavior of naional and inernaional economies, finance economiss
More informationTracking. Announcements
Tracking Tuesday, Nov 24 Krisen Grauman UT Ausin Announcemens Pse 5 ou onigh, due 12/4 Shorer assignmen Auo exension il 12/8 I will no hold office hours omorrow 5 6 pm due o Thanksgiving 1 Las ime: Moion
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 informationProbabilistic Robotics
Probabilisic Roboics Bayes Filer Implemenaions Gaussian filers Bayes Filer Reminder Predicion bel p u bel d Correcion bel η p z bel Gaussians : ~ π e p N p - Univariae / / : ~ μ μ μ e p Ν p d π Mulivariae
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 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 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 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 informationExponential Smoothing
Exponenial moohing Inroducion A simple mehod for forecasing. Does no require long series. Enables o decompose he series ino a rend and seasonal effecs. Paricularly useful mehod when here is a need o forecas
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 informationExponentially Weighted Moving Average (EWMA) Chart Based on Six Delta Initiatives
hps://doi.org/0.545/mjis.08.600 Exponenially Weighed Moving Average (EWMA) Char Based on Six Dela Iniiaives KALPESH S. TAILOR Deparmen of Saisics, M. K. Bhavnagar Universiy, Bhavnagar-36400 E-mail: kalpesh_lr@yahoo.co.in
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 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 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 informationProperties of Autocorrelated Processes Economics 30331
Properies of Auocorrelaed Processes Economics 3033 Bill Evans Fall 05 Suppose we have ime series daa series labeled as where =,,3, T (he final period) Some examples are he dail closing price of he S&500,
More informationObject tracking: Using HMMs to estimate the geographical location of fish
Objec racking: Using HMMs o esimae he geographical locaion of fish 02433 - Hidden Markov Models Marin Wæver Pedersen, Henrik Madsen Course week 13 MWP, compiled June 8, 2011 Objecive: Locae fish from agging
More informationRandom Walk with Anti-Correlated Steps
Random Walk wih Ani-Correlaed Seps John Noga Dirk Wagner 2 Absrac We conjecure he expeced value of random walks wih ani-correlaed seps o be exacly. We suppor his conjecure wih 2 plausibiliy argumens and
More informationNonstationarity-Integrated Models. Time Series Analysis Dr. Sevtap Kestel 1
Nonsaionariy-Inegraed Models Time Series Analysis Dr. Sevap Kesel 1 Diagnosic Checking Residual Analysis: Whie noise. P-P or Q-Q plos of he residuals follow a normal disribuion, he series is called a Gaussian
More informationChristos Papadimitriou & Luca Trevisan November 22, 2016
U.C. Bereley CS170: Algorihms Handou LN-11-22 Chrisos Papadimiriou & Luca Trevisan November 22, 2016 Sreaming algorihms In his lecure and he nex one we sudy memory-efficien algorihms ha process a sream
More informationEUROINDICATORS WORKING GROUP. A new method for assessing direct versus indirect adjustment
EUROINDICATOR WORKING GROUP 5 TH MEETING TH & TH JUNE 0 EUROTAT C4 DOC 330/ A new mehod for assessing direc versus indirec adjusmen ITEM 4.3 ON THE AGENDA OF THE MEETING OF THE WORKING GROUP ON EUROINDICATOR
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 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