# Univariate example
H_uni <- HermiteCoeff(Type = "Univariate", N = 5)
# Multivariate example
N <- 5; d <- 3
H_multi <- HermiteCoeff(Type = "Multivariate", N = N, d = d)
X <- c(1:3)
X3 <- kronecker(X, kronecker(X, X))
X5 <- kronecker(X3, kronecker(X, X))
Idv <- as.vector(diag(d)) # vector of variance matrix
# value of H5 at X is
vH5 <- H_multi[[1]] %*% X5 + H_multi[[2]] %*% kronecker(Idv, X3) +
H_multi[[3]] %*% kronecker(kronecker(Idv, Idv), X)
Run the code above in your browser using DataLab