dmvnorm berekenen over een raster

Om een bivariaat normale dichtheidsoppervlakte te visualiseren, bereken je de dichtheid over een dicht raster van x- en y-coördinaten, en gebruik je 3D-oppervlakteplotfuncties zoals persp() om het oppervlak te tonen.

Het 40 bij 40 dichte raster dat wordt gegenereerd met het commando mvals <- expand.grid(seq(-5, 10, length.out = 40), seq(-8, 4, length.out = 40)) is alvast voor je geladen.

Deze oefening maakt deel uit van de cursus

Multivariate kansverdelingen in R

Cursus bekijken

Oefeninstructies

Bereken de dichtheden van een normale verdeling met gemiddelde mu.sim en variantie-covariantiematrix sigma.sim op de rasterwaarden in het object mvals.
Gebruik de plotfunctie persp() om het dichtheidsoppervlak over het raster te visualiseren.

Praktische interactieve oefening

Probeer deze oefening eens door deze voorbeeldcode in te vullen.

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

Code bewerken en uitvoeren