Vẽ biểu đồ các cụm ước lượng
Đến thời điểm này, bạn đã phát hiện ra hai cụm trong data frame gaussian_sample. Trong bài tập này, bạn sẽ trực quan hóa cách các cụm được ước lượng ở vòng lặp 10 khớp với dữ liệu. Các vector means_iter10 và props_iter10 đã được lưu sẵn trong môi trường.
Để làm điều đó, bạn sẽ dùng hàm stat_function() của ggplot2, cho phép bạn chồng một hàm lên trên biểu đồ hiện có. Hàm bạn sẽ dùng là một hàm vẽ đường cong có sẵn tên fun_gaussian(), nhận các đối số là kỳ vọng (mean) và tỷ trọng (proportion) của phân phối Gaussian.
Bài tập này là một phần của khóa học
Các Mô hình Hỗn hợp trong R
Hướng dẫn bài tập
Vẽ biểu đồ histogram ở chế độ mật độ của biến x và thêm các đường cong ước lượng bằng cách dùng stat_function() kết hợp với hàm fun_gaussian().
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.
___ %>%
ggplot() + geom_histogram(aes(x = ___, y = ___), bins = 200) +
stat_function(geom = "line", fun = fun_gaussian,
args = list(mean = means_iter10[1], proportion = ___[1])) +
stat_function(geom = "line", fun = fun_gaussian,
args = list(mean = ___[2], proportion = props_iter10[2]))