Calcular la densidad de la normal multivariante
En muchas tareas estadísticas, como contrastes de hipótesis, clustering y cálculo de verosimilitud, necesitas calcular la densidad de una distribución normal multivariante especificada. En este ejercicio, usarás la función dmvnorm() para calcular densidades normales multivariantes con una media y una matriz de varianzas-covarianzas dadas en cada una de las observaciones de tu muestra generada previamente multnorm.sample.
La media y la matriz de varianzas-covarianzas están precargadas como los objetos mu.sim y sigma.sim.
Este ejercicio forma parte del curso
Distribuciones de probabilidad multivariantes en R
Instrucciones del ejercicio
- Usa
dmvnorm()para calcular las alturas de la densidad de las 100 muestras enmultnorm.samplepara una normal bivariante. - Usa
scatterplot3d()para trazar un diagrama de dispersión 3D de las alturas de la densidad en cada uno de los puntos de la muestra generada.
Ejercicio interactivo práctico
Prueba este ejercicio y completa el código de muestra.
# Calculate density
multnorm.dens <- dmvnorm(multnorm.sample, mean = ___, sigma = ___)
# Create scatter plot of density heights
___(cbind(___),
color="blue", pch="", type = "h",
xlab = "x", ylab = "y", zlab = "density")