Mô phỏng phân phối tiên nghiệm Beta

Giả sử bạn đang tranh cử một vị trí công quyền. Gọi \(p\) là mức ủng hộ thực sự của bạn, tức tỷ lệ cử tri dự định bỏ phiếu cho bạn. Dựa trên các khảo sát trước đây, mô hình tiên nghiệm của \(p\) được mô tả bằng phân phối Beta với các tham số shape là 45 và 55.

Bạn sẽ xấp xỉ tiên nghiệm Beta(45, 55) bằng cách lấy mẫu ngẫu nhiên từ hàm rbeta(). Hàm này nhận ba đối số: kích thước mẫu (n) và hai tham số shape (shape1,shape2). Tiếp theo, bạn sẽ dựng biểu đồ mật độ của các mẫu bằng ggplot(). Hàm này nhận hai đối số: bộ dữ liệu chứa các mẫu và, trong aes(), biến được vẽ trên trục x. Lớp biểu đồ mật độ được thêm bằng geom_density().

Bài tập này là một phần của khóa học

Mô hình hóa Bayesian với RJAGS

Xem khóa học

Hướng dẫn bài tập

Dùng rbeta() để lấy 10.000 mẫu từ Beta(45, 55). Gán kết quả cho prior_A.
Khung dữ liệu prior_sim đã bao gồm mẫu prior_A. Áp dụng ggplot() lên prior_sim để dựng biểu đồ mật độ cho các mẫu tiên nghiệm.

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.

# Sample 10000 draws from Beta(45,55) prior
prior_A <- rbeta(n = ___, shape1 = ___, shape2 = ___)

# Store the results in a data frame
prior_sim <- data.frame(prior_A)

# Construct a density plot of the prior sample
ggplot(prior_sim, aes(x = ___)) + 
    geom_density()

Chỉnh sửa và Chạy Mã

Bài tập này là một phần của khóa học

Mô hình hóa Bayesian với RJAGS

SkillTag.level.advancedSkillTag.label

4.9+

Bắt đầu khóa học miễn phí

Mô hình Bayesian kết hợp hiểu biết tiên nghiệm với hiểu biết từ dữ liệu quan sát để tạo thành hiểu biết hậu nghiệm, đã được cập nhật, về một tham số. Trong chương này, bạn sẽ ôn lại các khái niệm Bayesian trong bối cảnh mô hình Beta–Binomial nền tảng cho một tham số tỷ lệ. Bạn cũng sẽ học cách dùng gói rjags để định nghĩa, biên dịch và mô phỏng mô hình này trong R.

Exercise 1: Mô hình tiên nghiệm Exercise 2: Mô phỏng phân phối tiên nghiệm Beta

Bài tập hiện tại

Exercise 3: So sánh các prior Beta Exercise 4: Chọn prior nào?Exercise 5: Dữ liệu & hàm khả năng Exercise 6: Mô phỏng sự phụ thuộc của X vào p Exercise 7: Xấp xỉ hàm likelihood Exercise 8: Diễn giải hàm hợp lý Exercise 9: Mô hình hậu nghiệm Exercise 10: Định nghĩa, biên dịch và mô phỏng Exercise 11: Cập nhật hậu nghiệm Exercise 12: Ảnh hưởng của prior và dữ liệu lên posterior

Mô hình Bayesian Normal–Normal với hai tham số cung cấp nền tảng đơn giản cho các mô hình hồi quy Normal. Trong chương này, bạn sẽ xây dựng mô hình Normal–Normal và định nghĩa, biên dịch, cũng như mô phỏng mô hình này bằng rjags. Bạn cũng sẽ khám phá "ma thuật" của cơ chế chuỗi Markov đứng sau quá trình mô phỏng của rjags.

Exercise 1: Mô hình Normal-Normal Exercise 2: Prior Normal-Normal Exercise 3: Dữ liệu nghiên cứu giấc ngủ Exercise 4: Nhận định từ prior và dữ liệu Exercise 5: Mô phỏng Normal–Normal trong RJAGS Exercise 6: Định nghĩa, biên dịch và mô phỏng Normal–Normal Exercise 7: Nhận định hậu nghiệm về thiếu ngủ Exercise 8: Chuỗi Markov Exercise 9: Lưu trữ các chuỗi Markov Exercise 10: Biểu đồ vết (trace) của chuỗi Markov Exercise 11: Biểu đồ mật độ của chuỗi Markov Exercise 12: Chẩn đoán chuỗi Markov & khả năng tái lập Exercise 13: Nhiều chuỗi Exercise 14: Sai số chuẩn ngây thơ Exercise 15: Tái lập kết quả

Trong chương này, bạn sẽ mở rộng mô hình Normal–Normal thành một mô hình hồi quy Bayesian đơn giản. Trong bối cảnh đó, bạn sẽ khám phá cách dùng đầu ra mô phỏng từ rjags để thực hiện suy luận hậu nghiệm. Cụ thể, bạn sẽ xây dựng ước lượng hậu nghiệm cho các tham số hồi quy bằng trung bình hậu nghiệm và khoảng tin cậy (credible intervals), kiểm định giả thuyết bằng xác suất hậu nghiệm, và xây dựng phân phối dự đoán hậu nghiệm cho các quan sát mới.

Exercise 1: Một mô hình hồi quy Bayesian đơn giản Exercise 2: Tiên nghiệm cho hồi quy Exercise 3: Trực quan hóa prior của hồi quy Exercise 4: Dữ liệu cân nặng & chiều cao Exercise 5: Hồi quy Bayesian trong RJAGS Exercise 6: Định nghĩa, biên dịch và mô phỏng mô hình hồi quy Exercise 7: Chuỗi Markov trong hồi quy Exercise 8: Ước lượng hậu nghiệm & suy luận Exercise 9: Ước lượng điểm hậu nghiệm Exercise 10: Khoảng tin cậy hậu nghiệm Exercise 11: Xác suất hậu nghiệm Exercise 12: Dự đoán hậu nghiệm Exercise 13: Suy luận cho xu hướng hậu nghiệm Exercise 14: Tính toán dự đoán hậu nghiệm Exercise 15: Phân phối dự đoán hậu nghiệm

Ở chương cuối, bạn sẽ khái quát hóa mô hình hồi quy Normal đơn giản để áp dụng trong các bối cảnh rộng hơn. Bạn sẽ đưa biến dự báo phân loại vào mô hình, xây dựng một mô hình hồi quy đa biến với hai biến dự báo, và cuối cùng mở rộng phương pháp này sang các mô hình hồi quy Poisson đa biến cho biến đếm.

Exercise 1: Hồi quy Bayes với biến dự báo phân loại Exercise 2: Dữ liệu mẫu RailTrail Exercise 3: Mô phỏng RJAGS với biến phân loại Exercise 4: Diễn giải hệ số biến phân loại Exercise 5: Suy luận về lưu lượng theo ngày trong tuần Exercise 6: Hồi quy Bayes đa biến Exercise 7: Xem lại dữ liệu RailTrail Exercise 8: Mô phỏng RJAGS cho hồi quy đa biến Exercise 9: Diễn giải các tham số hồi quy đa biến Exercise 10: Suy luận hậu nghiệm cho hồi quy đa biến Exercise 11: Hồi quy Poisson theo phương pháp Bayes Exercise 12: Mô phỏng RJAGS cho hồi quy Poisson Exercise 13: Vẽ mô hình hồi quy Poisson Exercise 14: Suy luận cho tham số tốc độ Poisson Exercise 15: Dự đoán hậu nghiệm Poisson Exercise 16: Kết luận