CommencerCommencez gratuitement

Calculer dmvnorm sur une grille

Pour visualiser une surface de densité normale bivariée, il faut calculer la densité sur une grille dense de coordonnées x et y, puis utiliser des fonctions de tracé 3D comme persp() pour afficher la surface.

La grille dense 40 par 40 générée par la commande mvals <- expand.grid(seq(-5, 10, length.out = 40), seq(-8, 4, length.out = 40)) est préchargée pour vous.

Cet exercice fait partie du cours

<cours>Lois de probabilité multivariées en R</cours>
Voir le cours

Instructions de l’exercice

  • Calculez les densités d'une loi normale de moyenne mu.sim et de matrice de variance-covariance sigma.sim sur les valeurs de la grille contenues dans l'objet mvals.
  • Utilisez la fonction de tracé persp() pour visualiser la surface de densité sur la grille.

Exercice interactif pratique

Essayez cet exercice en complétant ce code d’exemple.

# Calculate density over the specified grid
mvds <- ___
matrix_mvds <-  matrix(___, nrow = 40)

# Create a perspective plot
persp(___, theta = 80, phi = 30, expand = 0.6, shade = 0.2, col = "lightblue", xlab = "x", ylab = "y", zlab = "dens")
Modifier et exécuter le code