Identification of Stochastic Systems Under Multiple Operating Conditions: The Vector Dependent FP ARX Parametrization
|
|
- Kathlyn Shaw
- 5 years ago
- Views:
Transcription
1 Identification of Stochastic Systems Under Multiple Operating Conditions: The Vector Dependent FP ARX Parametrization Fotis P Kopsaftopoulos and Spilios D Fassois Abstract The problem of identifying stochastic systems under multiple operating conditions, by using excitation response signals obtained from each condition, is addressed Each operating condition is characterized by several measurable variables forming a vector operating parameter The problem is tackled within a novel framework consisting of postulated Vector dependent Functionally Pooled ARX (VFP ARX) models, proper data pooling techniques, and statistical parameter estimation Least Squares (LS) and Maximum Likelihood (ML) estimation methods are developed Their strong consistency is established and their performance characteristics are assessed via a Monte Carlo study I INTRODUCTION In conventional system identification a mathematical model representing a system at a specific operating condition is identified based upon a single data record of excitation response signals Yet, in many applications, a system may operate under different operating conditions in different intervals of time, maintaining one such condition in each interval These operating conditions affect the system characteristics, and thus its dynamics Typical examples include physiological systems under different environmental conditions, mechanical systems under different load or lubrication conditions, systems under different configurations, hydraulic systems operating under different temperatures or fluid pressures, material and structures (civil mechanical aerospace) under different environmental conditions (such as temperature and humidity), and so on In such cases it is of interest to identify a global model describing the system under any operating condition, based upon excitation response data records available from each condition It could be, perhaps, argued that this may be handled by using conventional mathematical models and customary identification techniques that could artificially split the problem into a number of seemingly unrelated subproblems and derive a model based upon a single data record at a time Nevertheless, such a solution would be both awkward and statistically suboptimal Awkwardness has to do with the fact that a potentially large number of seemingly unrelated models (one per operating condition) would be obtained Statistical suboptimality has to do with the fact that the set of FP Kopsaftopoulos and SD Fassois are with the Stochastic Mechanical Systems & Automation (SMSA) Group, Department of Mechanical & Aeronautical Engineering, University of Patras, GR 6500, Patras, Greece Tel/fax: {fkopsaf,fassois}@mechupatrasgr, Internet: sms Corresponding author identified models would be of suboptimal accuracy This is due to two reasons The first is the violation of the principle of statistical parsimony (model economy) as a large number of models would be used for representing the system This would result in a large number of estimated parameters, and thus reduced accuracy The second is the ineffective use of the information available in the totality of the data records Indeed, not all available information would be extracted, as the interrelations among the different records would be ignored as a result of separating the problem into seemingly unrelated subproblems This work aims at the postulation of a proper framework and methods for effectively tackling the problem of identifying stochastic systems under multiple operating conditions This is to be based upon three important entities: 14th Mediterranean Conference on Control and Automation (MED 006) Ancona, Italy, June 006 (a) (b) (c) A novel, Functionally Pooled (FP), stochastic model structure that explicitly allows for system modelling under multiple operating conditions via a single mathematical representation This representation uses parameters that functionally depend upon the operating condition It also uses a stochastic structure that accounts for the statistical dependencies among the different data records Data pooling techniques (see [1) for combining and optimally treating (as one entity) the data obtained from the various experiments Statistical techniques for model estimation The resulting framework is referred to as a statistical Functional Pooling framework, and the corresponding models as stochastic Functionally Pooled (FP) models A schematic representation is provided in Fig 1 The only essential practical condition for using this framework and identifying global system models is that each operating condition corresponds to a specific value of a measurable variable, henceforth referred to as the operating parameter The case of a scalar operating parameter (for instance operating temperature) is treated in a companion paper [ The present paper focuses on the case of a vector operating parameter (consisting of two or more scalars, for instance operating temperature and humidity) It should be also noted that early versions of the Functional Pooling framework, including certain simple models and estimation methods, have been already applied to practical fault diagnosis problems with very promising results The interested reader is referred to [3, [4 for details
2 C5DAE F 7 ( & G G &, 0 /IH& 8 (, J0 (<0 K&<G 6, (,+ 8L/ 6H/ > + 7 & % & ' (*)& +,- / (,+ ' / / 1 +,-46, 7 (' ' & & 9:( ;& '<0 1 )&- 7 (, 0 & )> (' +, & (6/ 8?7 (' ' & & 9@( ;& '<7 ' (A/ /AB / & (,/!#"$ Fig 1 Schematic representation of the problem showing the operating points on the (k 1, k ) plane, an excitation response data set corresponding to one particular operating point, and the VFP-ARX model structure II THE DATA SET Excitation response data records from different operating points corresponding to various values of the (vector) operating parameter are used 1 : Z NM1M { x k, y k k [k 1 k T, with t 1,, N, k 1 {k 1 1,, k 1 M 1 }, k {k 1,, k M } } In this expression t designates normalized discrete time (the corresponding analog time being t T with T standing for the sampling period), k [k 1 k T the operating parameter (without loss of generality assumed to be twodimensional), and x k, y k the excitation and response signals corresponding to k N stands for the signal length (in samples) corresponding to each single experiment (each k) A total of M 1 M experiments (one for each element of k) are performed, with the complete series covering the required range of each scalar parameter, say [kmin 1, k1 max and [kmin, k max, via the discretizations k 1 k1, 1 k, 1, km 1 1 and k k1, k,, km Hence each experiment is characterized by a specific value of k, say k [ki 1, k j This vector is, for simplicity of notation, also designated as the duplet k i,j (ki 1, k j ) (the first subscript of k i,j designating the value of k 1 and the second that of k ) III THE VFP-ARX MODEL STRUCTURE The Vector dependent Functionally Pooled AutoRegressive with exogenous excitation (VFP-ARX) model structure postulated for treating the problem is of the form: na nb y k + a i (k) y k [t i b i (k) x k [t i+w k (1) i1 i0 1 Lower case/capital bold face symbols designate column-vector/matrix quantities, respectively a i (k) w k iid N ( 0, σ w(k) ) k R () a i,j G j (k), b i (k) j1 b i,j G j (k) (3) j1 E { w ki,j w km,n [t τ } γ w [k i,j, k m,n δ[τ (4) with na, nb designating the AutoRegressive (AR) and exogenous (X) orders, respectively, x k, y k the excitation and response signals, respectively, and w k the disturbance (innovations) signal that is a white (serially uncorrelated) zero-mean with variance σ w(k) and potentially crosscorrelated with its counterparts corresponding to different experiments The symbol E{ } designates statistical expectation, δ[τ the Kronecker delta (equal to unity for τ 0 and equal to zero for τ 0), N (, ) Gaussian distribution with the indicated mean and variance, and iid stands for identically independently distributed As (3) indicates, the AR and X parameters a i (k), b i (k) are modelled as explicit functions of the vector k belonging to a p-dimensional functional subspace spanned by the (mutually independent) functions G 1 (k), G (k),, G p (k) (functional basis) The functional basis consists of polynomials of two variables (vector polynomials) obtained as crossproducts from univariate polynomials (of the Chebyshev, Legendre, Jacobi and other families [5) The constants a i,j, b i,j designate the AR and X, respectively, coefficients of projection Defining: A[B, k na 1 + a i (k)b i, B[B, k nb b i (k)b i i1 i0 where A[B, k, B[B, k are the AutoRegressive (AR) and exogenous (X) polynomials in the backshift operator B ( B j u k u k [t j ), the VFP-ARX representation of (1) is
3 rewritten as: A[B, k y k B[B, k x k + w k (5) As already mentioned, the innovations sequences w k corresponding to different operating conditions may be contemporaneously correlated, that is E{w ki,j w ki,j } σ w([k i,j and E{w ki,j w km,n } γ w [k i,j, k m,n Defining the VFP-ARX model s cross-section innovations vector as: w [ w k1,1 w k1, w k1,m w km1,m T with covariance matrix: Γw E ww T σw[k 1,1 γ w[k 1,1, k M1,M 6 4 γ w[k M1,M, k 1,1 σw[k M1,M (6) (7) 3 then the covariance matrix corresponding to the time instants t 1,, N is given by: 7 5 Γw Γw I N (8) with designating Kronecker product [6, chap 7 In the case of cross-sectionally uncorrelated innovations sequences with different variances (σw[k 1,1 σw[k 1, σw[k M1,M, groupwise heteroscedasticity), the covariance matrix is given by: Γw σw[k 1,1 I N 0 0 σw[k M1,M I N (9) In the simpler case of cross-sectionally uncorrelated innovations sequences with equal variances (σ w[k 1,1 σ w[k 1, σ w[k M1,M σ w, groupwise homoscedasticity), the covariance matrix is given by Γw σ wi NM1M with I NM1M indicating the unity matrix The representation of equations (1) (4) is referred to as a VFP-ARX model of orders (n a, n b ) and functional subspace dimensionality p, or in short a VFP-ARX(n a, n b ) p model It is parameterized in terms of the parameter vector: θ [ α i,j bi,j γw [k i,j, k m,n T i, j, m, n (10) with γ w [k i,j, k i,j σ w[k i,j The VFP-ARX representation is assumed to satisfy the following conditions: A1 Stability condition The poles of the AR polynomial (see (5)) lie inside the unit circle for all operating parameters k A Irreducibility condition The polynomials A[B, k and B[B, k are coprime (have no common factors) k A3 The input signal x k is stationary, ergodic and persistently exciting with E { x ki,j w km,n } 0 i, j, m, n IV MODEL ESTIMATION A VFP-ARX model corresponding to the true system of (1) (4) may be expressed as: na nb y k + a i (k) y k [t i b i (k) x k [t i+e k (11) i1 a i (k) i0 e k iid N ( 0, σe(k) ) k R (1) a i,j G j (k), b i (k) b i,j G j (k) (13) j1 j1 E { e ki,j e km,n [t τ } γ e [k i,j, k m,n δ[τ (14) with e k designating the model s one-step-ahead prediction error or residual (corresponding to w k ) with variance σe(k) In the general case the model s one-step-ahead prediction error (residual) sequences e k may be contemporaneously correlated, that is E{e ki,j e ki,j } σe[k i,j and E{e ki,j e km,n } γ e [k i,j, k m,n, with the model residual cross-section vector defined as e [ e ek1,1 k T M1,M The cross-section vector covariance then is: Γe E ee T σe[k 1,1 γ e[k 1,1, k M1,M γ e[k M1,M, k 1,1 σe[k M1,M and the covariance matrix for the time instants t 1,, N is given as: Γe Γe I N The VFP-ARX model estimation problem may then be stated as follows: Given the excitation response data records select an element M( θ) from the VFP-ARX model set: M { M( θ) : A[B, k, θ y k B[B, k, θ x k +e k [t, θ γ w [k i,j, k m,n E{e ki,j [t, θe km,n [t, θ}, i, j, m, n } that best fits the measured data The model identification problem is usually distinguished into two subproblems: the parameter estimation subproblem and the model structure selection subproblem The present paper focuses on the parameter estimation part, while the model structure selection subproblem is treated in a forthcoming paper [7 A A Functionally Pooled Linear Regression Framework The VFP-ARX model (11) may be rewritten as: y k [ ϕ T k gt (k) θ+e k φ T k θ+e k (15) with: ϕ k g(k) θ h y k [t 1 y k [t na i T xk x k [t nb ˆG 1(k) G p(k) T h i a 1,1 a T na,p b0,1 b nb,p
4 Pooling together the expressions of the VFP-ARX model [Equation (15) corresponding to all operating parameters k (k 1,1, k 1,,, k M1,M ) considered in the experiments (cross-sectional pooling) yields: y k1,1 y km1,m φ T k 1,1 φ T k M1,M θ+ y Φ θ + e e k1,1 e km1,m Then, following substitution of the data for t 1,, N the following expression is obtained: with: y y[1 y[n, Φ y Φ θ + e (16) Φ[1 Φ[N, e B Least Squares (LS) Based Estimation Methods e[1 e[n Using the above linear regression framework the simplest possible approach for estimating the projection coefficient vector θ is based upon minimization of the Ordinary Least Squares criterion: J OLS (θ, Z NM1M ) 1 N e T e which leads to the Ordinary Least Squares (OLS) estimator: ˆθ OLS [ Φ T Φ 1[ Φ T y (17) A more appropriate criterion for the contemporaneously correlated residual case is (in view of the Gauss-Markov theorem [8) the Weighted Least Squares (WLS) criterion: J WLS (θ, Z NM1M ) 1 N e T Γ 1 w e 1 N et Γ 1 w e with Γw, Γw given by (7) and (8), respectively This leads to the Weighted Least Squares (WLS) estimator: ˆθ WLS [ Φ T Γ 1 w Φ 1[ Φ T Γ 1 w y (18) As the covariance matrix Γw is practically unavailable, it may be consistently estimated by using the Ordinary Least Squares (OLS) estimator, thus: Γ OLS w 1 N e[t, ˆθ OLS e T [t, ˆθ OLS with e[t, ˆθ OLS designating the residuals e for θ ˆθ OLS Then: Γ OLS OLS w Γ w I N The estimator in (18) is then expressed as: ˆθ WLS [ Φ T ( Γ OLS w ) 1 Φ 1[ Φ T ( Γ OLS w ) 1 y (19) while the final residual covariance matrix is estimated as: Γ WLS w 1 e[t, N ˆθ WLS e T [t, ˆθ WLS In the case of cross-sectionally uncorrelated residual sequences with different variances (σ e[k 1,1 σ e[k 1, σ e[k M1,M, groupwise heteroscedasticity) the residual covariance matrix Γw for all k has the same form as (9) As the variances are practically unavailable, they may be consistently estimated as [9: ˆσ e(k, ˆθ OLS ) 1 N e k [t, ˆθ OLS (0) for all k, with e k [t, ˆθ OLS designating the residual sequences obtained by applying OLS The ˆθ WLS estimator is then given by (19) The final residual variance is estimated as: ˆσ w(k) ˆσ e(k, ˆθ WLS ) 1 N e k [t, ˆθ WLS (1) In the simpler case of cross-sectionally uncorrelated residual sequences with equal variances (σ e[k 1,1 σ e[k 1, σ e[k M1,M σ e, groupwise homoscedasticity) the covariance matrix is Γw σ wi NM1M with I NM1M designating the unit matrix In this case the WLS estimator coincides with its OLS counterpart The residual variance is estimated by (1) C The Maximum Likelihood (ML) Estimation Method The complete parameter vector θ is estimated as: ˆ θ ML arg max L(θ, Γw/e) θ with L( ) the natural logarithm of the conditional likelihood function [10, [11 In the general case of normally distributed and contemporaneously correlated residuals e k k [10, p 198 we have: N L(θ, Γw/e[t 1,, e[t N ) ln p(e/θ, Γw) 1 e T Γ 1 w e NM 1M ln π N ln det{γ w} () with p( ) designating the Gaussian probability density function By setting: () becomes: Λ(θ) 1 N e[t, θe T [t, θ (3) L(θ, Γw/e) N T rλ(θ)γ 1 w N ln det{γ w} NM 1M ln π (4)
5 The first derivative of (4) with respect to Γw leads to: L(θ, Γw/e) Γw N Γ 1 w Λ(θ)Γ 1 w N Γ 1 w and equating it to zero yields Γw Λ(θ) As it may be shown, L(θ, Γw/e) is maximized with respect to Γw for Γw Λ(θ) and the maximum likelihood estimate of Λ(θ) is given by (3) for the optimum value of θ that has to be determined By replacing Γw with Λ(θ) in (4) yields: L(θ/e) NM 1M (ln π + 1) N ln det{λ(θ)} (5) Maximizing equation (5) with respect to θ leads to the ML estimator: ˆθ ML arg min det{λ(θ)} (6) θ and Γw Λ(ˆθ ML ) 1 N N e[t, ˆθ ML e T [t, ˆθ ML Notice that obtaining ˆθ requires the use of iterative opti- ML mization techniques [10 In the heteroscedastic case we have: ln det{λ(θ)} ln ( σe[k 1,1, θ σe[k M1,M, θ ) ln σ e[k 1,1, θ + + ln σ e[k M1,M, θ k 1 M 1 k M k 1 k1 1 k k1 ln σ e(k, θ) (7) Maximizing (5) with respect to θ leads to the optimal value of θ (as in (6)) and: ˆσ w(k) ˆσ e(k, ˆθ ML ) 1 N In the homoscedastic case we have: e k [t, ˆθ ML (8) ln det{λ(θ)} ln [ σ e(θ) M 1M M1 M ln σ e(θ) (9) and the final residual variance is given by (8) V CONSISTENCY ANALYSIS The consistency of the OLS, WLS and ML estimators of the previous section is examined For simplicity, the case of cross-sectionally uncorrelated innovations sequences with different variances (heteroscedastic case) is considered The estimated model is assumed to have the exact structure of the true system, with the latter and the excitation signals satisfying the assumptions A1, A and A3 of section III The proofs of the theorems are provided in [7 For the Least Squares (LS) estimators of the previous section we have the following theorem: Theorem 1: Least Squares estimator consistency Let θ o be the true projection coefficient vector, w k a white zero mean process with E{w k } σ w(k) for every operating point, and E{φ k φ k T } a nonsingular matrix Then: ˆθ LS N as θ o (N ) with as designating convergence in the almost sure sense [9, pp Note that, using the Kolmogorov theorem [9, p 3, it is easily seen that: as well ˆσ w(k, N) as σ w(k) (N ) Theorem : Maximum Likelihood estimator consistency Let θo [ θ T o γ w [k i,j, k m,n be the true parameter vector, w k a normally distributed zero mean white process with E{w k } σ w(k) for every operating point, and E{φ k φ T k } a nonsingular matrix Then: ˆ θ ML N as θ o (N ) VI MONTE CARLO STUDY The effectiveness of the OLS, WLS and ML estimators for VFP-ARX models is now examined via a Monte Carlo study It is noted that the ML estimator is initialized by the WLS estimates, and makes use of the Gauss Newton nonlinear optimization scheme [10 (maximum number of iterations 500; maximum number of function evaluations 5000; termination tolerance of the loss function 10 ; termination tolerance of the estimated parameters 10 1 ) The study is based upon a VFP ARX(4,) 9 model with zero delay (b 0 0 in the exogenous polynomial) and AR, X subspaces consisting of the cross-products of the first three (hence functional dimensionality p 9) shifted Chebyshev polynomials of the second kind [5 It includes 500 runs, in each one of which the first scalar operating parameter takes 16 values (ki 1 [1, 16) and the second scalar operating parameter takes 0 values (kj [1, 0) Thus, each run includes excitation response signals (of length equal to N 104 samples) from M 1 M 30 operating conditions Each response is corrupted by random noise at the 10% standard deviation level in accordance with the ARX structure expression (innovations standard deviation over the noise free response standard deviation equal to 010) The innovations sequences corresponding to different operating conditions are cross-sectionally uncorrelated, but characterized by different variances (groupwise heteroscedasticity) Some of the true system coefficients of projection (out of a total of 54) are indicated in the second column of Table I Monte Carlo partial estimation results by the Ordinary Least Squares (OLS), the Weighted Least Squares (WLS) and the Maximum Likelihood (ML) methods are presented in Table I (mean estimates ± standard deviations) Some of these results are pictorially depicted in Fig (AR/X coefficient of projection estimates, 95% confidence intervals) As may be readily observed, the results are all very accurate All three methods provide essentially unbiased estimates, with the WLS and ML methods expectedly providing better accuracy for the coefficients of projection (smaller standard deviations, thus narrower confidence intervals)
6 TABLE I INDICATIVE MONTE CARLO ESTIMATION RESULTS FOR THE VFP-ARX(4,) 9 MODEL (SELECTED PARAMETERS; 500 RUNS PER METHOD; MEAN ESTIMATE ± STANDARD DEVIATION) COEFF TRUE OLS ESTIMATE WLS ESTIMATE ML ESTIMATE a 1, ± ± ± a 1, ± ± ± a, ± ± ± a, ± ± ± a 3, ± ± ± a 3, ± ± ± a 4, ± ± ± 0000 a 4, ± ± ± b 1, ± ± ± b 1, ± ± ± b, ± ± ± b, ± ± ± Fig 3 finally depicts the true system s frequency response magnitude versus frequency and the first scalar operating parameter k 1 (for set k ) along with their mean estimated (by the OLS, WLS and ML methods) counterparts The agreement between each estimate and the true frequency response magnitude is excellent ACKNOWLEDGMENT The authors acknowledge the funding of this work by the European Social Fund (ESF), Operational Program for Education and Vocational Training II (EPEAEK II), and particularly the program PYTHAGORAS II of the Greek Ministry of Education REFERENCES [1 T Dielman, Pooled Cross-Sectional and Time Series Data Analysis Marcel Dekker, Inc, 1989 [ J Sakellariou and S Fassois, Externally Dependent Functional ARX (FARX) Models for the Identification of Dynamical Systems Under Multiple Operating Conditions, submitted for publication [3 J Sakellariou, K Petsounis, and S Fassois, A Functional Model Based Method for Fault Detection and Identification in Stochastic Structural Systems, in Proc of the Second European Workshop on Structural Health Monitoring (SHM 004), Munich, Germany, 004, pp [4 J Sakellariou and S Fassois, Fault Detection and Identification in a Scale Aircraft Skeleton Structure via a Functional Model Based Method, in Proc of the International Conference on Noise and Vibration Engineering (ISMA 004), Leuven, Belgium, 004, pp [5 C Dunkl and Y Xu, Orthogonal Polynomials of Several Variables Cambridge University Press, 001 [6 D Bernstein, Matrix Mathematics Princeton University Press, 005 [7 F Kopsaftopoulos and S Fassois, On the Estimation of Vector Dependent Functionally Pooled ARX (VFP-ARX) Models, under preparation for publication [8 W Greene, Econometric Analysis Prentice-Hall International, Inc, 1997 [9 H White, Asymptotic Theory for Econometricians, nd ed Academic Press, 001 [10 T Soderstrom and P Stoica, System Identification Prentice Hall, 1989 [11 J Mendel, Lessons in Estimation Theory for Signal Processing, Communications and Control Prentice Hall, 1995 α 1, α,4 α 3,6 b 1, OLS WLS ML α 1, α,9 α 4,5 b 1, OLS WLS ML α,3 α 3,3 b, α 4, OLS WLS ML Fig VFP-ARX(4,) 9 true values of the coefficients of projection (dashed lines) and Monte Carlo estimates (boxes indicate mean ±196 standard deviations) based upon the OLS, WLS and ML methods (500 runs per method) Fig 3 VFP-ARX(4,) 9 based frequency response magnitude versus frequency and k 1 (k is set to k9 ): (a) true system, (b) mean OLS estimate, (c) mean WLS estimate and (d) mean ML estimate (mean parameters over 500 runs)
Vector-dependent Functionally Pooled ARX Models for the Identification of Systems Under Multiple Operating Conditions
Preprints of the 16th IFAC Symposium on System Identification The International Federation of Automatic Control Vector-dependent Functionally Pooled ARX Models for the Identification of Systems Under Multiple
More informationVibration Based Health Monitoring for a Thin Aluminum Plate: Experimental Assessment of Several Statistical Time Series Methods
Vibration Based Health Monitoring for a Thin Aluminum Plate: Experimental Assessment of Several Statistical Time Series Methods Fotis P. Kopsaftopoulos and Spilios D. Fassois Stochastic Mechanical Systems
More informationNONLINEAR INTEGRAL MINIMUM VARIANCE-LIKE CONTROL WITH APPLICATION TO AN AIRCRAFT SYSTEM
NONLINEAR INTEGRAL MINIMUM VARIANCE-LIKE CONTROL WITH APPLICATION TO AN AIRCRAFT SYSTEM D.G. Dimogianopoulos, J.D. Hios and S.D. Fassois DEPARTMENT OF MECHANICAL & AERONAUTICAL ENGINEERING GR-26500 PATRAS,
More informationParametric Output Error Based Identification and Fault Detection in Structures Under Earthquake Excitation
Parametric Output Error Based Identification and Fault Detection in Structures Under Earthquake Excitation J.S. Sakellariou and S.D. Fassois Department of Mechanical & Aeronautical Engr. GR 265 Patras,
More informationthe Functional Model Based Method
Multi-Site Damage Localization via the Functional Model Based Method Christos S. Sakaris, John S. Sakellariou and Spilios D. Fassois Stochastic Mechanical Systems & Automation (SMSA) Laboratory Department
More informationStationary or Non-Stationary Random Excitation for Vibration-Based Structural Damage Detection? An exploratory study
6th International Symposium on NDT in Aerospace, 12-14th November 2014, Madrid, Spain - www.ndt.net/app.aerondt2014 More Info at Open Access Database www.ndt.net/?id=16938 Stationary or Non-Stationary
More informationMulti Channel Output Only Identification of an Extendable Arm Structure Under Random Excitation: A comparison of parametric methods
Multi Channel Output Only Identification of an Extendable Arm Structure Under Random Excitation: A comparison of parametric methods Minas Spiridonakos and Spilios Fassois Stochastic Mechanical Systems
More informationA. Poulimenos, M. Spiridonakos, and S. Fassois
PARAMETRIC TIME-DOMAIN METHODS FOR NON-STATIONARY RANDOM VIBRATION IDENTIFICATION AND ANALYSIS: AN OVERVIEW AND COMPARISON A. Poulimenos, M. Spiridonakos, and S. Fassois DEPARTMENT OF MECHANICAL & AERONAUTICAL
More informationOnboard Engine FDI in Autonomous Aircraft Using Stochastic Nonlinear Modelling of Flight Signal Dependencies
Onboard Engine FDI in Autonomous Aircraft Using Stochastic Nonlinear Modelling of Flight Signal Dependencies Dimitrios G. Dimogianopoulos, John D. Hios and Spilios D. Fassois Stochastic Mechanical Systems
More informationOutput Only Parametric Identification of a Scale Cable Stayed Bridge Structure: a comparison of vector AR and stochastic subspace methods
Output Only Parametric Identification of a Scale Cable Stayed Bridge Structure: a comparison of vector AR and stochastic subspace methods Fotis P. Kopsaftopoulos, Panagiotis G. Apostolellis and Spilios
More informationFunctionally Pooled Models for the Global Identification of Stochastic Systems Under Different Pseudo Static Operating Conditions
Functionally Pooled Models for the Global Identification of Stochastic Systems Under Different Pseudo Static Operating Conditions JS Sakellariou, and SD Fassois,2 Stochastic Mechanical Systems&Automation
More informationStationary or Non-Stationary Random Excitation for Vibration-Based Structural Damage Detection? An exploratory study
Stationary or Non-Stationary Random Excitation for Vibration-Based Structural Damage Detection? An exploratory study Andriana S. GEORGANTOPOULOU & Spilios D. FASSOIS Stochastic Mechanical Systems & Automation
More informationEconometrics II - EXAM Answer each question in separate sheets in three hours
Econometrics II - EXAM Answer each question in separate sheets in three hours. Let u and u be jointly Gaussian and independent of z in all the equations. a Investigate the identification of the following
More informationCONTROL SYSTEMS, ROBOTICS, AND AUTOMATION - Vol. V - Prediction Error Methods - Torsten Söderström
PREDICTIO ERROR METHODS Torsten Söderström Department of Systems and Control, Information Technology, Uppsala University, Uppsala, Sweden Keywords: prediction error method, optimal prediction, identifiability,
More informationExpressions for the covariance matrix of covariance data
Expressions for the covariance matrix of covariance data Torsten Söderström Division of Systems and Control, Department of Information Technology, Uppsala University, P O Box 337, SE-7505 Uppsala, Sweden
More informationOUTPUT-ONLY STATISTICAL TIME SERIES METHODS FOR STRUCTURAL HEALTH MONITORING: A COMPARATIVE STUDY
7th European Workshop on Structural Health Monitoring July 8-11, 2014. La Cité, Nantes, France More Info at Open Access Database www.ndt.net/?id=17198 OUTPUT-ONLY STATISTICAL TIME SERIES METHODS FOR STRUCTURAL
More informationScalar and Vector Time Series Methods for Vibration Based Damage Diagnosis in a Scale Aircraft Skeleton Structure
Scalar and Vector Time Series Methods for Vibration Based Damage Diagnosis in a Scale Aircraft Skeleton Structure Fotis P. Kopsaftopoulos and Spilios D. Fassois Stochastic Mechanical Systems & Automation
More informationT.-C.J. Aravanis, J.S. Sakellariou and S.D. Fassois
Vibration based fault detection under variable non-measurable, operating conditions via a stochastic Functional Model method and application to railway vehicle suspensions T.-C.J. Aravanis, J.S. Sakellariou
More informationMechanical Systems and Signal Processing
Mechanical Systems and Signal Processing 39 (213) 143 161 Contents lists available at SciVerse ScienceDirect Mechanical Systems and Signal Processing journal homepage: www.elsevier.com/locate/ymssp A functional
More informationEL1820 Modeling of Dynamical Systems
EL1820 Modeling of Dynamical Systems Lecture 9 - Parameter estimation in linear models Model structures Parameter estimation via prediction error minimization Properties of the estimate: bias and variance
More informationEECE Adaptive Control
EECE 574 - Adaptive Control Basics of System Identification Guy Dumont Department of Electrical and Computer Engineering University of British Columbia January 2010 Guy Dumont (UBC) EECE574 - Basics of
More information12. Prediction Error Methods (PEM)
12. Prediction Error Methods (PEM) EE531 (Semester II, 2010) description optimal prediction Kalman filter statistical results computational aspects 12-1 Description idea: determine the model parameter
More informationEcon 583 Final Exam Fall 2008
Econ 583 Final Exam Fall 2008 Eric Zivot December 11, 2008 Exam is due at 9:00 am in my office on Friday, December 12. 1 Maximum Likelihood Estimation and Asymptotic Theory Let X 1,...,X n be iid random
More informationStatistics 910, #5 1. Regression Methods
Statistics 910, #5 1 Overview Regression Methods 1. Idea: effects of dependence 2. Examples of estimation (in R) 3. Review of regression 4. Comparisons and relative efficiencies Idea Decomposition Well-known
More informationLinear Regression. Junhui Qian. October 27, 2014
Linear Regression Junhui Qian October 27, 2014 Outline The Model Estimation Ordinary Least Square Method of Moments Maximum Likelihood Estimation Properties of OLS Estimator Unbiasedness Consistency Efficiency
More informationOn Input Design for System Identification
On Input Design for System Identification Input Design Using Markov Chains CHIARA BRIGHENTI Masters Degree Project Stockholm, Sweden March 2009 XR-EE-RT 2009:002 Abstract When system identification methods
More informationFor more information about how to cite these materials visit
Author(s): Kerby Shedden, Ph.D., 2010 License: Unless otherwise noted, this material is made available under the terms of the Creative Commons Attribution Share Alike 3.0 License: http://creativecommons.org/licenses/by-sa/3.0/
More informationARIMA Modelling and Forecasting
ARIMA Modelling and Forecasting Economic time series often appear nonstationary, because of trends, seasonal patterns, cycles, etc. However, the differences may appear stationary. Δx t x t x t 1 (first
More informationLecture 2: ARMA(p,q) models (part 2)
Lecture 2: ARMA(p,q) models (part 2) Florian Pelgrin University of Lausanne, École des HEC Department of mathematics (IMEA-Nice) Sept. 2011 - Jan. 2012 Florian Pelgrin (HEC) Univariate time series Sept.
More informationNON-STATIONARY MECHANICAL VIBRATION MODELING AND ANALYSIS
NON-STATIONARY MECHANICAL VIBRATION MODELING AND ANALYSIS VIA FUNCTIONAL SERIES TARMA MODELS A.G. Poulimenos and S.D. Fassois DEPARTMENT OF MECHANICAL &AERONAUTICAL ENGINEERING GR-26500 PATRAS, GREECE
More informationMultivariate Regression
Multivariate Regression The so-called supervised learning problem is the following: we want to approximate the random variable Y with an appropriate function of the random variables X 1,..., X p with the
More informationBusiness Economics BUSINESS ECONOMICS. PAPER No. : 8, FUNDAMENTALS OF ECONOMETRICS MODULE No. : 3, GAUSS MARKOV THEOREM
Subject Business Economics Paper No and Title Module No and Title Module Tag 8, Fundamentals of Econometrics 3, The gauss Markov theorem BSE_P8_M3 1 TABLE OF CONTENTS 1. INTRODUCTION 2. ASSUMPTIONS OF
More informationNon-stationary functional series modeling and analysis of hardware reliability series: a comparative study using rail vehicle interfailure times
Reliability Engineering and System Safety 68 (2000) 169 183 www.elsevier.com/locate/ress Non-stationary functional series modeling and analysis of hardware reliability series: a comparative study using
More informationIdentification of ARX, OE, FIR models with the least squares method
Identification of ARX, OE, FIR models with the least squares method CHEM-E7145 Advanced Process Control Methods Lecture 2 Contents Identification of ARX model with the least squares minimizing the equation
More informationTime series methods for fault detection and identification in vibrating structures
Time series methods for fault detection and identification in vibrating structures By Spilios D. Fassois and John S. Sakellariou Stochastic Mechanical Systems (SMS) Group Department of Mechanical & Aeronautical
More informationA Primer on Asymptotics
A Primer on Asymptotics Eric Zivot Department of Economics University of Washington September 30, 2003 Revised: October 7, 2009 Introduction The two main concepts in asymptotic theory covered in these
More informationAssociation studies and regression
Association studies and regression CM226: Machine Learning for Bioinformatics. Fall 2016 Sriram Sankararaman Acknowledgments: Fei Sha, Ameet Talwalkar Association studies and regression 1 / 104 Administration
More informationTime Series Analysis. James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY
Time Series Analysis James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY PREFACE xiii 1 Difference Equations 1.1. First-Order Difference Equations 1 1.2. pth-order Difference Equations 7
More informationParametric Techniques Lecture 3
Parametric Techniques Lecture 3 Jason Corso SUNY at Buffalo 22 January 2009 J. Corso (SUNY at Buffalo) Parametric Techniques Lecture 3 22 January 2009 1 / 39 Introduction In Lecture 2, we learned how to
More informationIntroduction to Estimation Methods for Time Series models. Lecture 1
Introduction to Estimation Methods for Time Series models Lecture 1 Fulvio Corsi SNS Pisa Fulvio Corsi Introduction to Estimation () Methods for Time Series models Lecture 1 SNS Pisa 1 / 19 Estimation
More informationEstimation, Detection, and Identification
Estimation, Detection, and Identification Graduate Course on the CMU/Portugal ECE PhD Program Spring 2008/2009 Chapter 5 Best Linear Unbiased Estimators Instructor: Prof. Paulo Jorge Oliveira pjcro @ isr.ist.utl.pt
More informationParametric Techniques
Parametric Techniques Jason J. Corso SUNY at Buffalo J. Corso (SUNY at Buffalo) Parametric Techniques 1 / 39 Introduction When covering Bayesian Decision Theory, we assumed the full probabilistic structure
More informationSGN Advanced Signal Processing: Lecture 8 Parameter estimation for AR and MA models. Model order selection
SG 21006 Advanced Signal Processing: Lecture 8 Parameter estimation for AR and MA models. Model order selection Ioan Tabus Department of Signal Processing Tampere University of Technology Finland 1 / 28
More informationy k = ( ) x k + v k. w q wk i 0 0 wk
Four telling examples of Kalman Filters Example : Signal plus noise Measurement of a bandpass signal, center frequency.2 rad/sec buried in highpass noise. Dig out the quadrature part of the signal while
More informationLinear models. Linear models are computationally convenient and remain widely used in. applied econometric research
Linear models Linear models are computationally convenient and remain widely used in applied econometric research Our main focus in these lectures will be on single equation linear models of the form y
More informationLTI Systems, Additive Noise, and Order Estimation
LTI Systems, Additive oise, and Order Estimation Soosan Beheshti, Munther A. Dahleh Laboratory for Information and Decision Systems Department of Electrical Engineering and Computer Science Massachusetts
More informationMS&E 226: Small Data. Lecture 11: Maximum likelihood (v2) Ramesh Johari
MS&E 226: Small Data Lecture 11: Maximum likelihood (v2) Ramesh Johari ramesh.johari@stanford.edu 1 / 18 The likelihood function 2 / 18 Estimating the parameter This lecture develops the methodology behind
More informationTAKEHOME FINAL EXAM e iω e 2iω e iω e 2iω
ECO 513 Spring 2015 TAKEHOME FINAL EXAM (1) Suppose the univariate stochastic process y is ARMA(2,2) of the following form: y t = 1.6974y t 1.9604y t 2 + ε t 1.6628ε t 1 +.9216ε t 2, (1) where ε is i.i.d.
More informationEstimation theory. Parametric estimation. Properties of estimators. Minimum variance estimator. Cramer-Rao bound. Maximum likelihood estimators
Estimation theory Parametric estimation Properties of estimators Minimum variance estimator Cramer-Rao bound Maximum likelihood estimators Confidence intervals Bayesian estimation 1 Random Variables Let
More informationIdentification Methods for Structural Systems
Prof. Dr. Eleni Chatzi Lecture 13-29 May, 2013 Courtesy of Prof. S. Fassois & Dr. F. Kopsaftopoulos, SMSA Group, University of Patras Statistical methods for SHM courtesy of Prof. S. Fassois & Dr. F. Kopsaftopoulos,
More informationIf we want to analyze experimental or simulated data we might encounter the following tasks:
Chapter 1 Introduction If we want to analyze experimental or simulated data we might encounter the following tasks: Characterization of the source of the signal and diagnosis Studying dependencies Prediction
More informationCOMPLEX CONSTRAINED CRB AND ITS APPLICATION TO SEMI-BLIND MIMO AND OFDM CHANNEL ESTIMATION. Aditya K. Jagannatham and Bhaskar D.
COMPLEX CONSTRAINED CRB AND ITS APPLICATION TO SEMI-BLIND MIMO AND OFDM CHANNEL ESTIMATION Aditya K Jagannatham and Bhaskar D Rao University of California, SanDiego 9500 Gilman Drive, La Jolla, CA 92093-0407
More informationModel structure. Lecture Note #3 (Chap.6) Identification of time series model. ARMAX Models and Difference Equations
System Modeling and Identification Lecture ote #3 (Chap.6) CHBE 70 Korea University Prof. Dae Ryoo Yang Model structure ime series Multivariable time series x [ ] x x xm Multidimensional time series (temporal+spatial)
More informationThe Bootstrap: Theory and Applications. Biing-Shen Kuo National Chengchi University
The Bootstrap: Theory and Applications Biing-Shen Kuo National Chengchi University Motivation: Poor Asymptotic Approximation Most of statistical inference relies on asymptotic theory. Motivation: Poor
More informationEconometrics Summary Algebraic and Statistical Preliminaries
Econometrics Summary Algebraic and Statistical Preliminaries Elasticity: The point elasticity of Y with respect to L is given by α = ( Y/ L)/(Y/L). The arc elasticity is given by ( Y/ L)/(Y/L), when L
More informationCOMPARISON OF GMM WITH SECOND-ORDER LEAST SQUARES ESTIMATION IN NONLINEAR MODELS. Abstract
Far East J. Theo. Stat. 0() (006), 179-196 COMPARISON OF GMM WITH SECOND-ORDER LEAST SQUARES ESTIMATION IN NONLINEAR MODELS Department of Statistics University of Manitoba Winnipeg, Manitoba, Canada R3T
More informationLecture: Simultaneous Equation Model (Wooldridge s Book Chapter 16)
Lecture: Simultaneous Equation Model (Wooldridge s Book Chapter 16) 1 2 Model Consider a system of two regressions y 1 = β 1 y 2 + u 1 (1) y 2 = β 2 y 1 + u 2 (2) This is a simultaneous equation model
More informationFurther Results on Model Structure Validation for Closed Loop System Identification
Advances in Wireless Communications and etworks 7; 3(5: 57-66 http://www.sciencepublishinggroup.com/j/awcn doi:.648/j.awcn.735. Further esults on Model Structure Validation for Closed Loop System Identification
More informationTitle. Description. var intro Introduction to vector autoregressive models
Title var intro Introduction to vector autoregressive models Description Stata has a suite of commands for fitting, forecasting, interpreting, and performing inference on vector autoregressive (VAR) models
More informationLECTURE 2 LINEAR REGRESSION MODEL AND OLS
SEPTEMBER 29, 2014 LECTURE 2 LINEAR REGRESSION MODEL AND OLS Definitions A common question in econometrics is to study the effect of one group of variables X i, usually called the regressors, on another
More informationEconomics 583: Econometric Theory I A Primer on Asymptotics
Economics 583: Econometric Theory I A Primer on Asymptotics Eric Zivot January 14, 2013 The two main concepts in asymptotic theory that we will use are Consistency Asymptotic Normality Intuition consistency:
More informationMultivariate ARMA Processes
LECTURE 8 Multivariate ARMA Processes A vector y(t) of n elements is said to follow an n-variate ARMA process of orders p and q if it satisfies the equation (1) A 0 y(t) + A 1 y(t 1) + + A p y(t p) = M
More informationBSCB 6520: Statistical Computing Homework 4 Due: Friday, December 5
BSCB 6520: Statistical Computing Homework 4 Due: Friday, December 5 All computer code should be produced in R or Matlab and e-mailed to gjh27@cornell.edu (one file for the entire homework). Responses to
More informationLinear Regression and Its Applications
Linear Regression and Its Applications Predrag Radivojac October 13, 2014 Given a data set D = {(x i, y i )} n the objective is to learn the relationship between features and the target. We usually start
More information1 Teaching notes on structural VARs.
Bent E. Sørensen November 8, 2016 1 Teaching notes on structural VARs. 1.1 Vector MA models: 1.1.1 Probability theory The simplest to analyze, estimation is a different matter time series models are the
More informationNONUNIFORM SAMPLING FOR DETECTION OF ABRUPT CHANGES*
CIRCUITS SYSTEMS SIGNAL PROCESSING c Birkhäuser Boston (2003) VOL. 22, NO. 4,2003, PP. 395 404 NONUNIFORM SAMPLING FOR DETECTION OF ABRUPT CHANGES* Feza Kerestecioğlu 1,2 and Sezai Tokat 1,3 Abstract.
More informationElements of Multivariate Time Series Analysis
Gregory C. Reinsel Elements of Multivariate Time Series Analysis Second Edition With 14 Figures Springer Contents Preface to the Second Edition Preface to the First Edition vii ix 1. Vector Time Series
More informationArma-Arch Modeling Of The Returns Of First Bank Of Nigeria
Arma-Arch Modeling Of The Returns Of First Bank Of Nigeria Emmanuel Alphonsus Akpan Imoh Udo Moffat Department of Mathematics and Statistics University of Uyo, Nigeria Ntiedo Bassey Ekpo Department of
More informationSystem Identification, Lecture 4
System Identification, Lecture 4 Kristiaan Pelckmans (IT/UU, 2338) Course code: 1RT880, Report code: 61800 - Spring 2012 F, FRI Uppsala University, Information Technology 30 Januari 2012 SI-2012 K. Pelckmans
More information1 Teaching notes on structural VARs.
Bent E. Sørensen February 22, 2007 1 Teaching notes on structural VARs. 1.1 Vector MA models: 1.1.1 Probability theory The simplest (to analyze, estimation is a different matter) time series models are
More informationSystem Identification, Lecture 4
System Identification, Lecture 4 Kristiaan Pelckmans (IT/UU, 2338) Course code: 1RT880, Report code: 61800 - Spring 2016 F, FRI Uppsala University, Information Technology 13 April 2016 SI-2016 K. Pelckmans
More informationLessons learned from the theory and practice of. change detection. Introduction. Content. Simulated data - One change (Signal and spectral densities)
Lessons learned from the theory and practice of change detection Simulated data - One change (Signal and spectral densities) - Michèle Basseville - 4 6 8 4 6 8 IRISA / CNRS, Rennes, France michele.basseville@irisa.fr
More informationChapter 4 Neural Networks in System Identification
Chapter 4 Neural Networks in System Identification Gábor HORVÁTH Department of Measurement and Information Systems Budapest University of Technology and Economics Magyar tudósok körútja 2, 52 Budapest,
More informationTime Series Analysis. James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY
Time Series Analysis James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY & Contents PREFACE xiii 1 1.1. 1.2. Difference Equations First-Order Difference Equations 1 /?th-order Difference
More informationTime Series Analysis
Time Series Analysis hm@imm.dtu.dk Informatics and Mathematical Modelling Technical University of Denmark DK-2800 Kgs. Lyngby 1 Outline of the lecture Regression based methods, 1st part: Introduction (Sec.
More informationOn Moving Average Parameter Estimation
On Moving Average Parameter Estimation Niclas Sandgren and Petre Stoica Contact information: niclas.sandgren@it.uu.se, tel: +46 8 473392 Abstract Estimation of the autoregressive moving average (ARMA)
More information1. The OLS Estimator. 1.1 Population model and notation
1. The OLS Estimator OLS stands for Ordinary Least Squares. There are 6 assumptions ordinarily made, and the method of fitting a line through data is by least-squares. OLS is a common estimation methodology
More informationLecture 2: Univariate Time Series
Lecture 2: Univariate Time Series Analysis: Conditional and Unconditional Densities, Stationarity, ARMA Processes Prof. Massimo Guidolin 20192 Financial Econometrics Spring/Winter 2017 Overview Motivation:
More informationComputer Intensive Methods in Mathematical Statistics
Computer Intensive Methods in Mathematical Statistics Department of mathematics johawes@kth.se Lecture 16 Advanced topics in computational statistics 18 May 2017 Computer Intensive Methods (1) Plan of
More informationMGR-815. Notes for the MGR-815 course. 12 June School of Superior Technology. Professor Zbigniew Dziong
Modeling, Estimation and Control, for Telecommunication Networks Notes for the MGR-815 course 12 June 2010 School of Superior Technology Professor Zbigniew Dziong 1 Table of Contents Preface 5 1. Example
More informationPeter Hoff Linear and multilinear models April 3, GLS for multivariate regression 5. 3 Covariance estimation for the GLM 8
Contents 1 Linear model 1 2 GLS for multivariate regression 5 3 Covariance estimation for the GLM 8 4 Testing the GLH 11 A reference for some of this material can be found somewhere. 1 Linear model Recall
More informationA unified framework for EIV identification methods in the presence of mutually correlated noises
Preprints of the 19th World Congress The International Federation of Automatic Control A unified framework for EIV identification methods in the presence of mutually correlated noises Torsten Söderström
More informationBias Correction Methods for Dynamic Panel Data Models with Fixed Effects
MPRA Munich Personal RePEc Archive Bias Correction Methods for Dynamic Panel Data Models with Fixed Effects Mohamed R. Abonazel Department of Applied Statistics and Econometrics, Institute of Statistical
More informationOutline. Possible Reasons. Nature of Heteroscedasticity. Basic Econometrics in Transportation. Heteroscedasticity
1/25 Outline Basic Econometrics in Transportation Heteroscedasticity What is the nature of heteroscedasticity? What are its consequences? How does one detect it? What are the remedial measures? Amir Samimi
More informationThe Algebra of the Kronecker Product. Consider the matrix equation Y = AXB where
21 : CHAPTER Seemingly-Unrelated Regressions The Algebra of the Kronecker Product Consider the matrix equation Y = AXB where Y =[y kl ]; k =1,,r,l =1,,s, (1) X =[x ij ]; i =1,,m,j =1,,n, A=[a ki ]; k =1,,r,i=1,,m,
More informationVector Auto-Regressive Models
Vector Auto-Regressive Models Laurent Ferrara 1 1 University of Paris Nanterre M2 Oct. 2018 Overview of the presentation 1. Vector Auto-Regressions Definition Estimation Testing 2. Impulse responses functions
More informationEconomics 536 Lecture 7. Introduction to Specification Testing in Dynamic Econometric Models
University of Illinois Fall 2016 Department of Economics Roger Koenker Economics 536 Lecture 7 Introduction to Specification Testing in Dynamic Econometric Models In this lecture I want to briefly describe
More informationVAR Models and Applications
VAR Models and Applications Laurent Ferrara 1 1 University of Paris West M2 EIPMC Oct. 2016 Overview of the presentation 1. Vector Auto-Regressions Definition Estimation Testing 2. Impulse responses functions
More informationReview of Classical Least Squares. James L. Powell Department of Economics University of California, Berkeley
Review of Classical Least Squares James L. Powell Department of Economics University of California, Berkeley The Classical Linear Model The object of least squares regression methods is to model and estimate
More informationp(d θ ) l(θ ) 1.2 x x x
p(d θ ).2 x 0-7 0.8 x 0-7 0.4 x 0-7 l(θ ) -20-40 -60-80 -00 2 3 4 5 6 7 θ ˆ 2 3 4 5 6 7 θ ˆ 2 3 4 5 6 7 θ θ x FIGURE 3.. The top graph shows several training points in one dimension, known or assumed to
More informationIntroduction to Econometrics
Introduction to Econometrics T H I R D E D I T I O N Global Edition James H. Stock Harvard University Mark W. Watson Princeton University Boston Columbus Indianapolis New York San Francisco Upper Saddle
More informationBasic concepts in estimation
Basic concepts in estimation Random and nonrandom parameters Definitions of estimates ML Maimum Lielihood MAP Maimum A Posteriori LS Least Squares MMS Minimum Mean square rror Measures of quality of estimates
More informationSystem Identification
System Identification Lecture : Statistical properties of parameter estimators, Instrumental variable methods Roy Smith 8--8. 8--8. Statistical basis for estimation methods Parametrised models: G Gp, zq,
More informationBootstrap Approach to Comparison of Alternative Methods of Parameter Estimation of a Simultaneous Equation Model
Bootstrap Approach to Comparison of Alternative Methods of Parameter Estimation of a Simultaneous Equation Model Olubusoye, O. E., J. O. Olaomi, and O. O. Odetunde Abstract A bootstrap simulation approach
More informationExploring Granger Causality for Time series via Wald Test on Estimated Models with Guaranteed Stability
Exploring Granger Causality for Time series via Wald Test on Estimated Models with Guaranteed Stability Nuntanut Raksasri Jitkomut Songsiri Department of Electrical Engineering, Faculty of Engineering,
More informationAn estimate of the long-run covariance matrix, Ω, is necessary to calculate asymptotic
Chapter 6 ESTIMATION OF THE LONG-RUN COVARIANCE MATRIX An estimate of the long-run covariance matrix, Ω, is necessary to calculate asymptotic standard errors for the OLS and linear IV estimators presented
More informationAdaptive Dual Control
Adaptive Dual Control Björn Wittenmark Department of Automatic Control, Lund Institute of Technology Box 118, S-221 00 Lund, Sweden email: bjorn@control.lth.se Keywords: Dual control, stochastic control,
More information1 Estimation of Persistent Dynamic Panel Data. Motivation
1 Estimation of Persistent Dynamic Panel Data. Motivation Consider the following Dynamic Panel Data (DPD) model y it = y it 1 ρ + x it β + µ i + v it (1.1) with i = {1, 2,..., N} denoting the individual
More informationDiscrete time processes
Discrete time processes Predictions are difficult. Especially about the future Mark Twain. Florian Herzog 2013 Modeling observed data When we model observed (realized) data, we encounter usually the following
More informationMarcia Gumpertz and Sastry G. Pantula Department of Statistics North Carolina State University Raleigh, NC
A Simple Approach to Inference in Random Coefficient Models March 8, 1988 Marcia Gumpertz and Sastry G. Pantula Department of Statistics North Carolina State University Raleigh, NC 27695-8203 Key Words
More information