Nothing

```
rmvsn <- function(n, D, Mu, Sigma){
if(!is.matrix(D) || nrow(D) != ncol(D))
stop("'D' has to be a square matrix")
if(any(D[upper.tri(D)] != 0) | any(D[lower.tri(D)] != 0))
stop("'D' has to be a diagnonal matrix")
if(!is.vector(Mu) || length(Mu)!=nrow(D))
stop("'Mu' has to a vector with the length equal to the number of rows of 'D'.")
if(! isSymmetric(Sigma))
stop("'Sigma' has to be a symmetric matrix.")
if(! all(eigen(Sigma)$values > 0))
stop("'Sigma' has to be positive definite.")
if(nrow(Sigma) != nrow(D))
stop("The dimension of 'Sigma' does not match that of 'D'.")
p <- nrow(D)
Y <- matrix(0, n, p)
u <- runif(n*p, 1/2, 1)
z <- qnorm(u, 0, 1)
Z <- matrix(z, n, p)
for(i in 1:n){
Y[i,] <- mvrnorm(1, Mu+D%*%Z[i,], Sigma)
}
Y
}
```

**Any scripts or data that you put into this service are public.**

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.