xfit <- function(T, P, xmeans = rep(0, dim(P)[1])) {
T <- .matrix(T)
P <- .matrix(P)
.center(tcrossprod(T, P), -xmeans)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.