Fast algorithms that utilise the sparsity of Q exist (c-package ( Q 1. GMRFlib) Seattle February 3, 2009
|
|
- Maurice Boyd
- 6 years ago
- Views:
Transcription
1 GMRF:s Results References Some Theory for and Applications of Gaussian Markov Random Fields Johan Lindström 1 David Bolin 1 Finn Lindgren 1 Håvard Rue 2 1 Centre for Mathematical Sciences Lund University & 2 Department of Mathematical Sciences NTNU, Trondheim Seattle February 3, 2009 GMRF:s Results References Overview Gaussian Markov random fields Basics Approximating Matérn covariances INLA Fast estimation Results The Sahel Vegetation Precipitation Depth Data Gaussian Markov Random Fields (GMRF:s) A Gaussian Markov random field (GMRF) is a Gaussian random field with a Markov property. The neighbours Ni to a point si are the points that in some sense are close to si. The Gaussian random field x N ( Ñ,Q 1) has a joint distribution that satisfies p(xi {xj : j i}) = p(xi {xj : j Ni}). j / Ni xi xj { xk : k / {i,j} } Qi,j = 0. The density is p(x) = Q 1/2 exp (2Ô) n/2 ( 1 2 (x Ñ) Q(x Ñ) ) Fast algorithms that utilise the sparsity of Q exist (c-package GMRFlib) See Rue and Held (2005) for extensive details on GMRF:s. GMRF:s Simplified expressions Formulating the density through the precision matrix instead of through the covariance matrix simplifies several common expression for Gaussian models. Conditional expectation {xi xj,j I} N ( Ñ I Q 1 I,I Q I,J(xJ Ñ J ),Q 1 I,I ). Hirearchical modelling x θ N ( Bθ,Q 1), θ N ( Ñ,Q 1 ) θ, [ ] ( [BÑ ] [ x Q QB N, θ B Q B QB + Q θ ] 1 ) Ñ
2 GMRF:s Simplified expressions Formulating the density through the precision matrix instead of through the covariance matrix simplifies several common expression for Gaussian models. Conditional expectation {xi xj,j I} N ( Ñ I Q 1 I,I Q I,J(xJ Ñ J ),Q 1 I,I ). Hirearchical modelling x θ N ( Bθ,Q 1), θ N ( Ñ,Q 1 ) θ, [ ] ( [BÑ ] [ x Q QB N, θ B Q B QB + Q θ ] 1 ) GMRF:s Previous limitations How to choose or construct Q? A GMRF might be computationally effective but it is difficult in construct precision matrices that result in reasonable covariance functions for the underlying Gaussian fields. Various ad-hoc methods exist. A common solution is to use a small neighbourhood and let the precision between two points depend on the distance between the points. Rue and Tjelmeland (2002) created GMRF:s on rectangular grids in that approximate Gaussian fields with a wide class of covariance functions. Having the field defined only on a regular grid leads to issues with mapping the observations to the grid points. GMRF:s Previous limitations How to choose or construct Q? A GMRF might be computationally effective but it is difficult in construct precision matrices that result in reasonable covariance functions for the underlying Gaussian fields. Various ad-hoc methods exist. A common solution is to use a small neighbourhood and let the precision between two points depend on the distance between the points. Rue and Tjelmeland (2002) created GMRF:s on rectangular grids in that approximate Gaussian fields with a wide class of covariance functions. Having the field defined only on a regular grid leads to issues with mapping the observations to the grid points. GMRF:s Previous limitations How to choose or construct Q? A GMRF might be computationally effective but it is difficult in construct precision matrices that result in reasonable covariance functions for the underlying Gaussian fields. Various ad-hoc methods exist. A common solution is to use a small neighbourhood and let the precision between two points depend on the distance between the points. Rue and Tjelmeland (2002) created GMRF:s on rectangular grids in that approximate Gaussian fields with a wide class of covariance functions. Having the field defined only on a regular grid leads to issues with mapping the observations to the grid points. Ñ
3 Matérn covariances (Bertil Matérn, ) The Matérn covariance family on u R d : r(u,v) = C(x(u),x(v)) = 2 21 Ò (Ò) ( v u )Ò K Ò ( v u ) with scale (inverse range) > 0 and shape/smoothness Ò > 0, and K Ò a modified Bessel function. Fields with Matérn covariances are solutions to an SPDE Whittle (1954) based on the Laplacian, =, ( 2 ) /2 x(u) = Ø 2 (u), where (u) is spatial white noise, and = Ò + d/2. Parameter link: 2 = Ø 2 (Ò) SPDE issues Non-uniqueness: If x(u) is a solution to the SPDE for = 2, so is x(u) + c exp( e u), for any unit length vector e and any constant c. Non-stationarity: On a bounded domain, the SPDE solutions are non-stationary, unless conditioned on suitable boundary distributions. Practical solution to the non-uniqueness and non-stationarity: Zero-normal-derivative (Neumann) boundaries reduce the impact of the null-space solutions. Resulting covariance, for Ï = [0, L] R: C(x(u),x(v)) rm(u,v) + rm(u, v) + rm(u,2l v) = rm(0,v u) + rm(0,v + u) + rm(0,2l (v + u)). SPDE issues Non-uniqueness: If x(u) is a solution to the SPDE for = 2, so is x(u) + c exp( e u), for any unit length vector e and any constant c. Non-stationarity: On a bounded domain, the SPDE solutions are non-stationary, unless conditioned on suitable boundary distributions. Practical solution to the non-uniqueness and non-stationarity: Zero-normal-derivative (Neumann) boundaries reduce the impact of the null-space solutions. Resulting covariance, for Ï = [0, L] R: C(x(u),x(v)) rm(u,v) + rm(u, v) + rm(u,2l v) = rm(0,v u) + rm(0,v + u) + rm(0,2l (v + u)). The finite element method A stochastic weak formulation of the SPDE states that [ k, ( 2 ) /2 x ] k=1,...,n D = [ k, W ] k=1,...,n for each set of test functions { k}. We use N simple piecewise linear test functions, equal to the piecewise linear basis functions, x(u) = j k(u)wj, and compute the resulting distributions by explicitly calculating the expectation vector and precision matrix for w (= x). Ý For = 2, the weak formulation can be written Kw = [ ] i, 2 )Ýj w = D [ k, W ] i,j=1,...,n ( k=1,...,n ( ) 2Ò (4Ô) d/2
4 Construction of Q With the help of Green s first identity, Ci,j = Ýi, Ýj, Ki,j = Ýi, ( 2 )Ýj = 2 C i,j + Ýi, Ýj, for Neumann boundaries. Markovified Least Squares and Galerkin solutions: C = diag( Ýi, 1 ), ( optimal approximation) Q 1, = K, (Least Squares) Q 2, = KC 1 K, (Galerkin) Q, = KC 1 Q 2, C 1 K, = 3,4,... (Galerkin recursion) Neighbourhood radius equals. Simple closed-form expressions for the matrix elements. Approximating Matérn covariances (Lindgren and Rue, 2007) That fields with Matérn covariances are solutions to an SPDE has been used to construct GMRF:s that approximate fields with Matérn covariance for Z +. Well-defined SPDE on curved manifolds (e.g. a globe). Possible to introduce both anisotropy and non-stationarity. Spatially oscillating fields can be introduced via a complex version of the SPDE: (h1 + ih2 (H1 + ih2) )(x1(u) + ix2(u)) = 1(u) + i 2(u) The precision matrix of the approximating GMRF is found using the finite element method on a triangulation of irregularly spaced points. The resulting GMRF is defined on the points of the triangulation, making it suitable for modelling fields that are observed at irregular locations. Approximating Matérn covariances (Lindgren and Rue, 2007) That fields with Matérn covariances are solutions to an SPDE has been used to construct GMRF:s that approximate fields with Matérn covariance for Z +. Well-defined SPDE on curved manifolds (e.g. a globe). Possible to introduce both anisotropy and non-stationarity. Spatially oscillating fields can be introduced via a complex version of the SPDE: (h1 + ih2 (H1 + ih2) )(x1(u) + ix2(u)) = 1(u) + i 2(u) The precision matrix of the approximating GMRF is found using the finite element method on a triangulation of irregularly spaced points. The resulting GMRF is defined on the points of the triangulation, making it suitable for modelling fields that are observed at irregular locations. Approximating Matérn covariances (Lindgren and Rue, 2007) That fields with Matérn covariances are solutions to an SPDE has been used to construct GMRF:s that approximate fields with Matérn covariance for Z +. Well-defined SPDE on curved manifolds (e.g. a globe). Possible to introduce both anisotropy and non-stationarity. Spatially oscillating fields can be introduced via a complex version of the SPDE: (h1 + ih2 (H1 + ih2) )(x1(u) + ix2(u)) = 1(u) + i 2(u) The precision matrix of the approximating GMRF is found using the finite element method on a triangulation of irregularly spaced points. The resulting GMRF is defined on the points of the triangulation, making it suitable for modelling fields that are observed at irregular locations.
5 INLA Fast estimation (Rue and Martino, 2007) Assume an underlying GMRF with (possibly non-gaussian) point observations, i.e. x N ( Ñ( ),Q( ) 1) p(yi x, ) = p(yi xi, ). If the observations are non-gaussian the log-likelihood can be approximated with a Gaussian by a Taylor expansion of the observation density, p(yi xi, ). Do numerical optimisation of the log-likelihood, using the above approximation. Once the MAP-estimator is found, use the Hessian of the log-likelihood to determine the likely variations of the likelihood and do numerical integration. INLA Fast estimation (Rue and Martino, 2007) Assume an underlying GMRF with (possibly non-gaussian) point observations, i.e. x N ( Ñ( ),Q( ) 1) p(yi x, ) = p(yi xi, ). If the observations are non-gaussian the log-likelihood can be approximated with a Gaussian by a Taylor expansion of the observation density, p(yi xi, ). Do numerical optimisation of the log-likelihood, using the above approximation. Once the MAP-estimator is found, use the Hessian of the log-likelihood to determine the likely variations of the likelihood and do numerical integration. INLA Fast estimation (Rue and Martino, 2007) Assume an underlying GMRF with (possibly non-gaussian) point observations, i.e. x N ( Ñ( ),Q( ) 1) p(yi x, ) = p(yi xi, ). If the observations are non-gaussian the log-likelihood can be approximated with a Gaussian by a Taylor expansion of the observation density, p(yi xi, ). Do numerical optimisation of the log-likelihood, using the above approximation. Once the MAP-estimator is found, use the Hessian of the log-likelihood to determine the likely variations of the likelihood and do numerical integration. INLA Fast estimation (Rue and Martino, 2007) Assume an underlying GMRF with (possibly non-gaussian) point observations, i.e. x N ( Ñ( ),Q( ) 1) p(yi x, ) = p(yi xi, ). If the observations are non-gaussian the log-likelihood can be approximated with a Gaussian by a Taylor expansion of the observation density, p(yi xi, ). Do numerical optimisation of the log-likelihood, using the above approximation. Once the MAP-estimator is found, use the Hessian of the log-likelihood to determine the likely variations of the likelihood and do numerical integration.
6 INLA Fast estimation (cont.) The Good: Very fast compared to MCMC. Gives GOOD approximations of the posterior densities, p(xi Y, map), p(xi Y), p( i Y). The Bad: Limited to a few (5-10) hyperparameters,. Gives limited information about interaction between different parameters, p(xi,xj Y), p( i, j Y). The Ugly: Large memory requirments compared to MCMC of a GMRF. This is mainly due to the numerical integration. INLA Fast estimation (cont.) The Good: Very fast compared to MCMC. Gives GOOD approximations of the posterior densities, p(xi Y, map), p(xi Y), p( i Y). The Bad: Limited to a few (5-10) hyperparameters,. Gives limited information about interaction between different parameters, p(xi,xj Y), p( i, j Y). The Ugly: Large memory requirments compared to MCMC of a GMRF. This is mainly due to the numerical integration. INLA Fast estimation (cont.) The Good: Very fast compared to MCMC. Gives GOOD approximations of the posterior densities, p(xi Y, map), p(xi Y), p( i Y). The Bad: Limited to a few (5-10) hyperparameters,. Gives limited information about interaction between different parameters, p(xi,xj Y), p( i, j Y). The Ugly: Large memory requirments compared to MCMC of a GMRF. This is mainly due to the numerical integration. African Sahel The inspiration for the following applications was an article by Eklundh and Olsson (2003). A semi-arid region directly south of the Sahara desert. Starting in the late 1960 s, the area suffered droughts for over twenty years. Recent studies indicate a vegetation recovery.
7 Normalised Difference Vegetation Index Data taken from the NOAA/NASA Pathfinder AVHRR Land (PAL) data set. NDVI is calculated using satellite measurements of the reflectance from the Earth s surface. The PAL data set contains 36 measurements per year. We use aggregated yearly data. The question is Has the amount of vegetation increased? Figure: NDVI data for the Sahel 1983 Vegetation model (Bolin et al., 2008) Spatial model for measured NDVI: Y(si,t) = xt(si) + it, Yt = AtXt + Et. where xt are latent fields constrained to a sum of temporal trends, xt(si) = m j=1 fj(t) Kj(si) where fj are temporal trend functions with spatially varying regression coefficients Kj The observation matrices At determine which points are observed at each time point Regression model Without the trend constraint, a simple prior model for xt would be a Whittle field with = 0, with precision Q x. The trend restriction provides a natural prior for K = [K 1,...,K m] as K Ø N(0,(ØQ) 1 ) where Q = (F F) Q x, and F = [f1,...,fm] The residual variance is allowed to vary across space, Ë = diag( 2 i ) Collecting the data and observation matrices in a similar manner, Y = AK + E = diag(at)(f In)K + E Posterior distribution The posterior distribution for K given the data Y = [Y 0,...,Y T 1] and the parameters is Ë (K Y,, Ø) N(Ñ K,Q 1 K ), with K A Ë 1 Ñ K = Q 1 Q K = ØQ + A 1 Ë Y, The unknown parameters (Ø, 2 1,..., 2 n) can be estimated using an EM algorithm, exploiting the GMRF structure in the updating equations. A.
8 Linear trend estimation Slope (upper) and Intercept (lower) estimates. K 2 Estimate K 1 Estimate Trend significance Significance estimates for the slope of the linear trends OLS GMRF Significant Negative Non Significant Negative Non Significant Positive Significant Positive Precipitation (Lindström and Lindgren, 2008) Yearly precipitation, in metres, for the 443 stations that have reported measurements for Data from the Global Historical Climatology Network. Irregular triangulation We do not need a dense regular grid for modelling the precipitation. The GMRF construction allows for irregular triangulations, which reduces the computational burden
9 Spatio-temporal model Given a latent precipitation field x, transformed observations are assumed Gaussian: Y(si,t) x N(x(si,t), 2 ). We use an AR(1)-process in time, with spatially correlated noise and spatially varying mean: Xt+1 Ñ = a (Xt Ñ ) + t, S ), X1 N(Ñ, Q 1 S /(1 a2 )), = Bθ. Ñ ( t 2, Ø) N(0, Q 1 and we can write X N(1 Bθ, ( Q t Q s ) 1) Basis functions for the mean transformed precipitation latitude Left Latitude is modelled using a broken linear trend. Right In addition to the broken trend 15 B-spline surface basis functions are used. Hierarchical model graph Ø Ø Ñ θ Q θ Ø 2 a θ X 2 Y Directed acyclic graph describing the resulting hierarchical model. The jointly Gaussian fields X and θ can be integrated out. Markov chain Monte Carlo First we transform the parameters to obtain parameters valid on R: 2 = log( 2 ), ã = log(1 + a) log(1 a), 2 = 2 ), Õ = log(õ), log( Construct a proposal distribution by using a random walk proposal on the transformed variable space. However initial runs on a reduced datasets with coarser triangulation showed strong dependencies in the posterior distribution, leading us to a correlated proposal distribution: Ýnew N ( Ýold, Ë proposal).
10 Markov chain Monte Carlo First we transform the parameters to obtain parameters valid on R: 2 = log( 2 ), ã = log(1 + a) log(1 a), 2 = 2 ), Õ = log(õ), log( Construct a proposal distribution by using a random walk proposal on the transformed variable space. However initial runs on a reduced datasets with coarser triangulation showed strong dependencies in the posterior distribution, leading us to a correlated proposal distribution: Ýnew N ( Ýold, Ë proposal). Markov chain Monte Carlo (cont.) 2 a 2 a 2 2 ã Õ 2 ã 2 Two-dimensional histograms illustrating the dependence between the components of (Ý Y) before (left pane) and after (right pane) the transformation. Computational burden The dominating cost for each MCMC iteration is calculation of the Cholesky factorisation of Q. Inverting a full covariance matrix is O ( n 3), Given a spatial GMRF on a lattice with n points the Cholesky factor is O ( n 3/2), Given a spatio-temporal GMRF on a lattice with n points the Cholesky factor is O ( n 2), How bad is the additional burden for the temporal dependence? n Spatial Spatio-temporal s = s Computational burden The dominating cost for each MCMC iteration is calculation of the Cholesky factorisation of Q. Inverting a full covariance matrix is O ( n 3), Given a spatial GMRF on a lattice with n points the Cholesky factor is O ( n 3/2), Given a spatio-temporal GMRF on a lattice with n points the Cholesky factor is O ( n 2), How bad is the additional burden for the temporal dependence? n Spatial Spatio-temporal s = s Õ
11 Computational burden The dominating cost for each MCMC iteration is calculation of the Cholesky factorisation of Q. Inverting a full covariance matrix is O ( n 3), Given a spatial GMRF on a lattice with n points the Cholesky factor is O ( n 3/2), Given a spatio-temporal GMRF on a lattice with n points the Cholesky factor is O ( n 2), How bad is the additional burden for the temporal dependence? n Spatial Spatio-temporal s = s 46.8 s Kriging results from MCMC simulation of (Ø, 2, a, 2 ) Interpolated yearly precipitation, in metres, for Depth Data We have measurements of seabead depth for an area outside the Swedish coast. The data has been modelled using an underlying GMRF field with Matérn covariance and unknown expectation consisting of an arbitrary plane (e.g. Kriging). The roughly measurement locations have been trinagulated resulting in a grid containing points. Posterior densities for the three hyperparameters of the model and for the underlying field where estimated using INLA. Depth Data We have measurements of seabead depth for an area outside the Swedish coast. The data has been modelled using an underlying GMRF field with Matérn covariance and unknown expectation consisting of an arbitrary plane (e.g. Kriging). The roughly measurement locations have been trinagulated resulting in a grid containing points. Posterior densities for the three hyperparameters of the model and for the underlying field where estimated using INLA. Total estimation time on a Core2Duo desktop: 8 minutes.
12 Depth Data Results GMRF:s Results References Bibliography Bolin, D., Lindström, J., Eklundh, L., and Lindgren, F. (2008), Fast Estimation of Spatially Dependent Temporal Vegetation Trends using Gaussian Markov Random Fields, Submitted to Comput. Statist. and Data Anal. Eklundh, L. and Olsson, L. (2003), Vegetation index trends for the African Sahel , Geophys. Res. Lett., 30, Lindgren, F. and Rue, H. (2007), Explicit construction of GMRF approximations to generalised Matérn Fields on irregular grids, Tech. Rep. 12, Centre for Mathematical Sciences, Lund University, Lund, Sweden. Lindström, J. and Lindgren, F. (2008), Modeling Yearly cumulative Precipitation over the African Sahel using a Gaussian Markov Random Field, In preparation. Rue, H. and Held, L. (2005), Gaussian Markov Random Fields; Theory and Applications, vol. 104 of Monographs on Statistics and Applied Probability, Chapman & Hall/CRC. Rue, H. and Martino, S. (2007), Approximate Bayesian inference for hierarchical Gaussian Markov random field models, J. Statist. Plann. and Inference, 137, Rue, H. and Tjelmeland, H. (2002), Fitting Gaussian Markov Random Fields to Gaussian Fields, Scand. J. Statist., 29, Whittle, P. (1954), On Stationary Processes in the Plane, Biometrika, 41,
Data are collected along transect lines, with dense data along. Spatial modelling using GMRFs. 200m. Today? Spatial modelling using GMRFs
Centre for Mathematical Sciences Lund University Engineering geology Lund University Results A non-stationary extension The model Estimation Gaussian Markov random fields Basics Approximating Mate rn covariances
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 informationA GAUSSIAN MARKOV RANDOM FIELD MODEL FOR TOTAL YEARLY PRECIPITATION
A GAUSSIAN MARKOV RANDOM FIELD MODEL FOR TOTAL YEARLY PRECIPITATION OVER THE AFRICAN SAHEL JOHAN LINDSTRÖM AND FINN LINDGREN Preprints in Mathematical Sciences 2008:8 Centre for Mathematical Sciences Mathematical
More informationR-INLA. Sam Clifford Su-Yun Kang Jeff Hsieh. 30 August Bayesian Research and Analysis Group 1 / 14
1 / 14 R-INLA Sam Clifford Su-Yun Kang Jeff Hsieh Bayesian Research and Analysis Group 30 August 2012 What is R-INLA? R package for Bayesian computation Integrated Nested Laplace Approximation MCMC free
More informationMarkov random fields. The Markov property
Markov random fields The Markov property Discrete time: (X k X k!1,x k!2,... = (X k X k!1 A time symmetric version: (X k! X!k = (X k X k!1,x k+1 A more general version: Let A be a set of indices >k, B
More informationLatent Gaussian Processes and Stochastic Partial Differential Equations
Latent Gaussian Processes and Stochastic Partial Differential Equations Johan Lindström 1 1 Centre for Mathematical Sciences Lund University Lund 2016-02-04 Johan Lindström - johanl@maths.lth.se Gaussian
More informationSpatial Statistics with Image Analysis. Lecture L08. Computer exercise 3. Lecture 8. Johan Lindström. November 25, 2016
C3 Repetition Creating Q Spectral Non-grid Spatial Statistics with Image Analysis Lecture 8 Johan Lindström November 25, 216 Johan Lindström - johanl@maths.lth.se FMSN2/MASM25L8 1/39 Lecture L8 C3 Repetition
More informationComputation fundamentals of discrete GMRF representations of continuous domain spatial models
Computation fundamentals of discrete GMRF representations of continuous domain spatial models Finn Lindgren September 23 2015 v0.2.2 Abstract The fundamental formulas and algorithms for Bayesian spatial
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 informationarxiv: v4 [stat.me] 14 Sep 2015
Does non-stationary spatial data always require non-stationary random fields? Geir-Arne Fuglstad 1, Daniel Simpson 1, Finn Lindgren 2, and Håvard Rue 1 1 Department of Mathematical Sciences, NTNU, Norway
More informationGaussian Processes 1. Schedule
1 Schedule 17 Jan: Gaussian processes (Jo Eidsvik) 24 Jan: Hands-on project on Gaussian processes (Team effort, work in groups) 31 Jan: Latent Gaussian models and INLA (Jo Eidsvik) 7 Feb: Hands-on project
More informationBayesian spatial hierarchical modeling for temperature extremes
Bayesian spatial hierarchical modeling for temperature extremes Indriati Bisono Dr. Andrew Robinson Dr. Aloke Phatak Mathematics and Statistics Department The University of Melbourne Maths, Informatics
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 informationMultivariate Gaussian Random Fields with SPDEs
Multivariate Gaussian Random Fields with SPDEs Xiangping Hu Daniel Simpson, Finn Lindgren and Håvard Rue Department of Mathematics, University of Oslo PASI, 214 Outline The Matérn covariance function and
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 informationOverlapping block proposals for latent Gaussian Markov random fields
NORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET Overlapping block proposals for latent Gaussian Markov random fields by Ingelin Steinsland and Håvard Rue PREPRINT STATISTICS NO. 8/3 NORWEGIAN UNIVERSITY
More informationNearest Neighbor Gaussian Processes for Large Spatial Data
Nearest Neighbor Gaussian Processes for Large Spatial Data Abhi Datta 1, Sudipto Banerjee 2 and Andrew O. Finley 3 July 31, 2017 1 Department of Biostatistics, Bloomberg School of Public Health, Johns
More informationFast kriging of large data sets with Gaussian Markov random fields
Computational Statistics & Data Analysis 52 (2008) 233 2349 www.elsevier.com/locate/csda Fast kriging of large data sets with Gaussian Markov random fields Linda Hartman a,,, Ola Hössjer b,2 a Centre for
More informationLikelihood NIPS July 30, Gaussian Process Regression with Student-t. Likelihood. Jarno Vanhatalo, Pasi Jylanki and Aki Vehtari NIPS-2009
with with July 30, 2010 with 1 2 3 Representation Representation for Distribution Inference for the Augmented Model 4 Approximate Laplacian Approximation Introduction to Laplacian Approximation Laplacian
More informationSpatial Statistics with Image Analysis. Lecture L11. Home assignment 3. Lecture 11. Johan Lindström. December 5, 2016.
HA3 MRF:s Simulation Estimation Spatial Statistics with Image Analysis Lecture 11 Johan Lindström December 5, 2016 Lecture L11 Johan Lindström - johanl@maths.lth.se FMSN20/MASM25 L11 1/22 HA3 MRF:s Simulation
More informationarxiv: v1 [stat.co] 16 May 2011
Fast approximate inference with INLA: the past, the present and the future arxiv:1105.2982v1 [stat.co] 16 May 2011 Daniel Simpson, Finn Lindgren and Håvard Rue Department of Mathematical Sciences Norwegian
More informationMarkov Chain Monte Carlo methods
Markov Chain Monte Carlo methods By Oleg Makhnin 1 Introduction a b c M = d e f g h i 0 f(x)dx 1.1 Motivation 1.1.1 Just here Supresses numbering 1.1.2 After this 1.2 Literature 2 Method 2.1 New math As
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 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 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 informationarxiv: v1 [stat.ap] 18 Apr 2011
The Annals of Applied Statistics 2011, Vol. 5, No. 1, 523 550 DOI: 10.1214/10-AOAS383 c Institute of Mathematical Statistics, 2011 arxiv:1104.3436v1 [stat.ap] 18 Apr 2011 SPATIAL MODELS GENERATED BY NESTED
More informationRECONSTRUCTION OF NDVI DATA WITH NORMAL VARIANCE-MEAN MIXTURE NOISE USING STOCHASTIC GRADIENT METHODS
RECONSTRUCTION OF NDVI DATA WITH NORMAL VARIANCE-MEAN MIXTURE NOISE USING STOCHASTIC GRADIENT METHODS DAVID KRAFT Master s thesis 217:E1 Faculty of Engineering Centre for Mathematical Sciences Mathematical
More informationSpatial point processes in the modern world an
Spatial point processes in the modern world an interdisciplinary dialogue Janine Illian University of St Andrews, UK and NTNU Trondheim, Norway Bristol, October 2015 context statistical software past to
More informationAn explicit link between Gaussian fields and Gaussian Markov random fields: the stochastic partial differential equation approach
J. R. Statist. Soc. B (2011) 73, Part 4, pp. 423 498 An explicit link between Gaussian fields and Gaussian Markov random fields: the stochastic partial differential equation approach Finn Lindgren and
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 informationList of projects. FMS020F NAMS002 Statistical inference for partially observed stochastic processes, 2016
List of projects FMS020F NAMS002 Statistical inference for partially observed stochastic processes, 206 Work in groups of two (if this is absolutely not possible for some reason, please let the lecturers
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 informationComputer Vision Group Prof. Daniel Cremers. 4. Gaussian Processes - Regression
Group Prof. Daniel Cremers 4. Gaussian Processes - Regression Definition (Rep.) Definition: A Gaussian process is a collection of random variables, any finite number of which have a joint Gaussian distribution.
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 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 informationIMAGE MODELLING AND ESTIMATION FINN LINDGREN A STATISTICAL APPROACH CENTRUM SCIENTIARUM MATHEMATICARUM. Third edition, , Chapter 1 5
IMAGE MODELLING AND ESTIMATION A STATISTICAL APPROACH FINN LINDGREN Third edition, 2006-10-20, Chapter 1 5 Centre for Mathematical Sciences Mathematical Statistics CENTRUM SCIENTIARUM MATHEMATICARUM 1
More informationBayesian Inference: Principles and Practice 3. Sparse Bayesian Models and the Relevance Vector Machine
Bayesian Inference: Principles and Practice 3. Sparse Bayesian Models and the Relevance Vector Machine Mike Tipping Gaussian prior Marginal prior: single α Independent α Cambridge, UK Lecture 3: Overview
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 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 informationGaussian predictive process models for large spatial data sets.
Gaussian predictive process models for large spatial data sets. Sudipto Banerjee, Alan E. Gelfand, Andrew O. Finley, and Huiyan Sang Presenters: Halley Brantley and Chris Krut September 28, 2015 Overview
More informationDevelopment of Stochastic Artificial Neural Networks for Hydrological Prediction
Development of Stochastic Artificial Neural Networks for Hydrological Prediction G. B. Kingston, M. F. Lambert and H. R. Maier Centre for Applied Modelling in Water Engineering, School of Civil and Environmental
More informationTechnical Vignette 5: Understanding intrinsic Gaussian Markov random field spatial models, including intrinsic conditional autoregressive models
Technical Vignette 5: Understanding intrinsic Gaussian Markov random field spatial models, including intrinsic conditional autoregressive models Christopher Paciorek, Department of Statistics, University
More informationIndex. Geostatistics for Environmental Scientists, 2nd Edition R. Webster and M. A. Oliver 2007 John Wiley & Sons, Ltd. ISBN:
Index Akaike information criterion (AIC) 105, 290 analysis of variance 35, 44, 127 132 angular transformation 22 anisotropy 59, 99 affine or geometric 59, 100 101 anisotropy ratio 101 exploring and displaying
More informationInference and estimation in probabilistic time series models
1 Inference and estimation in probabilistic time series models David Barber, A Taylan Cemgil and Silvia Chiappa 11 Time series The term time series refers to data that can be represented as a sequence
More informationSpatial Statistics with Image Analysis. Lecture L02. Computer exercise 0 Daily Temperature. Lecture 2. Johan Lindström.
C Stochastic fields Covariance Spatial Statistics with Image Analysis Lecture 2 Johan Lindström November 4, 26 Lecture L2 Johan Lindström - johanl@maths.lth.se FMSN2/MASM2 L /2 C Stochastic fields Covariance
More informationOn Gaussian Process Models for High-Dimensional Geostatistical Datasets
On Gaussian Process Models for High-Dimensional Geostatistical Datasets Sudipto Banerjee Joint work with Abhirup Datta, Andrew O. Finley and Alan E. Gelfand University of California, Los Angeles, USA May
More informationA Framework for Daily Spatio-Temporal Stochastic Weather Simulation
A Framework for Daily Spatio-Temporal Stochastic Weather Simulation, Rick Katz, Balaji Rajagopalan Geophysical Statistics Project Institute for Mathematics Applied to Geosciences National Center for Atmospheric
More informationBayesian multiscale analysis for time series data
Computational Statistics & Data Analysis 5 (26) 79 73 www.elsevier.com/locate/csda Bayesian multiscale analysis for time series data Tor Arne ]igård a,,håvard Rue b, Fred Godtliebsen c a Department of
More informationMonte Carlo in Bayesian Statistics
Monte Carlo in Bayesian Statistics Matthew Thomas SAMBa - University of Bath m.l.thomas@bath.ac.uk December 4, 2014 Matthew Thomas (SAMBa) Monte Carlo in Bayesian Statistics December 4, 2014 1 / 16 Overview
More informationStatistics for extreme & sparse data
Statistics for extreme & sparse data University of Bath December 6, 2018 Plan 1 2 3 4 5 6 The Problem Climate Change = Bad! 4 key problems Volcanic eruptions/catastrophic event prediction. Windstorms
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 informationSpatial smoothing over complex domain
Spatial smoothing over complex domain Alessandro Ottavi NTNU Department of Mathematical Science August 16-20 2010 Outline Spatial smoothing Complex domain The SPDE approch Manifold Some tests Spatial smoothing
More informationHierarchical Nearest-Neighbor Gaussian Process Models for Large Geo-statistical Datasets
Hierarchical Nearest-Neighbor Gaussian Process Models for Large Geo-statistical Datasets Abhirup Datta 1 Sudipto Banerjee 1 Andrew O. Finley 2 Alan E. Gelfand 3 1 University of Minnesota, Minneapolis,
More informationGaussian Markov Random Fields: Theory and Applications
Gaussian Markov Random Fields: Theory and Applications Håvard Rue 1 1 Department of Mathematical Sciences NTNU, Norway September 26, 2008 Part I Definition and basic properties Outline I Introduction Why?
More informationComputer Vision Group Prof. Daniel Cremers. 9. Gaussian Processes - Regression
Group Prof. Daniel Cremers 9. Gaussian Processes - Regression Repetition: Regularized Regression Before, we solved for w using the pseudoinverse. But: we can kernelize this problem as well! First step:
More informationMarkov Chain Monte Carlo methods
Markov Chain Monte Carlo methods Tomas McKelvey and Lennart Svensson Signal Processing Group Department of Signals and Systems Chalmers University of Technology, Sweden November 26, 2012 Today s learning
More informationDAG models and Markov Chain Monte Carlo methods a short overview
DAG models and Markov Chain Monte Carlo methods a short overview Søren Højsgaard Institute of Genetics and Biotechnology University of Aarhus August 18, 2008 Printed: August 18, 2008 File: DAGMC-Lecture.tex
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 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 informationCovariance function estimation in Gaussian process regression
Covariance function estimation in Gaussian process regression François Bachoc Department of Statistics and Operations Research, University of Vienna WU Research Seminar - May 2015 François Bachoc Gaussian
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 informationNonstationary spatial process modeling Part II Paul D. Sampson --- Catherine Calder Univ of Washington --- Ohio State University
Nonstationary spatial process modeling Part II Paul D. Sampson --- Catherine Calder Univ of Washington --- Ohio State University this presentation derived from that presented at the Pan-American Advanced
More informationGRAPH SIGNAL PROCESSING: A STATISTICAL VIEWPOINT
GRAPH SIGNAL PROCESSING: A STATISTICAL VIEWPOINT Cha Zhang Joint work with Dinei Florêncio and Philip A. Chou Microsoft Research Outline Gaussian Markov Random Field Graph construction Graph transform
More informationUncertainty quantification for Wavefield Reconstruction Inversion
Uncertainty quantification for Wavefield Reconstruction Inversion Zhilong Fang *, Chia Ying Lee, Curt Da Silva *, Felix J. Herrmann *, and Rachel Kuske * Seismic Laboratory for Imaging and Modeling (SLIM),
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 informationMulti-resolution models for large data sets
Multi-resolution models for large data sets Douglas Nychka, National Center for Atmospheric Research National Science Foundation NORDSTAT, Umeå, June, 2012 Credits Steve Sain, NCAR Tia LeRud, UC Davis
More informationAn EM algorithm for Gaussian Markov Random Fields
An EM algorithm for Gaussian Markov Random Fields Will Penny, Wellcome Department of Imaging Neuroscience, University College, London WC1N 3BG. wpenny@fil.ion.ucl.ac.uk October 28, 2002 Abstract Lavine
More informationComputer Intensive Methods in Mathematical Statistics
Computer Intensive Methods in Mathematical Statistics Department of mathematics johawes@kth.se Lecture 16 Advanced topics in computational statistics 18 May 2017 Computer Intensive Methods (1) Plan of
More informationProbabilistic Graphical Models
Probabilistic Graphical Models Brown University CSCI 2950-P, Spring 2013 Prof. Erik Sudderth Lecture 13: Learning in Gaussian Graphical Models, Non-Gaussian Inference, Monte Carlo Methods Some figures
More informationDeblurring Jupiter (sampling in GLIP faster than regularized inversion) Colin Fox Richard A. Norton, J.
Deblurring Jupiter (sampling in GLIP faster than regularized inversion) Colin Fox fox@physics.otago.ac.nz Richard A. Norton, J. Andrés Christen Topics... Backstory (?) Sampling in linear-gaussian hierarchical
More informationA BAYESIAN SPATIO-TEMPORAL GEOSTATISTICAL MODEL WITH AN AUXILIARY LATTICE FOR LARGE DATASETS
Statistica Sinica 25 (2015), 61-79 doi:http://dx.doi.org/10.5705/ss.2013.085w A BAYESIAN SPATIO-TEMPORAL GEOSTATISTICAL MODEL WITH AN AUXILIARY LATTICE FOR LARGE DATASETS Ganggang Xu 1, Faming Liang 1
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 informationA STATISTICAL TECHNIQUE FOR MODELLING NON-STATIONARY SPATIAL PROCESSES
A STATISTICAL TECHNIQUE FOR MODELLING NON-STATIONARY SPATIAL PROCESSES JOHN STEPHENSON 1, CHRIS HOLMES, KERRY GALLAGHER 1 and ALEXANDRE PINTORE 1 Dept. Earth Science and Engineering, Imperial College,
More informationProbabilistic Graphical Models
2016 Robert Nowak Probabilistic Graphical Models 1 Introduction We have focused mainly on linear models for signals, in particular the subspace model x = Uθ, where U is a n k matrix and θ R k is a vector
More informationMCMC algorithms for fitting Bayesian models
MCMC algorithms for fitting Bayesian models p. 1/1 MCMC algorithms for fitting Bayesian models Sudipto Banerjee sudiptob@biostat.umn.edu University of Minnesota MCMC algorithms for fitting Bayesian models
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 informationGAUSSIAN PROCESS TRANSFORMS
GAUSSIAN PROCESS TRANSFORMS Philip A. Chou Ricardo L. de Queiroz Microsoft Research, Redmond, WA, USA pachou@microsoft.com) Computer Science Department, Universidade de Brasilia, Brasilia, Brazil queiroz@ieee.org)
More informationMarkov Chain Monte Carlo
Chapter 5 Markov Chain Monte Carlo MCMC is a kind of improvement of the Monte Carlo method By sampling from a Markov chain whose stationary distribution is the desired sampling distributuion, it is possible
More informationA Search and Jump Algorithm for Markov Chain Monte Carlo Sampling. Christopher Jennison. Adriana Ibrahim. Seminar at University of Kuwait
A Search and Jump Algorithm for Markov Chain Monte Carlo Sampling Christopher Jennison Department of Mathematical Sciences, University of Bath, UK http://people.bath.ac.uk/mascj Adriana Ibrahim Institute
More informationTheory of Stochastic Processes 8. Markov chain Monte Carlo
Theory of Stochastic Processes 8. Markov chain Monte Carlo Tomonari Sei sei@mist.i.u-tokyo.ac.jp Department of Mathematical Informatics, University of Tokyo June 8, 2017 http://www.stat.t.u-tokyo.ac.jp/~sei/lec.html
More informationConnections between score matching, contrastive divergence, and pseudolikelihood for continuous-valued variables. Revised submission to IEEE TNN
Connections between score matching, contrastive divergence, and pseudolikelihood for continuous-valued variables Revised submission to IEEE TNN Aapo Hyvärinen Dept of Computer Science and HIIT University
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 informationSUPPLEMENT TO MARKET ENTRY COSTS, PRODUCER HETEROGENEITY, AND EXPORT DYNAMICS (Econometrica, Vol. 75, No. 3, May 2007, )
Econometrica Supplementary Material SUPPLEMENT TO MARKET ENTRY COSTS, PRODUCER HETEROGENEITY, AND EXPORT DYNAMICS (Econometrica, Vol. 75, No. 3, May 2007, 653 710) BY SANGHAMITRA DAS, MARK ROBERTS, AND
More informationHow to solve the stochastic partial differential equation that gives a Matérn random field using the finite element method
arxiv:1803.03765v1 [stat.co] 10 Mar 2018 How to solve the stochastic partial differential equation that gives a Matérn random field using the finite element method Haakon Bakka bakka@r-inla.org March 13,
More informationA Bayesian Spatio-Temporal Geostatistical Model with an Auxiliary Lattice for Large Datasets
Statistica Sinica (2013): Preprint 1 A Bayesian Spatio-Temporal Geostatistical Model with an Auxiliary Lattice for Large Datasets Ganggang Xu 1, Faming Liang 1 and Marc G. Genton 2 1 Texas A&M University
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 informationMultivariate modelling and efficient estimation of Gaussian random fields with application to roller data
Multivariate modelling and efficient estimation of Gaussian random fields with application to roller data Reinhard Furrer, UZH PASI, Búzios, 14-06-25 NZZ.ch Motivation Microarray data: construct alternative
More informationAn explicit link between Gaussian fields and Gaussian Markov random fields: The SPDE approach
An explicit link between Gaussian fields and Gaussian Markov random fields: The SPDE approach Finn Lindgren and Håvard Rue Department of Mathematical Sciences Norwegian University of Science and Technology
More informationNon-stationary Gaussian models with physical barriers
Non-stationary Gaussian models with physical barriers Haakon Bakka; in collaboration with Jarno Vanhatalo, Janine Illian, Daniel Simpson and Håvard Rue King Abdullah University of Science and Technology
More informationOrganization. I MCMC discussion. I project talks. I Lecture.
Organization I MCMC discussion I project talks. I Lecture. Content I Uncertainty Propagation Overview I Forward-Backward with an Ensemble I Model Reduction (Intro) Uncertainty Propagation in Causal Systems
More informationat least 50 and preferably 100 observations should be available to build a proper model
III Box-Jenkins Methods 1. Pros and Cons of ARIMA Forecasting a) need for data at least 50 and preferably 100 observations should be available to build a proper model used most frequently for hourly or
More informationINLA for Spatial Statistics
INLA for Spatial Statistics 5. Log-Gaussian Co Processes Daniel Simpson Department of Mathematical Sciences University of Bath Outline spatial point processes fitting log Gaussian Co processes eamples
More informationInversion Base Height. Daggot Pressure Gradient Visibility (miles)
Stanford University June 2, 1998 Bayesian Backtting: 1 Bayesian Backtting Trevor Hastie Stanford University Rob Tibshirani University of Toronto Email: trevor@stat.stanford.edu Ftp: stat.stanford.edu:
More informationBayesian Dynamic Linear Modelling for. Complex Computer Models
Bayesian Dynamic Linear Modelling for Complex Computer Models Fei Liu, Liang Zhang, Mike West Abstract Computer models may have functional outputs. With no loss of generality, we assume that a single computer
More informationProbabilistic Graphical Models
Probabilistic Graphical Models Brown University CSCI 2950-P, Spring 2013 Prof. Erik Sudderth Lecture 12: Gaussian Belief Propagation, State Space Models and Kalman Filters Guest Kalman Filter Lecture by
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 informationDisease mapping with Gaussian processes
EUROHEIS2 Kuopio, Finland 17-18 August 2010 Aki Vehtari (former Helsinki University of Technology) Department of Biomedical Engineering and Computational Science (BECS) Acknowledgments Researchers - Jarno
More informationAn introduction to Sequential Monte Carlo
An introduction to Sequential Monte Carlo Thang Bui Jes Frellsen Department of Engineering University of Cambridge Research and Communication Club 6 February 2014 1 Sequential Monte Carlo (SMC) methods
More informationTowards a Bayesian model for Cyber Security
Towards a Bayesian model for Cyber Security Mark Briers (mbriers@turing.ac.uk) Joint work with Henry Clausen and Prof. Niall Adams (Imperial College London) 27 September 2017 The Alan Turing Institute
More informationChapter 4 Dynamic Bayesian Networks Fall Jin Gu, Michael Zhang
Chapter 4 Dynamic Bayesian Networks 2016 Fall Jin Gu, Michael Zhang Reviews: BN Representation Basic steps for BN representations Define variables Define the preliminary relations between variables Check
More information