多変量t分布の密度
この演習では、前の演習で生成して multt.sample オブジェクトに格納した200サンプルを使い、多変量t分布の密度を計算します。
この演習はコースの一部です
Rによる多変量確率分布
演習の手順
multt.sampleに含まれる200サンプルについて、位置パラメータmu.sim、分散共分散行列sigma.sim、自由度5の二変量t分布における密度の高さを計算します。- 200点それぞれの密度の高さを可視化するために、3次元散布図を作成します。
実践的なインタラクティブ演習
このサンプルコードを完成させて、この演習に挑戦してみましょう。
# Calculate densities
multt.dens <- dmvt(___)
# Plot 3D heights of densities
___(___(___, multt.dens),
color = "blue", pch = "", type = "h",
xlab = "x", ylab = "y", zlab = "density")