計算多變量常態的密度
在許多統計任務中，例如假設檢定、分群，以及概似值計算，你需要對指定的多變量常態分配計算其密度。在本練習中，你會使用 dmvnorm() 函式，針對先前產生的樣本 multnorm.sample 的每一筆觀測，在指定的平均向量與變異數-共變異數矩陣下，計算多變量常態的密度。
平均向量與變異數-共變異數矩陣已預先載入為物件 mu.sim 與 sigma.sim。
本練習屬於課程
R 的多變量機率分配
練習說明
- 使用
dmvnorm()，針對multnorm.sample中的 100 筆樣本，計算二變量常態的密度高度。 - 使用
scatterplot3d()，在每一個產生的樣本點上，繪製其密度高度的 3D 散佈圖。
動手互動練習
試著完成這個範例程式碼，體驗一下這個練習。
# Calculate density
multnorm.dens <- dmvnorm(multnorm.sample, mean = ___, sigma = ___)
# Create scatter plot of density heights
___(cbind(___),
color="blue", pch="", type = "h",
xlab = "x", ylab = "y", zlab = "density")