likfitBGCCM: Fits the bivariate Gaussian common component geostatistical...

In geoR: Analysis of Geostatistical Data

Fits the bivariate Gaussian common component geostatistical model

Description

Computes maximum likelihood estimates of the bivariate Gaussian common component geostatistical model.

Usage

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"), ...)

Arguments

geodata1	an object of the class geodata with the first variable.
geodata2	an object of the class geodata with the second variable.
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 cov.spatial for details.
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.
...	further arguments to be passed to optim or nlminb.

Value

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	results returned by optim or nlminb.
...	and other information related to the model fitting.

Warning

This is a new function and still in draft format and pretty much untested.

Author(s)

Paulo J. Ribeiro Jr. [e-mail],
Peter J. Diggle [e-mail].

See Also

optim, nlminb, varcovBGCCM, as.geodata, likfit.

Examples

# see http://www.leg.ufpr.br/geoR/tutorials/CCM.R

geoR documentation built on May 29, 2024, 1:36 a.m.