Learn R Programming

GpGp (version 0.5.0)

matern45_isotropic: Isotropic Matern covariance function, smoothness = 4.5

Description

From a matrix of locations and covariance parameters of the form (variance, range, nugget), return the square matrix of all pairwise covariances.

Usage

matern45_isotropic(covparms, locs)

Value

A matrix with n rows and n columns, with the i,j entry containing the covariance between observations at locs[i,] and locs[j,].

Arguments

covparms

A vector with covariance parameters in the form (variance, range, nugget)

locs

A matrix with n rows and d columns. Each row of locs is a point in R^d.

Parameterization

The covariance parameter vector is (variance, range, nugget) = \((\sigma^2,\alpha,\tau^2)\), and the covariance function is parameterized as $$ M(x,y) = \sigma^2 ( \sum_{j=0}^4 c_j || x - y ||^j/ \alpha^j ) exp( - || x - y ||/ \alpha )$$ where c_0 = 1, c_1 = 1, c_2 = 3/7, c_3 = 2/21, c_4 = 1/105. The nugget value \( \sigma^2 \tau^2 \) is added to the diagonal of the covariance matrix. NOTE: the nugget is \( \sigma^2 \tau^2 \), not \( \tau^2 \).