Tính mật độ của phân phối chuẩn đa biến
Với nhiều tác vụ thống kê như kiểm định giả thuyết, phân cụm và tính hợp lý (likelihood), bạn cần tính mật độ của một phân phối chuẩn đa biến đã cho. Trong bài này, bạn sẽ dùng hàm dmvnorm() để tính các giá trị mật độ chuẩn đa biến với vector kỳ vọng và ma trận phương sai–hiệp phương sai chỉ định, tại từng quan sát trong mẫu multnorm.sample mà bạn đã tạo trước đó.
Vector kỳ vọng và ma trận phương sai–hiệp phương sai đã được nạp sẵn cho bạn trong các đối tượng mu.sim và sigma.sim.
Bài tập này là một phần của khóa học
Phân phối Xác suất Đa biến trong R
Hướng dẫn bài tập
- Dùng
dmvnorm()để tính độ cao mật độ cho 100 mẫu trongmultnorm.sampleđối với phân phối chuẩn hai biến. - Dùng
scatterplot3d()để vẽ biểu đồ phân tán 3D của các độ cao mật độ tại từng điểm mẫu đã tạo.
Bài tập tương tác thực hành trực tiếp
Hãy thử làm bài tập này bằng cách hoàn thành đoạn mã mẫu này.
# 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")