Nothing
root.matrix <- function(X)
{
if((ncol(X) == 1L)&&(nrow(X) == 1L)) return(sqrt(X))
else
{
X.eigen <- eigen(X, symmetric=TRUE)
if(any(X.eigen$values < 0)) stop("matrix is not positive semidefinite")
sqomega <- sqrt(diag(X.eigen$values))
V <- X.eigen$vectors
V <- V %*% sqomega %*% t(V)
dimnames(V) <- dimnames(X)
return(V)
}
}
solveCrossprod <- function(X, method = c("qr", "chol", "solve")) {
switch(match.arg(method),
"qr" = chol2inv(qr.R(qr(X))),
"chol" = chol2inv(chol(crossprod(X))),
"solve" = solve(crossprod(X)))
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.