The Effects of Autocorrelated Noise and Biased HRF in fmri Analysis Error Rates
|
|
- Alaina Jefferson
- 5 years ago
- Views:
Transcription
1 The Effects of Autocorrelated Noise and Biased HRF in fmri Analysis Error Rates Ariana Anderson University of California, Los Angeles Departments of Psychiatry and Behavioral Sciences David Geffen School of Medicine July 29, / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
2 1 Background of fmri Imaging 2 Model-Driven Analysis Methods 3 Simulation 4 Conclusion 2 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
3 Mechanics of fmri The fmri Blood Oxygen Level Dependent (BOLD) signal is an indirect measure of neuronal activity captured during an fmri scan. Neuronal Activation Oxygenated Blood Influx Oxygen Concentration Increases Regional MR Signal Increases 3 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
4 Mechanics of fmri The fmri Blood Oxygen Level Dependent (BOLD) signal is an indirect measure of neuronal activity captured during an fmri scan. Analysis of the BOLD signal is performed under the assumption that neuronal activity coincides with increased blood flow. Neuronal Activation Oxygenated Blood Influx Oxygen Concentration Increases Regional MR Signal Increases 3 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
5 Mechanics of fmri Some regions with high CSF content will have a high-intensity signal under a T2-weighted scan. 4 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
6 Terminology Voxel: A 3-dimensional analogue of a pixel. A voxel is a region usually between 1 5mm 3 5 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
7 Terminology Voxel: A 3-dimensional analogue of a pixel. A voxel is a region usually between 1 5mm 3 Hemodynamic Response Function: A function modeling the change in blood flow due to some trigger. 5 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
8 6 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
9 General Linear Model Analysis Y = XB + ɛ, where Y = time course of voxels, X = design matrix and regressors, B = parameters, ɛ = error. Use parametric tests to infer how well a given voxel follows the time-course of a theoretical active voxel. 7 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
10 General Linear Model Analysis Y = XB + ɛ, where Y = time course of voxels, X = design matrix and regressors, B = parameters, ɛ = error. Use parametric tests to infer how well a given voxel follows the time-course of a theoretical active voxel. To get the design matrix for the task function, convolve the hemodynamic response function with task paradigm g f = b a f (τ)g(t τ)dτ 7 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
11 General Linear Model Analysis Y = XB + ɛ, where Y = time course of voxels, X = design matrix and regressors, B = parameters, ɛ = error. Use parametric tests to infer how well a given voxel follows the time-course of a theoretical active voxel. To get the design matrix for the task function, convolve the hemodynamic response function with task paradigm g f = b a f (τ)g(t τ)dτ This method allows one to account for confounding effects such as drift or respiration, by incorporating them into the design matrix. 7 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
12 8 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
13 HRF Variations Aging 9 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
14 HRF Variations Aging Disease (Alzheimers), Injury (Stroke) 9 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
15 HRF Variations Aging Disease (Alzheimers), Injury (Stroke) Pharmocological treatments 9 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
16 HRF Variations Aging Disease (Alzheimers), Injury (Stroke) Pharmocological treatments CO 2 Concentrations 9 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
17 10 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
18 Y = XB + ɛ, where Y = time course of voxels, X = design matrix and regressors, B = parameters, ɛ = error. The Gauss-Markov Theorem guarantees that this estimate is the Best Linear Unbiased Estimate (BLUE) if E(ɛ i ) = 0, V (ɛ i ) = σ 2 <, and cov(ɛ i, ɛ j ) = 0, for i j. 11 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
19 Y = XB + ɛ, where Y = time course of voxels, X = design matrix and regressors, B = parameters, ɛ = error. The Gauss-Markov Theorem guarantees that this estimate is the Best Linear Unbiased Estimate (BLUE) if E(ɛ i ) = 0, V (ɛ i ) = σ 2 <, and cov(ɛ i, ɛ j ) = 0, for i j. The final assumption is violated as the data are serially correlated over time both during rest and task periods (HRF). 11 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
20 Generalized Least Squares: Estimate ˆβ GLS s.t. ɛ N(0, σ 2 Σ) 12 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
21 Generalized Least Squares: Estimate ˆβ GLS s.t. ɛ N(0, σ 2 Σ) The GLS introduces a matrix P frequently used to decorrelate the residuals, and models the data as P 1 Y = P 1 X β + P 1 ɛ and minimizes the Mahalanobis distance to find β, such that ˆβ = argmin β (Y X β) Σ 1 (Y X β) 12 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
22 Generalized Least Squares: Estimate ˆβ GLS s.t. ɛ N(0, σ 2 Σ) The GLS introduces a matrix P frequently used to decorrelate the residuals, and models the data as P 1 Y = P 1 X β + P 1 ɛ and minimizes the Mahalanobis distance to find β, such that ˆβ = argmin β (Y X β) Σ 1 (Y X β) If P P = Σ, then Var(P 1 ɛ) = P 1 Var(ɛ)P 1 = σ 2 I 12 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
23 Generalized Least Squares: Estimate ˆβ GLS s.t. ɛ N(0, σ 2 Σ) The GLS introduces a matrix P frequently used to decorrelate the residuals, and models the data as P 1 Y = P 1 X β + P 1 ɛ and minimizes the Mahalanobis distance to find β, such that ˆβ = argmin β (Y X β) Σ 1 (Y X β) If P P = Σ, then Var(P 1 ɛ) = P 1 Var(ɛ)P 1 = σ 2 I This also guarantees that cov(ɛ i, ɛ j ) = 0 for i j and the conditions are met for this estimate of ˆβ GLS to be optimal as BLUE. In practice, P either imposed (coloring) or estimated (whitening), but the efficiency and accuracy of these approaches aren t guaranteed. 12 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
24 Autocorrelations in data result in serially correlated residuals, which in turn do not guarantee BLUE for ˆβ 13 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
25 Autocorrelations in data result in serially correlated residuals, which in turn do not guarantee BLUE for ˆβ HRF is assumed to be independent of covariates and location, which is not always true. 13 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
26 Simulating and measuring bias Simulate two time series of length 100 that differ by the addition of autocorrelated noise that is unrelated to the task being performed, and compare the parameter estimates ˆβ with those estimated from data that has uncorrelated (independent) additive white noise, 14 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
27 Simulating and measuring bias Simulate two time series of length 100 that differ by the addition of autocorrelated noise that is unrelated to the task being performed, and compare the parameter estimates ˆβ with those estimated from data that has uncorrelated (independent) additive white noise, y autocorrelated = β(h X ) + g E + ɛ 14 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
28 Simulating and measuring bias Simulate two time series of length 100 that differ by the addition of autocorrelated noise that is unrelated to the task being performed, and compare the parameter estimates ˆβ with those estimated from data that has uncorrelated (independent) additive white noise, y autocorrelated = β(h X ) + g E + ɛ y independent = β(h X ) + ɛ 14 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
29 Simulating and measuring bias Simulate two time series of length 100 that differ by the addition of autocorrelated noise that is unrelated to the task being performed, and compare the parameter estimates ˆβ with those estimated from data that has uncorrelated (independent) additive white noise, y autocorrelated = β(h X ) + g E + ɛ y independent = β(h X ) + ɛ 200,000 total realizations with randomization of β, X, g, E, h and ɛ 14 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
30 Parameter Settings X = 10 activations of length TR Unif [1, 5]. 15 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
31 Parameter Settings X = 10 activations of length TR Unif [1, 5]. β Unif [1, 5]. 15 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
32 Parameter Settings X = 10 activations of length TR Unif [1, 5]. β Unif [1, 5]. h is the HRF modeled by a double-gamma function dilated by δ Unif [0, 1]. 15 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
33 Parameter Settings X = 10 activations of length TR Unif [1, 5]. β Unif [1, 5]. h is the HRF modeled by a double-gamma function dilated by δ Unif [0, 1]. ɛ N(0, 1). 15 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
34 Parameter Settings X = 10 activations of length TR Unif [1, 5]. β Unif [1, 5]. h is the HRF modeled by a double-gamma function dilated by δ Unif [0, 1]. ɛ N(0, 1). Autocorrelation is introduced into y autocorrelated by convolving a random noise time series E = (η 1, η 2,..., η T ) with η N(0, 1), with an exponential decay function g(t) = r t where r Unif [0, 1]. 15 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
35 16 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
36 Bias Model Parameter HRF Known HRF Biased HRF Ignored y autocorrelated Mean y independent Mean y autocorrelated Variance y independent Variance Table: Mean and Variance of Bias, demonstrating that in the presence of autocorrelated noise, a biased HRF leads to drastic bias in the estimates of ˆβ. Autocorrelated noise is a known occurrence that is usually ignored in most theoretical models and analysis techniques. 17 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
37 Data HRF Task Length β ˆβ decay rate r δ HRF Bias y auto Known y auto Biased y auto Ignored y ind Known y ind Biased y ind Ignored Table: Correlations of Parameters with Bias = β ˆβ for yautocorrelated and y independent data. A higher correlation indicates that biased parameter assumptions lead to biased estimates of ˆβ, which in turn leads to even higher false positive and negatives in detecting active regions. 18 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
38 Typical Solution Include temporal derivative of the task function in design matrix: y autocorrelated = β(h X ) + g E + ɛ, where X = (x, x, confounders). 19 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
39 Model Sensitivity Specificity First Order (no derivatives) Second Order (with derivatives) Table: Derivative Adjustment within GLM Sensitivity: Probability of detection given activation occurred. 20 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
40 Model Sensitivity Specificity First Order (no derivatives) Second Order (with derivatives) Table: Derivative Adjustment within GLM Sensitivity: Probability of detection given activation occurred. Specificity: Probability of no detection given no activation occurred. 20 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
41 Positive correlation of the estimate ˆβ with bias shows that estimates that are large are more likely to be biased, implying a greater Type 1 error. 21 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
42 Positive correlation of the estimate ˆβ with bias shows that estimates that are large are more likely to be biased, implying a greater Type 1 error. The decay rate parameter r was positively correlated with bias, indicating that the longer the temporal autocorrelation present in the white noise, the greater the expected variance. 21 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
43 Conclusions Estimates become more biased and noisy both because of the addition of autocorrelated noise in the data and also because of mis-specification of the HRF. 22 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
44 Conclusions Estimates become more biased and noisy both because of the addition of autocorrelated noise in the data and also because of mis-specification of the HRF. When the HRF is ignored completely, the model doesn t exhibit changes in the estimates during the fitting of y autocorrelated data, but its bias is greater than for all other models and data. 22 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
45 Conclusions Estimates become more biased and noisy both because of the addition of autocorrelated noise in the data and also because of mis-specification of the HRF. When the HRF is ignored completely, the model doesn t exhibit changes in the estimates during the fitting of y autocorrelated data, but its bias is greater than for all other models and data. Models that include temporal derivatives perform worse than those without, in terms of sensitivity and specificity. 22 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
46 Conclusions Estimates become more biased and noisy both because of the addition of autocorrelated noise in the data and also because of mis-specification of the HRF. When the HRF is ignored completely, the model doesn t exhibit changes in the estimates during the fitting of y autocorrelated data, but its bias is greater than for all other models and data. Models that include temporal derivatives perform worse than those without, in terms of sensitivity and specificity. Suggestion: When temporal derivatives are modeled for, significance tests need to assess for both first order and derivative effects jointly. 22 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
47 Thank you! 23 / 23 Ariana Anderson Autocorrelation Induced Bias in fmri GLM Analysis
Advanced Data Modelling & Inference
Edinburgh 2015: biennial SPM course Advanced Data Modelling & Inference Cyril Pernet Centre for Clinical Brain Sciences (CCBS) Neuroimaging Sciences Modelling? Y = XB + e here is my model B might be biased
More informationInferential Brain Mapping
Inferential Brain Mapping Cyril Pernet Centre for Clinical Brain Sciences (CCBS) Neuroimaging Sciences Edinburgh Biennial SPM course 2019 @CyrilRPernet Overview Statistical inference is the process of
More informationMEA DISCUSSION PAPERS
Inference Problems under a Special Form of Heteroskedasticity Helmut Farbmacher, Heinrich Kögel 03-2015 MEA DISCUSSION PAPERS mea Amalienstr. 33_D-80799 Munich_Phone+49 89 38602-355_Fax +49 89 38602-390_www.mea.mpisoc.mpg.de
More informationInverse problems in functional brain imaging Identification of the hemodynamic response in fmri
Inverse problems in functional brain imaging Identification of the hemodynamic response in fmri Ph. Ciuciu1,2 philippe.ciuciu@cea.fr 1: CEA/NeuroSpin/LNAO May 7, 2010 www.lnao.fr 2: IFR49 GDR -ISIS Spring
More informationFunctional MRI Mapping Cognition
Outline Functional MRI Mapping Cognition Michael A. Yassa, B.A. Division of Psychiatric Neuro-imaging Psychiatry and Behavioral Sciences Johns Hopkins School of Medicine Why fmri? fmri - How it works Research
More informationPrediction of Successful Memory Encoding from fmri Data
Prediction of Successful Memory Encoding from fmri Data S.K. Balci 1, M.R. Sabuncu 1, J. Yoo 2, S.S. Ghosh 3, S. Whitfield-Gabrieli 2, J.D.E. Gabrieli 2 and P. Golland 1 1 CSAIL, MIT, Cambridge, MA, USA
More informationEdinburgh Imaging Academy online distance learning courses. Functional Imaging
Functional Imaging Semester 2 / Commences January 10 Credits Each Course is composed of Modules & Activities. Modules: BOLD Signal IMSc NI4R Experimental Design IMSc NI4R Pre-processing IMSc NI4R GLM IMSc
More informationTemporal preprocessing of fmri data
Temporal preprocessing of fmri data Blaise Frederick, Ph.D. McLean Hospital Brain Imaging Center Scope fmri data contains temporal noise and acquisition artifacts that complicate the interpretation of
More informationRepresentational similarity analysis
School of Psychology Representational similarity analysis Dr Ian Charest Representational similarity analysis representational dissimilarity matrices (RDMs) stimulus (e.g. images, sounds, other experimental
More informationCitation for published version (APA): Ebbes, P. (2004). Latent instrumental variables: a new approach to solve for endogeneity s.n.
University of Groningen Latent instrumental variables Ebbes, P. IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish to cite from it. Please check the document
More informationInstrumental Variables Estimation: An Introduction
Instrumental Variables Estimation: An Introduction Susan L. Ettner, Ph.D. Professor Division of General Internal Medicine and Health Services Research, UCLA The Problem The Problem Suppose you wish to
More informationModel-Based fmri Analysis. Will Alexander Dept. of Experimental Psychology Ghent University
Model-Based fmri Analysis Will Alexander Dept. of Experimental Psychology Ghent University Motivation Models (general) Why you ought to care Model-based fmri Models (specific) From model to analysis Extended
More informationInference with Difference-in-Differences Revisited
Inference with Difference-in-Differences Revisited M. Brewer, T- F. Crossley and R. Joyce Journal of Econometric Methods, 2018 presented by Federico Curci February 22nd, 2018 Brewer, Crossley and Joyce
More informationExperimental design. Experimental design. Experimental design. Guido van Wingen Department of Psychiatry Academic Medical Center
Experimental design Guido van Wingen Department of Psychiatry Academic Medical Center guidovanwingen@gmail.com Experimental design Experimental design Image time-series Spatial filter Design matrix Statistical
More informationTemporal preprocessing of fmri data
Temporal preprocessing of fmri data Blaise Frederick, Ph.D., Yunjie Tong, Ph.D. McLean Hospital Brain Imaging Center Scope This talk will summarize the sources and characteristics of unwanted temporal
More informationModel Specification and the Reliability of fmri Results: Implications for Longitudinal Neuroimaging Studies in Psychiatry
Model Specification and the Reliability of fmri Results: Implications for Longitudinal Neuroimaging Studies in Psychiatry Jay C. Fournier 1 *, Henry W. Chase 1, Jorge Almeida 2, Mary L. Phillips 1 1 Department
More informationStatistical Analysis of fmri Time-Series: A Critical Evaluation of the GLM Approach
Statistical Analysis of fmri Time-Series: A Critical Evaluation of the GLM Approach Martin M. Monti Princeton University, Dept. of Psychology April 18, 2006 NUMBER OF PAGES: 44 ABBREVIATED TITLE: Statistical
More informationDesign Optimization. Optimization?
Edinburgh 2015: biennial SPM course Design Optimization Cyril Pernet Centre for Clinical Brain Sciences (CCBS) Neuroimaging Sciences Optimization? Making a design that is the most favourable or desirable,
More informationLecture 14: Adjusting for between- and within-cluster covariates in the analysis of clustered data May 14, 2009
Measurement, Design, and Analytic Techniques in Mental Health and Behavioral Sciences p. 1/3 Measurement, Design, and Analytic Techniques in Mental Health and Behavioral Sciences Lecture 14: Adjusting
More informationContributions to Brain MRI Processing and Analysis
Contributions to Brain MRI Processing and Analysis Dissertation presented to the Department of Computer Science and Artificial Intelligence By María Teresa García Sebastián PhD Advisor: Prof. Manuel Graña
More informationComparing event-related and epoch analysis in blocked design fmri
Available online at www.sciencedirect.com R NeuroImage 18 (2003) 806 810 www.elsevier.com/locate/ynimg Technical Note Comparing event-related and epoch analysis in blocked design fmri Andrea Mechelli,
More informationDesign Optimization. SBIRC/SINAPSE University of Edinburgh
Design Optimization Cyril Pernet,, PhD SBIRC/SINAPSE University of Edinburgh fmri designs: overview Introduction: fmri noise fmri design types: blocked designs event-related designs mixed design adaptation
More informationSupplementary Information
Supplementary Information The neural correlates of subjective value during intertemporal choice Joseph W. Kable and Paul W. Glimcher a 10 0 b 10 0 10 1 10 1 Discount rate k 10 2 Discount rate k 10 2 10
More informationClassification and Statistical Analysis of Auditory FMRI Data Using Linear Discriminative Analysis and Quadratic Discriminative Analysis
International Journal of Innovative Research in Computer Science & Technology (IJIRCST) ISSN: 2347-5552, Volume-2, Issue-6, November-2014 Classification and Statistical Analysis of Auditory FMRI Data Using
More informationFREQUENCY DOMAIN HYBRID INDEPENDENT COMPONENT ANALYSIS OF FUNCTIONAL MAGNETIC RESONANCE IMAGING DATA
FREQUENCY DOMAIN HYBRID INDEPENDENT COMPONENT ANALYSIS OF FUNCTIONAL MAGNETIC RESONANCE IMAGING DATA J.D. Carew, V.M. Haughton, C.H. Moritz, B.P. Rogers, E.V. Nordheim, and M.E. Meyerand Departments of
More informationExperimental design of fmri studies
Experimental design of fmri studies Kerstin Preuschoff Computational Neuroscience Lab, EPFL LREN SPM Course Lausanne April 10, 2013 With many thanks for slides & images to: Rik Henson Christian Ruff Sandra
More informationHST.583 Functional Magnetic Resonance Imaging: Data Acquisition and Analysis Fall 2008
MIT OpenCourseWare http://ocw.mit.edu HST.583 Functional Magnetic Resonance Imaging: Data Acquisition and Analysis Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More information(b) empirical power. IV: blinded IV: unblinded Regr: blinded Regr: unblinded α. empirical power
Supplementary Information for: Using instrumental variables to disentangle treatment and placebo effects in blinded and unblinded randomized clinical trials influenced by unmeasured confounders by Elias
More informationProblem set 2: understanding ordinary least squares regressions
Problem set 2: understanding ordinary least squares regressions September 12, 2013 1 Introduction This problem set is meant to accompany the undergraduate econometrics video series on youtube; covering
More informationHierarchical Bayesian Analyses for Modeling BOLD Time Series Data
Computational Brain & Behavior (2018) 1:184 213 https://doi.org/10.1007/s42113-018-0013-5 Hierarchical Bayesian Analyses for Modeling BOLD Time Series Data M. Fiona Molloy 1 Giwon Bahg 1 Xiangrui Li 1
More informationGeorge B. Ploubidis. The role of sensitivity analysis in the estimation of causal pathways from observational data. Improving health worldwide
George B. Ploubidis The role of sensitivity analysis in the estimation of causal pathways from observational data Improving health worldwide www.lshtm.ac.uk Outline Sensitivity analysis Causal Mediation
More informationAnalysis of Hearing Loss Data using Correlated Data Analysis Techniques
Analysis of Hearing Loss Data using Correlated Data Analysis Techniques Ruth Penman and Gillian Heller, Department of Statistics, Macquarie University, Sydney, Australia. Correspondence: Ruth Penman, Department
More informationMultilevel IRT for group-level diagnosis. Chanho Park Daniel M. Bolt. University of Wisconsin-Madison
Group-Level Diagnosis 1 N.B. Please do not cite or distribute. Multilevel IRT for group-level diagnosis Chanho Park Daniel M. Bolt University of Wisconsin-Madison Paper presented at the annual meeting
More informationMethods for Addressing Selection Bias in Observational Studies
Methods for Addressing Selection Bias in Observational Studies Susan L. Ettner, Ph.D. Professor Division of General Internal Medicine and Health Services Research, UCLA What is Selection Bias? In the regression
More informationExample 7.2. Autocorrelation. Pilar González and Susan Orbe. Dpt. Applied Economics III (Econometrics and Statistics)
Example 7.2 Autocorrelation Pilar González and Susan Orbe Dpt. Applied Economics III (Econometrics and Statistics) Pilar González and Susan Orbe OCW 2014 Example 7.2. Autocorrelation 1 / 17 Questions.
More informationExperimental design for Cognitive fmri
Experimental design for Cognitive fmri Alexa Morcom Edinburgh SPM course 2017 Thanks to Rik Henson, Thomas Wolbers, Jody Culham, and the SPM authors for slides Overview Categorical designs Factorial designs
More informationCLASSICAL AND. MODERN REGRESSION WITH APPLICATIONS
- CLASSICAL AND. MODERN REGRESSION WITH APPLICATIONS SECOND EDITION Raymond H. Myers Virginia Polytechnic Institute and State university 1 ~l~~l~l~~~~~~~l!~ ~~~~~l~/ll~~ Donated by Duxbury o Thomson Learning,,
More informationROC Analysis of Statistical Methods Used in Functional MRI: Individual Subjects
NeuroImage 9, 311 329 (1999) Article ID nimg.1999.0402, available online at http://www.idealibrary.com on ROC Analysis of Statistical Methods Used in Functional MRI: Individual Subjects Pawel Skudlarski,
More informationGroup-Wise FMRI Activation Detection on Corresponding Cortical Landmarks
Group-Wise FMRI Activation Detection on Corresponding Cortical Landmarks Jinglei Lv 1,2, Dajiang Zhu 2, Xintao Hu 1, Xin Zhang 1,2, Tuo Zhang 1,2, Junwei Han 1, Lei Guo 1,2, and Tianming Liu 2 1 School
More informationMS&E 226: Small Data
MS&E 226: Small Data Lecture 10: Introduction to inference (v2) Ramesh Johari ramesh.johari@stanford.edu 1 / 17 What is inference? 2 / 17 Where did our data come from? Recall our sample is: Y, the vector
More informationfmri Acquisition: Temporal Effects
Functional MRI Data Acquisition: Temporal fmri Acquisition: Temporal Effects Session length Repetition time Fixed vs. distributed temporal sampling Sparse temporal sampling Noise source recording Prospective
More informationBayesian hierarchical modelling
Bayesian hierarchical modelling Matthew Schofield Department of Mathematics and Statistics, University of Otago Bayesian hierarchical modelling Slide 1 What is a statistical model? A statistical model:
More informationCausality from fmri?
Causality from fmri? Olivier David, PhD Brain Function and Neuromodulation, Joseph Fourier University Olivier.David@inserm.fr Grenoble Brain Connectivity Course Yes! Experiments (from 2003 on) Friston
More informationTutorial 3: MANOVA. Pekka Malo 30E00500 Quantitative Empirical Research Spring 2016
Tutorial 3: Pekka Malo 30E00500 Quantitative Empirical Research Spring 2016 Step 1: Research design Adequacy of sample size Choice of dependent variables Choice of independent variables (treatment effects)
More informationMendelian Randomization
Mendelian Randomization Drawback with observational studies Risk factor X Y Outcome Risk factor X? Y Outcome C (Unobserved) Confounders The power of genetics Intermediate phenotype (risk factor) Genetic
More informationEPSE 594: Meta-Analysis: Quantitative Research Synthesis
EPSE 594: Meta-Analysis: Quantitative Research Synthesis Ed Kroc University of British Columbia ed.kroc@ubc.ca March 28, 2019 Ed Kroc (UBC) EPSE 594 March 28, 2019 1 / 32 Last Time Publication bias Funnel
More informationST440/550: Applied Bayesian Statistics. (10) Frequentist Properties of Bayesian Methods
(10) Frequentist Properties of Bayesian Methods Calibrated Bayes So far we have discussed Bayesian methods as being separate from the frequentist approach However, in many cases methods with frequentist
More informationMarno Verbeek Erasmus University, the Netherlands. Cons. Pros
Marno Verbeek Erasmus University, the Netherlands Using linear regression to establish empirical relationships Linear regression is a powerful tool for estimating the relationship between one variable
More informationQuantitative Analysis of Arterial Spin Labeling FMRI Data Using a General Linear Model
Marquette University e-publications@marquette MSCS Faculty Research and Publications Mathematics, Statistics and Computer Science, Department of 9-1-2010 Quantitative Analysis of Arterial Spin Labeling
More informationInternal Model Calibration Using Overlapping Data
Internal Model Calibration Using Overlapping Data Ralph Frankland, James Sharpe, Gaurang Mehta and Rishi Bhatia members of the Extreme Events Working Party 23 November 2017 Agenda Problem Statement Overview
More informationPopulation Inference post Model Selection in Neuroscience
Population Inference post Model Selection in Neuroscience Genevera I. Allen Dobelman Family Junior Chair, Department of Statistics and Electrical and Computer Engineering, Rice University, Department of
More informationHebbian Plasticity for Improving Perceptual Decisions
Hebbian Plasticity for Improving Perceptual Decisions Tsung-Ren Huang Department of Psychology, National Taiwan University trhuang@ntu.edu.tw Abstract Shibata et al. reported that humans could learn to
More informationExperimental Design. Outline. Outline. A very simple experiment. Activation for movement versus rest
Experimental Design Kate Watkins Department of Experimental Psychology University of Oxford With thanks to: Heidi Johansen-Berg Joe Devlin Outline Choices for experimental paradigm Subtraction / hierarchical
More informationA NEW TRIAL DESIGN FULLY INTEGRATING BIOMARKER INFORMATION FOR THE EVALUATION OF TREATMENT-EFFECT MECHANISMS IN PERSONALISED MEDICINE
A NEW TRIAL DESIGN FULLY INTEGRATING BIOMARKER INFORMATION FOR THE EVALUATION OF TREATMENT-EFFECT MECHANISMS IN PERSONALISED MEDICINE Dr Richard Emsley Centre for Biostatistics, Institute of Population
More informationLec 02: Estimation & Hypothesis Testing in Animal Ecology
Lec 02: Estimation & Hypothesis Testing in Animal Ecology Parameter Estimation from Samples Samples We typically observe systems incompletely, i.e., we sample according to a designed protocol. We then
More informationUnit 1 Exploring and Understanding Data
Unit 1 Exploring and Understanding Data Area Principle Bar Chart Boxplot Conditional Distribution Dotplot Empirical Rule Five Number Summary Frequency Distribution Frequency Polygon Histogram Interquartile
More informationPERFUSION MRI CONTRAST BASED TECHNIQUES
PERFUSION MRI CONTRAST BASED TECHNIQUES by Kenny K Israni Mar 28, 2006 PERFUSION - MRI Dynamic Susceptibility contrast Dynamic Relaxivity contrast STEADY-STATE STATE TECHNIQUES Steady-state Susceptibility
More informationChanging expectations about speed alters perceived motion direction
Current Biology, in press Supplemental Information: Changing expectations about speed alters perceived motion direction Grigorios Sotiropoulos, Aaron R. Seitz, and Peggy Seriès Supplemental Data Detailed
More informationChallenges in Developing Learning Algorithms to Personalize mhealth Treatments
Challenges in Developing Learning Algorithms to Personalize mhealth Treatments JOOLHEALTH Bar-Fit Susan A Murphy 01.16.18 HeartSteps SARA Sense 2 Stop Continually Learning Mobile Health Intervention 1)
More informationSum of Neurally Distinct Stimulus- and Task-Related Components.
SUPPLEMENTARY MATERIAL for Cardoso et al. 22 The Neuroimaging Signal is a Linear Sum of Neurally Distinct Stimulus- and Task-Related Components. : Appendix: Homogeneous Linear ( Null ) and Modified Linear
More informationSupervised Learning Models of fmri Data for Inferring Brain Function and Predicting Behavior. Yunzhi Wang
Supervised Learning Models of fmri Data for Inferring Brain Function and Predicting Behavior by Yunzhi Wang A thesis submitted to the Graduate Faculty of Auburn University in partial fulfillment of the
More informationIs model fitting necessary for model-based fmri?
Is model fitting necessary for model-based fmri? Robert C. Wilson Princeton Neuroscience Institute Princeton University Princeton, NJ 85 rcw@princeton.edu Yael Niv Princeton Neuroscience Institute Princeton
More informationIntroduction to Computational Neuroscience
Introduction to Computational Neuroscience Lecture 10: Brain-Computer Interfaces Ilya Kuzovkin So Far Stimulus So Far So Far Stimulus What are the neuroimaging techniques you know about? Stimulus So Far
More informationOutlier Analysis. Lijun Zhang
Outlier Analysis Lijun Zhang zlj@nju.edu.cn http://cs.nju.edu.cn/zlj Outline Introduction Extreme Value Analysis Probabilistic Models Clustering for Outlier Detection Distance-Based Outlier Detection Density-Based
More informationPower-Based Connectivity. JL Sanguinetti
Power-Based Connectivity JL Sanguinetti Power-based connectivity Correlating time-frequency power between two electrodes across time or over trials Gives you flexibility for analysis: Test specific hypotheses
More informationMultiple Regression. James H. Steiger. Department of Psychology and Human Development Vanderbilt University
Multiple Regression James H. Steiger Department of Psychology and Human Development Vanderbilt University James H. Steiger (Vanderbilt University) Multiple Regression 1 / 19 Multiple Regression 1 The Multiple
More informationBayesian Inference. Thomas Nichols. With thanks Lee Harrison
Bayesian Inference Thomas Nichols With thanks Lee Harrison Attention to Motion Paradigm Results Attention No attention Büchel & Friston 1997, Cereb. Cortex Büchel et al. 1998, Brain - fixation only -
More informationOn Regression Analysis Using Bivariate Extreme Ranked Set Sampling
On Regression Analysis Using Bivariate Extreme Ranked Set Sampling Atsu S. S. Dorvlo atsu@squ.edu.om Walid Abu-Dayyeh walidan@squ.edu.om Obaid Alsaidy obaidalsaidy@gmail.com Abstract- Many forms of ranked
More informationMultilevel analysis quantifies variation in the experimental effect while optimizing power and preventing false positives
DOI 10.1186/s12868-015-0228-5 BMC Neuroscience RESEARCH ARTICLE Open Access Multilevel analysis quantifies variation in the experimental effect while optimizing power and preventing false positives Emmeke
More informationMethods for meta-analysis of individual participant data from Mendelian randomization studies with binary outcomes
Methods for meta-analysis of individual participant data from Mendelian randomization studies with binary outcomes Stephen Burgess Simon G. Thompson CRP CHD Genetics Collaboration May 24, 2012 Abstract
More informationSelected Topics in Biostatistics Seminar Series. Missing Data. Sponsored by: Center For Clinical Investigation and Cleveland CTSC
Selected Topics in Biostatistics Seminar Series Missing Data Sponsored by: Center For Clinical Investigation and Cleveland CTSC Brian Schmotzer, MS Biostatistician, CCI Statistical Sciences Core brian.schmotzer@case.edu
More informationPopulation-adjusted treatment comparisons Overview and recommendations from the NICE Decision Support Unit
Population-adjusted treatment comparisons Overview and recommendations from the NICE Decision Support Unit David M Phillippo, University of Bristol 3 Available from www.nicedsu.org.uk Outline Background
More informationHierarchical Bayesian Modeling of Individual Differences in Texture Discrimination
Hierarchical Bayesian Modeling of Individual Differences in Texture Discrimination Timothy N. Rubin (trubin@uci.edu) Michael D. Lee (mdlee@uci.edu) Charles F. Chubb (cchubb@uci.edu) Department of Cognitive
More informationPCA Enhanced Kalman Filter for ECG Denoising
IOSR Journal of Electronics & Communication Engineering (IOSR-JECE) ISSN(e) : 2278-1684 ISSN(p) : 2320-334X, PP 06-13 www.iosrjournals.org PCA Enhanced Kalman Filter for ECG Denoising Febina Ikbal 1, Prof.M.Mathurakani
More informationMSc Neuroimaging for Clinical & Cognitive Neuroscience
MSc Neuroimaging for Clinical & Cognitive Neuroscience School of Psychological Sciences Faculty of Medical & Human Sciences Module Information *Please note that this is a sample guide to modules. The exact
More informationNeuroImage 44 (2009) Contents lists available at ScienceDirect. NeuroImage. journal homepage:
NeuroImage 44 (2009) 893 905 Contents lists available at ScienceDirect NeuroImage journal homepage: www.elsevier.com/locate/ynimg The impact of global signal regression on resting state correlations: Are
More informationSensory Cue Integration
Sensory Cue Integration Summary by Byoung-Hee Kim Computer Science and Engineering (CSE) http://bi.snu.ac.kr/ Presentation Guideline Quiz on the gist of the chapter (5 min) Presenters: prepare one main
More informationFunctional connectivity in fmri
Functional connectivity in fmri Cyril Pernet, PhD Language and Categorization Laboratory, Brain Research Imaging Centre, University of Edinburgh Studying networks fmri can be used for studying both, functional
More informationBayesian and Frequentist Approaches
Bayesian and Frequentist Approaches G. Jogesh Babu Penn State University http://sites.stat.psu.edu/ babu http://astrostatistics.psu.edu All models are wrong But some are useful George E. P. Box (son-in-law
More informationExperimental design. Alexa Morcom Edinburgh SPM course Thanks to Rik Henson, Thomas Wolbers, Jody Culham, and the SPM authors for slides
Experimental design Alexa Morcom Edinburgh SPM course 2013 Thanks to Rik Henson, Thomas Wolbers, Jody Culham, and the SPM authors for slides Overview of SPM Image Image timeseries timeseries Design Design
More informationCarrying out an Empirical Project
Carrying out an Empirical Project Empirical Analysis & Style Hint Special program: Pre-training 1 Carrying out an Empirical Project 1. Posing a Question 2. Literature Review 3. Data Collection 4. Econometric
More informationGraphical Modeling Approaches for Estimating Brain Networks
Graphical Modeling Approaches for Estimating Brain Networks BIOS 516 Suprateek Kundu Department of Biostatistics Emory University. September 28, 2017 Introduction My research focuses on understanding how
More informationDaniel Bulte. Centre for Functional Magnetic Resonance Imaging of the Brain. University of Oxford
Daniel Bulte Centre for Functional Magnetic Resonance Imaging of the Brain University of Oxford Overview Signal Sources BOLD Contrast Mechanism of MR signal change FMRI Modelling Scan design details Factors
More informationHST 583 fmri DATA ANALYSIS AND ACQUISITION
HST 583 fmri DATA ANALYSIS AND ACQUISITION Neural Signal Processing for Functional Neuroimaging Neuroscience Statistics Research Laboratory Massachusetts General Hospital Harvard Medical School/MIT Division
More informationPsychology, 2010, 1: doi: /psych Published Online August 2010 (
Psychology, 2010, 1: 194-198 doi:10.4236/psych.2010.13026 Published Online August 2010 (http://www.scirp.org/journal/psych) Using Generalizability Theory to Evaluate the Applicability of a Serial Bayes
More informationPRINCIPLES OF EFFECTIVE MACHINE LEARNING APPLICATIONS IN REAL-WORLD EVIDENCE
PRINCIPLES OF EFFECTIVE MACHINE LEARNING APPLICATIONS IN REAL-WORLD EVIDENCE Prepared and Presented by: Gorana Capkun-Niggli, PhD, Global Head of Innovation, Health Economics and Outcomes Research, Novartis,
More informationThe SAGE Encyclopedia of Educational Research, Measurement, and Evaluation Multivariate Analysis of Variance
The SAGE Encyclopedia of Educational Research, Measurement, Multivariate Analysis of Variance Contributors: David W. Stockburger Edited by: Bruce B. Frey Book Title: Chapter Title: "Multivariate Analysis
More informationA Bayesian Perspective on Unmeasured Confounding in Large Administrative Databases
A Bayesian Perspective on Unmeasured Confounding in Large Administrative Databases Lawrence McCandless lmccandl@sfu.ca Faculty of Health Sciences, Simon Fraser University, Vancouver Canada Summer 2014
More informationThe Relative Performance of Full Information Maximum Likelihood Estimation for Missing Data in Structural Equation Models
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Educational Psychology Papers and Publications Educational Psychology, Department of 7-1-2001 The Relative Performance of
More informationExperimental Design and the Relative Sensitivity of BOLD and Perfusion fmri
NeuroImage 15, 488 500 (2002) doi:10.1006/nimg.2001.0990, available online at http://www.idealibrary.com on Experimental Design and the Relative Sensitivity of BOLD and Perfusion fmri G. K. Aguirre,*,1
More informationIntroduction to Multilevel Models for Longitudinal and Repeated Measures Data
Introduction to Multilevel Models for Longitudinal and Repeated Measures Data Today s Class: Features of longitudinal data Features of longitudinal models What can MLM do for you? What to expect in this
More informationNPTEL Project. Econometric Modelling. Module 14: Heteroscedasticity Problem. Module 16: Heteroscedasticity Problem. Vinod Gupta School of Management
1 P age NPTEL Project Econometric Modelling Vinod Gupta School of Management Module 14: Heteroscedasticity Problem Module 16: Heteroscedasticity Problem Rudra P. Pradhan Vinod Gupta School of Management
More informationEC352 Econometric Methods: Week 07
EC352 Econometric Methods: Week 07 Gordon Kemp Department of Economics, University of Essex 1 / 25 Outline Panel Data (continued) Random Eects Estimation and Clustering Dynamic Models Validity & Threats
More informationIdentifying Peer Influence Effects in Observational Social Network Data: An Evaluation of Propensity Score Methods
Identifying Peer Influence Effects in Observational Social Network Data: An Evaluation of Propensity Score Methods Dean Eckles Department of Communication Stanford University dean@deaneckles.com Abstract
More informationIntroduction to Observational Studies. Jane Pinelis
Introduction to Observational Studies Jane Pinelis 22 March 2018 Outline Motivating example Observational studies vs. randomized experiments Observational studies: basics Some adjustment strategies Matching
More informationChapter 1: Exploring Data
Chapter 1: Exploring Data Key Vocabulary:! individual! variable! frequency table! relative frequency table! distribution! pie chart! bar graph! two-way table! marginal distributions! conditional distributions!
More informationMaria-Athina Altzerinakou1, Xavier Paoletti2. 9 May, 2017
An adaptive design for the identification of the optimal dose using joint modelling of efficacy and toxicity in phase I/II clinical trials of molecularly targeted agents Maria-Athina Altzerinakou1, Xavier
More informationOn Test Scores (Part 2) How to Properly Use Test Scores in Secondary Analyses. Structural Equation Modeling Lecture #12 April 29, 2015
On Test Scores (Part 2) How to Properly Use Test Scores in Secondary Analyses Structural Equation Modeling Lecture #12 April 29, 2015 PRE 906, SEM: On Test Scores #2--The Proper Use of Scores Today s Class:
More information