Calcularea densității distribuției normale multivariante

Pentru multe sarcini statistice – cum ar fi testarea ipotezelor, clusterizarea și calculul funcției de verosimilitate – este necesar să calculezi densitatea unei distribuții normale multivariante specificate. În acest exercițiu, vei folosi funcția dmvnorm() pentru a calcula densitățile normale multivariante cu media și matricea de varianță-covarianță specificate, la fiecare dintre observațiile din eșantionul generat anterior, multnorm.sample.

Media și matricea de varianță-covarianță sunt preîncărcate ca obiectele mu.sim și sigma.sim.

Acest exercițiu face parte din cursul

Distribuții de probabilitate multivariate în R

Vezi cursul

Instrucțiuni pentru exercițiu

Folosește dmvnorm() pentru a calcula înălțimile densității pentru cele 100 de eșantioane din multnorm.sample, în cazul unei distribuții normale bivariante.
Folosește scatterplot3d() pentru a reprezenta grafic un scatterplot 3D al înălțimilor densității în fiecare dintre punctele de eșantionare generate.

Exercițiu interactiv practic

Încearcă acest exercițiu completând acest cod de exemplu.

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

Editează și rulează codul