グリッド上での dmvnorm の計算
二変量正規分布の密度曲面を可視化するには、x と y の座標を細かいグリッド上で計算し、persp() のような3次元曲面プロット関数を使って描画する必要があります。
mvals <- expand.grid(seq(-5, 10, length.out = 40), seq(-8, 4, length.out = 40)) というコマンドで生成された 40×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")