rmvnormp <- function(n,m,Q) { require(Matrix) Q=chol(Q) w=solve(t(Q),m) mu=solve(Q,w) z=rnorm(length(m)) v=solve(Q,z) beta=mu+v return(beta) }
Add the following code to your website.
REMOVE THIS Copy to clipboard
For more information on customizing the embed code, read Embedding Snippets.