data(mm)
xpx <- crossprod(mm)
str(xpx)
str(ch <- chol(xpx))
## Visualize the sparseness:
dq <- function(ch) paste('"',ch,'"', sep="") ## dQuote(<UTF-8>) gives bad plots
image(xpx, main=paste("crossprod(mm) : Sparse", dq(class(xpx))))
image(ch, main= paste("chol (crossprod(mm)) : Sparse", dq(class(ch))))
Run the code above in your browser using DataLab