Calcul de la densité normale multivariée
Pour de nombreuses tâches statistiques, comme les tests d'hypothèses, le regroupement (clustering) et le calcul de vraisemblance, vous devez calculer la densité d'une loi normale multivariée donnée. Dans cet exercice, vous utiliserez la fonction dmvnorm() pour calculer les densités normales multivariées, avec une moyenne et une matrice de variance-covariance spécifiées, à chacune des observations de votre échantillon multnorm.sample généré auparavant.
La moyenne et la matrice de variance-covariance sont déjà chargées pour vous dans les objets mu.sim et sigma.sim.
Cette activité fait partie du cours
Lois de probabilité multivariée en R
Instructions de l’exercice
- Utilisez
dmvnorm()pour calculer les hauteurs de densité des 100 échantillons dansmultnorm.samplepour une normale bivariée. - Utilisez
scatterplot3d()pour tracer un nuage de points 3D des hauteurs de densité à chacun des points d'échantillon générés.
Exercice interactif pratique
Essayez cet exercice en complétant ce code d’exemple.
# 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")