Design Efficiency

Total Page:16

File Type:pdf, Size:1020Kb

Design Efficiency This article was originally published in Brain Mapping: An Encyclopedic Reference, published by Elsevier, and the attached copy is provided by Elsevier for the author's benefit and for the benefit of the author's institution, for non-commercial research and educational use including without limitation use in instruction at your institution, sending it to specific colleagues who you know, and providing a copy to your institution’s administrator. All other uses, reproduction and distribution, including without limitation commercial reprints, selling or licensing copies or access, or posting on open internet sites, your personal or institution’s website or repository, are prohibited. For exceptions, permission may be sought for such use through Elsevier's permissions site at: http://www.elsevier.com/locate/permissionusematerial Henson R.N. (2015) Design Efficiency. In: Arthur W. Toga, editor. Brain Mapping: An Encyclopedic Reference, vol. 1, pp. 489-494. Academic Press: Elsevier. Author's personal copy Design Efficiency RN Henson, MRC Cognition and Brain Sciences Unit, Cambridge, UK ã 2015 Elsevier Inc. All rights reserved. Abbreviations GLM General linear model AR(p) Autoregressive model of order p HRF Hemodynamic response function BOLD Blood oxygenation level-dependent (signal ReML Restricted maximum likelihood normally measured with fMRI) SOA Stimulus-onset asynchrony (SOAmin ¼minimal DCT Discrete cosine transform SOA) FIR Finite impulse response (basis set) TR Interscan interval (repetition time) fMRI Functional magnetic resonance imaging ÀÁ À1 Formal Definition of Efficiency e ¼ 1= cXTX cT [3] which can be seen as inversely related to the denominator of The general linear model (GLM) normally used for mass- the T-statistic in eqn [2]. Thus, if we increase e, we also increase univariate statistical analysis of functional magnetic resonance  imaging (fMRI) data can be written for a single voxel as T. (For multiple contrasts, where c is an M P matrix of M contrasts, such as an F-contrast, we can define the average ÀÁ T À1 T 2 efficiency as 1/trace{c(X X) c )}.) y ¼ Xb þ e, e N 0, s Ce [1] Note that the scaling of e is arbitrary (depending on the scaling of the contrast, scaling of regressors, and number of where y is an NÂ1 column vector of the data time series scans), so the precise relationship between e and T is best sampled every TR for N scans, X is an NÂP design matrix in assumed only to be monotonic. Note also that this statement which the P columns are regressors for the time series of pre- assumes that the estimate of the error variance (s^2) is indepen- dicted experimental effects, b is a PÂ1 column vector of dent of the design (X), which may not always be true (see parameters for each regressor in X (whose values are estimated later). Given these assumptions, and that the contrasts are when fitting the model to the data), and e is NÂ1 vector of specified a priori, then to maximize the efficiency of our residual errors. The second expression in eqn [1] denotes that design, we simply need to vary X. We now consider how X is the residuals come from a zero-mean, multivariate normal defined for fMRI. (Gaussian) distribution with covariance Ce. Normally, the residuals are assumed to be drawn independently from the same distribution (white residuals), or if not, then the data and model are filtered, or prewhitened, by an estimate of the 2 HRF Convolution error covariance (see later). This means that Ce ¼s I, corre- sponding to an NÂN identity matrix (I) scaled by a single We can start by assuming that stimuli elicit brief bursts of variance term s2. neural activity, or events, which are modeled by delta functions Assuming white residuals, the parameters can be estimated by every time a stimulus is presented. Then, for the jth of N event minimizing the sum of squares of the residuals, to give the so- j types (conditions), the neural activity over time, or neural time called ordinary least squares (OLS)estimates,b^. The planned com- course, u (t), can be expressed as parisons we want to test with our experiment are a linear combi- j  ¼ ðÞ nation of these parameter estimates, specified by a 1 Pcontrast i XNi j ÀÁ c ¼ ½À vector, . For example, c 1 1 would test whether the ujðÞ¼t d t À Tji parameter estimate for the first of two regressors is greater than i¼1 the second. Significance can be assessed by a T-statistic, defined by ¼ ... d where Tji is a vector of i 1 Ni(j) onset times and is the Dirac delta function. With fMRI, we do not measure neural cb^ TðÞ¼df qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiÀÁ [2] activity directly, but rather the delayed and dispersed BOLD À1 cXTX cTs^2 impulse response, b(t), where t indexes poststimulus time (e.g., from 0 to 30 s). Given that b(t) may vary across voxels s^2 eTe T where is the error variance estimated by /df (where (and individuals), it can be modeled by linear combination of denotes the transpose of a matrix) and the degrees of freedom, t Nk hemodynamic response functions (HRFs), hk( ): df, are defined by NÀrank(X). The probability, p, of getting a ^ ¼ value of T or greater under the null hypothesis that cb ¼ 0, kXNk ðÞ¼t b ðÞt given the df, can then be calculated from Student’s T-distribu- bj kjhk ¼ tion, and the null hypothesis rejected if, for example, p<0.05. k 1 We are now in the position to define the efficiency of a where bkj are the parameters to be estimated for each HRF and contrast, e,as condition (and voxel). Brain Mapping: An Encyclopedic Reference http://dx.doi.org/10.1016/B978-0-12-397025-1.00321-3 489 Brain Mapping: An Encyclopedic Reference, (2015), vol. 1, pp. 489-494 Author's personal copy 490 INTRODUCTION TO METHODS AND MODELING | Design Efficiency Assuming that BOLD responses summate linearly (though of the BOLD impulse response, then a more general set such see later), the predicted BOLD time course over the experi- as an FIR is necessary. Using a canonical HRF would corre- ment, x(t), can then be expressed as the convolution of the spond to what Liu, Frank, Wong, and Buxton (2001) called neural time courses by the HRFs: detection power, while using an FIR would correspond to what they called estimation efficiency. This is important because ¼ ¼ ¼ ¼ ðÞ jXNj jXNj kXNk i XNi j ÀÁ the choice of HRF affects the optimal experimental design xtðÞ¼ ujðÞ t bðÞ¼t b hk t À Tji [4] kj (see later). j¼1 j¼1 k¼1 i¼1 resulting in a linearly separable equation that can be repre- ¼ sented by a design matrix X with P NjNk columns. At one extreme, we can assumed a fixed shape for the BOLD Filtering response by using a single canonical HRF (i.e., Nk ¼1). At the other extreme, we can make no assumptions about the shape So far, we have considered definition of the signal, x(t), but the of the BOLD response (up to a certain frequency limit) by other factor that affects the T-statistic in eqn [2] is the noise s^2 using a so-called finite impulse response (FIR) set (see Figure 1; variance, . fMRI is known to have a preponderance of low- for multiple basis functions, the contrasts become c INk , frequency noise, caused, for example, by scanner drift and by where c is a contrast across the Nj event types and INk is an biorhythms (e.g., pulse and respiration) that are aliased by Nk ÂNk identity matrix for the Nk basis functions). Normally, slower sample rates (1/TR). A common strategy therefore is to one is only interested in the magnitude of a BOLD response, high-pass filter the data. An example matrix, F, for implement- in which case a single canonical HRF is sufficient to estimate ing high-pass filtering within the GLM using a discrete cosine efficiency a priori (by assuming that a canonical HRF is a transform (DCT) set is shown in Figure 1. The reduction in sufficient approximation on average across voxels and individ- noise will improve sensitivity, as long as the filtering does not uals). If however one is interested in estimating the shape remove excessive signal too. Heuristics suggest that an Canonical HRF SOA t min hk( ) u1 (t) Amplitude (a.u.) and or 0 0 5 10 15 20 25 30 u (t) 2 Amplitude (a.u.) Post-stimulus time (s) Transition table 0 200 400 600 800 1000 FIR HRF set Time (s) t hk( ) x1 (t) X x2 (t) Amplitude (a.u.) 0 5 10 15 20 25 Post-stimulus time (s) Scans F c 1 0 e = 1/(c((KX)T (KX))-1 cT ) −1 Figure 1 Ingredients for efficiency: the minimal SOA, SOAmin, and stimulus transition table determine the neural time course, uj(t), which is convolved with HRFs of poststimulus time, hk(t), to create the design matrix, X. This, together with an a priori contrast, c, and any (high-pass) filter matrix, À1 K (here generated by K¼IN ÀFF ), then determines the efficiency, e. Brain Mapping: An Encyclopedic Reference, (2015), vol. 1, pp. 489-494 Author's personal copy INTRODUCTION TO METHODS AND MODELING | Design Efficiency 491 approximate inflection in the noise power spectrum typically baseline. For a canonical HRF, the efficiency of these two occurs at around 1/120 s, which is why it is inadvisable to have contrasts is plotted against SOA in Figure 2(a). As can be min designs with changes in signal slower than this (e.g., alternat- seen, the optimal SOA for the common effect is around 18 s, ing blocks of more than 60 s; see later).
Recommended publications
  • Lecture 13: Simple Linear Regression in Matrix Format
    11:55 Wednesday 14th October, 2015 See updates and corrections at http://www.stat.cmu.edu/~cshalizi/mreg/ Lecture 13: Simple Linear Regression in Matrix Format 36-401, Section B, Fall 2015 13 October 2015 Contents 1 Least Squares in Matrix Form 2 1.1 The Basic Matrices . .2 1.2 Mean Squared Error . .3 1.3 Minimizing the MSE . .4 2 Fitted Values and Residuals 5 2.1 Residuals . .7 2.2 Expectations and Covariances . .7 3 Sampling Distribution of Estimators 8 4 Derivatives with Respect to Vectors 9 4.1 Second Derivatives . 11 4.2 Maxima and Minima . 11 5 Expectations and Variances with Vectors and Matrices 12 6 Further Reading 13 1 2 So far, we have not used any notions, or notation, that goes beyond basic algebra and calculus (and probability). This has forced us to do a fair amount of book-keeping, as it were by hand. This is just about tolerable for the simple linear model, with one predictor variable. It will get intolerable if we have multiple predictor variables. Fortunately, a little application of linear algebra will let us abstract away from a lot of the book-keeping details, and make multiple linear regression hardly more complicated than the simple version1. These notes will not remind you of how matrix algebra works. However, they will review some results about calculus with matrices, and about expectations and variances with vectors and matrices. Throughout, bold-faced letters will denote matrices, as a as opposed to a scalar a. 1 Least Squares in Matrix Form Our data consists of n paired observations of the predictor variable X and the response variable Y , i.e., (x1; y1);::: (xn; yn).
    [Show full text]
  • Introducing the Game Design Matrix: a Step-By-Step Process for Creating Serious Games
    Air Force Institute of Technology AFIT Scholar Theses and Dissertations Student Graduate Works 3-2020 Introducing the Game Design Matrix: A Step-by-Step Process for Creating Serious Games Aaron J. Pendleton Follow this and additional works at: https://scholar.afit.edu/etd Part of the Educational Assessment, Evaluation, and Research Commons, Game Design Commons, and the Instructional Media Design Commons Recommended Citation Pendleton, Aaron J., "Introducing the Game Design Matrix: A Step-by-Step Process for Creating Serious Games" (2020). Theses and Dissertations. 4347. https://scholar.afit.edu/etd/4347 This Thesis is brought to you for free and open access by the Student Graduate Works at AFIT Scholar. It has been accepted for inclusion in Theses and Dissertations by an authorized administrator of AFIT Scholar. For more information, please contact [email protected]. INTRODUCING THE GAME DESIGN MATRIX: A STEP-BY-STEP PROCESS FOR CREATING SERIOUS GAMES THESIS Aaron J. Pendleton, Captain, USAF AFIT-ENG-MS-20-M-054 DEPARTMENT OF THE AIR FORCE AIR UNIVERSITY AIR FORCE INSTITUTE OF TECHNOLOGY Wright-Patterson Air Force Base, Ohio DISTRIBUTION STATEMENT A APPROVED FOR PUBLIC RELEASE; DISTRIBUTION UNLIMITED. The views expressed in this document are those of the author and do not reflect the official policy or position of the United States Air Force, the United States Department of Defense or the United States Government. This material is declared a work of the U.S. Government and is not subject to copyright protection in the United States. AFIT-ENG-MS-20-M-054 INTRODUCING THE GAME DESIGN MATRIX: A STEP-BY-STEP PROCESS FOR CREATING LEARNING OBJECTIVE BASED SERIOUS GAMES THESIS Presented to the Faculty Department of Electrical and Computer Engineering Graduate School of Engineering and Management Air Force Institute of Technology Air University Air Education and Training Command in Partial Fulfillment of the Requirements for the Degree of Master of Science in Cyberspace Operations Aaron J.
    [Show full text]
  • Stat 5102 Notes: Regression
    Stat 5102 Notes: Regression Charles J. Geyer April 27, 2007 In these notes we do not use the “upper case letter means random, lower case letter means nonrandom” convention. Lower case normal weight letters (like x and β) indicate scalars (real variables). Lowercase bold weight letters (like x and β) indicate vectors. Upper case bold weight letters (like X) indicate matrices. 1 The Model The general linear model has the form p X yi = βjxij + ei (1.1) j=1 where i indexes individuals and j indexes different predictor variables. Ex- plicit use of (1.1) makes theory impossibly messy. We rewrite it as a vector equation y = Xβ + e, (1.2) where y is a vector whose components are yi, where X is a matrix whose components are xij, where β is a vector whose components are βj, and where e is a vector whose components are ei. Note that y and e have dimension n, but β has dimension p. The matrix X is called the design matrix or model matrix and has dimension n × p. As always in regression theory, we treat the predictor variables as non- random. So X is a nonrandom matrix, β is a nonrandom vector of unknown parameters. The only random quantities in (1.2) are e and y. As always in regression theory the errors ei are independent and identi- cally distributed mean zero normal. This is written as a vector equation e ∼ Normal(0, σ2I), where σ2 is another unknown parameter (the error variance) and I is the identity matrix. This implies y ∼ Normal(µ, σ2I), 1 where µ = Xβ.
    [Show full text]
  • Uncertainty of the Design and Covariance Matrices in Linear Statistical Model*
    Acta Univ. Palacki. Olomuc., Fac. rer. nat., Mathematica 48 (2009) 61–71 Uncertainty of the design and covariance matrices in linear statistical model* Lubomír KUBÁČEK 1, Jaroslav MAREK 2 Department of Mathematical Analysis and Applications of Mathematics, Faculty of Science, Palacký University, tř. 17. listopadu 12, 771 46 Olomouc, Czech Republic 1e-mail: [email protected] 2e-mail: [email protected] (Received January 15, 2009) Abstract The aim of the paper is to determine an influence of uncertainties in design and covariance matrices on estimators in linear regression model. Key words: Linear statistical model, uncertainty, design matrix, covariance matrix. 2000 Mathematics Subject Classification: 62J05 1 Introduction Uncertainties in entries of design and covariance matrices influence the variance of estimators and cause their bias. A problem occurs mainly in a linearization of nonlinear regression models, where the design matrix is created by deriva- tives of some functions. The question is how precise must these derivatives be. Uncertainties of covariance matrices must be suppressed under some reasonable bound as well. The aim of the paper is to give the simple rules which enables us to decide how many ciphers an entry of the mentioned matrices must be consisted of. *Supported by the Council of Czech Government MSM 6 198 959 214. 61 62 Lubomír KUBÁČEK, Jaroslav MAREK 2 Symbols used In the following text a linear regression model (in more detail cf. [2]) is denoted as k Y ∼n (Fβ, Σ), β ∈ R , (1) where Y is an n-dimensional random vector with the mean value E(Y) equal to Fβ and with the covariance matrix Var(Y)=Σ.ThesymbolRk means the k-dimensional linear vector space.
    [Show full text]
  • The Concept of a Generalized Inverse for Matrices Was Introduced by Moore(1920)
    J. Japanese Soc. Comp. Statist. 2(1989), 1-7 THE MOORE-PENROSE INVERSE MATRIX POR THE BALANCED ANOVA MODELS Byung Chun Kim* and Ha Sik Sunwoo* ABSTRACT Since the design matrix of the balanced linear model with no interactions has special form, the general solution of the normal equations can be easily found. From the relationships between the minimum norm least squares solution and the Moore-Penrose inverse we can obtain the explicit form of the Moore-Penrose inverse X+ of the design matrix of the model y = XƒÀ + ƒÃ for the balanced model with no interaction. 1. Introduction The concept of a generalized inverse for matrices was introduced by Moore(1920). He developed it in the context of linear transformations from n-dimensional to m-dimensional vector space over a complex field with usual Euclidean norm. Unaware of Moore's work, Penrose(1955) established the existence and the uniqueness of the generalized inverse that satisfies the Moore's condition under L2 norm. It is commonly known as the Moore-Penrose inverse. In the early 1970's many statisticians, Albert (1972), Ben-Israel and Greville(1974), Graybill(1976),Lawson and Hanson(1974), Pringle and Raynor(1971), Rao and Mitra(1971), Rao(1975), Schmidt(1976), and Searle(1971, 1982), worked on this subject and yielded sev- eral texts that treat the mathematics and its application to statistics in some depth in detail. Especially, Shinnozaki et. al.(1972a, 1972b) have surveyed the general methods to find the Moore-Penrose inverse in two subject: direct methods and iterative methods. Also they tested their accuracy for some cases.
    [Show full text]
  • Kriging Models for Linear Networks and Non‐Euclidean Distances
    Received: 13 November 2017 | Accepted: 30 December 2017 DOI: 10.1111/2041-210X.12979 RESEARCH ARTICLE Kriging models for linear networks and non-Euclidean distances: Cautions and solutions Jay M. Ver Hoef Marine Mammal Laboratory, NOAA Fisheries Alaska Fisheries Science Center, Abstract Seattle, WA, USA 1. There are now many examples where ecological researchers used non-Euclidean Correspondence distance metrics in geostatistical models that were designed for Euclidean dis- Jay M. Ver Hoef tance, such as those used for kriging. This can lead to problems where predictions Email: [email protected] have negative variance estimates. Technically, this occurs because the spatial co- Handling Editor: Robert B. O’Hara variance matrix, which depends on the geostatistical models, is not guaranteed to be positive definite when non-Euclidean distance metrics are used. These are not permissible models, and should be avoided. 2. I give a quick review of kriging and illustrate the problem with several simulated examples, including locations on a circle, locations on a linear dichotomous net- work (such as might be used for streams), and locations on a linear trail or road network. I re-examine the linear network distance models from Ladle, Avgar, Wheatley, and Boyce (2017b, Methods in Ecology and Evolution, 8, 329) and show that they are not guaranteed to have a positive definite covariance matrix. 3. I introduce the reduced-rank method, also called a predictive-process model, for creating valid spatial covariance matrices with non-Euclidean distance metrics. It has an additional advantage of fast computation for large datasets. 4. I re-analysed the data of Ladle et al.
    [Show full text]
  • Linear Regression with Shuffled Data: Statistical and Computational
    3286 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 64, NO. 5, MAY 2018 Linear Regression With Shuffled Data: Statistical and Computational Limits of Permutation Recovery Ashwin Pananjady ,MartinJ.Wainwright,andThomasA.Courtade Abstract—Consider a noisy linear observation model with an permutation matrix, and w Rn is observation noise. We refer ∈ unknown permutation, based on observing y !∗ Ax∗ w, to the setting where w 0asthenoiseless case.Aswithlinear d = + where x∗ R is an unknown vector, !∗ is an unknown = ∈ n regression, there are two settings of interest, corresponding to n n permutation matrix, and w R is additive Gaussian noise.× We analyze the problem of∈ permutation recovery in a whether the design matrix is (i) deterministic (the fixed design random design setting in which the entries of matrix A are case), or (ii) random (the random design case). drawn independently from a standard Gaussian distribution and There are also two complementary problems of interest: establish sharp conditions on the signal-to-noise ratio, sample recovery of the unknown !∗,andrecoveryoftheunknownx∗. size n,anddimensiond under which !∗ is exactly and approx- In this paper, we focus on the former problem; the latter imately recoverable. On the computational front, we show that problem is also known as unlabelled sensing [5]. the maximum likelihood estimate of !∗ is NP-hard to compute for general d,whilealsoprovidingapolynomialtimealgorithm The observation model (1) is frequently encountered in when d 1. scenarios where there is uncertainty in the order in which = Index Terms—Correspondenceestimation,permutationrecov- measurements are taken. An illustrative example is that of ery, unlabeled sensing, information-theoretic bounds, random sampling in the presence of jitter [6], in which the uncertainty projections.
    [Show full text]
  • Week 7: Multiple Regression
    Week 7: Multiple Regression Brandon Stewart1 Princeton October 24, 26, 2016 1These slides are heavily influenced by Matt Blackwell, Adam Glynn, Jens Hainmueller and Danny Hidalgo. Stewart (Princeton) Week 7: Multiple Regression October 24, 26, 2016 1 / 145 Where We've Been and Where We're Going... Last Week I regression with two variables I omitted variables, multicollinearity, interactions This Week I Monday: F matrix form of linear regression I Wednesday: F hypothesis tests Next Week I break! I then ::: regression in social science Long Run I probability ! inference ! regression Questions? Stewart (Princeton) Week 7: Multiple Regression October 24, 26, 2016 2 / 145 1 Matrix Algebra Refresher 2 OLS in matrix form 3 OLS inference in matrix form 4 Inference via the Bootstrap 5 Some Technical Details 6 Fun With Weights 7 Appendix 8 Testing Hypotheses about Individual Coefficients 9 Testing Linear Hypotheses: A Simple Case 10 Testing Joint Significance 11 Testing Linear Hypotheses: The General Case 12 Fun With(out) Weights Stewart (Princeton) Week 7: Multiple Regression October 24, 26, 2016 3 / 145 Why Matrices and Vectors? Here's one way to write the full multiple regression model: yi = β0 + xi1β1 + xi2β2 + ··· + xiK βK + ui Notation is going to get needlessly messy as we add variables Matrices are clean, but they are like a foreign language You need to build intuitions over a long period of time (and they will return in Soc504) Reminder of Parameter Interpretation: β1 is the effect of a one-unit change in xi1 conditional on all other xik . We are going to review the key points quite quickly just to refresh the basics.
    [Show full text]
  • Stat 714 Linear Statistical Models
    STAT 714 LINEAR STATISTICAL MODELS Fall, 2010 Lecture Notes Joshua M. Tebbs Department of Statistics The University of South Carolina TABLE OF CONTENTS STAT 714, J. TEBBS Contents 1 Examples of the General Linear Model 1 2 The Linear Least Squares Problem 13 2.1 Least squares estimation . 13 2.2 Geometric considerations . 15 2.3 Reparameterization . 25 3 Estimability and Least Squares Estimators 28 3.1 Introduction . 28 3.2 Estimability . 28 3.2.1 One-way ANOVA . 35 3.2.2 Two-way crossed ANOVA with no interaction . 37 3.2.3 Two-way crossed ANOVA with interaction . 39 3.3 Reparameterization . 40 3.4 Forcing least squares solutions using linear constraints . 46 4 The Gauss-Markov Model 54 4.1 Introduction . 54 4.2 The Gauss-Markov Theorem . 55 4.3 Estimation of σ2 in the GM model . 57 4.4 Implications of model selection . 60 4.4.1 Underfitting (Misspecification) . 60 4.4.2 Overfitting . 61 4.5 The Aitken model and generalized least squares . 63 5 Distributional Theory 68 5.1 Introduction . 68 i TABLE OF CONTENTS STAT 714, J. TEBBS 5.2 Multivariate normal distribution . 69 5.2.1 Probability density function . 69 5.2.2 Moment generating functions . 70 5.2.3 Properties . 72 5.2.4 Less-than-full-rank normal distributions . 73 5.2.5 Independence results . 74 5.2.6 Conditional distributions . 76 5.3 Noncentral χ2 distribution . 77 5.4 Noncentral F distribution . 79 5.5 Distributions of quadratic forms . 81 5.6 Independence of quadratic forms . 85 5.7 Cochran's Theorem .
    [Show full text]
  • Fitting Differential Equations to Functional Data: Principal Differential Analysis
    19 Fitting differential equations to functional data: Principal differential analysis 19.1 Introduction Now that we have fastened a belt of tools around our waists for tinker- ing with differential equations, we return to the problems introduced in Chapter 17 ready to get down to some serious work. Using a differential equation as a modelling object involves concepts drawn from both the functional linear model and from principal com- ponents analysis. A differential equation can certainly capture the shape features in both the curve and its derivative for a single functional datum such as the oil refinery observation shown in Figure 1.4. But because the set of solutions to a differential equation is an entire function space, it can also model variation across observations when N>1. In this sense, it also has the flavor of principal components analysis where we find a subspace of functions able to capture the dominant modes of variation in the data. We have, then, a question of emphasis or perspective. On one hand, the data analyst may want to capture important features of the dynamics of a single observation, and thus look within the m-dimensional space of solutions of an estimated equation to find that which gives the best account of the data. On the other hand, the goal may be to see how much functional variation can be explained across multiple realizations of a process. Thus, linear modelling and variance decomposition merge into one analysis in this environment. We introduce a new term here: principal differential analysis means the fitting of a differential equation to noisy data so as to capture either the 328 19.
    [Show full text]
  • Linear Regression in Matrix Form
    Statistics 512: Applied Linear Models Topic 3 Topic Overview This topic will cover • thinking in terms of matrices • regression on multiple predictor variables • case study: CS majors • Text Example (KNNL 236) Chapter 5: Linear Regression in Matrix Form The SLR Model in Scalar Form iid 2 Yi = β0 + β1Xi + i where i ∼ N(0,σ ) Consider now writing an equation for each observation: Y1 = β0 + β1X1 + 1 Y2 = β0 + β1X2 + 2 . Yn = β0 + β1Xn + n TheSLRModelinMatrixForm Y1 β0 + β1X1 1 Y2 β0 + β1X2 2 . = . + . . . . Yn β0 + β1Xn n Y X 1 1 1 1 Y X 2 1 2 β0 2 . = . + . . . β1 . Yn 1 Xn n (I will try to use bold symbols for matrices. At first, I will also indicate the dimensions as a subscript to the symbol.) 1 • X is called the design matrix. • β is the vector of parameters • is the error vector • Y is the response vector The Design Matrix 1 X1 1 X2 Xn×2 = . . 1 Xn Vector of Parameters β0 β2×1 = β1 Vector of Error Terms 1 2 n×1 = . . n Vector of Responses Y1 Y2 Yn×1 = . . Yn Thus, Y = Xβ + Yn×1 = Xn×2β2×1 + n×1 2 Variance-Covariance Matrix In general, for any set of variables U1,U2,... ,Un,theirvariance-covariance matrix is defined to be 2 σ {U1} σ{U1,U2} ··· σ{U1,Un} . σ{U ,U } σ2{U } ... σ2{ } 2 1 2 . U = . . .. .. σ{Un−1,Un} 2 σ{Un,U1} ··· σ{Un,Un−1} σ {Un} 2 where σ {Ui} is the variance of Ui,andσ{Ui,Uj} is the covariance of Ui and Uj.
    [Show full text]
  • Some Matrix Results
    APPENDIX Some matrix results The main matrix results required for this book are given in this Appendix. Most of them are stated without proof as they are relatively straightforward and can be found in most books on matrix algebra; see for instance Searle (1966), Basilevsky (1983) and Graybill (1983). The properties of circulant matrices are considered in more detail as they have particular relevance in the book and are perhaps not widely known. A detailed account of circulant matrices can be found in Davis (1979). A.1 Notation An n x m matrix A = ((aij)) has n rows and m columns with aij being the element in the ith row andjth column. A' will denote the transpose of A. A is symmetrie if A = A'. If m = 1 then the matrix will be called a eolumn veetor and will usually be denoted by a lower ca se letter, a say. The transpose of a will be a row veetor a'. A diagonal matrix whose diagonal elements are those ofthe vector a, and whose ofT-diagonal elements are zero, will be denoted by a6• The inverse of a6 will be denoted by a- 6• More generally, if each element in a is raised to power p then the corresponding diagonal matrix will be denoted by a P6• Some special matrices and vectors are: In' the n x 1 vector with every element unity In = 1~, the n x n identity matrix Jnxm = I nl;", the n x m matrix with every element unity J n = J nxn Kn = (l/n)Jn the n x n matrix with every element n -1.
    [Show full text]