绘制二维高斯的均值和等高线。

Question

假设我有一个均值和方差-协方差矩阵为一定的高斯数

mu <- c(1, 2)
Sigma <- matrix(c(2, 0.5, 
                  0.5, 1), nrow=2)

我想绘制一个平均值的点（在(1, 2)），并在其周围绘制一个对应于一个标准差的椭圆。我应该怎样做呢？

Answer 1

该 ellipse 包是几个包中的一个，它有绘制椭圆的功能。安装这个包，然后加载它。你必须计算标准差数的概率水平，因为函数不直接接受标准差。

library(ellipse)
plot(mu[1], mu[2], xlim=c(-2, 4), ylim=c(0, 4))
lvl <- diff(pnorm(c(-1, 1)))
lines(ellipse(Sigma, scale=c(1.4, 1), centre=c(1, 2), level=lvl))