gp-class | R Documentation |
gp
classThis is an S4 class definition for gp
in the GaSP
package.
formula
an object of formula
class that specifies regressors; see formula
for details.
output
a numerical vector including observations or outputs in a GaSP
input
a matrix including inputs in a GaSP
param
a list including values for regression parameters, correlation parameters, and nugget variance parameter. The specification of param should depend on the covariance model.
The regression parameters are denoted by coeff. Default value is \mathbf{0}
.
The marginal variance or partial sill is denoted by sig2. Default value is 1.
The nugget variance parameter is denoted by nugget for all covariance models. Default value is 0.
For the Confluent Hypergeometric class, range is used to denote the range parameter \beta
.
tail is used to denote the tail decay parameter \alpha
. nu is used to denote the
smoothness parameter \nu
.
For the generalized Cauchy class, range is used to denote the range parameter \phi
.
tail is used to denote the tail decay parameter \alpha
. nu is used to denote the
smoothness parameter \nu
.
For the Matérn class, range is used to denote the range parameter \phi
.
nu is used to denote the smoothness parameter \nu
. When \nu=0.5
, the
Matérn class corresponds to the exponential covariance.
For the powered-exponential class, range is used to denote the range parameter \phi
.
nu is used to denote the smoothness parameter. When \nu=2
, the powered-exponential class
corresponds to the Gaussian covariance.
cov.model
a list of two strings: family, form, where family indicates the family of covariance functions including the Confluent Hypergeometric class, the Matérn class, the Cauchy class, the powered-exponential class. form indicates the specific form of covariance structures including the isotropic form, tensor form, automatic relevance determination form.
The Confluent Hypergeometric correlation function is given by
C(h) = \frac{\Gamma(\nu+\alpha)}{\Gamma(\nu)}
\mathcal{U}\left(\alpha, 1-\nu, \left(\frac{h}{\beta}\right)^2\right),
where \alpha
is the tail decay parameter. \beta
is the range parameter.
\nu
is the smoothness parameter. \mathcal{U}(\cdot)
is the confluent hypergeometric
function of the second kind. Note that this parameterization of the CH covariance is different from the one in Ma and Bhadra (2023). For details about this covariance,
see Ma and Bhadra (2023; \Sexpr[results=rd]{tools:::Rd_expr_doi("10.1080/01621459.2022.2027775")}).
The generalized Cauchy covariance is given by
C(h) = \left\{ 1 + \left( \frac{h}{\phi} \right)^{\nu}
\right\}^{-\alpha/\nu},
where \phi
is the range parameter. \alpha
is the tail decay parameter.
\nu
is the smoothness parameter with default value at 2.
The Matérn correlation function is given by
C(h)=\frac{2^{1-\nu}}{\Gamma(\nu)} \left( \frac{h}{\phi} \right)^{\nu}
\mathcal{K}_{\nu}\left( \frac{h}{\phi} \right),
where \phi
is the range parameter. \nu
is the smoothness parameter.
\mathcal{K}_{\nu}(\cdot)
is the modified Bessel function of the second kind of order \nu
.
The exponential correlation function is given by
C(h)=\exp(-h/\phi),
where \phi
is the range parameter. This is the Matérn correlation with \nu=0.5
.
The Matérn correlation with \nu=1.5
.
The Matérn correlation with \nu=2.5
.
The powered-exponential correlation function is given by
C(h)=\exp\left\{-\left(\frac{h}{\phi}\right)^{\nu}\right\},
where \phi
is the range parameter. \nu
is the smoothness parameter.
The Gaussian correlation function is given by
C(h)=\exp\left(-\frac{h^2}{\phi^2}\right),
where \phi
is the range parameter.
This indicates the isotropic form of covariance functions. That is,
C(\mathbf{h}) = C^0(\|\mathbf{h}\|; \boldsymbol \theta),
where \| \mathbf{h}\|
denotes the
Euclidean distance or the great circle distance for data on sphere. C^0(\cdot)
denotes
any isotropic covariance family specified in family.
This indicates the tensor product of correlation functions. That is,
C(\mathbf{h}) = \prod_{i=1}^d C^0(|h_i|; \boldsymbol \theta_i),
where d
is the dimension of input space. h_i
is the distance along the i
th input dimension. This type of covariance structure has been often used in Gaussian process emulation for computer experiments.
This indicates the automatic relevance determination form. That is,
C(\mathbf{h}) = C^0\left(\sqrt{\sum_{i=1}^d\frac{h_i^2}{\phi^2_i}}; \boldsymbol \theta \right),
where \phi_i
denotes the range parameter along the i
th input dimension.
smooth.est
a logical value. If it is TRUE
, the smoothness parameter
will be estimated; otherwise the smoothness is not estimated.
dtype
a string indicating the type of distance:
Euclidean distance is used. This is the default choice.
Great circle distance is used for data on sphere.
loglik
a numerical value containing the log-likelihood with current
gp
object.
mcmc
a list containing MCMC samples if available.
prior
a list containing tuning parameters in prior distribution. This is used only if a Bayes estimation method with informative priors is used.
proposal
a list containing tuning parameters in proposal distribution. This is used only if a Bayes estimation method is used.
info
a list containing the maximum distance in the input space. It should be a vector if isotropic covariance is used, otherwise it is vector of maximum distances along each input dimension
Pulong Ma mpulong@gmail.com
GPBayes-package, GaSP
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.