r - 如何用扩展椭圆绘制双变量正态分布

Question

如何使用扩展椭圆绘制二元正态分布并添加 5%、25%、50%、75% 和图中有 95% 的标签？谢谢!

Answer 1

您可以使用名为 mvtnorm 的 R 包创建等高线图。

假设您正在尝试绘制一个双变量正态分布，其中 mu_x = 1 且 mu_y = 1，方差矩阵为 c(2,1,1,1)。为 x、y、z 生成 100 个观测值。您可以为这种情况创建等高线图:

library(mvtnorm)
x.points <- seq(-3,3,length.out=100)
y.points <- x.points
z <- matrix(0,nrow=100,ncol=100)
mu <- c(1,1)
sigma <- matrix(c(2,1,1,1),nrow=2)
for (i in 1:100) {
   for (j in 1:100) {
    z[i,j] <- dmvnorm(c(x.points[i],y.points[j]),
    mean=mu,sigma=sigma)
   }
}
contour(x.points,y.points,z)