Beräkna tätheten för multivariat normalfördelning

Många statistiska uppgifter – som hypotestestning, klustring och likelihoodberäkning – kräver att du beräknar tätheten för en given multivariat normalfördelning. I den här övningen använder du funktionen dmvnorm() för att beräkna multivariata normaltätheter med angivet medelvärde och angiven varians-kovariansmatris för varje observation i ditt tidigare genererade sampel multnorm.sample.

Medelvärdet och varians-kovariansmatrisen är fördefinierade som objekten mu.sim och sigma.sim.

Den här övningen är en del av kursen

Multivariata sannolikhetsfördelningar i R

Visa kurs

Övningsinstruktioner

Använd dmvnorm() för att beräkna densitetshöjderna för de 100 samplen i multnorm.sample för en bivariat normalfördelning.
Använd scatterplot3d() för att plotta ett 3D-spridningsdiagram över densitetshöjderna i varje genererad sampelpunkt.

Interaktiv övning med praktiskt arbete

Testa den här övningen genom att slutföra den här exempelkoden.

# 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")

Redigera och kör kod