Matérn covariance function

In statistics, the Matérn covariance (named after the Swedish forestry statistician Bertil Matérn[1]) is a covariance function used in spatial statistics, geostatistics, machine learning, image analysis, and other applications of multivariate statistical analysis on metric spaces. It is commonly used to define the statistical covariance between measurements made at two points that are d units distant from each other. Since the covariance only depends on distances between points, it is stationary. If the distance is Euclidean distance, the Matérn covariance is also isotropic.

The Matérn covariance between two points separated by d distance units is given by [2]


C(d) = \sigma^2\frac{2^{1-\nu}}{\Gamma(\nu)}\Bigg(\sqrt{2\nu}\frac{d}{\rho}\Bigg)^\nu K_\nu\Bigg(\sqrt{2\nu}\frac{d}{\rho}\Bigg),

where Γ is the gamma function, Kν is the modified Bessel function of the second kind, and ρ and ν are non-negative parameters of the covariance.

A Gaussian process with Matérn covariance has sample paths that are \lceil \nu-1 \rceil times differentiable.[3] As \nu\rightarrow\infty, the Matérn covariance converges to the squared exponential covariance function


C(d) = \sigma^2\exp(-d^2/ 2\rho^2).

When  \nu = 1/2 , the Matérn covariance is identical to the exponential covariance function. In fact,[3]

C(d) = \sigma^2 \exp \Bigg(-\frac{d}{\rho} \Bigg) \quad \quad \nu= \tfrac{1}{2},
C(d) = \sigma^2 \Bigg(1 + \frac{ \sqrt{3}d }{\rho} \Bigg) \exp \Bigg(-\frac{\sqrt{3}d}{\rho} \Bigg) \quad \quad \nu= \tfrac{3}{2},
C(d) = \sigma^2 \Bigg(1 + \frac{ \sqrt{5}d }{\rho} +\frac{ 5d^2}{3 \rho^2 }   \Bigg) \exp \Bigg(-\frac{\sqrt{5}d}{\rho} \Bigg) \quad \quad \nu= \tfrac{5}{2}.

See also

References

  1. Minasny, B.; McBratney, A. B. (2005). "The Matérn function as a general model for soil variograms". Geoderma 128 (3–4): 192–207. doi:10.1016/j.geoderma.2005.04.003.
  2. Rasmussen, Carl Edward (2006) Gaussian Processes for Machine Learning
  3. 1 2 Rasmussen, Carl Edward (2006) Gaussian Processes Covariance Functions and Classification. Presentation at Gaussian Processes in Practice
This article is issued from Wikipedia - version of the Friday, February 12, 2016. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.