loglik.GRF: Log-likelihood

Description Usage Arguments See Also

Description

Same as geoR's loglik.GRF but taking into account non-Euclidean distances if pertinent.

Usage

1
2
3
4
loglik.GRF(geodata, coords = geodata$coords, data = geodata$data,
  obj.model = NULL, cov.model = "exp", cov.pars, nugget = 0,
  kappa = 0.5, lambda = 1, psiR = 1, psiA = 0, trend = "cte",
  method.lik = "ML", compute.dists = TRUE, realisations = NULL)

Arguments

geodata

a list containing elements coords and data as described next. Typically an object of the class "geodata" - a geoR data-set. If not provided the arguments coords and data must be provided instead.

coords

an n x 2 matrix, each row containing Euclidean coordinates of the n data locations. By default it takes the element coords of the argument geodata.

data

a vector with data values. By default it takes the element data of the argument geodata.

obj.model

a object of the class variomodel with a fitted model. Tipically an output of likfit or variofit.

cov.model

a string specifying the model for the correlation function. For further details see documentation for cov.spatial.

cov.pars

a vector with 2 elements with values of the covariance parameters sigma^2 (partial sill) and phi (range parameter).

nugget

value of the nugget parameter. Defaults to 0.

kappa

value of the smoothness parameter. Defaults to 0.5.

lambda

value of the Box-Cox tranformation parameter. Defaults to 1.

psiR

value of the anisotropy ratio parameter. Defaults to 1, corresponding to isotropy.

psiA

value (in radians) of the anisotropy rotation parameter. Defaults to zero.

trend

specifies the mean part of the model. The options are: "cte" (constant mean), "1st" (a first order polynomial on the coordinates), "2nd" (a second order polynomial on the coordinates), or a formula of the type ~X where X is a matrix with the covariates (external trend). Defaults to "cte".

method.lik

options are "ML" for likelihood and "REML" for restricted likelihood. Defaults to "ML".

compute.dists

for internal use with other function. Don't change the default unless you know what you are doing.

realisations

optional. A vector indicating replication number for each data. For more details see as.geodata.

See Also

loglik.GRF


famuvie/geoRcb documentation built on May 16, 2019, 10:04 a.m.