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 completando este código de ejemplo.
# 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")