그리드에서 dmvnorm 계산하기
이변량 정규 밀도 곡면을 시각화하려면, 촘촘한 x, y 좌표 그리드 전체에 걸쳐 밀도를 계산한 뒤 persp()와 같은 3D 곡면 플롯 함수로 시각화해야 합니다.
mvals <- expand.grid(seq(-5, 10, length.out = 40), seq(-8, 4, length.out = 40)) 명령으로 생성된 40×40 크기의 촘촘한 그리드가 미리 불러와져 있습니다.
이 연습은 강의의 일부입니다
R로 배우는 다변량 확률 분포
연습 안내
- 평균
mu.sim과 분산-공분산 행렬sigma.sim을 갖는 정규분포의 밀도를mvals객체의 그리드 값에 대해 계산하세요. persp()플롯 함수를 사용하여 그리드 위의 밀도 곡면을 시각화하세요.
실습형 인터랙티브 연습
이 예제를 이 샘플 코드를 완성하여 풀어보세요.
# 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")