Menghitung densitas normal multivariat
Untuk banyak tugas statistik, seperti pengujian hipotesis, pengelompokan (clustering), dan perhitungan likelihood, Anda perlu menghitung densitas dari suatu distribusi normal multivariat tertentu. Pada latihan ini, Anda akan menggunakan fungsi dmvnorm() untuk menghitung densitas normal multivariat dengan mean dan matriks varians-kovarians yang ditentukan pada setiap observasi dari sampel multnorm.sample yang sebelumnya Anda hasilkan.
Mean dan matriks varians-kovarians sudah dimuat untuk Anda sebagai objek mu.sim dan sigma.sim.
Latihan ini merupakan bagian dari kursus
Distribusi Probabilitas Multivariat di R
Instruksi latihan
- Gunakan
dmvnorm()untuk menghitung ketinggian densitas dari 100 sampel dalammultnorm.sampleuntuk normal bivariat. - Gunakan
scatterplot3d()untuk membuat scatterplot 3D dari ketinggian densitas pada setiap titik sampel yang dihasilkan.
Latihan interaktif langsung praktik
Cobalah latihan ini dengan melengkapi kode contoh ini.
# 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")