TDtT <- function(A, ...){ ch <- chol(A) dd <- diag(ch) return(list(T = t(drop0(zapsmall(ch / dd, ...))), D = Diagonal(nrow(A), dd^2))) }
Add the following code to your website.
REMOVE THIS Copy to clipboard
For more information on customizing the embed code, read Embedding Snippets.