Review from last time Time Series Analysis, Fall 2007 Professor Anna Mikusheva Paul Schrimpf, scribe October 23, 2007.

Size: px
Start display at page:

Download "Review from last time Time Series Analysis, Fall 2007 Professor Anna Mikusheva Paul Schrimpf, scribe October 23, 2007."

Transcription

1 Review fro last tie 4384 ie Series Analsis, Fall 007 Professor Anna Mikusheva Paul Schrif, scribe October 3, 007 Lecture 3 Unit Roots Review fro last tie Let t be a rando walk t = ρ t + ɛ t, ρ = where ɛ t is a artingale difference sequence (E[ɛ t I t] = 0, with E[ɛ It ] σ as and E ɛ 4 t t < K < hen {ɛ t } satisf a functional central liit theore, ξ (τ = [τ ] ɛt t= σ ( e showed last tie that: σ (td (t σ (t 0 σ (t t ɛ t 3/ t t 0 0 and our OLS estiator and t-statistic have non-standard distributions: d ( ρˆ ρ d t his is quite different fro the case with ρ < If x t = ρx t + ɛ t, ρ < then, xt ɛ t x 3/ t x t σ N(0, 4 ρ Ex t = 0 Ex t = σ ρ and the OLS estiate and t-stat converge to noral distributions: (ρˆ ρ N(0, ( ρ t N(0, Cite as: Anna Mikusheva, course aterials for 4384 ie Series Analsis, Fall 007 MI OenCourseare (htt://ocwitedu, Massachusetts Institute of echnolog Downloaded on [DD Month YYYY]

2 Adding a constant Adding a constant Suose we add a constant to our OLS regression of t his is equivalent to running OLS on deeaned t Let t = t ȳ hen, ρˆ = t ɛ t ( t Consider the nuerator: t ɛ t = ( t ȳɛ t = t ɛ t ȳɛ ȳɛ = 3/ t ɛ t e know that t σ, and 3/ we know that t ɛ t σ d Cobining this, we have: ( t ɛ t σ (sd (s ( ɛt σ ( Also, fro before (t e can think of this as the integral of a deeaned Brownian otion, ( (sd (s ( (t = (s (t d (s = (sd (s he liiting distributions of all statistics (ρˆ, t, etc would change In ost cases, the change is onl in relacing b one also can include a linear trend in the regression, then the liiting distribution would deend on a detrended Brownian otion Liiting Distribution Let s R return to the case without a constant he liiting distribution of the OLS esitate is (ρˆ R d his distribution is skewed and shifted to the left If we include a constant, the distribution is even ore shifted If we also include a trend, the distribution shifts et ore For exale, the 5%-tile without a constant is -05, with a constant is -69, and with a trend is -5 he 975%-tiles are 6, 04, and -8, resectivel hus, estiates of ρ have bias of order his bias is quite large in sall sales he distribution of the t-statistic is also shifted to the left and skewed, but less so than the distribution of the OLS estiate he 5%-tiles are -3, -3, -366 and of the 975%-tiles are 6, 03, and -066 for no constant, with a constant, and with a trend, resectivel Allowing Auto-correlation So far, we have assued that ɛ t are not correlated his assution is too strong for eirical work Let, t = t + v t where v t is a zero-ean stationar rocess with an MA reresentation, v t = c j ɛ j j=0 Cite as: Anna Mikusheva, course aterials for 4384 ie Series Analsis, Fall 007 MI OenCourseare (htt://ocwitedu, Massachusetts Institute of echnolog Downloaded on [DD Month YYYY]

3 Augented Dicke Fuller 3 where ɛ j is white noise Now we need to look at the liiting distributions of each of the ters we looked at before ( v t, t v t, t, etc v t is a ter we alread encountered in HAC estiation: where ω is the long-ter variance, vt N(0, ω t= ω = γj = σ c( ω j= For linearl deendent rocess, v t, we have the following central liit theore: ξ (τ = [ τ] vt t= ω ( his can be roven b using the Beveridge-Nelson decoosition All the other ters converge to things siilar to the uncorrelated case, excet with σ relaced b ω For exale, An excetion is: = ξ (t/ ω (t/ 3/ t ω (sds t v t = his leads to an extra constant in the distribution of ρˆ: t v (ω ( σ v ω ω d + σ v (ρˆ d ω σ v ω so it is iossible to just look u critical values in a table Phillis & Perron (988 suggested corrected statistics: ωˆ σˆv d (ρˆ + where ωˆ t is an estiate of the long-run variance, sa Newe-est, and σˆv is the variance of the residuals Augented Dicke Fuller Another aroach is the augented Dicke-Fuller test Suose, t = a j t j + ɛ t j= Cite as: Anna Mikusheva, course aterials for 4384 ie Series Analsis, Fall 007 MI OenCourseare (htt://ocwitedu, Massachusetts Institute of echnolog Downloaded on [DD Month YYYY]

4 Augented Dicke Fuller 4 where z = is a root, Suose we factor a(l as a j z j =0 j= with b(l = j j= β jl e can write the odel as his suggests estiating: = a j a j L j = ( Lb(L j t = β j t + ɛ t i= t = ρ t + β j t + ɛ t j i= and testing whether ρ = to see if we have a unit root his is another wa of allowing auto-correlation because we can write the odel as: a(l t =ɛ t b(l t =ɛ t t = t + v t t =b(l ɛ t where v t = b(l ɛ t is auto-correlated he nice thing about the augeented Dicke-Fuller test is that the coefficients on t j each converge to a noral distribution at rate, and the coefficient on t converges to a non-standard distribution without nuisance araeters o see this, let x t = [ t, t,, t + ] and let θ = [ρ, β,, β ] Consider: θˆ θ = (X X (X ɛ e need to noralize this so that it converges Our noralizing atrix will be Q = 0 Multiling b Q gives: ˆ Q(θ θ =(Q (X X Q (Q X ɛ he denoinator is: Q (X X Q = t t 3/ t t 3/ t X X Cite as: Anna Mikusheva, course aterials for 4384 ie Series Analsis, Fall 007 MI OenCourseare (htt://ocwitedu, Massachusetts Institute of echnolog Downloaded on [DD Month YYYY]

5 Other ests 5 Note that t 3/ t j 0 Also, we know that t ω, so ω 0 Q (X X Q 0 EX X Siilarl, ( Q X ɛ = ( t ɛ t σω d t j ɛ t N(0, E[X X]σ hus, or ωˆ (ρˆ σˆ σ d (ρˆ ω (ρˆ = d βˆ j Other ests Sargan-Bhargrava (983 look at: For a non-stationar rocess, the variance of t grows with tie, so we could when ɛ t is an ds, this converges to a distribution, converges to 0 in robabilit t t P t P t ( For a stationar rocess, this Range (Mandelbrot and coauthors You could also look at the range of t, this should blow u for non-stationar rocesses for a rando walk, we have: est Coarison (ax t in t su (λ t λ inf (λ λ he Sargan-Bhargrava (983 test has otial astotic ower against local alternatives, but has bad size distortions in sall sales, eseciall if errors are negativel autocorrelated he augented Dicke-Fuller test with the BIC used to choose lag-length has overall good size in finite sales Cite as: Anna Mikusheva, course aterials for 4384 ie Series Analsis, Fall 007 MI OenCourseare (htt://ocwitedu, Massachusetts Institute of echnolog Downloaded on [DD Month YYYY]

13.2 Fully Polynomial Randomized Approximation Scheme for Permanent of Random 0-1 Matrices

13.2 Fully Polynomial Randomized Approximation Scheme for Permanent of Random 0-1 Matrices CS71 Randoness & Coputation Spring 018 Instructor: Alistair Sinclair Lecture 13: February 7 Disclaier: These notes have not been subjected to the usual scrutiny accorded to foral publications. They ay

More information

Change-point detection for recursive Bayesian geoacoustic inversion

Change-point detection for recursive Bayesian geoacoustic inversion Change-oint detection for recursive Baesian geoacoustic inversion Bien Aik Tan Peter Gerstoft Caglar Yardi and Willia S. Hodgkiss Universit of California San Diego htt://www.l.ucsd.edu/eole/btan/ Overview

More information

[95/95] APPROACH FOR DESIGN LIMITS ANALYSIS IN VVER. Shishkov L., Tsyganov S. Russian Research Centre Kurchatov Institute Russian Federation, Moscow

[95/95] APPROACH FOR DESIGN LIMITS ANALYSIS IN VVER. Shishkov L., Tsyganov S. Russian Research Centre Kurchatov Institute Russian Federation, Moscow [95/95] APPROACH FOR DESIGN LIMITS ANALYSIS IN VVER Shishkov L., Tsyganov S. Russian Research Centre Kurchatov Institute Russian Federation, Moscow ABSTRACT The aer discusses a well-known condition [95%/95%],

More information

Binomial and Poisson Probability Distributions

Binomial and Poisson Probability Distributions Binoial and Poisson Probability Distributions There are a few discrete robability distributions that cro u any ties in hysics alications, e.g. QM, SM. Here we consider TWO iortant and related cases, the

More information

Some simple continued fraction expansions for an in nite product Part 1. Peter Bala, January ax 4n+3 1 ax 4n+1. (a; x) =

Some simple continued fraction expansions for an in nite product Part 1. Peter Bala, January ax 4n+3 1 ax 4n+1. (a; x) = Soe sile continued fraction exansions for an in nite roduct Part. Introduction The in nite roduct Peter Bala, January 3 (a; x) = Y ax 4n+3 ax 4n+ converges for arbitrary colex a rovided jxj

More information

A Simple Regression Problem

A Simple Regression Problem A Siple Regression Proble R. M. Castro March 23, 2 In this brief note a siple regression proble will be introduced, illustrating clearly the bias-variance tradeoff. Let Y i f(x i ) + W i, i,..., n, where

More information

1 Brownian motion and the Langevin equation

1 Brownian motion and the Langevin equation Figure 1: The robust appearance of Robert Brown (1773 1858) 1 Brownian otion and the Langevin equation In 1827, while exaining pollen grains and the spores of osses suspended in water under a icroscope,

More information

Modi ed Local Whittle Estimator for Long Memory Processes in the Presence of Low Frequency (and Other) Contaminations

Modi ed Local Whittle Estimator for Long Memory Processes in the Presence of Low Frequency (and Other) Contaminations Modi ed Local Whittle Estiator for Long Meory Processes in the Presence of Low Frequency (and Other Containations Jie Hou y Boston University Pierre Perron z Boston University March 5, 203; Revised: January

More information

Questions and Answers on Unit Roots, Cointegration, VARs and VECMs

Questions and Answers on Unit Roots, Cointegration, VARs and VECMs Questions and Answers on Unit Roots, Cointegration, VARs and VECMs L. Magee Winter, 2012 1. Let ɛ t, t = 1,..., T be a series of independent draws from a N[0,1] distribution. Let w t, t = 1,..., T, be

More information

Time Series Analysis Fall 2008

Time Series Analysis Fall 2008 MIT OpenCourseWare http://ocw.mit.edu 14.384 Time Series Analysis Fall 008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Introduction 1 14.384 Time

More information

TEST OF HOMOGENEITY OF PARALLEL SAMPLES FROM LOGNORMAL POPULATIONS WITH UNEQUAL VARIANCES

TEST OF HOMOGENEITY OF PARALLEL SAMPLES FROM LOGNORMAL POPULATIONS WITH UNEQUAL VARIANCES TEST OF HOMOGENEITY OF PARALLEL SAMPLES FROM LOGNORMAL POPULATIONS WITH UNEQUAL VARIANCES S. E. Ahed, R. J. Tokins and A. I. Volodin Departent of Matheatics and Statistics University of Regina Regina,

More information

1 Introduction to Generalized Least Squares

1 Introduction to Generalized Least Squares ECONOMICS 7344, Spring 2017 Bent E. Sørensen April 12, 2017 1 Introduction to Generalized Least Squares Consider the model Y = Xβ + ɛ, where the N K matrix of regressors X is fixed, independent of the

More information

Quadratic Reciprocity. As in the previous notes, we consider the Legendre Symbol, defined by

Quadratic Reciprocity. As in the previous notes, we consider the Legendre Symbol, defined by Math 0 Sring 01 Quadratic Recirocity As in the revious notes we consider the Legendre Sybol defined by $ ˆa & 0 if a 1 if a is a quadratic residue odulo. % 1 if a is a quadratic non residue We also had

More information

Random Process Review

Random Process Review Rando Process Review Consider a rando process t, and take k saples. For siplicity, we will set k. However it should ean any nuber of saples. t () t x t, t, t We have a rando vector t, t, t. If we find

More information

Block designs and statistics

Block designs and statistics Bloc designs and statistics Notes for Math 447 May 3, 2011 The ain paraeters of a bloc design are nuber of varieties v, bloc size, nuber of blocs b. A design is built on a set of v eleents. Each eleent

More information

Risk & Safety in Engineering. Dr. Jochen Köhler

Risk & Safety in Engineering. Dr. Jochen Köhler Risk & afety in Engineering Dr. Jochen Köhler Contents of Today's Lecture Introduction to Classical Reliability Theory tructural Reliability The fundaental case afety argin Introduction to Classical Reliability

More information

Plasma-Wall Interaction: Sheath and Pre-sheath

Plasma-Wall Interaction: Sheath and Pre-sheath Plasa-Wall Interaction: Sheath and Pre-sheath Under ost conditions, a very thin negative sheath appears in the vicinity of walls, due to accuulation of electrons on the wall. This is in turn necessitated

More information

1 Proof of learning bounds

1 Proof of learning bounds COS 511: Theoretical Machine Learning Lecturer: Rob Schapire Lecture #4 Scribe: Akshay Mittal February 13, 2013 1 Proof of learning bounds For intuition of the following theore, suppose there exists a

More information

5. Dimensional Analysis. 5.1 Dimensions and units

5. Dimensional Analysis. 5.1 Dimensions and units 5. Diensional Analysis In engineering the alication of fluid echanics in designs ake uch of the use of eirical results fro a lot of exerients. This data is often difficult to resent in a readable for.

More information

EE5900 Spring Lecture 4 IC interconnect modeling methods Zhuo Feng

EE5900 Spring Lecture 4 IC interconnect modeling methods Zhuo Feng EE59 Spring Parallel LSI AD Algoriths Lecture I interconnect odeling ethods Zhuo Feng. Z. Feng MTU EE59 So far we ve considered only tie doain analyses We ll soon see that it is soeties preferable to odel

More information

Variance Decomposition

Variance Decomposition Variance Decomposition 1 14.384 Time Series Analysis, Fall 2007 Recitation by Paul Schrimpf Supplementary to lectures given by Anna Mikusheva October 5, 2007 Recitation 5 Variance Decomposition Suppose

More information

Machine Learning Basics: Estimators, Bias and Variance

Machine Learning Basics: Estimators, Bias and Variance Machine Learning Basics: Estiators, Bias and Variance Sargur N. srihari@cedar.buffalo.edu This is part of lecture slides on Deep Learning: http://www.cedar.buffalo.edu/~srihari/cse676 1 Topics in Basics

More information

The Schrödinger Equation and the Scale Principle

The Schrödinger Equation and the Scale Principle Te Scrödinger Equation and te Scale Princile RODOLFO A. FRINO Jul 014 Electronics Engineer Degree fro te National Universit of Mar del Plata - Argentina rodolfo_frino@aoo.co.ar Earlier tis ear (Ma) I wrote

More information

CSE525: Randomized Algorithms and Probabilistic Analysis May 16, Lecture 13

CSE525: Randomized Algorithms and Probabilistic Analysis May 16, Lecture 13 CSE55: Randoied Algoriths and obabilistic Analysis May 6, Lecture Lecturer: Anna Karlin Scribe: Noah Siegel, Jonathan Shi Rando walks and Markov chains This lecture discusses Markov chains, which capture

More information

1. (2.5.1) So, the number of moles, n, contained in a sample of any substance is equal N n, (2.5.2)

1. (2.5.1) So, the number of moles, n, contained in a sample of any substance is equal N n, (2.5.2) Lecture.5. Ideal gas law We have already discussed general rinciles of classical therodynaics. Classical therodynaics is a acroscoic science which describes hysical systes by eans of acroscoic variables,

More information

ESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS. A Thesis. Presented to. The Faculty of the Department of Mathematics

ESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS. A Thesis. Presented to. The Faculty of the Department of Mathematics ESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS A Thesis Presented to The Faculty of the Departent of Matheatics San Jose State University In Partial Fulfillent of the Requireents

More information

arxiv: v1 [math.ds] 19 Jun 2012

arxiv: v1 [math.ds] 19 Jun 2012 Rates in the strong invariance rincile for ergodic autoorhiss of the torus Jérôe Dedecker a, Florence Merlevède b and Françoise Pène c 1 a Université Paris Descartes, Sorbonne Paris Cité, Laboratoire MAP5

More information

Estimating Parameters for a Gaussian pdf

Estimating Parameters for a Gaussian pdf Pattern Recognition and achine Learning Jaes L. Crowley ENSIAG 3 IS First Seester 00/0 Lesson 5 7 Noveber 00 Contents Estiating Paraeters for a Gaussian pdf Notation... The Pattern Recognition Proble...3

More information

W-BASED VS LATENT VARIABLES SPATIAL AUTOREGRESSIVE MODELS: EVIDENCE FROM MONTE CARLO SIMULATIONS

W-BASED VS LATENT VARIABLES SPATIAL AUTOREGRESSIVE MODELS: EVIDENCE FROM MONTE CARLO SIMULATIONS W-BASED VS LATENT VARIABLES SPATIAL AUTOREGRESSIVE MODELS: EVIDENCE FROM MONTE CARLO SIMULATIONS. Introduction When it coes to applying econoetric odels to analyze georeferenced data, researchers are well

More information

arxiv: v2 [math.st] 11 Dec 2018

arxiv: v2 [math.st] 11 Dec 2018 esting for high-diensional network paraeters in auto-regressive odels arxiv:803659v [aths] Dec 08 Lili Zheng and Garvesh Raskutti Abstract High-diensional auto-regressive odels provide a natural way to

More information

Keywords: Estimator, Bias, Mean-squared error, normality, generalized Pareto distribution

Keywords: Estimator, Bias, Mean-squared error, normality, generalized Pareto distribution Testing approxiate norality of an estiator using the estiated MSE and bias with an application to the shape paraeter of the generalized Pareto distribution J. Martin van Zyl Abstract In this work the norality

More information

Online Supplement. and. Pradeep K. Chintagunta Graduate School of Business University of Chicago

Online Supplement. and. Pradeep K. Chintagunta Graduate School of Business University of Chicago Online Suppleent easuring Cross-Category Price Effects with Aggregate Store Data Inseong Song ksong@usthk Departent of arketing, HKUST Hong Kong University of Science and Technology and Pradeep K Chintagunta

More information

P (t) = P (t = 0) + F t Conclusion: If we wait long enough, the velocity of an electron will diverge, which is obviously impossible and wrong.

P (t) = P (t = 0) + F t Conclusion: If we wait long enough, the velocity of an electron will diverge, which is obviously impossible and wrong. 4 Phys520.nb 2 Drude theory ~ Chapter in textbook 2.. The relaxation tie approxiation Here we treat electrons as a free ideal gas (classical) 2... Totally ignore interactions/scatterings Under a static

More information

SUPPORTING INFORMATION FOR. Mass Spectrometrically-Detected Statistical Aspects of Ligand Populations in Mixed Monolayer Au 25 L 18 Nanoparticles

SUPPORTING INFORMATION FOR. Mass Spectrometrically-Detected Statistical Aspects of Ligand Populations in Mixed Monolayer Au 25 L 18 Nanoparticles SUPPORTIG IFORMATIO FOR Mass Sectroetrically-Detected Statistical Asects of Lig Poulations in Mixed Monolayer Au 25 L 8 anoarticles Aala Dass,,a Kennedy Holt, Joseh F. Parer, Stehen W. Feldberg, Royce

More information

A method to determine relative stroke detection efficiencies from multiplicity distributions

A method to determine relative stroke detection efficiencies from multiplicity distributions A ethod to deterine relative stroke detection eiciencies ro ultiplicity distributions Schulz W. and Cuins K. 2. Austrian Lightning Detection and Inoration Syste (ALDIS), Kahlenberger Str.2A, 90 Vienna,

More information

Thus, the formula weight per C-atom is: 1*12 (g/mol) +1.85*1(g/mol)+0.55*16 (g/mol)+0.13*14(g/mol)=24.5 g/mol.

Thus, the formula weight per C-atom is: 1*12 (g/mol) +1.85*1(g/mol)+0.55*16 (g/mol)+0.13*14(g/mol)=24.5 g/mol. Proble. a. C: H: O: N47.6%/: 7.%/: (-47.6%-7.%-7.%-.%)/6: 7.%/4:.85:.55:. Therefore, the eleental coposition for the ash-free bioass is CH.85 O.55 N.. Thus, the forula weight per C-ato is: * (g/ol).85*(g/ol).55*6

More information

Pseudo-marginal Metropolis-Hastings: a simple explanation and (partial) review of theory

Pseudo-marginal Metropolis-Hastings: a simple explanation and (partial) review of theory Pseudo-arginal Metropolis-Hastings: a siple explanation and (partial) review of theory Chris Sherlock Motivation Iagine a stochastic process V which arises fro soe distribution with density p(v θ ). Iagine

More information

General Physical Chemistry I

General Physical Chemistry I General Physical Cheistry I Lecture 12 Aleksey Kocherzhenko Aril 2, 2015" Last tie " Gibbs free energy" In order to analyze the sontaneity of cheical reactions, we need to calculate the entroy changes

More information

Introduction to Robotics (CS223A) (Winter 2006/2007) Homework #5 solutions

Introduction to Robotics (CS223A) (Winter 2006/2007) Homework #5 solutions Introduction to Robotics (CS3A) Handout (Winter 6/7) Hoework #5 solutions. (a) Derive a forula that transfors an inertia tensor given in soe frae {C} into a new frae {A}. The frae {A} can differ fro frae

More information

Extension of CSRSM for the Parametric Study of the Face Stability of Pressurized Tunnels

Extension of CSRSM for the Parametric Study of the Face Stability of Pressurized Tunnels Extension of CSRSM for the Paraetric Study of the Face Stability of Pressurized Tunnels Guilhe Mollon 1, Daniel Dias 2, and Abdul-Haid Soubra 3, M.ASCE 1 LGCIE, INSA Lyon, Université de Lyon, Doaine scientifique

More information

Random Process Examples 1/23

Random Process Examples 1/23 ando Process Eaples /3 E. #: D-T White Noise Let ] be a sequence of V s where each V ] in the sequence is uncorrelated with all the others: E{ ] ] } 0 for This DEFINES a DT White Noise Also called Uncorrelated

More information

Pattern Recognition and Machine Learning. Artificial Neural networks

Pattern Recognition and Machine Learning. Artificial Neural networks Pattern Recognition and Machine Learning Jaes L. Crowley ENSIMAG 3 - MMIS Fall Seester 2016 Lessons 7 14 Dec 2016 Outline Artificial Neural networks Notation...2 1. Introduction...3... 3 The Artificial

More information

Vibrations: Two Degrees of Freedom Systems - Wilberforce Pendulum and Bode Plots

Vibrations: Two Degrees of Freedom Systems - Wilberforce Pendulum and Bode Plots .003J/.053J Dynaics and Control, Spring 007 Professor Peacock 5/4/007 Lecture 3 Vibrations: Two Degrees of Freedo Systes - Wilberforce Pendulu and Bode Plots Wilberforce Pendulu Figure : Wilberforce Pendulu.

More information

More Empirical Process Theory

More Empirical Process Theory More Empirical Process heory 4.384 ime Series Analysis, Fall 2008 Recitation by Paul Schrimpf Supplementary to lectures given by Anna Mikusheva October 24, 2008 Recitation 8 More Empirical Process heory

More information

Bayesian inference for stochastic differential mixed effects models - initial steps

Bayesian inference for stochastic differential mixed effects models - initial steps Bayesian inference for stochastic differential ixed effects odels - initial steps Gavin Whitaker 2nd May 2012 Supervisors: RJB and AG Outline Mixed Effects Stochastic Differential Equations (SDEs) Bayesian

More information

Automated Frequency Domain Decomposition for Operational Modal Analysis

Automated Frequency Domain Decomposition for Operational Modal Analysis Autoated Frequency Doain Decoposition for Operational Modal Analysis Rune Brincker Departent of Civil Engineering, University of Aalborg, Sohngaardsholsvej 57, DK-9000 Aalborg, Denark Palle Andersen Structural

More information

STOPPING SIMULATED PATHS EARLY

STOPPING SIMULATED PATHS EARLY Proceedings of the 2 Winter Siulation Conference B.A.Peters,J.S.Sith,D.J.Medeiros,andM.W.Rohrer,eds. STOPPING SIMULATED PATHS EARLY Paul Glasseran Graduate School of Business Colubia University New Yor,

More information

AN ACCURACY OF ASYMPTOTIC FORMULAS IN CALCULATIONS OF A RANDOM NETWORK RELIABILATY. Tsitsiashvili G.Sh., Losev A.S.

AN ACCURACY OF ASYMPTOTIC FORMULAS IN CALCULATIONS OF A RANDOM NETWORK RELIABILATY. Tsitsiashvili G.Sh., Losev A.S. Tsitsiashvili G, Losev A AN ACCUACY OF ASYMPTOTIC FOMULAS IN CALCULATIONS OF A ANDOM NETWOK ELIABILATY &ATA # 3 (Vol) 009, Seteber AN ACCUACY OF ASYMPTOTIC FOMULAS IN CALCULATIONS OF A ANDOM NETWOK ELIABILATY

More information

Course Notes for EE227C (Spring 2018): Convex Optimization and Approximation

Course Notes for EE227C (Spring 2018): Convex Optimization and Approximation Course Notes for EE7C (Spring 018: Convex Optiization and Approxiation Instructor: Moritz Hardt Eail: hardt+ee7c@berkeley.edu Graduate Instructor: Max Sichowitz Eail: sichow+ee7c@berkeley.edu October 15,

More information

Principles of Optimal Control Spring 2008

Principles of Optimal Control Spring 2008 MIT OpenCourseWare http://ocw.it.edu 16.323 Principles of Optial Control Spring 2008 For inforation about citing these aterials or our Ters of Use, visit: http://ocw.it.edu/ters. 16.323 Lecture 10 Singular

More information

Some Proofs: This section provides proofs of some theoretical results in section 3.

Some Proofs: This section provides proofs of some theoretical results in section 3. Testing Jups via False Discovery Rate Control Yu-Min Yen. Institute of Econoics, Acadeia Sinica, Taipei, Taiwan. E-ail: YMYEN@econ.sinica.edu.tw. SUPPLEMENTARY MATERIALS Suppleentary Materials contain

More information

Problem Set 2 Due Sept, 21

Problem Set 2 Due Sept, 21 EE6: Rando Processes in Sstes Lecturer: Jean C. Walrand Proble Set Due Sept, Fall 6 GSI: Assane Guee This proble set essentiall reviews notions of conditional epectation, conditional distribution, and

More information

Generalized Method of Moments (GMM) Estimation

Generalized Method of Moments (GMM) Estimation Econometrics 2 Fall 2004 Generalized Method of Moments (GMM) Estimation Heino Bohn Nielsen of29 Outline of the Lecture () Introduction. (2) Moment conditions and methods of moments (MM) estimation. Ordinary

More information

OPENING THE NEURAL NETWORK BLACKBOX: AN ALGORITHM FOR EXTRACTING RULES FROM FUNCTION APPROXIMATING ARTIFICIAL NEURAL NETWORKS

OPENING THE NEURAL NETWORK BLACKBOX: AN ALGORITHM FOR EXTRACTING RULES FROM FUNCTION APPROXIMATING ARTIFICIAL NEURAL NETWORKS OENING THE NEURAL NETWORK BLACKBOX: AN ALGORITHM FOR EXTRACTING RULES FROM FUNCTION AROXIMATING ARTIFICIAL NEURAL NETWORKS Rud Setiono Wee Kheng Leow National Universit of Singaore Singaore Jaes Y.L. Thong

More information

PERIODIC STEADY STATE ANALYSIS, EFFECTIVE VALUE,

PERIODIC STEADY STATE ANALYSIS, EFFECTIVE VALUE, PERIODIC SEADY SAE ANALYSIS, EFFECIVE VALUE, DISORSION FACOR, POWER OF PERIODIC CURRENS t + Effective value of current (general definition) IRMS i () t dt Root Mean Square, in Czech boo denoted I he value

More information

Seismic Analysis of Structures by TK Dutta, Civil Department, IIT Delhi, New Delhi.

Seismic Analysis of Structures by TK Dutta, Civil Department, IIT Delhi, New Delhi. Seisic Analysis of Structures by K Dutta, Civil Departent, II Delhi, New Delhi. Module 5: Response Spectru Method of Analysis Exercise Probles : 5.8. or the stick odel of a building shear frae shown in

More information

Questions and Answers on Heteroskedasticity, Autocorrelation and Generalized Least Squares

Questions and Answers on Heteroskedasticity, Autocorrelation and Generalized Least Squares Questions and Answers on Heteroskedasticity, Autocorrelation and Generalized Least Squares L Magee Fall, 2008 1 Consider a regression model y = Xβ +ɛ, where it is assumed that E(ɛ X) = 0 and E(ɛɛ X) =

More information

Classical systems in equilibrium

Classical systems in equilibrium 35 Classical systes in equilibriu Ideal gas Distinguishable particles Here we assue that every particle can be labeled by an index i... and distinguished fro any other particle by its label if not by any

More information

The proofs of Theorem 1-3 are along the lines of Wied and Galeano (2013).

The proofs of Theorem 1-3 are along the lines of Wied and Galeano (2013). A Appendix: Proofs The proofs of Theore 1-3 are along the lines of Wied and Galeano (2013) Proof of Theore 1 Let D[d 1, d 2 ] be the space of càdlàg functions on the interval [d 1, d 2 ] equipped with

More information

Optimal Jackknife for Discrete Time and Continuous Time Unit Root Models

Optimal Jackknife for Discrete Time and Continuous Time Unit Root Models Optial Jackknife for Discrete Tie and Continuous Tie Unit Root Models Ye Chen and Jun Yu Singapore Manageent University January 6, Abstract Maxiu likelihood estiation of the persistence paraeter in the

More information

Symmetrization and Rademacher Averages

Symmetrization and Rademacher Averages Stat 928: Statistical Learning Theory Lecture: Syetrization and Radeacher Averages Instructor: Sha Kakade Radeacher Averages Recall that we are interested in bounding the difference between epirical and

More information

arxiv: v1 [math.pr] 17 May 2009

arxiv: v1 [math.pr] 17 May 2009 A strong law of large nubers for artingale arrays Yves F. Atchadé arxiv:0905.2761v1 [ath.pr] 17 May 2009 March 2009 Abstract: We prove a artingale triangular array generalization of the Chow-Birnbau- Marshall

More information

Ufuk Demirci* and Feza Kerestecioglu**

Ufuk Demirci* and Feza Kerestecioglu** 1 INDIRECT ADAPTIVE CONTROL OF MISSILES Ufuk Deirci* and Feza Kerestecioglu** *Turkish Navy Guided Missile Test Station, Beykoz, Istanbul, TURKEY **Departent of Electrical and Electronics Engineering,

More information

arxiv: v4 [math.st] 9 Aug 2017

arxiv: v4 [math.st] 9 Aug 2017 PARALLELIZING SPECTRAL ALGORITHMS FOR KERNEL LEARNING GILLES BLANCHARD AND NICOLE MÜCKE arxiv:161007487v4 [athst] 9 Aug 2017 Abstract We consider a distributed learning aroach in suervised learning for

More information

Statistical Techniques II

Statistical Techniques II Statistical Techniques II EST705 Regression with atrix Algebra 06a_atrix SLR atrix Algebra We will not be doing our regressions with matrix algebra, except that the computer does employ matrices. In fact,

More information

A GENERAL THEORY OF PARTICLE FILTERS IN HIDDEN MARKOV MODELS AND SOME APPLICATIONS. By Hock Peng Chan National University of Singapore and

A GENERAL THEORY OF PARTICLE FILTERS IN HIDDEN MARKOV MODELS AND SOME APPLICATIONS. By Hock Peng Chan National University of Singapore and Subitted to the Annals of Statistics A GENERAL THEORY OF PARTICLE FILTERS IN HIDDEN MARKOV MODELS AND SOME APPLICATIONS By Hock Peng Chan National University of Singaore and By Tze Leung Lai Stanford University

More information

CS Lecture 13. More Maximum Likelihood

CS Lecture 13. More Maximum Likelihood CS 6347 Lecture 13 More Maxiu Likelihood Recap Last tie: Introduction to axiu likelihood estiation MLE for Bayesian networks Optial CPTs correspond to epirical counts Today: MLE for CRFs 2 Maxiu Likelihood

More information

Linear Regression with Time Series Data

Linear Regression with Time Series Data u n i v e r s i t y o f c o p e n h a g e n d e p a r t m e n t o f e c o n o m i c s Econometrics II Linear Regression with Time Series Data Morten Nyboe Tabor u n i v e r s i t y o f c o p e n h a g

More information

Generalized eigenfunctions and a Borel Theorem on the Sierpinski Gasket.

Generalized eigenfunctions and a Borel Theorem on the Sierpinski Gasket. Generalized eigenfunctions and a Borel Theore on the Sierpinski Gasket. Kasso A. Okoudjou, Luke G. Rogers, and Robert S. Strichartz May 26, 2006 1 Introduction There is a well developed theory (see [5,

More information

Solutions of some selected problems of Homework 4

Solutions of some selected problems of Homework 4 Solutions of soe selected probles of Hoework 4 Sangchul Lee May 7, 2018 Proble 1 Let there be light A professor has two light bulbs in his garage. When both are burned out, they are replaced, and the next

More information

A Sequential Tracking Filter without Requirement of Measurement Decorrelation

A Sequential Tracking Filter without Requirement of Measurement Decorrelation A Sequential racing Filter without Requient of Measuent Decorlation Gongjian Zhou, Changjun Yu, aifan Quan Deartent of Electronic Engineering Harbin Institute of echnology Harbin, China zhougj@hit.edu.cn

More information

Lecture 3: October 2, 2017

Lecture 3: October 2, 2017 Inforation and Coding Theory Autun 2017 Lecturer: Madhur Tulsiani Lecture 3: October 2, 2017 1 Shearer s lea and alications In the revious lecture, we saw the following stateent of Shearer s lea. Lea 1.1

More information

3.3 Variational Characterization of Singular Values

3.3 Variational Characterization of Singular Values 3.3. Variational Characterization of Singular Values 61 3.3 Variational Characterization of Singular Values Since the singular values are square roots of the eigenvalues of the Heritian atrices A A and

More information

Chapter 8 Markov Chains and Some Applications ( 馬哥夫鏈 )

Chapter 8 Markov Chains and Some Applications ( 馬哥夫鏈 ) Chater 8 arkov Chains and oe Alications ( 馬哥夫鏈 Consider a sequence of rando variables,,, and suose that the set of ossible values of these rando variables is {,,,, }, which is called the state sace. It

More information

Lecture 21. Interior Point Methods Setup and Algorithm

Lecture 21. Interior Point Methods Setup and Algorithm Lecture 21 Interior Point Methods In 1984, Kararkar introduced a new weakly polynoial tie algorith for solving LPs [Kar84a], [Kar84b]. His algorith was theoretically faster than the ellipsoid ethod and

More information

The Distribution of the Covariance Matrix for a Subset of Elliptical Distributions with Extension to Two Kurtosis Parameters

The Distribution of the Covariance Matrix for a Subset of Elliptical Distributions with Extension to Two Kurtosis Parameters journal of ultivariate analysis 58, 96106 (1996) article no. 0041 The Distribution of the Covariance Matrix for a Subset of Elliptical Distributions with Extension to Two Kurtosis Paraeters H. S. Steyn

More information

Order Recursion Introduction Order versus Time Updates Matrix Inversion by Partitioning Lemma Levinson Algorithm Interpretations Examples

Order Recursion Introduction Order versus Time Updates Matrix Inversion by Partitioning Lemma Levinson Algorithm Interpretations Examples Order Recursion Introduction Order versus Tie Updates Matrix Inversion by Partitioning Lea Levinson Algorith Interpretations Exaples Introduction Rc d There are any ways to solve the noral equations Solutions

More information

Week 5 Quantitative Analysis of Financial Markets Characterizing Cycles

Week 5 Quantitative Analysis of Financial Markets Characterizing Cycles Week 5 Quantitative Analysis of Financial Markets Characterizing Cycles Christopher Ting http://www.mysmu.edu/faculty/christophert/ Christopher Ting : christopherting@smu.edu.sg : 6828 0364 : LKCSB 5036

More information

Slide10. Haykin Chapter 8: Principal Components Analysis. Motivation. Principal Component Analysis: Variance Probe

Slide10. Haykin Chapter 8: Principal Components Analysis. Motivation. Principal Component Analysis: Variance Probe Slide10 Motivation Haykin Chapter 8: Principal Coponents Analysis 1.6 1.4 1.2 1 0.8 cloud.dat 0.6 CPSC 636-600 Instructor: Yoonsuck Choe Spring 2015 0.4 0.2 0 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 How can we

More information

Topic 5a Introduction to Curve Fitting & Linear Regression

Topic 5a Introduction to Curve Fitting & Linear Regression /7/08 Course Instructor Dr. Rayond C. Rup Oice: A 337 Phone: (95) 747 6958 E ail: rcrup@utep.edu opic 5a Introduction to Curve Fitting & Linear Regression EE 4386/530 Coputational ethods in EE Outline

More information

Broadband Synthetic Aperture Matched Field Geoacoustic Inversion

Broadband Synthetic Aperture Matched Field Geoacoustic Inversion Broadband Snthetic Aerture Matched Field Geoacoustic Inversion PhD Candidate: Bien Aik Tan htt://www.l.ucsd.edu/eole/btan/ PhD Coittee: Prof. Willia Hodgkiss Chair Prof. Peter Gerstoft Co-chair Prof. Willia

More information

Model Mis-specification

Model Mis-specification Model Mis-specification Carlo Favero Favero () Model Mis-specification 1 / 28 Model Mis-specification Each specification can be interpreted of the result of a reduction process, what happens if the reduction

More information

Parallelizing Spectrally Regularized Kernel Algorithms

Parallelizing Spectrally Regularized Kernel Algorithms Journal of Machine Learning Research 19 (2018) 1-29 Subitted 11/16; Revised 8/18; Published 8/18 Parallelizing Sectrally Regularized Kernel Algoriths Nicole Mücke nicole.uecke@atheatik.uni-stuttgart.de

More information

Fundamentals of Astrodynamics and Applications 3 rd Ed

Fundamentals of Astrodynamics and Applications 3 rd Ed Fundaentals of Astrodynaics and Alications 3 rd Ed Errata June 0, 0 This listing is an on-going docuent of corrections and clarifications encountered in the book. I areciate any coents and questions you

More information

Numerical Method for Obtaining a Predictive Estimator for the Geometric Distribution

Numerical Method for Obtaining a Predictive Estimator for the Geometric Distribution British Journal of Matheatics & Couter Science 19(5): 1-13, 2016; Article no.bjmcs.29941 ISSN: 2231-0851 SCIENCEDOMAIN international www.sciencedoain.org Nuerical Method for Obtaining a Predictive Estiator

More information

Approximation by Piecewise Constants on Convex Partitions

Approximation by Piecewise Constants on Convex Partitions Aroxiation by Piecewise Constants on Convex Partitions Oleg Davydov Noveber 4, 2011 Abstract We show that the saturation order of iecewise constant aroxiation in L nor on convex artitions with N cells

More information

Learnability of Gaussians with flexible variances

Learnability of Gaussians with flexible variances Learnability of Gaussians with flexible variances Ding-Xuan Zhou City University of Hong Kong E-ail: azhou@cityu.edu.hk Supported in part by Research Grants Council of Hong Kong Start October 20, 2007

More information

COS 424: Interacting with Data. Written Exercises

COS 424: Interacting with Data. Written Exercises COS 424: Interacting with Data Hoework #4 Spring 2007 Regression Due: Wednesday, April 18 Written Exercises See the course website for iportant inforation about collaboration and late policies, as well

More information

Quantum algorithms (CO 781, Winter 2008) Prof. Andrew Childs, University of Waterloo LECTURE 15: Unstructured search and spatial search

Quantum algorithms (CO 781, Winter 2008) Prof. Andrew Childs, University of Waterloo LECTURE 15: Unstructured search and spatial search Quantu algoriths (CO 781, Winter 2008) Prof Andrew Childs, University of Waterloo LECTURE 15: Unstructured search and spatial search ow we begin to discuss applications of quantu walks to search algoriths

More information

Functions: Review of Algebra and Trigonometry

Functions: Review of Algebra and Trigonometry Sec. and. Functions: Review o Algebra and Trigonoetry A. Functions and Relations DEFN Relation: A set o ordered pairs. (,y) (doain, range) DEFN Function: A correspondence ro one set (the doain) to anther

More information

The degree of a typical vertex in generalized random intersection graph models

The degree of a typical vertex in generalized random intersection graph models Discrete Matheatics 306 006 15 165 www.elsevier.co/locate/disc The degree of a typical vertex in generalized rando intersection graph odels Jerzy Jaworski a, Michał Karoński a, Dudley Stark b a Departent

More information

Proc. of the IEEE/OES Seventh Working Conference on Current Measurement Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES

Proc. of the IEEE/OES Seventh Working Conference on Current Measurement Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES Proc. of the IEEE/OES Seventh Working Conference on Current Measureent Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES Belinda Lipa Codar Ocean Sensors 15 La Sandra Way, Portola Valley, CA 98 blipa@pogo.co

More information

Non-Stationary Time Series and Unit Root Testing

Non-Stationary Time Series and Unit Root Testing Econometrics II Non-Stationary Time Series and Unit Root Testing Morten Nyboe Tabor Course Outline: Non-Stationary Time Series and Unit Root Testing 1 Stationarity and Deviation from Stationarity Trend-Stationarity

More information

CHAPTER 2 THERMODYNAMICS

CHAPTER 2 THERMODYNAMICS CHAPER 2 HERMODYNAMICS 2.1 INRODUCION herodynaics is the study of the behavior of systes of atter under the action of external fields such as teerature and ressure. It is used in articular to describe

More information

1 Bounding the Margin

1 Bounding the Margin COS 511: Theoretical Machine Learning Lecturer: Rob Schapire Lecture #12 Scribe: Jian Min Si March 14, 2013 1 Bounding the Margin We are continuing the proof of a bound on the generalization error of AdaBoost

More information

Analyzing Simulation Results

Analyzing Simulation Results Analyzing Siulation Results Dr. John Mellor-Cruey Departent of Coputer Science Rice University johnc@cs.rice.edu COMP 528 Lecture 20 31 March 2005 Topics for Today Model verification Model validation Transient

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.436J/15.085J Fall 2008 Lecture 11 10/15/2008 ABSTRACT INTEGRATION I

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.436J/15.085J Fall 2008 Lecture 11 10/15/2008 ABSTRACT INTEGRATION I MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.436J/15.085J Fall 2008 Lecture 11 10/15/2008 ABSTRACT INTEGRATION I Contents 1. Preliinaries 2. The ain result 3. The Rieann integral 4. The integral of a nonnegative

More information

Heteroskedasticity and Autocorrelation Consistent Standard Errors

Heteroskedasticity and Autocorrelation Consistent Standard Errors NBER Summer Institute Minicourse What s New in Econometrics: ime Series Lecture 9 July 6, 008 Heteroskedasticity and Autocorrelation Consistent Standard Errors Lecture 9, July, 008 Outline. What are HAC

More information

2 Q 10. Likewise, in case of multiple particles, the corresponding density in 2 must be averaged over all

2 Q 10. Likewise, in case of multiple particles, the corresponding density in 2 must be averaged over all Lecture 6 Introduction to kinetic theory of plasa waves Introduction to kinetic theory So far we have been odeling plasa dynaics using fluid equations. The assuption has been that the pressure can be either

More information

Multi-Dimensional Hegselmann-Krause Dynamics

Multi-Dimensional Hegselmann-Krause Dynamics Multi-Diensional Hegselann-Krause Dynaics A. Nedić Industrial and Enterprise Systes Engineering Dept. University of Illinois Urbana, IL 680 angelia@illinois.edu B. Touri Coordinated Science Laboratory

More information