在網格上計算 dmvnorm
若要視覺化二變量常態的密度曲面，需要在密集的 x 與 y 座標網格上計算密度，並使用像 persp() 這類的 3D 曲面繪圖函式來呈現。
以指令 mvals <- expand.grid(seq(-5, 10, length.out = 40), seq(-8, 4, length.out = 40)) 產生的 40×40 密集網格已為你預先載入。
本練習屬於課程
R 的多變量機率分配
練習說明
- 在物件
mvals的網格值上，計算平均為mu.sim、變異—共變異矩陣為sigma.sim的常態分配之密度。 - 使用
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")