Среднее изображение и контур 2D гауссовского

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