Calcolare dmvnorm su una griglia
Per visualizzare una superficie di densità normale bivariata, devi calcolare la densità su una griglia fitta di coordinate x e y e usare funzioni di grafico 3D come persp() per visualizzare la superficie.
La griglia fitta 40 per 40 generata dal comando mvals <- expand.grid(seq(-5, 10, length.out = 40), seq(-8, 4, length.out = 40)) è già caricata per te.
Questo esercizio fa parte del corso
Distribuzioni di probabilità multivariate in R
Istruzioni dell'esercizio
- Calcola le densità di una normale con media
mu.sime matrice varianza-covarianzasigma.simsui valori della griglia nell'oggettomvals. - Usa la funzione di grafico
persp()per visualizzare la superficie di densità sulla griglia.
Esercizio pratico interattivo
Prova a risolvere questo esercizio completando il codice di esempio.
# 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")