cocons (version 0.1.5)

getCovMatrix: Covariance matrix for "coco" class

Description

Compute the covariance matrix of coco.object.

Usage

getCovMatrix(coco.object, type = 'global', index = NULL)

Value

(matrix or S4) a n x n covariance matrix (for 'dense' coco objects) or a S4 spam object (for 'sparse' coco objects).

Arguments

coco.object

(S4) a fitted coco() object.

type

(character) whether 'global' to retrieve the regular covariance matrix, or 'local' to retrieve global covariance. based on the local aspects of a specific location (not implemented yet).

index

(integer) index to perform local covariance matrix (not implemented yet).

Author

Federico Blasi

Examples

if (FALSE) {
model.list <- list('mean' = 0,
                   'std.dev' = formula( ~ 1 + cov_x + cov_y),
                   'scale' = formula( ~ 1 + cov_x + cov_y),
                   'aniso' = 0,
                   'tilt' = 0,
                   'smooth' = 3/2,
                   'nugget' = -Inf)
                   
coco_object <- coco(type = 'dense',
                    data = holes[[1]][1:100,],
                    locs = as.matrix(holes[[1]][1:100,1:2]),
                    z = holes[[1]][1:100,]$z,
                    model.list = model.list)
                    
optim_coco <- cocoOptim(coco_object,
boundaries = getBoundaries(coco_object,
lower.value = -3, 3))

getCovMatrix(optim_coco)

}