liu.seSök publikationer i DiVA
Ändra sökning
RefereraExporteraLänk till posten
Permanent länk

Direktlänk
Referera
Referensformat
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • oxford
  • Annat format
Fler format
Språk
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Annat språk
Fler språk
Utmatningsformat
  • html
  • text
  • asciidoc
  • rtf
Scalable Bayesian spatial analysis with Gaussian Markov random fields
Linköpings universitet, Institutionen för datavetenskap, Statistik och maskininlärning. Linköpings universitet, Filosofiska fakulteten.ORCID-id: 0000-0002-5115-5657
2020 (Engelska)Doktorsavhandling, sammanläggning (Övrigt vetenskapligt)Alternativ titel
Skalbar Bayesiansk spatial analys med Gaussiska Markov-fält (Svenska)
Abstract [en]

Accurate statistical analysis of spatial data is important in many applications. Failing to properly account for spatial autocorrelation may often lead to false conclusions. At the same time, the ever-increasing sizes of spatial datasets pose a great computational challenge, as many standard methods for spatial analysis are limited to a few thousand data points.

In this thesis, we explore how Gaussian Markov random fields (GMRFs) can be used for scalable analysis of spatial data. GMRFs are closely connected to the commonly used Gaussian processes, but have sparsity properties that make them computationally cheap both in time and memory. The Bayesian framework enables a GMRF to be used as a spatial prior, comprising the assumption of smooth variation over space, and gives a principled way to estimate the parameters and propagate uncertainty.

We develop new algorithms that enable applying GMRF priors in 3D to the brain activity inherent in functional magnetic resonance imaging (fMRI) data, with millions of observations. We show that our methods are both faster and more accurate than previous work. A method for approximating selected elements of the inverse precision matrix (i.e. the covariance matrix) is also proposed, which is important for evaluating the posterior uncertainty. In addition, we establish a link between GMRFs and deep convolutional neural networks, which have been successfully used in countless machine learning tasks for images, resulting in a deep GMRF model. Finally, we show how GMRFs can be used in real-time robotic search and rescue operations, for modeling the spatial distribution of injured persons.

Abstract [sv]

Tillförlitlig statistisk analys av spatiala data är viktigt inom många tillämpningar. Om inte korrekt hänsyn tas till spatial autokorrelation kan det ofta leda till felaktiga slutsatser. Samtidigt ökar ständigt storleken på de spatiala datamaterialen vilket utgör en stor beräkningsmässig utmaning, eftersom många standardmetoder för spatial analys är begränsade till några tusental datapunkter.

I denna avhandling utforskar vi hur Gaussiska Markov-fält (eng: Gaussian Markov random fields, GMRF) kan användas för mer skalbara analyser av spatiala data. GMRF-modeller är nära besläktade med de ofta använda Gaussiska processerna, men har gleshetsegenskaper som gör dem beräkningsmässigt effektiva både vad gäller tids- och minnesåtgång. Det Bayesianska synsättet gör det möjligt att använda GMRF som en spatial prior som innefattar antagandet om långsam spatial variation och ger ett principiellt tillvägagångssätt för att skatta parametrar och propagera osäkerhet.

Vi utvecklar nya algoritmer som gör det möjligt att använda GMRF-priors i 3D för den hjärnaktivitet som indirekt kan observeras i hjärnbilder framtagna med tekniken fMRI, som innehåller milliontals datapunkter. Vi visar att våra metoder är både snabbare och mer korrekta än tidigare forskning. En metod för att approximera utvalda element i den inversa precisionsmatrisen (dvs. kovariansmatrisen) framförs också, vilket är viktigt för att kunna evaluera osäkerheten i posteriorn. Vidare gör vi en koppling mellan GMRF och djupa neurala faltningsnätverk, som har använts framgångsrikt för mängder av bildrelaterade problem inom maskininlärning, vilket mynnar ut i en djup GMRF-modell. Slutligen visar vi hur GMRF kan användas i realtid av autonoma drönare för räddningsinsatser i katastrofområden för att modellera den spatiala fördelningen av skadade personer.

Ort, förlag, år, upplaga, sidor
Linköping: Linköping University Electronic Press, 2020. , s. 53
Serie
Linköping Studies in Arts and Sciences, ISSN 0282-9800 ; 790Linköping Studies in Statistics, ISSN 1651-1700 ; 15
Nyckelord [en]
Spatial statistics, Bayesian statistics, Gaussian Markov random fields, fMRI, Machine learning
Nyckelord [sv]
Spatial statistik, Bayesiansk statistik, Gaussiska Markov-fält, fMRI, Maskininlärning
Nationell ämneskategori
Sannolikhetsteori och statistik
Identifikatorer
URN: urn:nbn:se:liu:diva-165872DOI: 10.3384/diss.diva-165872ISBN: 9789179298180 (tryckt)OAI: oai:DiVA.org:liu-165872DiVA, id: diva2:1433819
Disputation
2020-09-18, Ada Lovelace, B Building, Campus Valla, Linköping, 10:15 (Engelska)
Opponent
Handledare
Tillgänglig från: 2020-08-17 Skapad: 2020-06-01 Senast uppdaterad: 2020-11-24Bibliografiskt granskad
Delarbeten
1. Fast Bayesian whole-brain fMRI analysis with spatial 3D priors
Öppna denna publikation i ny flik eller fönster >>Fast Bayesian whole-brain fMRI analysis with spatial 3D priors
2017 (Engelska)Ingår i: NeuroImage, ISSN 1053-8119, E-ISSN 1095-9572, Vol. 146, s. 211-225Artikel i tidskrift (Refereegranskat) Published
Abstract [en]

Spatial whole-brain Bayesian modeling of task-related functional magnetic resonance imaging (fMRI) is a great computational challenge. Most of the currently proposed methods therefore do inference in subregions of the brain separately or do approximate inference without comparison to the true posterior distribution. A popular such method, which is now the standard method for Bayesian single subject analysis in the SPM software, is introduced in Penny et al. (2005b). The method processes the data slice-by-slice and uses an approximate variational Bayes (VB) estimation algorithm that enforces posterior independence between activity coefficients in different voxels. We introduce a fast and practical Markov chain Monte Carlo (MCMC) scheme for exact inference in the same model, both slice-wise and for the whole brain using a 3D prior on activity coefficients. The algorithm exploits sparsity and uses modern techniques for efficient sampling from high-dimensional Gaussian distributions, leading to speed-ups without which MCMC would not be a practical option. Using MCMC, we are for the first time able to evaluate the approximate VB posterior against the exact MCMC posterior, and show that VB can lead to spurious activation. In addition, we develop an improved VB method that drops the assumption of independent voxels a posteriori. This algorithm is shown to be much faster than both MCMC and the original VB for large datasets, with negligible error compared to the MCMC posterior.

Ort, förlag, år, upplaga, sidor
Elsevier, 2017
Nyckelord
fMRI, Spatial priors, Variational Bayes, Markov chain Monte Carlo, Gaussian Markov random fields, General linear model
Nationell ämneskategori
Medicinteknik
Identifikatorer
urn:nbn:se:liu:diva-132945 (URN)10.1016/j.neuroimage.2016.11.040 (DOI)000394560700019 ()27876654 (PubMedID)2-s2.0-84999622353 (Scopus ID)
Anmärkning

Funding agencies: Swedish Research Council (Vetenskapsradet) [20135229]; Knut and Alice Wallenberg Foundation [KAW 20012.0067]

Tillgänglig från: 2016-12-05 Skapad: 2016-12-05 Senast uppdaterad: 2020-06-29Bibliografiskt granskad
2. Efficient Covariance Approximations for Large Sparse Precision Matrices
Öppna denna publikation i ny flik eller fönster >>Efficient Covariance Approximations for Large Sparse Precision Matrices
2018 (Engelska)Ingår i: Journal of Computational And Graphical Statistics, ISSN 1061-8600, E-ISSN 1537-2715, Vol. 27, nr 4, s. 898-909Artikel i tidskrift (Refereegranskat) Published
Abstract [en]

The use of sparse precision (inverse covariance) matrices has become popular because they allow for efficient algorithms for joint inference in high-dimensional models. Many applications require the computation of certain elements of the covariance matrix, such as the marginal variances, which may be nontrivial to obtain when the dimension is large. This article introduces a fast Rao-Blackwellized Monte Carlo sampling-based method for efficiently approximating selected elements of the covariance matrix. The variance and confidence bounds of the approximations can be precisely estimated without additional computational costs. Furthermore, a method that iterates over subdomains is introduced, and is shown to additionally reduce the approximation errors to practically negligible levels in an application on functional magnetic resonance imaging data. Both methods have low memory requirements, which is typically the bottleneck for competing direct methods.

Ort, förlag, år, upplaga, sidor
AMER STATISTICAL ASSOC, 2018
Nyckelord
Gaussian Markov random fields; Selected inversion; Sparse precision matrix; Spatial analysis; Stochastic approximation
Nationell ämneskategori
Sannolikhetsteori och statistik
Identifikatorer
urn:nbn:se:liu:diva-153715 (URN)10.1080/10618600.2018.1473782 (DOI)000453029500018 ()
Anmärkning

Funding Agencies|Swedish Research Council (Vetenskapsradet) [2013-5229, 2016-04187]; European Unions Horizon 2020 Programme for Research and Innovation [640171]

Tillgänglig från: 2019-01-07 Skapad: 2019-01-07 Senast uppdaterad: 2020-06-29
3. Deep Gaussian Markov Random Fields
Öppna denna publikation i ny flik eller fönster >>Deep Gaussian Markov Random Fields
2020 (Engelska)Ingår i: Proceedings of the 37th International Conference on Machine Learning / [ed] Hal Daumé III, Aarti Singh, PMLR , 2020, Vol. 119, s. 8916-8926Konferensbidrag, Publicerat paper (Refereegranskat)
Abstract [en]

Gaussian Markov random fields (GMRFs) are probabilistic graphical models widely used in spatial statistics and related fields to model dependencies over spatial structures. We establish a formal connection between GMRFs and convolutional neural networks (CNNs). Common GMRFs are special cases of a generative model where the inverse mapping from data to latent variables is given by a 1-layer linear CNN. This connection allows us to generalize GMRFs to multi-layer CNN architectures, effectively increasing the order of the corresponding GMRF in a way which has favorable computational scaling. We describe how well-established tools, such as autodiff and variational inference, can be used for simple and efficient inference and learning of the deep GMRF. We demonstrate the flexibility of the proposed model and show that it outperforms the state-of-the-art on a dataset of satellite temperatures, in terms of prediction and predictive uncertainty.

Ort, förlag, år, upplaga, sidor
PMLR, 2020
Serie
Proceedings of Machine Learning Research, ISSN 2640-3498 ; 119
Nationell ämneskategori
Sannolikhetsteori och statistik
Identifikatorer
urn:nbn:se:liu:diva-171582 (URN)
Konferens
International Conference on Machine Learning, 13-18 July 2020, Virtual
Tillgänglig från: 2020-11-24 Skapad: 2020-11-24 Senast uppdaterad: 2020-11-24
4. Real-Time Robotic Search using Structural Spatial Point Processes
Öppna denna publikation i ny flik eller fönster >>Real-Time Robotic Search using Structural Spatial Point Processes
Visa övriga...
2020 (Engelska)Ingår i: 35TH UNCERTAINTY IN ARTIFICIAL INTELLIGENCE CONFERENCE (UAI 2019), Association For Uncertainty in Artificial Intelligence (AUAI) , 2020, Vol. 115, s. 995-1005Konferensbidrag, Publicerat paper (Refereegranskat)
Abstract [en]

Aerial robots hold great potential for aiding Search and Rescue (SAR) efforts over large areas, such as during natural disasters. Traditional approaches typically search an area exhaustively, thereby ignoring that the density of victims varies based on predictable factors, such as the terrain, population density and the type of disaster. We present a probabilistic model to automate SAR planning, with explicit minimization of the expected time to discovery. The proposed model is a spatial point process with three interacting spatial fields for i) the point patterns of persons in the area, ii) the probability of detecting persons and iii) the probability of injury. This structure allows inclusion of informative priors from e.g. geographic or cell phone traffic data, while falling back to latent Gaussian processes when priors are missing or inaccurate. To solve this problem in real-time, we propose a combination of fast approximate inference using Integrated Nested Laplace Approximation (INLA), and a novel Monte Carlo tree search tailored to the problem. Experiments using data simulated from real world Geographic Information System (GIS) maps show that the framework outperforms competing approaches, finding many more injured in the crucial first hours.

Ort, förlag, år, upplaga, sidor
Association For Uncertainty in Artificial Intelligence (AUAI), 2020
Serie
Proceedings of Machine Learning Research (PMLR), E-ISSN 2640-3498 ; 115
Nationell ämneskategori
Data- och informationsvetenskap
Identifikatorer
urn:nbn:se:liu:diva-159698 (URN)000722423500092 ()2-s2.0-85084016675 (Scopus ID)
Konferens
Proceedings of the 35th Conference on Uncertainty in Artificial Intelligence (UAI 2019), Tel Aviv, Israel, July 22-25, 2019
Anmärkning

Funding: Wallenberg AI, Autonomous Systems and Software Program (WASP); WASP Autonomous Research Arenas - Knut and Alice Wallenberg Foundation; Swedish Foundation for Strategic Research (SSF)Swedish Foundation for Strategic Research; ELLIIT Excellence Center at Link opingLund for Information Technology

Tillgänglig från: 2019-08-19 Skapad: 2019-08-19 Senast uppdaterad: 2023-04-05Bibliografiskt granskad

Open Access i DiVA

fulltext(6923 kB)1765 nedladdningar
Filinformation
Filnamn FULLTEXT01.pdfFilstorlek 6923 kBChecksumma SHA-512
0ef4fa9c69c5bc784b5d08e517fcf86cd91999e22fe4cf049bb477baf197c24ef78320b727a8be63d61f4d1273d5bd48dc39f9bebdf1412212610cf4693a533d
Typ fulltextMimetyp application/pdf
Beställ online >>

Övriga länkar

Förlagets fulltext

Person

Sidén, Per

Sök vidare i DiVA

Av författaren/redaktören
Sidén, Per
Av organisationen
Statistik och maskininlärningFilosofiska fakulteten
Sannolikhetsteori och statistik

Sök vidare utanför DiVA

GoogleGoogle Scholar
Totalt: 1778 nedladdningar
Antalet nedladdningar är summan av nedladdningar för alla fulltexter. Det kan inkludera t.ex tidigare versioner som nu inte längre är tillgängliga.

doi
isbn
urn-nbn

Altmetricpoäng

doi
isbn
urn-nbn
Totalt: 3439 träffar
RefereraExporteraLänk till posten
Permanent länk

Direktlänk
Referera
Referensformat
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • oxford
  • Annat format
Fler format
Språk
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Annat språk
Fler språk
Utmatningsformat
  • html
  • text
  • asciidoc
  • rtf