powered by
mvnradiusdensity(cov, radii)
temp <- mvngen(c(0, 0), diag(2), 5) R <- chol(solve(diag(2))) # trivial here, as Cholesky R of I(-1) is I temp <- temp <!-- %*% t(R) %*% R %*% t(temp) --> mvnradiusdensity(diag(2), diag(temp))
Run the code above in your browser using DataLab