二维高斯的 Plot 均值和轮廓

Question

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

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

我想 plot 一个点作为平均值（在 (1, 2) 处）及其周围我想 plot 一个对应于一个标准偏差的椭圆。 我该怎么做呢？

Answer 1

ellipse package 只是几个具有绘制椭圆功能的软件包之一。 安装 package 然后加载。 您必须计算标准偏差数量的概率水平，因为 function 不直接接受标准偏差：

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))