José Alberto Mauricio. Instituto Complutense de Análisis Económico. Universidad Complutense de Madrid. Campus de Somosaguas Madrid - SPAIN
|
|
- Stella Francis
- 5 years ago
- Views:
Transcription
1 A CORRECTED ALGORITHM FOR COMPUTING THE THEORETICAL AUTOCOVARIANCE MATRICES OF A VECTOR ARMA MODEL José Alberto Mauricio Instituto Comlutense de Análisis Económico Universidad Comlutense de Madrid Camus de Somosaguas Madrid - SPAIN Abstract: The algorithm of Kohn and Ansley (1982) is reconsidered here, in order to correct several imlementation errors concerning the construction of the linear equations that must be solved for comuting the theoretical autocovariance matrices of a vector ARMA model. This note resents a concise descrition of the corrected algorithm. Resumen: En esta nota se corrigen algunos errores del algoritmo de Kohn y Ansley (1982), que tienen que ver con la construcción de un sistema de ecuaciones lineales ara calcular las matrices de autocovarianzas teóricas de un modelo ARMA multivariante. March 1995
2 1. INTRODUCTION Consider the stationary vector ARMA(,q) model: Φ(B) w t = Θ(B)a t, (1) where Φ(B) =I m Φ 1 B Φ 2 B 2 Φ B, Θ(B) =I m Θ 1 B Θ 2 B 2 Θ q B q, Φ i (i = 1,..., ) and Θ i (i = 1,..., q) are m m arameter matrices, I m is the identity matrix of order m, µ is the m 1 mean vector, w t = w t µ is the m 1 vector of mean-corrected observations at time t, B is the backshift oerator (B k z t = z t k ), and the a t 's are m 1 random vectors identically and indeendently distributed, with E[a t ]=0, E[a t a T t ]=Σ (symmetric and ositive definite) and E[a t at t+k ]=0 (k /= 0). The theoretical autocovariance matrices of the stationary rocess w t generated by (1) are defined by Γ k = E[ w t w T t+k ]=E[ w t k w T t ]=Γ T k, and can be shown to be given by the following exression (see, for examle, Jenkins and Alavi 1981): Γ k = Γ k 1 Φ T Γ k Φ T + Λ k Λ k 1 Θ T 1 Λ k q Θ T q (k 0), (2) where Λ i = E[ w t a T t +i ]=E[ w t i at t ] is the theoretical cross-covariance matrix of order i between w t and a t. The comutational roblem arises in the calculation of the Γ k 's, since, although the Λ i 's need to be obtained too, their calculation is straightforward (see Section 2). 1
3 Both the theoretical autocovariance and cross-covariance matrices are needed, for examle, to comute the objective function during exact maximum likelihood estimation of (1) (see, among others, Hall and Nicholls 1980; Shea 1989; and Mauricio 1995). Other alications include the generation of indeendent realizations from a vector ARMA model (see, for examle, Shea 1988), which is articularly useful in Monte Carlo studies. Thus, as Kohn and Ansley (1982) oint out, it is imortant to have a fast method of obtaining those matrices, since they are recalculated many times during both estimation and simulation runs. The method of Kohn and Ansley (1982) rovides a fast means of comuting the autocovariance matrices. In articular, it is faster than the methods of Hall and Nicholls (1980) and Ansley (1980). However, the imlementation in the algorithmic section of Kohn and Ansley (1982, ) is incorrect. In articular, the guidelines on the construction of the linear equations, that must be solved in order to comute the first autocovariance matrices, are wrong. A corrected version of that algorithm, based on the theoretical background of Section 2, is resented in Section 3. The construction of both the matrix and the right-hand-side vector of the system of linear equations for the first autocovariance matrices, is described in detail. Some further comments in Section 4 comlete this note. 2. THEORETICAL BACKGROUND In order to motivate the algorithm of Section 3, a concise summary of the method roosed by Kohn and Ansley (1982), using the notation introduced in Section 1, is now given. First, ostmultilying (1) by a T t +i, taking exectations, and noting that E[ w t i a T t ]=E[a t at t +i ]=0 for i > 0, it is straightforward to show that Λ i = Θ i Σ i + Φ j Λ j i (i 1), (3) where Φ j = 0 if j >, Λ 0 = Σ, and Λ i = 0 for i > 0. Then, writing out (2) for k = 0 and for k 1: 2
4 where Γ 0 = Γ i Φ T i q + Λ 0 Λ j Θ T j, Γ k = W k + Γ k j Φ T j (k 1), q W k = Λ k j Θ T j (k 1). j=k (4) (5) (6) Noting that Γ i = Γ T i and using (5), with k = i, in (4), it turns out that where Γ 0 = W 0 + Φ i Γ i j Φ T j, W 0 = Σ (B + B T q )+ Θ j ΣΘ T j, B = q Φ i Λ i j Θ T j. j=i (7) (8) (9) Finally, since Γ k = Γ T k, it is clear that, when (5) and (6) are written for k = 1,..., 1, only Γ 0 through Γ 1 (unknown) and Λ 0 through Λ q+1 (known) aear in equations (5) through (9). Thus, the following system of linear equations may be secified in order to calculate the first unknown theoretical autocovariance matrices: Γ 0 Φ i Γ 0 Φ T i k 1 Γ k Γ i Φ T k i 1 k i [ Φ i+j Γ i Φ T j + Φ j Γ T i Φ T i+j ]=W 0 Γ T i Φ T k+i = W k i=0 (k = 1,..., 1) (10) (11) Together, exressions (10) and (11) form a system of m 2 linear equations with m 2 unknowns (the m 2 elements of each of the matrices Γ k, k = 0, 1,..., 1). But since Γ 0 is symmetric, it contains only m(m+1)/2 different elements. Thus, letting Γ 0 and W 0 be the diagonals and uer triangles of Γ 0 and W 0, resectively, the m(m+1)/2 + m 2 ( 1) unknowns in (10) and (11) are the (unique) solution x of the linear 3
5 system Ax = b, where x = vec( Γ 0, Γ 1,..., Γ 1 ), b = vec( W 0, W 1,..., W 1 ) and A is the matrix of coefficients, which is constructed as described in the next section. 3. THE CORRECTED ALGORITHM In order to concentrate on the more relevant asects of the algorithm (the construction of A and b), it is assumed that Λ 0 and Λ 1 through Λ q+1, as well as W 0, have already been comuted using (3), (8) and (9). Now, it is shown how to construct A and b row by row. First, initialize A and b to zero; then, execute the following two stes: Ste 1. Comute the first m(m+1)/2 rows: FOR j =1TOm FOR i =1TOj row = j( j 1)/2 + i Ste 1.1. Comute the first m(m+1)/2 columns within row: FOR l =1TOm FOR k =1TOl col = l(l 1)/2 + k IF k = l ELSE A(row,col) = Φ r (i,k)φ r ( j,l) r=1 A(row,col) = [ Φ r (i,k)φ r ( j,l) +Φ r (i,l)φ r ( j,k) ] r=1 4
6 Ste 1.2. Comute the remaining m 2 ( 1) columns within row: FOR s =1TO 1 FOR l =1TOm FOR k =1TOm col = m(m+1)/2 + m 2 (s 1) + m(l 1) + k s A(row,col) = [ Φ r+s (i,k)φ r ( j,l) +Φ r+s ( j,k)φ r (i,l) ] r=1 Ste 1.3. Set u diagonal of A and right-hand-side b: A(row,row) =1+A(row,row) b(row) = W 0 (i,j) Ste 2. Comute the remaining m 2 ( 1) rows: FOR s =1TO 1 FOR i =1TOm FOR j =1TOm row = m(m+1)/2 + m 2 (s 1) + m(i 1) + j Ste 2.1. Comute the first m(m+1)/2 columns within row: FOR l =1TOm IF l j col = j( j 1)/2 + l ELSE col = l(l 1)/2 + j A(row,col) = Φ s (i,l) 5
7 Ste 2.2. Comute the remaining m 2 ( 1) columns within row: FOR r =1TO 1 FOR l =1TOm col = m(m+1)/2 + m 2 (r 1) + m( j 1) + l IF r+s A(row,col) = Φ r+s (i,l) IF s > r col = m(m+1)/2 + m 2 (r 1) + m(l 1) + j A(row,col) =A(row,col) Φ s r (i,l) Ste 2.3. Set u diagonal of A and right-hand-side b: A(row,row) =1+A(row,row) FOR h = s TO q m b(row) =b(row) Λ s h ( j,k)θ h (i,k) k=1 END. Once A and b are constructed, the solution x = A 1 b may be obtained using standard numerical linear algebra routines (see, for examle, Moler 1972; and Press et al. 1992). Then, the elements of the theoretical autocovariance matrices can be extracted from x as follows: Γ 0 (i,j) x j( j 1)/2+i (i = 1,..., m; j = i,..., m) ; Γ k (i,j) x m(m+1)/2+m 2 (k 1)+m( j 1)+i (i = 1,..., m; j = 1,..., m ),k = 1,..., 1. For k, the autocovariances can be comuted recursively using (5), (6) and (3). 6
8 4. CONCLUDING REMARKS Note that, excet for the first m(m+1)/2 coefficients in each of the first m(m+1)/2 rows of A (ste 1.1 above), the imlementation resented here differs from that of Kohn and Ansley (1982), which is incorrect. To note this, comuter rograms were written by the author, following the imlementation in Kohn and Ansley (1982) and the one described in Section 3. The theoretical autocovariance matrices, calculated through the imlementation in Kohn and Ansley (1982) for a wide range of vector ARMA models, always showed substantial differences with the results obtained through the imlementation described here. Furthermore, the latter results always agreed with those obtained through standard rocedures, such as that of Hall and Nicholls (1980) as imlemented in Shea (1989). REFERENCES Ansley, C. F. (1980), Comutation of the Theoretical Autocovariance Function for a Vector ARMA Process, Journal of Statistical Comutation and Simulation, 12, Hall, A. D. and Nicholls, D. F. (1980), The Evaluation of Exact Maximum Likelihood Estimates for VARMA Models, Journal of Statistical Comutation and Simulation, 10, Jenkins, G. M. and Alavi, A. S. (1981), Some Asects of Modeling and Forecasting Multivariate Time Series, Journal of Time Series Analysis, 2, Kohn, R. and Ansley, C. F. (1982), A Note on Obtaining the Theoretical Autocovariances of an ARMA Process, Journal of Statistical Comutation and Simulation, 15, Mauricio, J. A. (1995), Exact Maximum Likelihood Estimation of Stationary Vector ARMA Models, Journal of the American Statistical Association, 90, Moler, C. B. (1972), Algorithm 423: Linear Equation Solver, Communications of the Association for Comuting Machinery, 15,
9 Press, W. H., Teukolsky, S. A., Vetterling, W. T., and Flannery, B. P (1992), Numerical Recies in C (2nd edition), Cambridge: Cambridge University Press. Shea, B. L. (1988), A Note on the Generation of Indeendent Realizations of a Vector Autoregressive Moving-Average Process, Journal of Time Series Analysis, 9, Shea, B. L. (1989), Algorithm AS 242: The Exact Likelihood of a Vector Autoregressive Moving-Average Model, Alied Statistics, 38,
ALTERNATIVE SOLUTION TO THE QUARTIC EQUATION by Farid A. Chouery 1, P.E. 2006, All rights reserved
ALTERNATIVE SOLUTION TO THE QUARTIC EQUATION b Farid A. Chouer, P.E. 006, All rights reserved Abstract A new method to obtain a closed form solution of the fourth order olnomial equation is roosed in this
More informationNumerical Linear Algebra
Numerical Linear Algebra Numerous alications in statistics, articularly in the fitting of linear models. Notation and conventions: Elements of a matrix A are denoted by a ij, where i indexes the rows and
More informationarxiv: v1 [physics.data-an] 26 Oct 2012
Constraints on Yield Parameters in Extended Maximum Likelihood Fits Till Moritz Karbach a, Maximilian Schlu b a TU Dortmund, Germany, moritz.karbach@cern.ch b TU Dortmund, Germany, maximilian.schlu@cern.ch
More informationUse of Transformations and the Repeated Statement in PROC GLM in SAS Ed Stanek
Use of Transformations and the Reeated Statement in PROC GLM in SAS Ed Stanek Introduction We describe how the Reeated Statement in PROC GLM in SAS transforms the data to rovide tests of hyotheses of interest.
More informationResiduals in Time Series Models
Residuals in Time Series Models José Alberto Mauricio Universidad Complutense de Madrid, Facultad de Económicas, Campus de Somosaguas, 83 Madrid, Spain. (E-mail: jamauri@ccee.ucm.es.) Summary: Three types
More informationUsing the Divergence Information Criterion for the Determination of the Order of an Autoregressive Process
Using the Divergence Information Criterion for the Determination of the Order of an Autoregressive Process P. Mantalos a1, K. Mattheou b, A. Karagrigoriou b a.deartment of Statistics University of Lund
More informationOn Line Parameter Estimation of Electric Systems using the Bacterial Foraging Algorithm
On Line Parameter Estimation of Electric Systems using the Bacterial Foraging Algorithm Gabriel Noriega, José Restreo, Víctor Guzmán, Maribel Giménez and José Aller Universidad Simón Bolívar Valle de Sartenejas,
More informationEstimation of the large covariance matrix with two-step monotone missing data
Estimation of the large covariance matrix with two-ste monotone missing data Masashi Hyodo, Nobumichi Shutoh 2, Takashi Seo, and Tatjana Pavlenko 3 Deartment of Mathematical Information Science, Tokyo
More informationSystem Reliability Estimation and Confidence Regions from Subsystem and Full System Tests
009 American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June 0-, 009 FrB4. System Reliability Estimation and Confidence Regions from Subsystem and Full System Tests James C. Sall Abstract
More informationApproximating min-max k-clustering
Aroximating min-max k-clustering Asaf Levin July 24, 2007 Abstract We consider the roblems of set artitioning into k clusters with minimum total cost and minimum of the maximum cost of a cluster. The cost
More informationFinite Mixture EFA in Mplus
Finite Mixture EFA in Mlus November 16, 2007 In this document we describe the Mixture EFA model estimated in Mlus. Four tyes of deendent variables are ossible in this model: normally distributed, ordered
More informationTowards understanding the Lorenz curve using the Uniform distribution. Chris J. Stephens. Newcastle City Council, Newcastle upon Tyne, UK
Towards understanding the Lorenz curve using the Uniform distribution Chris J. Stehens Newcastle City Council, Newcastle uon Tyne, UK (For the Gini-Lorenz Conference, University of Siena, Italy, May 2005)
More informationLinear diophantine equations for discrete tomography
Journal of X-Ray Science and Technology 10 001 59 66 59 IOS Press Linear diohantine euations for discrete tomograhy Yangbo Ye a,gewang b and Jiehua Zhu a a Deartment of Mathematics, The University of Iowa,
More informationEstimating Time-Series Models
Estimating ime-series Models he Box-Jenkins methodology for tting a model to a scalar time series fx t g consists of ve stes:. Decide on the order of di erencing d that is needed to roduce a stationary
More informationCointegrated VARIMA models: specification and. simulation
Cointegrated VARIMA models: specification and simulation José L. Gallego and Carlos Díaz Universidad de Cantabria. Abstract In this note we show how specify cointegrated vector autoregressive-moving average
More informationAveraging sums of powers of integers and Faulhaber polynomials
Annales Mathematicae et Informaticae 42 (20. 0 htt://ami.ektf.hu Averaging sums of owers of integers and Faulhaber olynomials José Luis Cereceda a a Distrito Telefónica Madrid Sain jl.cereceda@movistar.es
More informationLower Confidence Bound for Process-Yield Index S pk with Autocorrelated Process Data
Quality Technology & Quantitative Management Vol. 1, No.,. 51-65, 15 QTQM IAQM 15 Lower onfidence Bound for Process-Yield Index with Autocorrelated Process Data Fu-Kwun Wang * and Yeneneh Tamirat Deartment
More informationarxiv: v2 [quant-ph] 2 Aug 2012
Qcomiler: quantum comilation with CSD method Y. G. Chen a, J. B. Wang a, a School of Physics, The University of Western Australia, Crawley WA 6009 arxiv:208.094v2 [quant-h] 2 Aug 202 Abstract In this aer,
More informationVerifying Two Conjectures on Generalized Elite Primes
1 2 3 47 6 23 11 Journal of Integer Sequences, Vol. 12 (2009), Article 09.4.7 Verifying Two Conjectures on Generalized Elite Primes Xiaoqin Li 1 Mathematics Deartment Anhui Normal University Wuhu 241000,
More informationDetection Algorithm of Particle Contamination in Reticle Images with Continuous Wavelet Transform
Detection Algorithm of Particle Contamination in Reticle Images with Continuous Wavelet Transform Chaoquan Chen and Guoing Qiu School of Comuter Science and IT Jubilee Camus, University of Nottingham Nottingham
More informationState Estimation with ARMarkov Models
Deartment of Mechanical and Aerosace Engineering Technical Reort No. 3046, October 1998. Princeton University, Princeton, NJ. State Estimation with ARMarkov Models Ryoung K. Lim 1 Columbia University,
More informationA Comparison between Biased and Unbiased Estimators in Ordinary Least Squares Regression
Journal of Modern Alied Statistical Methods Volume Issue Article 7 --03 A Comarison between Biased and Unbiased Estimators in Ordinary Least Squares Regression Ghadban Khalaf King Khalid University, Saudi
More informationMODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL
Technical Sciences and Alied Mathematics MODELING THE RELIABILITY OF CISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL Cezar VASILESCU Regional Deartment of Defense Resources Management
More information1. INTRODUCTION. Fn 2 = F j F j+1 (1.1)
CERTAIN CLASSES OF FINITE SUMS THAT INVOLVE GENERALIZED FIBONACCI AND LUCAS NUMBERS The beautiful identity R.S. Melham Deartment of Mathematical Sciences, University of Technology, Sydney PO Box 23, Broadway,
More informationRANDOM WALKS AND PERCOLATION: AN ANALYSIS OF CURRENT RESEARCH ON MODELING NATURAL PROCESSES
RANDOM WALKS AND PERCOLATION: AN ANALYSIS OF CURRENT RESEARCH ON MODELING NATURAL PROCESSES AARON ZWIEBACH Abstract. In this aer we will analyze research that has been recently done in the field of discrete
More informationConstruction of High-Girth QC-LDPC Codes
MITSUBISHI ELECTRIC RESEARCH LABORATORIES htt://wwwmerlcom Construction of High-Girth QC-LDPC Codes Yige Wang, Jonathan Yedidia, Stark Draer TR2008-061 Setember 2008 Abstract We describe a hill-climbing
More informationASYMPTOTIC RESULTS OF A HIGH DIMENSIONAL MANOVA TEST AND POWER COMPARISON WHEN THE DIMENSION IS LARGE COMPARED TO THE SAMPLE SIZE
J Jaan Statist Soc Vol 34 No 2004 9 26 ASYMPTOTIC RESULTS OF A HIGH DIMENSIONAL MANOVA TEST AND POWER COMPARISON WHEN THE DIMENSION IS LARGE COMPARED TO THE SAMPLE SIZE Yasunori Fujikoshi*, Tetsuto Himeno
More informationHENSEL S LEMMA KEITH CONRAD
HENSEL S LEMMA KEITH CONRAD 1. Introduction In the -adic integers, congruences are aroximations: for a and b in Z, a b mod n is the same as a b 1/ n. Turning information modulo one ower of into similar
More informationMATHEMATICAL MODELLING OF THE WIRELESS COMMUNICATION NETWORK
Comuter Modelling and ew Technologies, 5, Vol.9, o., 3-39 Transort and Telecommunication Institute, Lomonosov, LV-9, Riga, Latvia MATHEMATICAL MODELLIG OF THE WIRELESS COMMUICATIO ETWORK M. KOPEETSK Deartment
More information4. Score normalization technical details We now discuss the technical details of the score normalization method.
SMT SCORING SYSTEM This document describes the scoring system for the Stanford Math Tournament We begin by giving an overview of the changes to scoring and a non-technical descrition of the scoring rules
More informationThe Properties of Pure Diagonal Bilinear Models
American Journal of Mathematics and Statistics 016, 6(4): 139-144 DOI: 10.593/j.ajms.0160604.01 The roerties of ure Diagonal Bilinear Models C. O. Omekara Deartment of Statistics, Michael Okara University
More informationON THE SET a x + b g x (mod p) 1 Introduction
PORTUGALIAE MATHEMATICA Vol 59 Fasc 00 Nova Série ON THE SET a x + b g x (mod ) Cristian Cobeli, Marian Vâjâitu and Alexandru Zaharescu Abstract: Given nonzero integers a, b we rove an asymtotic result
More informationPrinciples of Computed Tomography (CT)
Page 298 Princiles of Comuted Tomograhy (CT) The theoretical foundation of CT dates back to Johann Radon, a mathematician from Vienna who derived a method in 1907 for rojecting a 2-D object along arallel
More informationA generalization of Amdahl's law and relative conditions of parallelism
A generalization of Amdahl's law and relative conditions of arallelism Author: Gianluca Argentini, New Technologies and Models, Riello Grou, Legnago (VR), Italy. E-mail: gianluca.argentini@riellogrou.com
More informationDeriving Indicator Direct and Cross Variograms from a Normal Scores Variogram Model (bigaus-full) David F. Machuca Mory and Clayton V.
Deriving ndicator Direct and Cross Variograms from a Normal Scores Variogram Model (bigaus-full) David F. Machuca Mory and Clayton V. Deutsch Centre for Comutational Geostatistics Deartment of Civil &
More informationMATH 2710: NOTES FOR ANALYSIS
MATH 270: NOTES FOR ANALYSIS The main ideas we will learn from analysis center around the idea of a limit. Limits occurs in several settings. We will start with finite limits of sequences, then cover infinite
More informationMathematics as the Language of Physics.
Mathematics as the Language of Physics. J. Dunning-Davies Deartment of Physics University of Hull Hull HU6 7RX England. Email: j.dunning-davies@hull.ac.uk Abstract. Courses in mathematical methods for
More informationChater Matrix Norms and Singular Value Decomosition Introduction In this lecture, we introduce the notion of a norm for matrices The singular value de
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A Dahleh George Verghese Deartment of Electrical Engineering and Comuter Science Massachuasetts Institute of Technology c Chater Matrix Norms
More informationAn Ant Colony Optimization Approach to the Probabilistic Traveling Salesman Problem
An Ant Colony Otimization Aroach to the Probabilistic Traveling Salesman Problem Leonora Bianchi 1, Luca Maria Gambardella 1, and Marco Dorigo 2 1 IDSIA, Strada Cantonale Galleria 2, CH-6928 Manno, Switzerland
More informationUncorrelated Multilinear Principal Component Analysis for Unsupervised Multilinear Subspace Learning
TNN-2009-P-1186.R2 1 Uncorrelated Multilinear Princial Comonent Analysis for Unsuervised Multilinear Subsace Learning Haiing Lu, K. N. Plataniotis and A. N. Venetsanooulos The Edward S. Rogers Sr. Deartment
More informationFE FORMULATIONS FOR PLASTICITY
G These slides are designed based on the book: Finite Elements in Plasticity Theory and Practice, D.R.J. Owen and E. Hinton, 1970, Pineridge Press Ltd., Swansea, UK. 1 Course Content: A INTRODUCTION AND
More informationsubstantial literature on emirical likelihood indicating that it is widely viewed as a desirable and natural aroach to statistical inference in a vari
Condence tubes for multile quantile lots via emirical likelihood John H.J. Einmahl Eindhoven University of Technology Ian W. McKeague Florida State University May 7, 998 Abstract The nonarametric emirical
More informationSCHUR S LEMMA AND BEST CONSTANTS IN WEIGHTED NORM INEQUALITIES. Gord Sinnamon The University of Western Ontario. December 27, 2003
SCHUR S LEMMA AND BEST CONSTANTS IN WEIGHTED NORM INEQUALITIES Gord Sinnamon The University of Western Ontario December 27, 23 Abstract. Strong forms of Schur s Lemma and its converse are roved for mas
More information22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 13: PF Design II The Coil Solver
.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 13: PF Design II The Coil Solver Introduction 1. Let us assume that we have successfully solved the Grad Shafranov equation for a fixed boundary
More informationPublished: 14 October 2013
Electronic Journal of Alied Statistical Analysis EJASA, Electron. J. A. Stat. Anal. htt://siba-ese.unisalento.it/index.h/ejasa/index e-issn: 27-5948 DOI: 1.1285/i275948v6n213 Estimation of Parameters of
More informationCHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules
CHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules. Introduction: The is widely used in industry to monitor the number of fraction nonconforming units. A nonconforming unit is
More informationBayesian inference & Markov chain Monte Carlo. Note 1: Many slides for this lecture were kindly provided by Paul Lewis and Mark Holder
Bayesian inference & Markov chain Monte Carlo Note 1: Many slides for this lecture were kindly rovided by Paul Lewis and Mark Holder Note 2: Paul Lewis has written nice software for demonstrating Markov
More informationA recipe for an unpredictable random number generator
Condensed Matter Physics 006, Vol. 9, No (46,. 367 37 A recie for an unredictable random number generator M.A.García-Ñustes, L.Trujillo, J.A.González Centro de Física, Instituto Venezolano de Investigaciones
More informationRound-off Errors and Computer Arithmetic - (1.2)
Round-off Errors and Comuter Arithmetic - (.). Round-off Errors: Round-off errors is roduced when a calculator or comuter is used to erform real number calculations. That is because the arithmetic erformed
More informationHotelling s Two- Sample T 2
Chater 600 Hotelling s Two- Samle T Introduction This module calculates ower for the Hotelling s two-grou, T-squared (T) test statistic. Hotelling s T is an extension of the univariate two-samle t-test
More informationON THE LEAST SIGNIFICANT p ADIC DIGITS OF CERTAIN LUCAS NUMBERS
#A13 INTEGERS 14 (014) ON THE LEAST SIGNIFICANT ADIC DIGITS OF CERTAIN LUCAS NUMBERS Tamás Lengyel Deartment of Mathematics, Occidental College, Los Angeles, California lengyel@oxy.edu Received: 6/13/13,
More information2. Sample representativeness. That means some type of probability/random sampling.
1 Neuendorf Cluster Analysis Assumes: 1. Actually, any level of measurement (nominal, ordinal, interval/ratio) is accetable for certain tyes of clustering. The tyical methods, though, require metric (I/R)
More informationLilian Markenzon 1, Nair Maria Maia de Abreu 2* and Luciana Lee 3
Pesquisa Oeracional (2013) 33(1): 123-132 2013 Brazilian Oerations Research Society Printed version ISSN 0101-7438 / Online version ISSN 1678-5142 www.scielo.br/oe SOME RESULTS ABOUT THE CONNECTIVITY OF
More informationAR PROCESSES AND SOURCES CAN BE RECONSTRUCTED FROM. Radu Balan, Alexander Jourjine, Justinian Rosca. Siemens Corporation Research
AR PROCESSES AND SOURCES CAN BE RECONSTRUCTED FROM DEGENERATE MIXTURES Radu Balan, Alexander Jourjine, Justinian Rosca Siemens Cororation Research 7 College Road East Princeton, NJ 8 fradu,jourjine,roscag@scr.siemens.com
More informationNode-voltage method using virtual current sources technique for special cases
Node-oltage method using irtual current sources technique for secial cases George E. Chatzarakis and Marina D. Tortoreli Electrical and Electronics Engineering Deartments, School of Pedagogical and Technological
More informationarxiv: v2 [stat.me] 3 Nov 2014
onarametric Stein-tye Shrinkage Covariance Matrix Estimators in High-Dimensional Settings Anestis Touloumis Cancer Research UK Cambridge Institute University of Cambridge Cambridge CB2 0RE, U.K. Anestis.Touloumis@cruk.cam.ac.uk
More informationChapter 10. Supplemental Text Material
Chater 1. Sulemental Tet Material S1-1. The Covariance Matri of the Regression Coefficients In Section 1-3 of the tetbook, we show that the least squares estimator of β in the linear regression model y=
More informationElliptic Curves and Cryptography
Ellitic Curves and Crytograhy Background in Ellitic Curves We'll now turn to the fascinating theory of ellitic curves. For simlicity, we'll restrict our discussion to ellitic curves over Z, where is a
More information8 STOCHASTIC PROCESSES
8 STOCHASTIC PROCESSES The word stochastic is derived from the Greek στoχαστικoς, meaning to aim at a target. Stochastic rocesses involve state which changes in a random way. A Markov rocess is a articular
More informationSupplement to the paper Accurate distributions of Mallows C p and its unbiased modifications with applications to shrinkage estimation
To aear in Economic Review (Otaru University of Commerce Vol. No..- 017. Sulement to the aer Accurate distributions of Mallows C its unbiased modifications with alications to shrinkage estimation Haruhiko
More informationCMSC 425: Lecture 4 Geometry and Geometric Programming
CMSC 425: Lecture 4 Geometry and Geometric Programming Geometry for Game Programming and Grahics: For the next few lectures, we will discuss some of the basic elements of geometry. There are many areas
More information216 S. Chandrasearan and I.C.F. Isen Our results dier from those of Sun [14] in two asects: we assume that comuted eigenvalues or singular values are
Numer. Math. 68: 215{223 (1994) Numerische Mathemati c Sringer-Verlag 1994 Electronic Edition Bacward errors for eigenvalue and singular value decomositions? S. Chandrasearan??, I.C.F. Isen??? Deartment
More informationOptimization of Gear Design and Manufacture. Vilmos SIMON *
7 International Conference on Mechanical and Mechatronics Engineering (ICMME 7) ISBN: 978--6595-44- timization of Gear Design and Manufacture Vilmos SIMN * Budaest Universit of Technolog and Economics,
More informationOn the Maximum and Minimum of Multivariate Pareto Random Variables
ISSN 1392 124X (rint) ISSN 2335 884X (online) INFORMATION TECHNOLOGY AND CONTROL 213 Vol.42 No.3 On the Maximum Minimum of Multivariate Pareto Rom Variables Jungsoo Woo Deartment of Statistics Yeungnam
More informationElementary Analysis in Q p
Elementary Analysis in Q Hannah Hutter, May Szedlák, Phili Wirth November 17, 2011 This reort follows very closely the book of Svetlana Katok 1. 1 Sequences and Series In this section we will see some
More informationLower bound solutions for bearing capacity of jointed rock
Comuters and Geotechnics 31 (2004) 23 36 www.elsevier.com/locate/comgeo Lower bound solutions for bearing caacity of jointed rock D.J. Sutcliffe a, H.S. Yu b, *, S.W. Sloan c a Deartment of Civil, Surveying
More informationPerformance of lag length selection criteria in three different situations
MPRA Munich Personal RePEc Archive Performance of lag length selection criteria in three different situations Zahid Asghar and Irum Abid Quaid-i-Azam University, Islamabad Aril 2007 Online at htts://mra.ub.uni-muenchen.de/40042/
More informationAN EXTENDED YULE-WALKER METHOD FOR ESTIMATING A VECTOR AUTOREGRESSIVE MODEL WITH MIXED-FREQUENCY DATA *
AN EXTENDED YULE-WALKER METHOD FOR ESTIMATING A VECTOR AUTOREGRESSIVE MODEL WITH MIXED-FREQUENCY DATA * Baoline Chen Deartment of Economics Rutgers University-Camden Camden, NJ 0802 e-mail:baoline@crabrutgersedu
More information2-D Analysis for Iterative Learning Controller for Discrete-Time Systems With Variable Initial Conditions Yong FANG 1, and Tommy W. S.
-D Analysis for Iterative Learning Controller for Discrete-ime Systems With Variable Initial Conditions Yong FANG, and ommy W. S. Chow Abstract In this aer, an iterative learning controller alying to linear
More informationSpectral Analysis by Stationary Time Series Modeling
Chater 6 Sectral Analysis by Stationary Time Series Modeling Choosing a arametric model among all the existing models is by itself a difficult roblem. Generally, this is a riori information about the signal
More information2. Multivariate ARMA
2. Multivariate ARMA JEM 140: Quantitative Multivariate Finance IES, Charles University, Prague Summer 2018 JEM 140 () 2. Multivariate ARMA Summer 2018 1 / 19 Multivariate AR I Let r t = (r 1t,..., r kt
More informationThe power performance of fixed-t panel unit root tests allowing for structural breaks in their deterministic components
ATHES UIVERSITY OF ECOOMICS AD BUSIESS DEPARTMET OF ECOOMICS WORKIG PAPER SERIES 23-203 The ower erformance of fixed-t anel unit root tests allowing for structural breaks in their deterministic comonents
More informationProbability Estimates for Multi-class Classification by Pairwise Coupling
Probability Estimates for Multi-class Classification by Pairwise Couling Ting-Fan Wu Chih-Jen Lin Deartment of Comuter Science National Taiwan University Taiei 06, Taiwan Ruby C. Weng Deartment of Statistics
More informationIntroduction to Probability and Statistics
Introduction to Probability and Statistics Chater 8 Ammar M. Sarhan, asarhan@mathstat.dal.ca Deartment of Mathematics and Statistics, Dalhousie University Fall Semester 28 Chater 8 Tests of Hyotheses Based
More informationChapter 6. Phillip Hall - Room 537, Huxley
Chater 6 6 Partial Derivatives.................................................... 72 6. Higher order artial derivatives...................................... 73 6.2 Matrix of artial derivatives.........................................74
More informationRESOLUTIONS OF THREE-ROWED SKEW- AND ALMOST SKEW-SHAPES IN CHARACTERISTIC ZERO
RESOLUTIONS OF THREE-ROWED SKEW- AND ALMOST SKEW-SHAPES IN CHARACTERISTIC ZERO MARIA ARTALE AND DAVID A. BUCHSBAUM Abstract. We find an exlicit descrition of the terms and boundary mas for the three-rowed
More informationDistributed Rule-Based Inference in the Presence of Redundant Information
istribution Statement : roved for ublic release; distribution is unlimited. istributed Rule-ased Inference in the Presence of Redundant Information June 8, 004 William J. Farrell III Lockheed Martin dvanced
More informationAnalysis of some entrance probabilities for killed birth-death processes
Analysis of some entrance robabilities for killed birth-death rocesses Master s Thesis O.J.G. van der Velde Suervisor: Dr. F.M. Sieksma July 5, 207 Mathematical Institute, Leiden University Contents Introduction
More informationCalculation of MTTF values with Markov Models for Safety Instrumented Systems
7th WEA International Conference on APPLIE COMPUTE CIENCE, Venice, Italy, November -3, 7 3 Calculation of MTTF values with Markov Models for afety Instrumented ystems BÖCÖK J., UGLJEA E., MACHMU. University
More informationMulti-Operation Multi-Machine Scheduling
Multi-Oeration Multi-Machine Scheduling Weizhen Mao he College of William and Mary, Williamsburg VA 3185, USA Abstract. In the multi-oeration scheduling that arises in industrial engineering, each job
More informationOn Wald-Type Optimal Stopping for Brownian Motion
J Al Probab Vol 34, No 1, 1997, (66-73) Prerint Ser No 1, 1994, Math Inst Aarhus On Wald-Tye Otimal Stoing for Brownian Motion S RAVRSN and PSKIR The solution is resented to all otimal stoing roblems of
More informationCombining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO)
Combining Logistic Regression with Kriging for Maing the Risk of Occurrence of Unexloded Ordnance (UXO) H. Saito (), P. Goovaerts (), S. A. McKenna (2) Environmental and Water Resources Engineering, Deartment
More informationMULTIVARIATE SHEWHART QUALITY CONTROL FOR STANDARD DEVIATION
MULTIVARIATE SHEWHART QUALITY CONTROL FOR STANDARD DEVIATION M. Jabbari Nooghabi, Deartment of Statistics, Faculty of Mathematical Sciences, Ferdowsi University of Mashhad, Mashhad-Iran. and H. Jabbari
More informationEvaluating Circuit Reliability Under Probabilistic Gate-Level Fault Models
Evaluating Circuit Reliability Under Probabilistic Gate-Level Fault Models Ketan N. Patel, Igor L. Markov and John P. Hayes University of Michigan, Ann Arbor 48109-2122 {knatel,imarkov,jhayes}@eecs.umich.edu
More informationControllability and Resiliency Analysis in Heat Exchanger Networks
609 A ublication of CHEMICAL ENGINEERING RANSACIONS VOL. 6, 07 Guest Editors: Petar S Varbanov, Rongxin Su, Hon Loong Lam, Xia Liu, Jiří J Klemeš Coyright 07, AIDIC Servizi S.r.l. ISBN 978-88-95608-5-8;
More informationRobust Performance Design of PID Controllers with Inverse Multiplicative Uncertainty
American Control Conference on O'Farrell Street San Francisco CA USA June 9 - July Robust Performance Design of PID Controllers with Inverse Multilicative Uncertainty Tooran Emami John M Watkins Senior
More informationTests for Two Proportions in a Stratified Design (Cochran/Mantel-Haenszel Test)
Chater 225 Tests for Two Proortions in a Stratified Design (Cochran/Mantel-Haenszel Test) Introduction In a stratified design, the subects are selected from two or more strata which are formed from imortant
More informationON THE DEVELOPMENT OF PARAMETER-ROBUST PRECONDITIONERS AND COMMUTATOR ARGUMENTS FOR SOLVING STOKES CONTROL PROBLEMS
Electronic Transactions on Numerical Analysis. Volume 44,. 53 72, 25. Coyright c 25,. ISSN 68 963. ETNA ON THE DEVELOPMENT OF PARAMETER-ROBUST PRECONDITIONERS AND COMMUTATOR ARGUMENTS FOR SOLVING STOKES
More informationTRACES OF SCHUR AND KRONECKER PRODUCTS FOR BLOCK MATRICES
Khayyam J. Math. DOI:10.22034/kjm.2019.84207 TRACES OF SCHUR AND KRONECKER PRODUCTS FOR BLOCK MATRICES ISMAEL GARCÍA-BAYONA Communicated by A.M. Peralta Abstract. In this aer, we define two new Schur and
More informationModel checking, verification of CTL. One must verify or expel... doubts, and convert them into the certainty of YES [Thomas Carlyle]
Chater 5 Model checking, verification of CTL One must verify or exel... doubts, and convert them into the certainty of YES or NO. [Thomas Carlyle] 5. The verification setting Page 66 We introduce linear
More informationMULTIVARIATE STATISTICAL PROCESS OF HOTELLING S T CONTROL CHARTS PROCEDURES WITH INDUSTRIAL APPLICATION
Journal of Statistics: Advances in heory and Alications Volume 8, Number, 07, Pages -44 Available at htt://scientificadvances.co.in DOI: htt://dx.doi.org/0.864/jsata_700868 MULIVARIAE SAISICAL PROCESS
More informationThe Recursive Fitting of Multivariate. Complex Subset ARX Models
lied Mathematical Sciences, Vol. 1, 2007, no. 23, 1129-1143 The Recursive Fitting of Multivariate Comlex Subset RX Models Jack Penm School of Finance and lied Statistics NU College of Business & conomics
More informationDimensional perturbation theory for Regge poles
Dimensional erturbation theory for Regge oles Timothy C. Germann Deartment of Chemistry, University of California, Berkeley, California 94720 Sabre Kais Deartment of Chemistry, Purdue University, West
More information19th Bay Area Mathematical Olympiad. Problems and Solutions. February 28, 2017
th Bay Area Mathematical Olymiad February, 07 Problems and Solutions BAMO- and BAMO- are each 5-question essay-roof exams, for middle- and high-school students, resectively. The roblems in each exam are
More informationStochastic processes: basic notions
Stochastic rocesses: basic notions Jean-Marie Dufour McGill University First version: March 2002 Revised: Setember 2002, Aril 2004, Setember 2004, January 2005, July 2011, May 2016, July 2016 This version:
More informationFuzzy Methods. Additions to Chapter 5: Fuzzy Arithmetic. Michael Hanss.
Fuzzy Methods Additions to Chater 5: Fuzzy Arithmetic Michael Hanss Part I: A short review of the Institute of Engineering Comutational Mechanics University of Stuttgart Germany Examle : q = f( ) = 2 2
More informationBooth School of Business, University of Chicago Business 41914, Spring Quarter 2013, Mr. Ruey S. Tsay. Midterm
Booth School of Business, University of Chicago Business 41914, Spring Quarter 2013, Mr. Ruey S. Tsay Midterm Chicago Booth Honor Code: I pledge my honor that I have not violated the Honor Code during
More informationEstimating function analysis for a class of Tweedie regression models
Title Estimating function analysis for a class of Tweedie regression models Author Wagner Hugo Bonat Deartamento de Estatística - DEST, Laboratório de Estatística e Geoinformação - LEG, Universidade Federal
More informationConvex Optimization methods for Computing Channel Capacity
Convex Otimization methods for Comuting Channel Caacity Abhishek Sinha Laboratory for Information and Decision Systems (LIDS), MIT sinhaa@mit.edu May 15, 2014 We consider a classical comutational roblem
More informationOn a Markov Game with Incomplete Information
On a Markov Game with Incomlete Information Johannes Hörner, Dinah Rosenberg y, Eilon Solan z and Nicolas Vieille x{ January 24, 26 Abstract We consider an examle of a Markov game with lack of information
More information