Obliczanie dmvnorm na siatce punktów
Aby zwizualizować powierzchnię gęstości dwuwymiarowego rozkładu normalnego, trzeba obliczyć gęstość na gęstej siatce współrzędnych x i y, a następnie użyć funkcji do rysowania powierzchni 3D, takich jak persp().
Gęsta siatka 40 na 40 punktów, wygenerowana poleceniem mvals <- expand.grid(seq(-5, 10, length.out = 40), seq(-8, 4, length.out = 40)), jest już wczytana.
To ćwiczenie jest częścią kursu
Wielowymiarowe rozkłady prawdopodobieństwa w R
Instrukcje do ćwiczenia
- Oblicz gęstości rozkładu normalnego ze średnią
mu.simi macierzą wariancji-kowariancjisigma.simdla wartości siatki zapisanych w obiekciemvals. - Użyj funkcji
persp(), aby zwizualizować powierzchnię gęstości na siatce.
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
# 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")