# NOT RUN {
# Simulation and density evaluation for p = 2
Lambda <- diag(c(5, 1))
n <- 1e3
x <- r_ACG(n = n, Lambda = Lambda)
col <- viridisLite::viridis(n)
r <- runif(n, 0.95, 1.05) # Radius perturbation to improve visualization
plot(r * x, pch = 16, col = col[rank(d_ACG(x = x, Lambda = Lambda))])
# Simulation and density evaluation for p = 3
Lambda <- rbind(c(5, 1, 0.5),
c(1, 2, 1),
c(0.5, 1, 1))
x <- r_ACG(n = n, Lambda = Lambda)
rgl::plot3d(x, col = col[rank(d_ACG(x = x, Lambda = Lambda))], size = 5)
# }
Run the code above in your browser using DataLab