多變量 t 分配的密度
在這個練習中,你要使用前一個練習所產生並儲存在 multt.sample 物件中的 200 個樣本,來計算多變量 t 分配的密度。
本練習屬於課程
R 的多變量機率分配
練習說明
- 對
multt.sample中的 200 個樣本,計算其在雙變量 t 分配下的密度高度。此分配的定位向量為mu.sim,變異—共變異矩陣為sigma.sim,自由度為 5。 - 建立 3D 散佈圖,以視覺化這 200 個點各自的密度高度。
動手互動練習
試著完成這個範例程式碼,體驗一下這個練習。
# Calculate densities
multt.dens <- dmvt(___)
# Plot 3D heights of densities
___(___(___, multt.dens),
color = "blue", pch = "", type = "h",
xlab = "x", ylab = "y", zlab = "density")