assist (version 3.1.9)

rk.prod: Calculate product of reproducing kernels

Description

Return a matix as product of reproducing kernels

Usage

rk.prod(x, ...)

Value

a matrix as the product of reproducing kernels. If one argument is a vector, a kron

kernel is constructed first.

Arguments

x

a matrix evaluating a reproducing kernel, or a vector.

...

optional lists of matrices evaluating reproducing kernels or vectors. All matrics must have the same dimensions. All vectors must have the same length. The length of each vector must equal to the column and row numbers of each matrix.

Author

Chunlei Ke and Yuedong Wang

Details

The product of reproducing kernels is agian a reproducing kernel. In SS ANOVA, product of reproduing kernels is often used to model interaction spline terms.

References

Gu, C. and Wahba, G. (1993a). Smoothing Spline ANOVA with component-wise Bayesian confidence intervals. Journal of Computational and Graphical Statistics 55, 353--368.

Gu, C. and Wahba, G. (1993b). Semiparametric analysis of variance with tensor product thin plate splines. JRSS B 55, 353--368.

See Also

kron, ssr

Examples

if (FALSE) {
x1<- 1:10/10
x2<- runif(10)
rk.prod(cubic(x1), periodic(x2))
}