我有一个散点图，我想绘制一个椭圆边界来识别异常值：
library(MASS)
mu <- c(7, 5);
sigma <- matrix(c(1, 1/2, 1/2, 1), 2);
nu <- 4
n <- 4e3
y <- t(t(mvrnorm(n, rep(0, length(mu)), sigma) / sqrt(nu / rchisq(n, nu))) + mu)
plot(y[,1],y[,2])
我可以绘制分位数线：
abline(v=quantile(y[,1],c(0.05,0.95)))
abline(h=quantile(y[,2],c(0.05,0.95)))
但是如何使用Copulas改善这一点？ 谢谢你的帮助。