dmvnorm über einem Gitter berechnen
Um eine bivariate Normaldichte als Oberfläche darzustellen, musst du die Dichte über einem dichten Gitter von x- und y-Koordinaten berechnen und 3D-Oberflächenfunktionen wie persp() zur Visualisierung verwenden.
Das 40-mal-40-Punkte-Gitter aus dem Befehl mvals <- expand.grid(seq(-5, 10, length.out = 40), seq(-8, 4, length.out = 40)) ist bereits für dich geladen.
Diese Übung ist Teil des Kurses
<Kurs>Multivariate Wahrscheinlichkeitsverteilungen in R</Kurs>Übungsanweisungen
- Berechne die Dichten einer Normalverteilung mit Mittelwert
mu.simund Varianz-Kovarianz-Matrixsigma.simauf den Gitterwerten im Objektmvals. - Verwende die Plot-Funktion
persp(), um die Dichteoberfläche über dem Gitter zu visualisieren.
Interaktive praktische Übung
Versuche dich an dieser Übung, indem du diesen Beispielcode vervollständigst.
# 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")