Disease mapping with Gaussian processes
|
|
- Rodger Ray
- 5 years ago
- Views:
Transcription
1 EUROHEIS2 Kuopio, Finland August 2010 Aki Vehtari (former Helsinki University of Technology) Department of Biomedical Engineering and Computational Science (BECS)
2 Acknowledgments Researchers - Jarno Vanhatalo - Jouni Hartikainen - Ville Pietiläinen - Jaakko Riihimäki - Simo Särkkä Collaborators and data - The National Institute for Health and Welfare - Finnish Cancer Registry - Statistics Finland
3 Example: alcohol related diseases in Finland Collaboration: The National Institute for Health and Welfare Data: Statistics Finland About inhabited 5km 5km cells in Finland - many cells with no inhabited neighbors In about died due to alcohol diseases - expected death count less than one per cell
4 Outline Related methods and previous work What is a Gaussian process (GP) Flexibility of GP covariance functions Inference Non-Gaussian observation models Computational challenges and some solutions Examples - spatial variation of alcohol-related diseases - spatio-temporal prediction of cancer incidences
5 Spatio-temporal smoothing Assuming nearby-regions have similar latent effect Popular choices - latent Gaussian Markov random field models (e.g. CAR, BYM) - latent Gaussian process models
6 Previous work in GP Kriging since 1950 s Bayesian approach for GPs since 1970 s (Empirical) Bayesian GPs in spatial epidemiology since 1980 s Popularity growing due to algorithm development and increase in the computer speed
7 Gaussian process Infinite-dimensional generalization of multivariate Gaussian Can be used to set a prior on function space Univariate Gaussian Bivariate Gaussian Samples from GP f N(µ, σ 2 ) f N(µ, Σ) f GP(µ, Σ)
8 Spatial smoothing Observations Y X y N(f, σ 2 )
9 Spatial smoothing Observations Y X1 X2 X3 X 1) Uncertainty related to f
10 Spatial smoothing Observations Y X1 X2 X3 X 1) Uncertainty related to f 2) Uncertainty related to y
11 Spatial smoothing Observations Y X1 X2 X3 X Let s focus on f f values given similar x values seem to correlate
12 Spatial smoothing Observations Covariance k(x 2,X j ) Y X1 X2 X3 X Describe correlation with covariance function k(x i, x j )
13 Spatial smoothing Observations Covariance k(x 2,X j ) Y X1 X2 X3 X Describe correlation with covariance function k(x i, x j ) e.g. k(x i, x j ) = φ exp ( (x i x j ) 2 l 2 )
14 Spatial smoothing Observations Covariance k(x 2,X j ) Y F(X 1 ) X1 X2 X3 X F(X 2 ) Function values given X 1 and X 2 are far in X space and do not correlate
15 Spatial smoothing Observations Covariance k(x 2,X j ) Y F(X 3 ) X1 X2 X3 X F(X 2 ) Function values given similar X 2 and X 3 correlate
16 Gaussian process Joint distribution for the observations y and a new f [ ] ( [ y K (X, X) + σ N 0, 2 ]) I K (X, X ) K (X, X) K (X, X ) f
17 Gaussian process Joint distribution for the observations y and a new f [ ] ( [ y K (X, X) + σ N 0, 2 ]) I K (X, X ) K (X, X) K (X, X ) f Since a finite number of random variables, computations as for a multivariate Gaussian E[f ] = K (X, X)[K (X, X) + σ 2 I] 1 y cov[f ] = K (X, X ) K (X, X)[K (X, X) + σ 2 I] 1 K (X, X )
18 Gaussian process regression Observations f=gp(x,x obs,y obs,θ) Uncertainty in f Uncertainty in y Y X E[f ] = K (X, X)[K (X, X) + σ 2 I] 1 y cov[f ] = K (X, X ) K (X, X)[K (X, X) + σ 2 I] 1 K (X, X ) Note: linear equations but nonlinear regression due to weighting by covariance function
19 Flexibility of GP covariance functions Squared exponential Matern(5/2) Exponential Two lengthscales
20 Spatial smoothing Correlation defined, e.g., using Euclidian distance Spatial correlation decreases as distance increases
21 Flexibility of GP covariance functions Examples of spatial surfaces sampled from different GP priors
22 Flexibility of GP covariance functions Non-linearities and interactions automatically for large number of covariates 0.75 readmission probability days of preceding care age
23 Computational challenges Matrix inversion O(n 3 ) Non-Gaussian observation models Estimation of the covariance function parameters
24 Inference given a covariance matrix Given a Gaussian observation model, a covariance function K and its parameters θ: analytic solution with matrix algebra E[f θ] = K (X, X)[K (X, X) + σ 2 I] 1 y cov[f θ] = K (X, X ) K (X, X)[K (X, X) + σ 2 I] 1 K (X, X ) Requires the inversion of a n n matrix, which is generally an O(n 3 ) operation
25 Covariance matrix approximations Matrix inversion O(n 3 ) With reduced rank or projection approximations O(m 2 n) - m n - useful for global spatial phenomena Compact support covariance functions - naturally sparse covariance matrix if the correlation length is short - useful for local spatial phenomena Combined approach for good global and local performance
26 Non-Gaussian observation models With a latent variable approach E.g. count observations y Poisson(α exp(f )) f Uncertainty in f 20 Observations µ=α exp(f) Uncertainty in µ F 0 Y 10 X 0 X
27 Non-Gaussian observation models With a latent variable approach E.g. count observations y Poisson(α exp(f )) f Uncertainty in f 20 Observations µ=α exp(f) Uncertainty in µ Uncertainty in y F 0 Y 10 X 0 X
28 Inference given a non-gaussian observation model Integration over the latent variables f i required No analytic solution - Monte Carlo methods have mixing and convergence problems when n is large New coming of analytic approximations - Laplace approximation ie. a Gaussian at the mode - expectation propagation (EP) ie. a Gaussian which minimizes KL-divergence from the true distribution - variational Bayes ie. a Gaussian which minimizes KL-divergence from the approximation
29 Integration over the latent space In our spatial and spatio-temporal experiments - Expectation propagation about as good as MCMC, but much faster - Laplace approximation almost as good as EP, but faster In others experiments - variational not as good as EP
30 Inference given unknown covariance function parameters θ Estimation of θ (empirical Bayes) - optimize θ given the marginal likelihood p(d θ) (from the properties of a multivariate Gaussian) - nonlinear optimization wrt. θ
31 Inference given unknown covariance function parameters θ Estimation of θ (empirical Bayes) - optimize θ given the marginal likelihood p(d θ) (from the properties of a multivariate Gaussian) - nonlinear optimization wrt. θ Full Bayes - integrate over the posterior of θ p(θ D) p(d θ)p(θ) - no analytic solution Monte Carlo, grid, or CCD approximation reasonable performance when latent values have been marginalized with analytic approximation
32 Integration over the covariance function parameters θ In our experiments - with large n optimization of θ as good as integration integration over the latent space more important If integration needed - if computation of the marginal density is slow, the integration need to be made with as few density evaluations as possible - eg. grid based when few hyperparameters
33 Computational solutions Matrix inversion O(n 3 ) sparse approximations O(m 2 n) - reduced rank approximations compact support covariance functions - naturally sparse covariance matrix Integration over n dimensional latent space MCMC, Laplace, Expectation propagation, Variational Integration over hyperparameters Monte Carlo, grid evaluation, etc.
34 Improvement in computational speed Basic GP with MCMC vs. speeded-up GP 2006: resolution 20km (900 cells), analysis 2-3 days 2009: resolution 5km ( cells), analysis 20 min
35 Example: alcohol related diseases in Finland Collaboration: The National Institute for Health and Welfare Data: Statistics Finland About inhabited 5km 5km cells in Finland - many cells with no inhabited neighbors In about died due to alcohol diseases - expected death count less than one per cell
36 Example: alcohol related diseases in Finland Sex-age-education standardized expected death counts used to compute the raw risk Risk smoothed using GP with long and short length scale and negative-binomial observation model
37 Example: alcohol related diseases in Finland Sex-age-education standardized expected death counts used to compute the raw risk Risk smoothed using GP with long and short length scale and negative-binomial observation model Is the relative risk higher in the population centers?
38 Example: alcohol related diseases in Finland The smoothed relative risks vs. the population density
39 Example: alcohol related diseases in Finland The smoothed relative risks vs. the population density Add population density as explanatory variable
40 Example: alcohol related diseases in Finland Population density and spatial variation explain the variation in the risk Population density effect Spatial effect
41 Example: alcohol related diseases in Finland Spatial and population density explain the variation in the risk Combined effect
42 Example: alcohol related diseases in Finland With Google Maps With Google Earth
43 Example: alcohol related diseases in Finland Zoomed near Kuopio
44 Example: alcohol related diseases in Finland Jarno Vanhatalo, Ville Pietiläinen and Aki Vehtari (2010). Approximate inference for disease mapping with sparse Gaussian processes. Statistics in Medicine, 29(15): Jarno Vanhatalo, Pia Mäkelä, ja Aki Vehtari (2010). Alkoholikuolleisuuden alueelliset erot Suomessa 2000-luvun alussa. Yhteiskuntapolitiikka, 75(3): Online supplement research/mm/finnwell/alcoholfinland/ GPstuff toolbox In Finnish Newspapers: Helsingin Sanomat, Pohjalainen, Lapin Kansa
45 Compared to CAR In spatial epidemiology CAR (or BYM) is most used model Correlation defined conditionally based on a neighborhood structure discrete definition + major computational speed-up if a precision matrix is sparse due to small neighborhoods - describes only local correlation - neighborhood definition may be difficult for irregularly spaced data and high dimensional data
46 Example: alcohol related diseases in Finland Comparison to CAR Compared to CAR computed with INLA software (uses fast Integrated Nested Laplace Approximation) - CAR model lacks long range correlation part - CAR model has much higher variance, especially for cells having no or few inhabited neighbors - GP has a better predictive performance Variance CAR CAR Mean GP GP
47 Example: lung cancer women County incidences and background population for years years, 431 counties observations Data: Finnish Cancer Registry Model: GP with temporal + spatial + spatiotemporal component
48 Example: lung cancer women (a) Temporal (b) Spatial
49 Example: lung cancer women (a) 1953 (d) 1983 (b) 1963 (c) (e) 1993 (f) 2003
50 Work in progress Spatio-temporal GPs can be written as linear stochastic partial differential equations Reduces computational complexity from O(n 3 T 3 ) to O(n 3 T ), i.e. method scales linearily in T n limited as for spatial GP - few thousand with no sparse aproximations - more than ten thousand with sparse approximations Current tests with over million spatio-temporal points Makes it easier to specify non-stationary temporal dynamics, which are necessary, for example, when performing future predictions
51 Example: skin melanoma Prediction with a non-stationary spatio-temporal GP Ten year predictions (mean + 95 % CI) and observed counts in one city y Time
52 Example: skin melanoma Prediction with a non-stationary spatio-temporal GP Ten year predictions (mean + 95 % CI) and observed counts in one city y Time Age-period-cohort model in progress
53 Other applications GPs are not just for spatial statistics Length of stay in rehabilitation after hip fracture surgery Probability of readmission Clustering users of health care services Genome wide analysis Marked point processes Industrial processes
Disease mapping with Gaussian processes
Liverpool, UK, 4 5 November 3 Aki Vehtari Department of Biomedical Engineering and Computational Science (BECS) Outline Example: Alcohol related deaths in Finland Spatial priors and benefits of GP prior
More informationLecture: Gaussian Process Regression. STAT 6474 Instructor: Hongxiao Zhu
Lecture: Gaussian Process Regression STAT 6474 Instructor: Hongxiao Zhu Motivation Reference: Marc Deisenroth s tutorial on Robot Learning. 2 Fast Learning for Autonomous Robots with Gaussian Processes
More informationSummary STK 4150/9150
STK4150 - Intro 1 Summary STK 4150/9150 Odd Kolbjørnsen May 22 2017 Scope You are expected to know and be able to use basic concepts introduced in the book. You knowledge is expected to be larger than
More informationBayesian Hierarchical Models
Bayesian Hierarchical Models Gavin Shaddick, Millie Green, Matthew Thomas University of Bath 6 th - 9 th December 2016 1/ 34 APPLICATIONS OF BAYESIAN HIERARCHICAL MODELS 2/ 34 OUTLINE Spatial epidemiology
More informationA short introduction to INLA and R-INLA
A short introduction to INLA and R-INLA Integrated Nested Laplace Approximation Thomas Opitz, BioSP, INRA Avignon Workshop: Theory and practice of INLA and SPDE November 7, 2018 2/21 Plan for this talk
More informationGAUSSIAN PROCESS REGRESSION
GAUSSIAN PROCESS REGRESSION CSE 515T Spring 2015 1. BACKGROUND The kernel trick again... The Kernel Trick Consider again the linear regression model: y(x) = φ(x) w + ε, with prior p(w) = N (w; 0, Σ). The
More informationChapter 4 - Fundamentals of spatial processes Lecture notes
Chapter 4 - Fundamentals of spatial processes Lecture notes Geir Storvik January 21, 2013 STK4150 - Intro 2 Spatial processes Typically correlation between nearby sites Mostly positive correlation Negative
More informationIntegrated Non-Factorized Variational Inference
Integrated Non-Factorized Variational Inference Shaobo Han, Xuejun Liao and Lawrence Carin Duke University February 27, 2014 S. Han et al. Integrated Non-Factorized Variational Inference February 27, 2014
More informationNon-Gaussian likelihoods for Gaussian Processes
Non-Gaussian likelihoods for Gaussian Processes Alan Saul University of Sheffield Outline Motivation Laplace approximation KL method Expectation Propagation Comparing approximations GP regression Model
More informationBeyond MCMC in fitting complex Bayesian models: The INLA method
Beyond MCMC in fitting complex Bayesian models: The INLA method Valeska Andreozzi Centre of Statistics and Applications of Lisbon University (valeska.andreozzi at fc.ul.pt) European Congress of Epidemiology
More informationIntroduction to Gaussian Processes
Introduction to Gaussian Processes Iain Murray murray@cs.toronto.edu CSC255, Introduction to Machine Learning, Fall 28 Dept. Computer Science, University of Toronto The problem Learn scalar function of
More informationModel Selection for Gaussian Processes
Institute for Adaptive and Neural Computation School of Informatics,, UK December 26 Outline GP basics Model selection: covariance functions and parameterizations Criteria for model selection Marginal
More informationGaussian Processes. Le Song. Machine Learning II: Advanced Topics CSE 8803ML, Spring 2012
Gaussian Processes Le Song Machine Learning II: Advanced Topics CSE 8803ML, Spring 01 Pictorial view of embedding distribution Transform the entire distribution to expected features Feature space Feature
More informationStatistical Ecology with Gaussian Processes
Statistical Ecology with Gaussian Processes 5.6.2017 Machine Learning Coffee Seminar, University of Helsinki Jarno Vanhatalo Department of Mathematics and Statistics and Department of Biosciences University
More informationDensity Estimation. Seungjin Choi
Density Estimation Seungjin Choi Department of Computer Science and Engineering Pohang University of Science and Technology 77 Cheongam-ro, Nam-gu, Pohang 37673, Korea seungjin@postech.ac.kr http://mlg.postech.ac.kr/
More informationBayesian Regression Linear and Logistic Regression
When we want more than point estimates Bayesian Regression Linear and Logistic Regression Nicole Beckage Ordinary Least Squares Regression and Lasso Regression return only point estimates But what if we
More informationCPSC 540: Machine Learning
CPSC 540: Machine Learning MCMC and Non-Parametric Bayes Mark Schmidt University of British Columbia Winter 2016 Admin I went through project proposals: Some of you got a message on Piazza. No news is
More informationPart 1: Expectation Propagation
Chalmers Machine Learning Summer School Approximate message passing and biomedicine Part 1: Expectation Propagation Tom Heskes Machine Learning Group, Institute for Computing and Information Sciences Radboud
More informationNonparametric Bayesian Methods (Gaussian Processes)
[70240413 Statistical Machine Learning, Spring, 2015] Nonparametric Bayesian Methods (Gaussian Processes) Jun Zhu dcszj@mail.tsinghua.edu.cn http://bigml.cs.tsinghua.edu.cn/~jun State Key Lab of Intelligent
More informationMarkov Chain Monte Carlo Algorithms for Gaussian Processes
Markov Chain Monte Carlo Algorithms for Gaussian Processes Michalis K. Titsias, Neil Lawrence and Magnus Rattray School of Computer Science University of Manchester June 8 Outline Gaussian Processes Sampling
More informationBayesian Learning. HT2015: SC4 Statistical Data Mining and Machine Learning. Maximum Likelihood Principle. The Bayesian Learning Framework
HT5: SC4 Statistical Data Mining and Machine Learning Dino Sejdinovic Department of Statistics Oxford http://www.stats.ox.ac.uk/~sejdinov/sdmml.html Maximum Likelihood Principle A generative model for
More informationExpectation Propagation for Approximate Bayesian Inference
Expectation Propagation for Approximate Bayesian Inference José Miguel Hernández Lobato Universidad Autónoma de Madrid, Computer Science Department February 5, 2007 1/ 24 Bayesian Inference Inference Given
More informationSTA 4273H: Statistical Machine Learning
STA 4273H: Statistical Machine Learning Russ Salakhutdinov Department of Statistics! rsalakhu@utstat.toronto.edu! http://www.utstat.utoronto.ca/~rsalakhu/ Sidney Smith Hall, Room 6002 Lecture 7 Approximate
More informationBayesian SAE using Complex Survey Data Lecture 4A: Hierarchical Spatial Bayes Modeling
Bayesian SAE using Complex Survey Data Lecture 4A: Hierarchical Spatial Bayes Modeling Jon Wakefield Departments of Statistics and Biostatistics University of Washington 1 / 37 Lecture Content Motivation
More informationLinear regression example Simple linear regression: f(x) = ϕ(x)t w w ~ N(0, ) The mean and covariance are given by E[f(x)] = ϕ(x)e[w] = 0.
Gaussian Processes Gaussian Process Stochastic process: basically, a set of random variables. may be infinite. usually related in some way. Gaussian process: each variable has a Gaussian distribution every
More informationGaussian Process Regression Model in Spatial Logistic Regression
Journal of Physics: Conference Series PAPER OPEN ACCESS Gaussian Process Regression Model in Spatial Logistic Regression To cite this article: A Sofro and A Oktaviarina 018 J. Phys.: Conf. Ser. 947 01005
More informationA general mixed model approach for spatio-temporal regression data
A general mixed model approach for spatio-temporal regression data Thomas Kneib, Ludwig Fahrmeir & Stefan Lang Department of Statistics, Ludwig-Maximilians-University Munich 1. Spatio-temporal regression
More informationState Space Representation of Gaussian Processes
State Space Representation of Gaussian Processes Simo Särkkä Department of Biomedical Engineering and Computational Science (BECS) Aalto University, Espoo, Finland June 12th, 2013 Simo Särkkä (Aalto University)
More informationProbabilistic Graphical Models Lecture 20: Gaussian Processes
Probabilistic Graphical Models Lecture 20: Gaussian Processes Andrew Gordon Wilson www.cs.cmu.edu/~andrewgw Carnegie Mellon University March 30, 2015 1 / 53 What is Machine Learning? Machine learning algorithms
More informationCSci 8980: Advanced Topics in Graphical Models Gaussian Processes
CSci 8980: Advanced Topics in Graphical Models Gaussian Processes Instructor: Arindam Banerjee November 15, 2007 Gaussian Processes Outline Gaussian Processes Outline Parametric Bayesian Regression Gaussian
More informationModelling geoadditive survival data
Modelling geoadditive survival data Thomas Kneib & Ludwig Fahrmeir Department of Statistics, Ludwig-Maximilians-University Munich 1. Leukemia survival data 2. Structured hazard regression 3. Mixed model
More informationAnalysing geoadditive regression data: a mixed model approach
Analysing geoadditive regression data: a mixed model approach Institut für Statistik, Ludwig-Maximilians-Universität München Joint work with Ludwig Fahrmeir & Stefan Lang 25.11.2005 Spatio-temporal regression
More informationSTA414/2104 Statistical Methods for Machine Learning II
STA414/2104 Statistical Methods for Machine Learning II Murat A. Erdogdu & David Duvenaud Department of Computer Science Department of Statistical Sciences Lecture 3 Slide credits: Russ Salakhutdinov Announcements
More informationGaussian processes for spatial modelling in environmental health: parameterizing for flexibility vs. computational efficiency
Gaussian processes for spatial modelling in environmental health: parameterizing for flexibility vs. computational efficiency Chris Paciorek March 11, 2005 Department of Biostatistics Harvard School of
More informationIntroduction to Gaussian Processes
Introduction to Gaussian Processes Iain Murray School of Informatics, University of Edinburgh The problem Learn scalar function of vector values f(x).5.5 f(x) y i.5.2.4.6.8 x f 5 5.5 x x 2.5 We have (possibly
More informationNonparmeteric Bayes & Gaussian Processes. Baback Moghaddam Machine Learning Group
Nonparmeteric Bayes & Gaussian Processes Baback Moghaddam baback@jpl.nasa.gov Machine Learning Group Outline Bayesian Inference Hierarchical Models Model Selection Parametric vs. Nonparametric Gaussian
More informationLog Gaussian Cox Processes. Chi Group Meeting February 23, 2016
Log Gaussian Cox Processes Chi Group Meeting February 23, 2016 Outline Typical motivating application Introduction to LGCP model Brief overview of inference Applications in my work just getting started
More informationLecture 2: From Linear Regression to Kalman Filter and Beyond
Lecture 2: From Linear Regression to Kalman Filter and Beyond January 18, 2017 Contents 1 Batch and Recursive Estimation 2 Towards Bayesian Filtering 3 Kalman Filter and Bayesian Filtering and Smoothing
More informationPractical Bayesian Optimization of Machine Learning. Learning Algorithms
Practical Bayesian Optimization of Machine Learning Algorithms CS 294 University of California, Berkeley Tuesday, April 20, 2016 Motivation Machine Learning Algorithms (MLA s) have hyperparameters that
More informationGaussian processes. Chuong B. Do (updated by Honglak Lee) November 22, 2008
Gaussian processes Chuong B Do (updated by Honglak Lee) November 22, 2008 Many of the classical machine learning algorithms that we talked about during the first half of this course fit the following pattern:
More information(5) Multi-parameter models - Gibbs sampling. ST440/540: Applied Bayesian Analysis
Summarizing a posterior Given the data and prior the posterior is determined Summarizing the posterior gives parameter estimates, intervals, and hypothesis tests Most of these computations are integrals
More informationNonparameteric Regression:
Nonparameteric Regression: Nadaraya-Watson Kernel Regression & Gaussian Process Regression Seungjin Choi Department of Computer Science and Engineering Pohang University of Science and Technology 77 Cheongam-ro,
More informationLarge-scale Ordinal Collaborative Filtering
Large-scale Ordinal Collaborative Filtering Ulrich Paquet, Blaise Thomson, and Ole Winther Microsoft Research Cambridge, University of Cambridge, Technical University of Denmark ulripa@microsoft.com,brmt2@cam.ac.uk,owi@imm.dtu.dk
More informationApproaches for Multiple Disease Mapping: MCAR and SANOVA
Approaches for Multiple Disease Mapping: MCAR and SANOVA Dipankar Bandyopadhyay Division of Biostatistics, University of Minnesota SPH April 22, 2015 1 Adapted from Sudipto Banerjee s notes SANOVA vs MCAR
More informationRonald Christensen. University of New Mexico. Albuquerque, New Mexico. Wesley Johnson. University of California, Irvine. Irvine, California
Texts in Statistical Science Bayesian Ideas and Data Analysis An Introduction for Scientists and Statisticians Ronald Christensen University of New Mexico Albuquerque, New Mexico Wesley Johnson University
More informationLecture 2: From Linear Regression to Kalman Filter and Beyond
Lecture 2: From Linear Regression to Kalman Filter and Beyond Department of Biomedical Engineering and Computational Science Aalto University January 26, 2012 Contents 1 Batch and Recursive Estimation
More informationCommunity Health Needs Assessment through Spatial Regression Modeling
Community Health Needs Assessment through Spatial Regression Modeling Glen D. Johnson, PhD CUNY School of Public Health glen.johnson@lehman.cuny.edu Objectives: Assess community needs with respect to particular
More informationSTA 4273H: Sta-s-cal Machine Learning
STA 4273H: Sta-s-cal Machine Learning Russ Salakhutdinov Department of Computer Science! Department of Statistical Sciences! rsalakhu@cs.toronto.edu! h0p://www.cs.utoronto.ca/~rsalakhu/ Lecture 2 In our
More informationPumps, Maps and Pea Soup: Spatio-temporal methods in environmental epidemiology
Pumps, Maps and Pea Soup: Spatio-temporal methods in environmental epidemiology Gavin Shaddick Department of Mathematical Sciences University of Bath 2012-13 van Eeden lecture Thanks Constance van Eeden
More informationBAYESIAN METHODS FOR VARIABLE SELECTION WITH APPLICATIONS TO HIGH-DIMENSIONAL DATA
BAYESIAN METHODS FOR VARIABLE SELECTION WITH APPLICATIONS TO HIGH-DIMENSIONAL DATA Intro: Course Outline and Brief Intro to Marina Vannucci Rice University, USA PASI-CIMAT 04/28-30/2010 Marina Vannucci
More informationChapter 4 - Spatial processes R packages and software Lecture notes
TK4150 - Intro 1 Chapter 4 - Spatial processes R packages and software Lecture notes Odd Kolbjørnsen and Geir Storvik February 13, 2017 STK4150 - Intro 2 Last time General set up for Kriging type problems
More informationThe Bayesian approach to inverse problems
The Bayesian approach to inverse problems Youssef Marzouk Department of Aeronautics and Astronautics Center for Computational Engineering Massachusetts Institute of Technology ymarz@mit.edu, http://uqgroup.mit.edu
More informationPrinciples of Bayesian Inference
Principles of Bayesian Inference Sudipto Banerjee and Andrew O. Finley 2 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A. 2 Department of Forestry & Department
More informationPareto Smoothed Importance Sampling
Pareto moothed Importance ampling Aki Vehtari Andrew Gelman Jonah Gabry October 7 Abstract Importance weighting is a convenient general way to adjust for draws from the wrong distribution, but the resulting
More informationBayesian Inference and MCMC
Bayesian Inference and MCMC Aryan Arbabi Partly based on MCMC slides from CSC412 Fall 2018 1 / 18 Bayesian Inference - Motivation Consider we have a data set D = {x 1,..., x n }. E.g each x i can be the
More informationHow to build an automatic statistician
How to build an automatic statistician James Robert Lloyd 1, David Duvenaud 1, Roger Grosse 2, Joshua Tenenbaum 2, Zoubin Ghahramani 1 1: Department of Engineering, University of Cambridge, UK 2: Massachusetts
More informationADVANCED MACHINE LEARNING ADVANCED MACHINE LEARNING. Non-linear regression techniques Part - II
1 Non-linear regression techniques Part - II Regression Algorithms in this Course Support Vector Machine Relevance Vector Machine Support vector regression Boosting random projections Relevance vector
More informationPrinciples of Bayesian Inference
Principles of Bayesian Inference Sudipto Banerjee 1 and Andrew O. Finley 2 1 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A. 2 Department of Forestry & Department
More informationWrapped Gaussian processes: a short review and some new results
Wrapped Gaussian processes: a short review and some new results Giovanna Jona Lasinio 1, Gianluca Mastrantonio 2 and Alan Gelfand 3 1-Università Sapienza di Roma 2- Università RomaTRE 3- Duke University
More informationState Space Gaussian Processes with Non-Gaussian Likelihoods
State Space Gaussian Processes with Non-Gaussian Likelihoods Hannes Nickisch 1 Arno Solin 2 Alexander Grigorievskiy 2,3 1 Philips Research, 2 Aalto University, 3 Silo.AI ICML2018 July 13, 2018 Outline
More informationPattern Recognition and Machine Learning. Bishop Chapter 2: Probability Distributions
Pattern Recognition and Machine Learning Chapter 2: Probability Distributions Cécile Amblard Alex Kläser Jakob Verbeek October 11, 27 Probability Distributions: General Density Estimation: given a finite
More informationA Process over all Stationary Covariance Kernels
A Process over all Stationary Covariance Kernels Andrew Gordon Wilson June 9, 0 Abstract I define a process over all stationary covariance kernels. I show how one might be able to perform inference that
More informationApproximate inference for disease mapping with sparse Gaussian processes
Research Article Received 1 March 9, Accepted 4 February 1 Published online 9 March 1 in Wiley Interscience (www.interscience.wiley.com) DOI: 1.1/sim.3895 Approximate inference for disease mapping with
More informationarxiv: v1 [stat.ml] 22 Apr 2014
Approximate Inference for Nonstationary Heteroscedastic Gaussian process Regression arxiv:1404.5443v1 [stat.ml] 22 Apr 2014 Abstract Ville Tolvanen Department of Biomedical Engineering and Computational
More informationBayesian Modeling with Gaussian Processes using the GPstuff Toolbox
Bayesian Modeling with Gaussian Processes using the GPstuff Toolbox Jarno Vanhatalo Department of Environmental Sciences University of Helsinki P.O. Box 65, FI-00014 Helsinki, Finland Jaakko Riihimäki
More informationPMR Learning as Inference
Outline PMR Learning as Inference Probabilistic Modelling and Reasoning Amos Storkey Modelling 2 The Exponential Family 3 Bayesian Sets School of Informatics, University of Edinburgh Amos Storkey PMR Learning
More informationStatistícal Methods for Spatial Data Analysis
Texts in Statistícal Science Statistícal Methods for Spatial Data Analysis V- Oliver Schabenberger Carol A. Gotway PCT CHAPMAN & K Contents Preface xv 1 Introduction 1 1.1 The Need for Spatial Analysis
More informationPrinciples of Bayesian Inference
Principles of Bayesian Inference Sudipto Banerjee University of Minnesota July 20th, 2008 1 Bayesian Principles Classical statistics: model parameters are fixed and unknown. A Bayesian thinks of parameters
More informationGaussian Process Regression with Censored Data Using Expectation Propagation
Sixth European Workshop on Probabilistic Graphical Models, Granada, Spain, 01 Gaussian Process Regression with Censored Data Using Expectation Propagation Perry Groot, Peter Lucas Radboud University Nijmegen,
More informationBayesian Linear Regression
Bayesian Linear Regression Sudipto Banerjee 1 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A. September 15, 2010 1 Linear regression models: a Bayesian perspective
More informationProbabilistic machine learning group, Aalto University Bayesian theory and methods, approximative integration, model
Aki Vehtari, Aalto University, Finland Probabilistic machine learning group, Aalto University http://research.cs.aalto.fi/pml/ Bayesian theory and methods, approximative integration, model assessment and
More informationLecture 13 Fundamentals of Bayesian Inference
Lecture 13 Fundamentals of Bayesian Inference Dennis Sun Stats 253 August 11, 2014 Outline of Lecture 1 Bayesian Models 2 Modeling Correlations Using Bayes 3 The Universal Algorithm 4 BUGS 5 Wrapping Up
More informationIntroduction to Systems Analysis and Decision Making Prepared by: Jakub Tomczak
Introduction to Systems Analysis and Decision Making Prepared by: Jakub Tomczak 1 Introduction. Random variables During the course we are interested in reasoning about considered phenomenon. In other words,
More information20: Gaussian Processes
10-708: Probabilistic Graphical Models 10-708, Spring 2016 20: Gaussian Processes Lecturer: Andrew Gordon Wilson Scribes: Sai Ganesh Bandiatmakuri 1 Discussion about ML Here we discuss an introduction
More informationEfficient Likelihood-Free Inference
Efficient Likelihood-Free Inference Michael Gutmann http://homepages.inf.ed.ac.uk/mgutmann Institute for Adaptive and Neural Computation School of Informatics, University of Edinburgh 8th November 2017
More informationGaussian Process Regression
Gaussian Process Regression 4F1 Pattern Recognition, 21 Carl Edward Rasmussen Department of Engineering, University of Cambridge November 11th - 16th, 21 Rasmussen (Engineering, Cambridge) Gaussian Process
More informationBuilding an Automatic Statistician
Building an Automatic Statistician Zoubin Ghahramani Department of Engineering University of Cambridge zoubin@eng.cam.ac.uk http://learning.eng.cam.ac.uk/zoubin/ ALT-Discovery Science Conference, October
More informationBayesian Methods for Machine Learning
Bayesian Methods for Machine Learning CS 584: Big Data Analytics Material adapted from Radford Neal s tutorial (http://ftp.cs.utoronto.ca/pub/radford/bayes-tut.pdf), Zoubin Ghahramni (http://hunch.net/~coms-4771/zoubin_ghahramani_bayesian_learning.pdf),
More informationDynamic System Identification using HDMR-Bayesian Technique
Dynamic System Identification using HDMR-Bayesian Technique *Shereena O A 1) and Dr. B N Rao 2) 1), 2) Department of Civil Engineering, IIT Madras, Chennai 600036, Tamil Nadu, India 1) ce14d020@smail.iitm.ac.in
More informationLecture 6: Bayesian Inference in SDE Models
Lecture 6: Bayesian Inference in SDE Models Bayesian Filtering and Smoothing Point of View Simo Särkkä Aalto University Simo Särkkä (Aalto) Lecture 6: Bayesian Inference in SDEs 1 / 45 Contents 1 SDEs
More informationFast approximations for the Expected Value of Partial Perfect Information using R-INLA
Fast approximations for the Expected Value of Partial Perfect Information using R-INLA Anna Heath 1 1 Department of Statistical Science, University College London 22 May 2015 Outline 1 Health Economic
More informationReliability Monitoring Using Log Gaussian Process Regression
COPYRIGHT 013, M. Modarres Reliability Monitoring Using Log Gaussian Process Regression Martin Wayne Mohammad Modarres PSA 013 Center for Risk and Reliability University of Maryland Department of Mechanical
More informationFaculty of Health Sciences. Regression models. Counts, Poisson regression, Lene Theil Skovgaard. Dept. of Biostatistics
Faculty of Health Sciences Regression models Counts, Poisson regression, 27-5-2013 Lene Theil Skovgaard Dept. of Biostatistics 1 / 36 Count outcome PKA & LTS, Sect. 7.2 Poisson regression The Binomial
More informationModels for spatial data (cont d) Types of spatial data. Types of spatial data (cont d) Hierarchical models for spatial data
Hierarchical models for spatial data Based on the book by Banerjee, Carlin and Gelfand Hierarchical Modeling and Analysis for Spatial Data, 2004. We focus on Chapters 1, 2 and 5. Geo-referenced data arise
More informationCreating Non-Gaussian Processes from Gaussian Processes by the Log-Sum-Exp Approach. Radford M. Neal, 28 February 2005
Creating Non-Gaussian Processes from Gaussian Processes by the Log-Sum-Exp Approach Radford M. Neal, 28 February 2005 A Very Brief Review of Gaussian Processes A Gaussian process is a distribution over
More informationFast Likelihood-Free Inference via Bayesian Optimization
Fast Likelihood-Free Inference via Bayesian Optimization Michael Gutmann https://sites.google.com/site/michaelgutmann University of Helsinki Aalto University Helsinki Institute for Information Technology
More informationSTA414/2104. Lecture 11: Gaussian Processes. Department of Statistics
STA414/2104 Lecture 11: Gaussian Processes Department of Statistics www.utstat.utoronto.ca Delivered by Mark Ebden with thanks to Russ Salakhutdinov Outline Gaussian Processes Exam review Course evaluations
More informationIntroduction to Probabilistic Machine Learning
Introduction to Probabilistic Machine Learning Piyush Rai Dept. of CSE, IIT Kanpur (Mini-course 1) Nov 03, 2015 Piyush Rai (IIT Kanpur) Introduction to Probabilistic Machine Learning 1 Machine Learning
More informationSTA 4273H: Statistical Machine Learning
STA 4273H: Statistical Machine Learning Russ Salakhutdinov Department of Statistics! rsalakhu@utstat.toronto.edu! http://www.utstat.utoronto.ca/~rsalakhu/ Sidney Smith Hall, Room 6002 Lecture 3 Linear
More informationGWAS V: Gaussian processes
GWAS V: Gaussian processes Dr. Oliver Stegle Christoh Lippert Prof. Dr. Karsten Borgwardt Max-Planck-Institutes Tübingen, Germany Tübingen Summer 2011 Oliver Stegle GWAS V: Gaussian processes Summer 2011
More informationSTAT 518 Intro Student Presentation
STAT 518 Intro Student Presentation Wen Wei Loh April 11, 2013 Title of paper Radford M. Neal [1999] Bayesian Statistics, 6: 475-501, 1999 What the paper is about Regression and Classification Flexible
More informationData Integration Model for Air Quality: A Hierarchical Approach to the Global Estimation of Exposures to Ambient Air Pollution
Data Integration Model for Air Quality: A Hierarchical Approach to the Global Estimation of Exposures to Ambient Air Pollution Matthew Thomas 9 th January 07 / 0 OUTLINE Introduction Previous methods for
More informationPrinciples of Bayesian Inference
Principles of Bayesian Inference Sudipto Banerjee 1 and Andrew O. Finley 2 1 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A. 2 Department of Forestry & Department
More informationLecture 7 and 8: Markov Chain Monte Carlo
Lecture 7 and 8: Markov Chain Monte Carlo 4F13: Machine Learning Zoubin Ghahramani and Carl Edward Rasmussen Department of Engineering University of Cambridge http://mlg.eng.cam.ac.uk/teaching/4f13/ Ghahramani
More informationHilbert Space Methods for Reduced-Rank Gaussian Process Regression
Hilbert Space Methods for Reduced-Rank Gaussian Process Regression Arno Solin and Simo Särkkä Aalto University, Finland Workshop on Gaussian Process Approximation Copenhagen, Denmark, May 2015 Solin &
More informationVCMC: Variational Consensus Monte Carlo
VCMC: Variational Consensus Monte Carlo Maxim Rabinovich, Elaine Angelino, Michael I. Jordan Berkeley Vision and Learning Center September 22, 2015 probabilistic models! sky fog bridge water grass object
More informationGaussian with mean ( µ ) and standard deviation ( σ)
Slide from Pieter Abbeel Gaussian with mean ( µ ) and standard deviation ( σ) 10/6/16 CSE-571: Robotics X ~ N( µ, σ ) Y ~ N( aµ + b, a σ ) Y = ax + b + + + + 1 1 1 1 1 1 1 1 1 1, ~ ) ( ) ( ), ( ~ ), (
More informationContents. Part I: Fundamentals of Bayesian Inference 1
Contents Preface xiii Part I: Fundamentals of Bayesian Inference 1 1 Probability and inference 3 1.1 The three steps of Bayesian data analysis 3 1.2 General notation for statistical inference 4 1.3 Bayesian
More informationSpatial Statistics with Image Analysis. Outline. A Statistical Approach. Johan Lindström 1. Lund October 6, 2016
Spatial Statistics Spatial Examples More Spatial Statistics with Image Analysis Johan Lindström 1 1 Mathematical Statistics Centre for Mathematical Sciences Lund University Lund October 6, 2016 Johan Lindström
More informationPareto Smoothed Importance Sampling
Pareto Smoothed Importance Sampling Aki Vehtari Andrew Gelman Jonah Gabry 22 September 216 Abstract Importance weighting is a convenient general way to adjust for draws from the wrong distribution, but
More information