Calculando dmvnorm em uma grade
Para visualizar uma superfície de densidade normal bivariada, é preciso calcular a densidade em uma grade densa de coordenadas x e y e usar funções de plotagem 3D como persp() para visualizar a superfície.
A grade densa de 40 por 40 gerada com o comando mvals <- expand.grid(seq(-5, 10, length.out = 40), seq(-8, 4, length.out = 40)) já está pré-carregada para você.
Este exercicio faz parte do curso
Distribuições de Probabilidade Multivariadas em R
Instruções do exercicio
- Calcule as densidades de uma normal com média
mu.sime matriz de variância-covariânciasigma.simnos valores da grade no objetomvals. - Use a função de plotagem
persp()para visualizar a superfície de densidade sobre a grade.
exercicio interativo prático
Tente este exercicio completando este código de exemplo.
# 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")