Fits the bivariate Gaussian common component geostatistical model
Computes maximum likelihood estimates of the bivariate Gaussian common component geostatistical model.
likfitBGCCM(geodata1, geodata2, ini.sigmasq, ini.phi,
cov0.model="matern", cov1.model="matern", cov2.model="matern",
kappa0=0.5, kappa1=0.5, kappa2=0.5,
fc.min = c("optim", "nlminb"), ...)geodata1 |
an object of the class |
geodata2 |
an object of the class |
ini.sigmasq |
optional, a vector with initial values for the variance parameters. If not provided default values are used. |
ini.phi |
optional, a vector with initial values for the correlation range parameters. If not provided default values are used. |
cov0.model, cov1.model, cov2.model |
covariance model for each of
the processes. See |
kappa0, kappa1, kappa2 |
extra parameter for some covariance models. |
fc.min |
a string indication which function should be used to minimise the negative of the log-likelihood. |
... |
A list with model fitting information to which
the class BGCCM is assigned.
mu |
a 2 elements vector with mean estimates. |
sigmasq |
a 4 elements vector with variance estimates. |
phi |
a 3 elements vector with estimated correlation parameters values. |
loglik |
a scalar. Maximised value of the log-likelihood. |
optim |
|
... |
and other information related to the model fitting. |
This is a new function and still in draft format and pretty much untested.
Paulo J. Ribeiro Jr. paulojus@leg.ufpr.br,
Peter J. Diggle p.diggle@lancaster.ac.uk.
# see http://www.leg.ufpr.br/geoR/tutorials/CCM.R
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.