2
$\begingroup$

Suppose we want to model a stochastic process using Gaussian processes. We have data on $z$ (dependent variable) at some spatial points $(x_{i},y_{i})$. If our dataset is large then calculating the covariance matrix of the Gaussian process becomes difficult and that’s why approximation methods like using SPDE+FEA (Stochastic Partial Differential Equations + Finite Element Analysis) are suggested. Lindgren et al. (2011) in his paper explains how we can approximate a Gaussian process with Matern covariance function. He suggests that the solution to the below SPDE

$\kappa^{2} f-\Delta f=\epsilon / \tau$

is the Gaussian process we are looking for and the SPDE can be solved using a finite element method by using mesh and basis functions over the mesh to approximate the solution. I tried reading his paper and many other resources but only got some vague ideas about the approach and the big missing piece is that how data comes into play? I don’t see where the data is being used in the finite element method to solve the SPDE. Can someone clarify what’s happening under the hood with minimum use of sophisticated calculus and advanced Linear algebra?

$\endgroup$

1 Answer 1

1
$\begingroup$

What you're trying to do has been published last year in this paper (10.1016/j.cma.2020.113533 "The statistical finite element method (statFEM) for coherent synthesis of observation data and model predictions", a preprint: https://arxiv.org/abs/1905.06391) by Girolami et al. The statFEM approach works by applying Bayesian statistics to incorporate data: The stochastic FEM, modeled as a GP as you stated, serves as a Bayesian prior. Using the likelihood of the data, also modeled Gaussian, you are able to create a posterior GP which can be considered an updated FEM solution given your data. Unfortunately there is no code snippet or a single equation I can show you because the whole approach is rather lengthy. You can nevertheless have a look at my github where I have a working solver. https://github.com/herluc/Masterarbeit/blob/master/MA_LucasHermann/Python/statFEMSolver.py

$\endgroup$
2
  • $\begingroup$ Hi, welcome to CV. Please add full ref of your paper in case link dies in the future $\endgroup$
    – Antoine
    Jan 13, 2022 at 16:34
  • $\begingroup$ Hey, sure, done! $\endgroup$ Jan 13, 2022 at 17:39

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.