getDesignMatrix {cocons}

Create an efficient design matrix based on a list of aspect models

Description

Creates a unique design matrix based on model specification for each of the different potentially spatially varying aspects.

Usage

getDesignMatrix(model.list, data)

Arguments

model.list

a list of formulas, one for each source of nonstationarity, specifying the models.

data

a data.frame

Value

a list() with two elements: a design matrix of dimension (n x p), and a par.pos object, indexing columns of the design matrix refered to each aspect models.

Author(s)

Federico Blasi

Examples

model.list <- list(   "mean" = 0,
                      "std.dev" = as.formula(" ~ 1 + lati_s * long_s"),
                      "scale" = as.formula(" ~ 1 + elev_s"),
                      "aniso" = as.formula(" ~ 1 + elev_s"),
                      "tilt" = as.formula(" ~ 1 + elev_s"),
                      "smooth" = as.formula(" ~ 1"),
                      "nugget" = -Inf)