Partial Directed Coherence: Some Estimation Issues Baccalá LA 1, Sameshima K 2
|
|
- Laurel Spencer
- 5 years ago
- Views:
Transcription
1 Partial Directed Coherence: Some Estimation Issues Baccalá LA 1, Sameshima K 2 1 Telecommunications and Control Eng. Dept., Escola Politécnica and 2 Discipline of Medical Informatics & Functional Neurosurgery Laboratory, School of Medicine University of São Paulo baccala@lcs.poli.usp.br ksameshi@usp.br Summary This paper addresses some estimation performance issues of the Partial Directed Coherence (PDC), a newly introduced quantity related to Granger Causality, that is applicable to the functional and structural inference in neurobiology by processing many simultaneously measured signals. Exhaustive Monte Carlo simulations of simple models of interest show the asymptotic validity of previously used frequency independent thresholds in rejecting the hypothesis of lack of unidirectional connection between structures. When only short stationary time-series are available, we show that frequency dependent thresholds for PDC are more appropriate. Finally our results suggest that PDC converges faster to the correct structural inferences compared to classical Granger Causality tests. Introduction Recently, we proposed a new multiple-time-series based method to investigate structural relationships between simultaneous measurements taken from multiple electrodes (Baccalá and Sameshima, 21). The new method was named Partial Directed Coherence because, being an evolution of the well known frequency domain concept of partial coherence, it also permits determining the direction of information flow between each member of a pair of neural structures by subtracting the interactions and possible common influences due to other remaining yet simultaneously observed time series. The ability of PDC to disclose the direction of information flow is intimately related to the fundamental notion of Granger causality (Granger, 1969) which, roughly, refers to how much the predictability of the present of a time-series is enhanced by the exclusive knowledge of the past of another time series (Baccalá and Sameshima, 21). As a matter of fact, zero PDC from structure A to structure B at all frequencies is synonimous with lack of Granger causality from A to B. In turn this allows one to infer whether some possibly existing physical link from structure A to B is inactive during a given span of signal observation. As such, the value of the functional connectivity information PDC can provide is self evident in that it offers a glimpse on active signal pathways within the brain and deducible from the simultaneous measurement of multiple sensors. While a number of alternatives exist to gauge Granger causality (Lütkepohl, 1993, Baccalá et al., 1998), a number of open questions concerning reliable PDC inference remain. So much so that, as a conservative measure, to test the null
2 hypothesis of zero PDC in our previously published illustrative examples (Baccalá and Sameshima, 21,21b), we used the arbitrary threshold of.1 beyond which we deemed PDC significantly non-zero. This value was adopted because PDC reduces to the so called Directed Coherence proposed by Saito and Harashima (1983) when just two time series (normalized to unit power) are simultaneously analyzed; the value.1 itself was derived from extensive simulations reported by Schnider et al. (1989) for pairs of simultaneously analyzed time series. (Note that Directed Coherence was later generalized to dealing with more than two simultaneously processed time series (Baccalá et al., 1998; Baccalá and Sameshima, 1998) and is closely related to the estimator proposed by Kaminski and Blinowska (1991) known as Directed Transfer Function (see Baccalá and Sameshima, 21, for further comparisons)). As such, by contrast with Granger Causality tests (GCT) whose large sample behaviour is well known (Lütkepohl, 1993), more research is necessary to assess PDC s statistical effectiveness and this paper represents a step towards that goal. Here, rather than pursue a detailed mathematical investigation of the asymptotic statistics, we chose to investigate and illustrate PDC s behaviour via extensive Monte Carlo simulations of some scenarios of interest from a very simple model comprising just a pair of time series to a more complicated scenario involving more structures and varying complexities in dynamical behaviour. With neuroscience applications in mind, performing Monte Carlo studies is justifiable because many functional connectivity issues among brain structures, as in behaviour related studies, take place in small sample contexts, where large sample statistics are just a rough guide. Background and Methods In this paper, three scenarios of interest were simulated; their data were obtained by from multiple time series derived from autoregressive models. The first example covered just two time series with unidirectional connection only from x 1 (k) to x 2 (k) described by: x 1 (k) x 2 (k) =.75.5 x 1 (k 1) x 2 (k 1) + w 1(k) w 2 (k) (Model A) where we allowed the dynamics of x 2 (k) to vary through the parameter. Note that lack of Granger causality from x 2 (k) to x 1 (k) came from the zero entry in the matrix. The w i (k) are uncorrelated zero mean gaussian innovations processes. To measure the effect of a more complex time structure in the series we repeated PDC estimation after pre-processing Model A signals through identical moving average filters: y i (k) = 1, (Model A ) 2 (x i(k) +x i (k 1)) i = 1, 2 Remark: Though, as one can easily show, the theoretical values of PDC remain unaffected in going from Model A to A, estimability changes severely. Finally, we consider a scenario, borrowed from Baccalá and Sameshima (21b) involving a signal source and a signal sink in which information travels via
3 multiple pathways. The network comprises N = 7 structures and has two disconnected portions that are the object of our investigation here through x 1 (k)=.95 2 x 1 (k 1).925x 1 (k 2)+w 1 (k) x 2 (k)=.5x 1 (k 1)+w 2 (k) x 3 (k)=+.4x 1 (k 4).4x 2 (k 2)+w 3 (k) x 4 (k)=.5x 3 (k 1) x 4 (k 1) x 5 (k 1)+w 4 (k) (Model B) x 5 (k)=.25 2 x 4 (k 1) x 5 (k 1)+w 5 (k) x 6 (k)=.95 2 x 6 (k 1).925x 6 (k 2)+w 6 (k) x 7 (k)=.1x 6 (k 2)+w 7 (k) (see also Baccalá and Sameshima (21b) for a directed graph connectivity structure and a figure of the overall theoretical PDC). Here we concentrated on the interaction behaviour between x 1 (k) and x 6 (k) that represent identical but unlinked oscillators. Each model was simulated using random white zero mean unit variance gaussian uncorrelated innovations w i (k) and various signal lengths n. Simulations were repeated N r times to obtain the sample behaviour of PDC (from each x j (k) to each x i (k)) as defined in Baccalá and Sameshima (21) via ij (f) = A ij (f) N A & kj (f)a kj (f) i=1 with p 1 aij (r)e j2 fr, if i = j A ij (f) = r=1 p aij (r)e j2 fr, otherwise r=1 which was computed through the Nutall-Strand algorithm (Marple Jr.,1985) fit of the model: x 1 (k) x N (k) p = r=1 a 11 (r) a 12 (r) a 1N (r) a ij (r) a N1 (r) a NN (r) x 1 (k r) x N (k r) + w 1 (k) w N (k). Model orders p were estimated using Akaike s Information criterion (AIC), except for Model A where other criteria (Schwartz s (SC) and Hannan-Quinn s (HQ)) were also investigated (Lütkepohl, 1993). Model order was limited to p = 2 for practical computational reasons. In all cases, the corresponding Granger Causality tests (GCT) (see Baccalá et al (1998) for details) were performed at 5% confidence levels as comparison benchmarks. We also used the conservative PDC based criterion: ij (f) 2 >.1
4 (called the Spectral Causality Criterion (SCC)) to look for the consistency of its previous use as evidence of existing connectivity. Results The average behaviour of N r = 256 PDC estimates from x 2 to x 1 as a function of sample size n is depicted in Fig. 1 for Model A ( =.75). For smaller sample sizes the graph reveals the frequency dependence of PDC under the null hypothesis of no connection. In comparing GCT to SCC, the rough independence of GCT on sample size becomes evident (Fig. 2) and contrasts with SCC whose number of false positives decreases as n grows. Obviously this failure of SCC is a by-product of the inadequacy of a frequency independent threshold for small samples sizes. Note also that very few instances of false positive SCC inferences were observed for sufficiently large sample sizes. 1 π 12 (f) n=32 n=128 n= f Fig. 1 - Average behaviour of the estimated squared PDC, which is theoretically zero, as a function of the number of samples for the non-existing connection in Model A ( =.75). Note the explicit dependence on frequency (represented in normalized scale from to.5 of the Nyquist sampling frequency), Percentage of False Positives SCC Samples GCT Fig. 2 - Comparison as a function of signal length n between a canonical Granger Causality test (GCT) and SCC, i. e. assuming PDC connectivity detection based on a frequency independent decision threshold for Model A with. =.75
5 To appreciate the effect of variation in the temporal dynamics, we changed in Model A and computed the corresponding percentage of false positives for n = 256 and N r = 256. GCT results were fairly uniform whereas SCC depended more strongly on the underlying signal dynamics (Fig. 3) having its effectiveness decreased for larger values. Percentage of False Positives SCC GCT θ Fig. 3 - The effect on GCT and SSC of variations in Model A s dynamics through changes in (n=256 and Nr=256). Only SCC was affected by not taking account of PDC s small sample dependence on frequency. For Model A (Fig. 4) comparable PDC performance was only achieved for larger sample sizes. Also, GCT computation produced unacceptably high ratios of false positives that increased with sample size. Of the model order criteria, HQ provided the best PDC performance. Note the intriguing reversal in model order criterion performance in Fig. 4 depending on whether PDC or GCT was used for causalilty detection. ( N r = 64). Percentage of False Positives HQ SC HQ SC AIC AIC Samples Fig. 4 - Impact of different model order estimators on GCT (dashed) and SCC (solid) as connectivity estimators as a function of signal length in samples for Model A. For the next scenario, our focus was the behaviour of possible causality inference errors between the structurally disconnected oscillators represented by x 1 (k) and x 6 (k) in Model B. It is easy to see (Fig. 5) that on average PDC estimates
6 are larger around that common frequency. As for the previous models, SCC s percentages of false positives decrease monotonically and eventually fall bellow GCT s percentages for a sample size of n = 512 points. By contrast GCT s percentages, as for Model A are fairly constant throughout the sample length range (Fig. 5), yet remain well above 5% ( N r = 256)..9 Percentage of False Positives GCT SCC Samples Fig. 5 - Comparison between GCT and SCC in falsely detecting connections from x 6 (k) to x 1 (k) as a function of signal duration in samples for Model B based on 256 independent simulations 5 π 16 (f) Increasing sample size f x 6 (k) x 1 (k) Fig. 6 - Average squared PDC from to shows that most of the partial directed coherence inference errors reflected in the SCC test of Fig. 5 come from the frequency common to both oscillators (N r = 256). The frequency scale is normalized from to.5 of the Nyquist sampling frequency. Discussion The first emerging pattern is that PDC and its frequency independent threshold based criterion (SCC) are good estimators of the unidirectional connectivity for a sufficiently long signal sample. For small samples sizes, PDC dependends on frequency; for Model A, under the assumption of an estimated correct p = 1 model order, it is easy to explicitly quantify this dependence by computing
7 12 (f) 2 = 2, Re{ } cos(2 f) where { results from imprecise estimates of a 12 (1) due to the finite sample duration. This expression also predicts the dependence of PDC errors due to changes in the inherent time series dynamics (Fig. 3) as the denominator approaches zero for large. In fact if is negative, PDC increases for high frequencies (as in Fig. 1); conversely PDC estimates are increased for small frequencies if is positive. From this, one also concludes that a sufficiently long sample for correct SCC-based connectivity detection is a concept that depends on the inherent temporal dynamics of the series under study via parameters like. When sufficiently long signals were used, negligible error rates in directional connectivity inference resulted from the appplication of SCC. The same level of performance was not witnessed for GCT. In the computations for Model A, at least 2 to 4 times as many samples (Fig. 4) were necessary to achieve the same rate of correct conclusions as for Model A which has the same underlying directional connectivity structure as Model A. The common moving average filter applied to both series complicated the dynamics so that it no longer conformed to a low order autoregressive dynamic structure. In this case, the best model order estimator to produce the most correct connectivity inferences based on SCC with fewer time samples was Hannan-Quinn s. The exactly reverse (and as whole generally unsatisfactory) effect was observed for GCT and merits further investigation. The intent behind Model B s illustration was to show that though many more N = 7 time series were analysed simultaneously, the effect on causality inference was less pronounced than for Model A. This leads to the conjecture that reliable PDC-based connectivity inference from short time samples through multivariate autoregressions is possible only in so far as the signals are parsimoniously represented by MAR models, being less affected by issues of how many series are simultaneously processed. Though beyond the scope of this paper, our investigations show that bootstraping can be used to introduce frequency dependent thresholds to PDC connectivity estimation that further decrease the number of time samples necessary for correct inference. Conclusion Partial Directed Coherence affords better means of correct direct unidirectional connectivity inference when compared to a more traditional Granger Causality Test. PDC alse achieves negligibly few incorrect connectivity decisions if the sample size is sufficiently long even if a frequency independent threshold is used. What constitutes a sufficiently long sample depends on the temporal dynamics of the signals analysed. Supported by Grants from FAPESP (99/ to LAB, 96/ to KS) and CNPq (31273/96- to LAB).
8 References Baccalá LA, Sameshima K, Ballester G, Valle AC, Timo-Iaria C: Studying the interaction between brain structures via directed coherence and Granger causality. App. Sig Process. 1998; 5:4-48. Baccalá LA, Sameshima K: Directed Coherence: a tool for exploring functional interactions among brain Structures. In Methods for Neuronal Ensemble Recordings. Nicolelis MAL (ed) CRC Press, Boca Raton pp Baccalá LA, Sameshima K: Partial directed coherence: a new concept in neural structure determination. Biol. Cybern. 21; 84: Baccalá LA, Sameshima K: Overcoming the limitations of correlation analysis for many simultaneously processed neural structures. In Progress in Brain Research, vol. 13, Nicolelis MAL (ed). 21b. Elsevier Science. (in press). Granger CWJ: Investigating causal relations by econometric models and cross-spectral methods. Econometrica. 1969; 37: Kaminski MJ, Blinowska KJ: A new method of the description of the information flow in the brain structures. Biol. Cybern. 1991; 65: Lütkepohl, H: Introduction to Multiple Time Series Analysis. 2nd ed Springer-Verlag, Berlin. Marple Jr SL: Digital Spectral Analysis Prentice Hall, Englewood Cliffs. Saito Y and Harashima H: Tracking of information within multichannel EEG record - causal analysis in EEG. In N. Yamaguchi and K. Fujisawa (Eds.), Recent Advances in EEG and EMG Data Processing Elsevier, Amsterdam, pp Sameshima K, Baccalá LA: Using partial directed coherence to describe neuronal ensemble interactions. J. Neurosci. Methods, 1999; 94: Schnider SM, Kwong RH, Lenz FA, Kwan HC: Detection of feedback in the central nervous system using system identification techniques. Biol. Cybern. 1989; 6:
PERFORMANCE STUDY OF CAUSALITY MEASURES
PERFORMANCE STUDY OF CAUSALITY MEASURES T. Bořil, P. Sovka Department of Circuit Theory, Faculty of Electrical Engineering, Czech Technical University in Prague Abstract Analysis of dynamic relations in
More information9 Graphical modelling of dynamic relationships in multivariate time series
9 Graphical modelling of dynamic relationships in multivariate time series Michael Eichler Institut für Angewandte Mathematik Universität Heidelberg Germany SUMMARY The identification and analysis of interactions
More informationCausal Influence: Advances in Neurosignal Analysis
Critical Reviews in Biomedical Engineering, 33(4):347 430 (2005) Causal Influence: Advances in Neurosignal Analysis Maciej Kamiński¹ & Hualou Liang² 1 Department of Biomedical Physics, Institute of Experimental
More informationSpectral Interdependency Methods
Spectral Interdependency Methods Mukesh Dhamala* Department of Physics and Astronomy, Neuroscience Institute, Georgia State University, Atlanta, Georgia, USA Synonyms Coherence and Granger causality spectral
More informationVolume 30, Issue 1. The relationship between the F-test and the Schwarz criterion: Implications for Granger-causality tests
Volume 30, Issue 1 The relationship between the F-test and the Schwarz criterion: Implications for Granger-causality tests Erdal Atukeren ETH Zurich - KOF Swiss Economic Institute Abstract In applied research,
More informationThe Role of "Leads" in the Dynamic Title of Cointegrating Regression Models. Author(s) Hayakawa, Kazuhiko; Kurozumi, Eiji
he Role of "Leads" in the Dynamic itle of Cointegrating Regression Models Author(s) Hayakawa, Kazuhiko; Kurozumi, Eiji Citation Issue 2006-12 Date ype echnical Report ext Version publisher URL http://hdl.handle.net/10086/13599
More informationOn the evaluation of information flow in multivariate systems by the directed transfer function
Biol Cybern (6) 94: 469 48 DOI 1.17/s4-6-6-z ORIGINAL PAPER Michael Eichler On the evaluation of information flow in multivariate systems by the directed transfer function Received: 18 October 5 / Accepted
More informationA TIME SERIES PARADOX: UNIT ROOT TESTS PERFORM POORLY WHEN DATA ARE COINTEGRATED
A TIME SERIES PARADOX: UNIT ROOT TESTS PERFORM POORLY WHEN DATA ARE COINTEGRATED by W. Robert Reed Department of Economics and Finance University of Canterbury, New Zealand Email: bob.reed@canterbury.ac.nz
More informationEconometría 2: Análisis de series de Tiempo
Econometría 2: Análisis de series de Tiempo Karoll GOMEZ kgomezp@unal.edu.co http://karollgomez.wordpress.com Segundo semestre 2016 IX. Vector Time Series Models VARMA Models A. 1. Motivation: The vector
More informationTopic 4 Unit Roots. Gerald P. Dwyer. February Clemson University
Topic 4 Unit Roots Gerald P. Dwyer Clemson University February 2016 Outline 1 Unit Roots Introduction Trend and Difference Stationary Autocorrelations of Series That Have Deterministic or Stochastic Trends
More informationBrief Sketch of Solutions: Tutorial 3. 3) unit root tests
Brief Sketch of Solutions: Tutorial 3 3) unit root tests.5.4.4.3.3.2.2.1.1.. -.1 -.1 -.2 -.2 -.3 -.3 -.4 -.4 21 22 23 24 25 26 -.5 21 22 23 24 25 26.8.2.4. -.4 - -.8 - - -.12 21 22 23 24 25 26 -.2 21 22
More informationModel selection criteria Λ
Model selection criteria Λ Jean-Marie Dufour y Université de Montréal First version: March 1991 Revised: July 1998 This version: April 7, 2002 Compiled: April 7, 2002, 4:10pm Λ This work was supported
More informationA Multivariate Time-Frequency Based Phase Synchrony Measure for Quantifying Functional Connectivity in the Brain
A Multivariate Time-Frequency Based Phase Synchrony Measure for Quantifying Functional Connectivity in the Brain Dr. Ali Yener Mutlu Department of Electrical and Electronics Engineering, Izmir Katip Celebi
More informationBootstrapping the Grainger Causality Test With Integrated Data
Bootstrapping the Grainger Causality Test With Integrated Data Richard Ti n University of Reading July 26, 2006 Abstract A Monte-carlo experiment is conducted to investigate the small sample performance
More informationQuality Control Using Inferential Statistics In Weibull Based Reliability Analyses S. F. Duffy 1 and A. Parikh 2
Quality Control Using Inferential Statistics In Weibull Based Reliability Analyses S. F. Duffy 1 and A. Parikh 2 1 Cleveland State University 2 N & R Engineering www.inl.gov ASTM Symposium on Graphite
More informationFinancial Econometrics
Financial Econometrics Multivariate Time Series Analysis: VAR Gerald P. Dwyer Trinity College, Dublin January 2013 GPD (TCD) VAR 01/13 1 / 25 Structural equations Suppose have simultaneous system for supply
More information1 Introduction. 2 AIC versus SBIC. Erik Swanson Cori Saviano Li Zha Final Project
Erik Swanson Cori Saviano Li Zha Final Project 1 Introduction In analyzing time series data, we are posed with the question of how past events influences the current situation. In order to determine this,
More informationInference in VARs with Conditional Heteroskedasticity of Unknown Form
Inference in VARs with Conditional Heteroskedasticity of Unknown Form Ralf Brüggemann a Carsten Jentsch b Carsten Trenkler c University of Konstanz University of Mannheim University of Mannheim IAB Nuremberg
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 informationSelection Criteria Based on Monte Carlo Simulation and Cross Validation in Mixed Models
Selection Criteria Based on Monte Carlo Simulation and Cross Validation in Mixed Models Junfeng Shang Bowling Green State University, USA Abstract In the mixed modeling framework, Monte Carlo simulation
More informationAR-order estimation by testing sets using the Modified Information Criterion
AR-order estimation by testing sets using the Modified Information Criterion Rudy Moddemeijer 14th March 2006 Abstract The Modified Information Criterion (MIC) is an Akaike-like criterion which allows
More informationTesting composite hypotheses applied to AR-model order estimation; the Akaike-criterion revised
MODDEMEIJER: TESTING COMPOSITE HYPOTHESES; THE AKAIKE-CRITERION REVISED 1 Testing composite hypotheses applied to AR-model order estimation; the Akaike-criterion revised Rudy Moddemeijer Abstract Akaike
More informationEEG- Signal Processing
Fatemeh Hadaeghi EEG- Signal Processing Lecture Notes for BSP, Chapter 5 Master Program Data Engineering 1 5 Introduction The complex patterns of neural activity, both in presence and absence of external
More informationVector Autoregressive Model. Vector Autoregressions II. Estimation of Vector Autoregressions II. Estimation of Vector Autoregressions I.
Vector Autoregressive Model Vector Autoregressions II Empirical Macroeconomics - Lect 2 Dr. Ana Beatriz Galvao Queen Mary University of London January 2012 A VAR(p) model of the m 1 vector of time series
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 informationOld and New Methods in the Age of Big Data
Old and New Methods in the Age of Big Data Zoltán Somogyvári1,2 1Department of Theory Wigner Research Center for Physics of the Hungarian Academy of Sciences 2National Institute for Clinical Neuroscience
More informationLATVIAN GDP: TIME SERIES FORECASTING USING VECTOR AUTO REGRESSION
LATVIAN GDP: TIME SERIES FORECASTING USING VECTOR AUTO REGRESSION BEZRUCKO Aleksandrs, (LV) Abstract: The target goal of this work is to develop a methodology of forecasting Latvian GDP using ARMA (AutoRegressive-Moving-Average)
More informationUNIFORMLY MOST POWERFUL CYCLIC PERMUTATION INVARIANT DETECTION FOR DISCRETE-TIME SIGNALS
UNIFORMLY MOST POWERFUL CYCLIC PERMUTATION INVARIANT DETECTION FOR DISCRETE-TIME SIGNALS F. C. Nicolls and G. de Jager Department of Electrical Engineering, University of Cape Town Rondebosch 77, South
More informationEvaluating directed connectivity of BOLD signals applying the multivariate
Evaluating directed connectivity of BOLD signals applying the multivariate autoregressive model and relative power contribution Okito Yamashita, Norihiro Sadato 2,4, Tomohisa Okada 3, Tohru Ozaki Institute
More informationResponse surface models for the Elliott, Rothenberg, Stock DF-GLS unit-root test
Response surface models for the Elliott, Rothenberg, Stock DF-GLS unit-root test Christopher F Baum Jesús Otero Stata Conference, Baltimore, July 2017 Baum, Otero (BC, U. del Rosario) DF-GLS response surfaces
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 informationOn Autoregressive Order Selection Criteria
On Autoregressive Order Selection Criteria Venus Khim-Sen Liew Faculty of Economics and Management, Universiti Putra Malaysia, 43400 UPM, Serdang, Malaysia This version: 1 March 2004. Abstract This study
More informationIntroductory Econometrics
Based on the textbook by Wooldridge: : A Modern Approach Robert M. Kunst robert.kunst@univie.ac.at University of Vienna and Institute for Advanced Studies Vienna November 23, 2013 Outline Introduction
More informationAnalysis of the AIC Statistic for Optimal Detection of Small Changes in Dynamic Systems
Analysis of the AIC Statistic for Optimal Detection of Small Changes in Dynamic Systems Jeremy S. Conner and Dale E. Seborg Department of Chemical Engineering University of California, Santa Barbara, CA
More informationThe Influence of Additive Outliers on the Performance of Information Criteria to Detect Nonlinearity
The Influence of Additive Outliers on the Performance of Information Criteria to Detect Nonlinearity Saskia Rinke 1 Leibniz University Hannover Abstract In this paper the performance of information criteria
More informationISSN Article. Simulation Study of Direct Causality Measures in Multivariate Time Series
Entropy 2013, 15, 2635-2661; doi:10.3390/e15072635 OPEN ACCESS entropy ISSN 1099-4300 www.mdpi.com/journal/entropy Article Simulation Study of Direct Causality Measures in Multivariate Time Series Angeliki
More informationPolitical Science 236 Hypothesis Testing: Review and Bootstrapping
Political Science 236 Hypothesis Testing: Review and Bootstrapping Rocío Titiunik Fall 2007 1 Hypothesis Testing Definition 1.1 Hypothesis. A hypothesis is a statement about a population parameter The
More informationAn indicator for the number of clusters using a linear map to simplex structure
An indicator for the number of clusters using a linear map to simplex structure Marcus Weber, Wasinee Rungsarityotin, and Alexander Schliep Zuse Institute Berlin ZIB Takustraße 7, D-495 Berlin, Germany
More informationLeast Absolute Value vs. Least Squares Estimation and Inference Procedures in Regression Models with Asymmetric Error Distributions
Journal of Modern Applied Statistical Methods Volume 8 Issue 1 Article 13 5-1-2009 Least Absolute Value vs. Least Squares Estimation and Inference Procedures in Regression Models with Asymmetric Error
More informationOn Monitoring Shift in the Mean Processes with. Vector Autoregressive Residual Control Charts of. Individual Observation
Applied Mathematical Sciences, Vol. 8, 14, no. 7, 3491-3499 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/.12988/ams.14.44298 On Monitoring Shift in the Mean Processes with Vector Autoregressive Residual
More informationDETECTING PROCESS STATE CHANGES BY NONLINEAR BLIND SOURCE SEPARATION. Alexandre Iline, Harri Valpola and Erkki Oja
DETECTING PROCESS STATE CHANGES BY NONLINEAR BLIND SOURCE SEPARATION Alexandre Iline, Harri Valpola and Erkki Oja Laboratory of Computer and Information Science Helsinki University of Technology P.O.Box
More informationTime Series Forecasting: A Tool for Out - Sample Model Selection and Evaluation
AMERICAN JOURNAL OF SCIENTIFIC AND INDUSTRIAL RESEARCH 214, Science Huβ, http://www.scihub.org/ajsir ISSN: 2153-649X, doi:1.5251/ajsir.214.5.6.185.194 Time Series Forecasting: A Tool for Out - Sample Model
More informationTesting Error Correction in Panel data
University of Vienna, Dept. of Economics Master in Economics Vienna 2010 The Model (1) Westerlund (2007) consider the following DGP: y it = φ 1i + φ 2i t + z it (1) x it = x it 1 + υ it (2) where the stochastic
More informationRV Coefficient and Congruence Coefficient
RV Coefficient and Congruence Coefficient Hervé Abdi 1 1 Overview The congruence coefficient was first introduced by Burt (1948) under the name of unadjusted correlation as a measure of the similarity
More informationDarmstadt Discussion Papers in Economics
Darmstadt Discussion Papers in Economics The Effect of Linear Time Trends on Cointegration Testing in Single Equations Uwe Hassler Nr. 111 Arbeitspapiere des Instituts für Volkswirtschaftslehre Technische
More informationA Bootstrap Test for Causality with Endogenous Lag Length Choice. - theory and application in finance
CESIS Electronic Working Paper Series Paper No. 223 A Bootstrap Test for Causality with Endogenous Lag Length Choice - theory and application in finance R. Scott Hacker and Abdulnasser Hatemi-J April 200
More informationNonstationary time series models
13 November, 2009 Goals Trends in economic data. Alternative models of time series trends: deterministic trend, and stochastic trend. Comparison of deterministic and stochastic trend models The statistical
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 informationSTATS 200: Introduction to Statistical Inference. Lecture 29: Course review
STATS 200: Introduction to Statistical Inference Lecture 29: Course review Course review We started in Lecture 1 with a fundamental assumption: Data is a realization of a random process. The goal throughout
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 informationTIME-FREQUENCY AND TOPOGRAPHICAL ANALYSIS OF SIGNALS
XI Conference "Medical Informatics & Technologies" - 2006 Katarzyna. J. BLINOWSKA * Time-frequency methods, adaptive approximation, matching pursuit, multivariate autoregressive model, directed transfer
More informationTime Series: Theory and Methods
Peter J. Brockwell Richard A. Davis Time Series: Theory and Methods Second Edition With 124 Illustrations Springer Contents Preface to the Second Edition Preface to the First Edition vn ix CHAPTER 1 Stationary
More informationObservations and objectivity in statistics
PML seminar UvA 2010 Observations and objectivity in statistics Jan-Willem Romeijn Faculty of Philosophy University of Groningen Observing theory Observation is never independent of implicit, or explicit,
More informationMathematical Models of Dynamic Behavior of Individual Neural Networks of Central Nervous System
Mathematical Models of Dynamic Behavior of Individual Neural Networks of Central Nervous System Dimitra Despoina Pagania 1, Adam Adamopoulos 1,2 and Spiridon D. Likothanassis 1 1 Pattern Recognition Laboratory,
More informationReliability of Acceptance Criteria in Nonlinear Response History Analysis of Tall Buildings
Reliability of Acceptance Criteria in Nonlinear Response History Analysis of Tall Buildings M.M. Talaat, PhD, PE Senior Staff - Simpson Gumpertz & Heger Inc Adjunct Assistant Professor - Cairo University
More informationMinimum Message Length Autoregressive Model Order Selection
Minimum Message Length Autoregressive Model Order Selection Leigh J. Fitzgibbon School of Computer Science and Software Engineering, Monash University Clayton, Victoria 38, Australia leighf@csse.monash.edu.au
More informationResponse surface models for the Elliott, Rothenberg, Stock DF-GLS unit-root test
Response surface models for the Elliott, Rothenberg, Stock DF-GLS unit-root test Christopher F Baum Jesús Otero UK Stata Users Group Meetings, London, September 2017 Baum, Otero (BC, U. del Rosario) DF-GLS
More informationSelf Adaptive Particle Filter
Self Adaptive Particle Filter Alvaro Soto Pontificia Universidad Catolica de Chile Department of Computer Science Vicuna Mackenna 4860 (143), Santiago 22, Chile asoto@ing.puc.cl Abstract The particle filter
More informationSRMR in Mplus. Tihomir Asparouhov and Bengt Muthén. May 2, 2018
SRMR in Mplus Tihomir Asparouhov and Bengt Muthén May 2, 2018 1 Introduction In this note we describe the Mplus implementation of the SRMR standardized root mean squared residual) fit index for the models
More informationPramod K. Varshney. EECS Department, Syracuse University This research was sponsored by ARO grant W911NF
Pramod K. Varshney EECS Department, Syracuse University varshney@syr.edu This research was sponsored by ARO grant W911NF-09-1-0244 2 Overview of Distributed Inference U i s may be 1. Local decisions 2.
More informationResampling techniques for statistical modeling
Resampling techniques for statistical modeling Gianluca Bontempi Département d Informatique Boulevard de Triomphe - CP 212 http://www.ulb.ac.be/di Resampling techniques p.1/33 Beyond the empirical error
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 informationEcon 423 Lecture Notes: Additional Topics in Time Series 1
Econ 423 Lecture Notes: Additional Topics in Time Series 1 John C. Chao April 25, 2017 1 These notes are based in large part on Chapter 16 of Stock and Watson (2011). They are for instructional purposes
More informationMathematical Models of Dynamic Behavior of Individual Neural Networks of Central Nervous System
Mathematical Models of Dynamic Behavior of Individual Neural Networks of Central Nervous System Dimitra-Despoina Pagania,*, Adam Adamopoulos,2, and Spiridon D. Likothanassis Pattern Recognition Laboratory,
More informationNew Machine Learning Methods for Neuroimaging
New Machine Learning Methods for Neuroimaging Gatsby Computational Neuroscience Unit University College London, UK Dept of Computer Science University of Helsinki, Finland Outline Resting-state networks
More informationThe General Linear Model. Guillaume Flandin Wellcome Trust Centre for Neuroimaging University College London
The General Linear Model Guillaume Flandin Wellcome Trust Centre for Neuroimaging University College London SPM Course Lausanne, April 2012 Image time-series Spatial filter Design matrix Statistical Parametric
More informationAlgorithmisches Lernen/Machine Learning
Algorithmisches Lernen/Machine Learning Part 1: Stefan Wermter Introduction Connectionist Learning (e.g. Neural Networks) Decision-Trees, Genetic Algorithms Part 2: Norman Hendrich Support-Vector Machines
More informationVector autoregressions, VAR
1 / 45 Vector autoregressions, VAR Chapter 2 Financial Econometrics Michael Hauser WS17/18 2 / 45 Content Cross-correlations VAR model in standard/reduced form Properties of VAR(1), VAR(p) Structural VAR,
More informationHeteroskedasticity-Robust Inference in Finite Samples
Heteroskedasticity-Robust Inference in Finite Samples Jerry Hausman and Christopher Palmer Massachusetts Institute of Technology December 011 Abstract Since the advent of heteroskedasticity-robust standard
More informationOrder Selection for Vector Autoregressive Models
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 51, NO. 2, FEBRUARY 2003 427 Order Selection for Vector Autoregressive Models Stijn de Waele and Piet M. T. Broersen Abstract Order-selection criteria for vector
More informationEfficiently merging symbolic rules into integrated rules
Efficiently merging symbolic rules into integrated rules Jim Prentzas a, Ioannis Hatzilygeroudis b a Democritus University of Thrace, School of Education Sciences Department of Education Sciences in Pre-School
More informationQED. Queen s Economics Department Working Paper No Hypothesis Testing for Arbitrary Bounds. Jeffrey Penney Queen s University
QED Queen s Economics Department Working Paper No. 1319 Hypothesis Testing for Arbitrary Bounds Jeffrey Penney Queen s University Department of Economics Queen s University 94 University Avenue Kingston,
More informationTesting for spectral Granger causality
The Stata Journal (2015) 15, Number 4, pp. 1157 1166 Testing for spectral Granger causality Hüseyin Tastan Department of Economics Yildiz Technical University Istanbul, Turkey tastan@yildiz.edu.tr Abstract.
More informationLM threshold unit root tests
Lee, J., Strazicich, M.C., & Chul Yu, B. (2011). LM Threshold Unit Root Tests. Economics Letters, 110(2): 113-116 (Feb 2011). Published by Elsevier (ISSN: 0165-1765). http://0- dx.doi.org.wncln.wncln.org/10.1016/j.econlet.2010.10.014
More informationAPPLIED TIME SERIES ECONOMETRICS
APPLIED TIME SERIES ECONOMETRICS Edited by HELMUT LÜTKEPOHL European University Institute, Florence MARKUS KRÄTZIG Humboldt University, Berlin CAMBRIDGE UNIVERSITY PRESS Contents Preface Notation and Abbreviations
More informationError Reporting Recommendations: A Report of the Standards and Criteria Committee
Error Reporting Recommendations: A Report of the Standards and Criteria Committee Adopted by the IXS Standards and Criteria Committee July 26, 2000 1. Introduction The development of the field of x-ray
More informationImproved Statistical TRE Model When Using a Reference Frame
Improved Statistical TRE Model When Using a Reference Frame Andrew D. Wiles and Terry M. Peters Dept. of Medical Biophysics, The University of Western Ontario and Imaging Research Laboratories, Robarts
More informationImproved characterization of neural and behavioral response. state-space framework
Improved characterization of neural and behavioral response properties using point-process state-space framework Anna Alexandra Dreyer Harvard-MIT Division of Health Sciences and Technology Speech and
More informationBAYESIAN ESTIMATION OF LINEAR STATISTICAL MODEL BIAS
BAYESIAN ESTIMATION OF LINEAR STATISTICAL MODEL BIAS Andrew A. Neath 1 and Joseph E. Cavanaugh 1 Department of Mathematics and Statistics, Southern Illinois University, Edwardsville, Illinois 606, USA
More informationQuality Control Using Inferential Statistics in Weibull-based Reliability Analyses
GRAPHITE TESTING FOR NUCLEAR APPLICATIONS: THE SIGNIFICANCE OF TEST SPECIMEN VOLUME AND GEOMETRY AND THE STATISTICAL SIGNIFICANCE OF TEST SPECIMEN POPULATION 1 STP 1578, 2014 / available online at www.astm.org
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 information2 T. Schneider and A. Neumaier m-dimensional random vectors with zero mean and covariance matrix C. The m- dimensional parameter vector w of intercept
Algorithm: ARfit A Matlab Package for Estimation and Spectral Decomposition of Multivariate Autoregressive Models Tapio Schneider Princeton University and Arnold Neumaier Universitat Wien ARfit is a collection
More informationFACTOR ANALYSIS AND MULTIDIMENSIONAL SCALING
FACTOR ANALYSIS AND MULTIDIMENSIONAL SCALING Vishwanath Mantha Department for Electrical and Computer Engineering Mississippi State University, Mississippi State, MS 39762 mantha@isip.msstate.edu ABSTRACT
More informationProbability. We will now begin to explore issues of uncertainty and randomness and how they affect our view of nature.
Probability We will now begin to explore issues of uncertainty and randomness and how they affect our view of nature. We will explore in lab the differences between accuracy and precision, and the role
More informationPerformance Evaluation and Comparison
Outline Hong Chang Institute of Computing Technology, Chinese Academy of Sciences Machine Learning Methods (Fall 2012) Outline Outline I 1 Introduction 2 Cross Validation and Resampling 3 Interval Estimation
More informationVolume 03, Issue 6. Comparison of Panel Cointegration Tests
Volume 03, Issue 6 Comparison of Panel Cointegration Tests Deniz Dilan Karaman Örsal Humboldt University Berlin Abstract The main aim of this paper is to compare the size and size-adjusted power properties
More informationTesting Random Effects in Two-Way Spatial Panel Data Models
Testing Random Effects in Two-Way Spatial Panel Data Models Nicolas Debarsy May 27, 2010 Abstract This paper proposes an alternative testing procedure to the Hausman test statistic to help the applied
More informationSequential Procedure for Testing Hypothesis about Mean of Latent Gaussian Process
Applied Mathematical Sciences, Vol. 4, 2010, no. 62, 3083-3093 Sequential Procedure for Testing Hypothesis about Mean of Latent Gaussian Process Julia Bondarenko Helmut-Schmidt University Hamburg University
More informationUsing LSTM Network in Face Classification Problems
Using LS Network in Face Classification Problems Débora C. Corrêa, Denis H. P. Salvadeo, Alexandre L.. Levada, José H. Saito, Nelson D. A. ascarenhas, Jander oreira Departamento de Computação, Universidade
More information401 Review. 6. Power analysis for one/two-sample hypothesis tests and for correlation analysis.
401 Review Major topics of the course 1. Univariate analysis 2. Bivariate analysis 3. Simple linear regression 4. Linear algebra 5. Multiple regression analysis Major analysis methods 1. Graphical analysis
More informationA Review of the E 3 Algorithm: Near-Optimal Reinforcement Learning in Polynomial Time
A Review of the E 3 Algorithm: Near-Optimal Reinforcement Learning in Polynomial Time April 16, 2016 Abstract In this exposition we study the E 3 algorithm proposed by Kearns and Singh for reinforcement
More informationApproximating Bayesian Posterior Means Using Multivariate Gaussian Quadrature
Approximating Bayesian Posterior Means Using Multivariate Gaussian Quadrature John A.L. Cranfield Paul V. Preckel Songquan Liu Presented at Western Agricultural Economics Association 1997 Annual Meeting
More informationSimulating Properties of the Likelihood Ratio Test for a Unit Root in an Explosive Second Order Autoregression
Simulating Properties of the Likelihood Ratio est for a Unit Root in an Explosive Second Order Autoregression Bent Nielsen Nuffield College, University of Oxford J James Reade St Cross College, University
More informationComputing FMRI Activations: Coefficients and t-statistics by Detrending and Multiple Regression
Computing FMRI Activations: Coefficients and t-statistics by Detrending and Multiple Regression Daniel B. Rowe and Steven W. Morgan Division of Biostatistics Medical College of Wisconsin Technical Report
More informationLinear Combinations of Optic Flow Vectors for Estimating Self-Motion a Real-World Test of a Neural Model
Linear Combinations of Optic Flow Vectors for Estimating Self-Motion a Real-World Test of a Neural Model Matthias O. Franz MPI für biologische Kybernetik Spemannstr. 38 D-72076 Tübingen, Germany mof@tuebingen.mpg.de
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 informationThe Selection of Weighting Functions For Linear Arrays Using Different Techniques
The Selection of Weighting Functions For Linear Arrays Using Different Techniques Dr. S. Ravishankar H. V. Kumaraswamy B. D. Satish 17 Apr 2007 Rajarshi Shahu College of Engineering National Conference
More informationRidge Regression Revisited
Ridge Regression Revisited Paul M.C. de Boer Christian M. Hafner Econometric Institute Report EI 2005-29 In general ridge (GR) regression p ridge parameters have to be determined, whereas simple ridge
More informationApproximation of Average Run Length of Moving Sum Algorithms Using Multivariate Probabilities
Syracuse University SURFACE Electrical Engineering and Computer Science College of Engineering and Computer Science 3-1-2010 Approximation of Average Run Length of Moving Sum Algorithms Using Multivariate
More informationFORECASTING YIELD PER HECTARE OF RICE IN ANDHRA PRADESH
International Journal of Mathematics and Computer Applications Research (IJMCAR) ISSN 49-6955 Vol. 3, Issue 1, Mar 013, 9-14 TJPRC Pvt. Ltd. FORECASTING YIELD PER HECTARE OF RICE IN ANDHRA PRADESH R. RAMAKRISHNA
More information