Beräkna dmvnorm över ett rutnät

För att visualisera en bivariat normalfördelnings täthetssyta behöver du beräkna tätheten över ett tätt rutnät av x- och y-koordinater och sedan använda 3D-ytplotsfunktioner som persp() för att visualisera ytan.

Det 40 × 40 täta rutnät som genereras med kommandot mvals <- expand.grid(seq(-5, 10, length.out = 40), seq(-8, 4, length.out = 40)) är redan inläst åt dig.

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

Multivariata sannolikhetsfördelningar i R

Visa kurs

Övningsinstruktioner

Beräkna tätheterna för en normalfördelning med medelvärdet mu.sim och varians-kovariansmatrisen sigma.sim för rutnätsvärdena i objektet mvals.
Använd plotfunktionen persp() för att visualisera täthetssytan över rutnätet.

Interaktiv övning med praktiskt arbete

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

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

Redigera och kör kod