Výpočet dmvnorm na mřížce

Abychom mohli vizualizovat povrch hustoty bivariátního normálního rozdělení, je potřeba vypočítat hustotu na husté mřížce souřadnic x a y a k vizualizaci povrchu použít funkce pro 3D vykreslování, jako je persp().

Hustá mřížka 40 × 40 vygenerovaná příkazem mvals <- expand.grid(seq(-5, 10, length.out = 40), seq(-8, 4, length.out = 40)) je už předem načtená.

Toto cvičení je součástí kurzu

Multivariate Probability Distributions in R

Zobrazit kurz

Pokyny k cvičení

Vypočítej hustoty normálního rozdělení se střední hodnotou mu.sim a varianční-kovarianční maticí sigma.sim pro hodnoty mřížky uložené v objektu mvals.
Pomocí funkce persp() vizualizuj povrch hustoty nad touto mřížkou.

Interaktivní cvičení na vyzkoušení si v praxi

Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.

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

Upravit a spustit kód