Tự triển khai thuật toán Metropolis-Hastings (MH) bằng Mata Implementing the Metropolis-Hastings (MH) algorithm from scratch in Mata Từ việc sử dụng đến việc xây dựng công cụ Ở bài học trước, chúng ta đã thấy lệnh bayesmh của Stata có thể thực hiện phân tích Bayes một cách kỳ diệu như thế nào. Nhưng điều gì thực sự xảy ra đằng sau lệnh đó? Làm thế nào Stata có thể rút ra hàng ngàn mẫu từ một phân phối hậu nghiệm phức tạp mà chúng ta không thể giải quyết bằng toán học thông thường? Câu trả lời nằm ở các thuật toán (Markov chain Monte Carlo – MCMC), và một trong những thuật toán nền tảng và phổ biến nhất chính là Metropolis-Hastings (MH). Trong bài học này, chúng ta sẽ mở “hộp đen” đó ra. Thay vì chỉ sử dụng một lệnh có sẵn, chúng ta sẽ tự tay xây dựng thuật toán MH từ đầu bằng cách sử dụng Mata, ngôn ngữ lập trình ma trận mạnh mẽ được tích hợp trong Stata. Việc này …