Обчислення dmvnorm на сітці
Щоб візуалізувати поверхню двовимірної нормально розподіленої щільності, потрібно обчислити щільність на густій сітці координат x та y і скористатися 3D-функціями побудови графіків, як-от 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")