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

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.

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.

Author

Paulo J. Ribeiro Jr. paulojus@leg.ufpr.br,
Peter J. Diggle p.diggle@lancaster.ac.uk.

Warning

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

Examples

Run this code

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

Run the code above in your browser using DataLab