Вычисление dmvnorm на сетке значений
Чтобы визуализировать поверхность плотности двумерного нормального распределения, необходимо вычислить плотность на плотной сетке координат x и y, а затем использовать функции построения трёхмерных поверхностей, например persp().
Сетка размером 40 на 40, созданная командой mvals <- expand.grid(seq(-5, 10, length.out = 40), seq(-8, 4, length.out = 40)), уже загружена для вас.
Это упражнение является частью курса
Многомерные вероятностные распределения в R
Инструкции к упражнению
- Вычислите значения плотности нормального распределения с вектором средних
mu.simи матрицей дисперсий-ковариацийsigma.simдля всех точек сетки в объектеmvals. - Используйте функцию
persp(), чтобы визуализировать поверхность плотности на этой сетке.
Интерактивное практическое упражнение
Попробуйте выполнить это упражнение, дополнив этот пример кода.
# Calculate density over the specified grid
mvds <- ___
matrix_mvds <- matrix(___, nrow = 40)
# Create a perspective plot
persp(___, theta = 80, phi = 30, expand = 0.6, shade = 0.2, col = "lightblue", xlab = "x", ylab = "y", zlab = "dens")