Description Usage Arguments Details Value Author(s) References See Also Examples
This function computes the value of the log-likelihood for a Gaussian random field.
1 2 3 4 5 |
geodata |
a list containing elements |
coords |
an n x 2 matrix, each row containing Euclidean
coordinates of the n data locations. By default it takes the
element |
data |
a vector with data values. By default it takes the
element |
obj.model |
a object of the class |
cov.model |
a string specifying the model for the correlation
function. For further details see
documentation for |
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:
|
method.lik |
options are |
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 |
The expression log-likelihood is:
l(theta) = -(n/2) * log(2 * pi) - 0.5 * log|V| - 0.5 * (y - F b)' V^{-1} (y - F b),
where n is the size of the data vector y, b is the mean (vector) parameter with dimention p, V is the covariance matrix and F is the matrix with the values of the covariates (a vector of 1's if the mean is constant.
The expression restricted log-likelihood is:
rl(θ) = -((n-p)/2) * log (2 * pi) + 0.5 * log |F'F| - 0.5 * log |V| - 0.5 * log |F'VF| - 0.5 * (y - Fb)' V^(-1) (y - Fb).
The numerical value of the log-likelihood.
Paulo Justiniano Ribeiro Jr. paulojus@leg.ufpr.br,
Peter J. Diggle p.diggle@lancaster.ac.uk.
Further information on the package geoR can be found at:
http://www.leg.ufpr.br/geoR.
likfit
for likelihood-based parameter estimation.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 | loglik.GRF(s100, cov.pars=c(0.8, .25), nugget=0.2)
loglik.GRF(s100, cov.pars=c(0.8, .25), nugget=0.2, met="RML")
## Computing the likelihood of a model fitted by ML
s100.ml <- likfit(s100, ini=c(1, .5))
s100.ml
s100.ml$loglik
loglik.GRF(s100, obj=s100.ml)
## Computing the likelihood of a variogram fitted model
s100.v <- variog(s100, max.dist=1)
s100.vf <- variofit(s100.v, ini=c(1, .5))
s100.vf
loglik.GRF(s100, obj=s100.vf)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.