tensor.prod.model.matrix: Utility functions for constructing tensor product smooths

Description Usage Arguments Details Value Author(s) References See Also Examples

View source: R/mgcv.R

Description

Produce model matrices or penalty matrices for a tensor product smooth from the model matrices or penalty matrices for the marginal bases of the smooth.

Usage

1

Arguments

X

a list of model matrices for the marginal bases of a smooth

Details

If X[[1]], X[[2]] ... X[[m]] are the model matrices of the marginal bases of a tensor product smooth then the ith row of the model matrix for the whole tensor product smooth is given by X[[1]][i,]%x%X[[2]][i,]%x% ... X[[m]][i,], where %x% is the Kronecker product. Of course the routine operates column-wise, not row-wise!

Value

Either a single model matrix for a tensor product smooth, or a list of penalty terms for a tensor product smooth.

Author(s)

Simon N. Wood simon.wood@r-project.org

References

Wood, S.N. (2006) “Low Rank Scale Invariant Tensor Product Smooths for Generalized Additive Mixed Models”. Biometrics 62(4):1025-1036

See Also

te, smooth.construct.tensor.smooth.spec

Examples

1
2
X <- list(matrix(1:4,2,2),matrix(5:10,2,3))
tensor.prod.model.matrix(X)

Example output

Categorical Regression Splines (version 0.15-31)
[vignette("crs_faq") provides answers to frequently asked questions]
Warning messages:
1: In rgl.init(initValue, onlyNULL) : RGL: unable to open X11 display
2: 'rgl_init' failed, running with rgl.useNULL = TRUE 
3: .onUnload failed in unloadNamespace() for 'rgl', details:
  call: fun(...)
  error: object 'rgl_quit' not found 
     [,1] [,2] [,3] [,4] [,5] [,6]
[1,]    5    7    9   15   21   27
[2,]   12   16   20   24   32   40

crs documentation built on Feb. 2, 2021, 5:13 p.m.