Biểu đồ vết (trace) của chuỗi Markov

Một trace plot giúp bạn trực quan hóa diễn biến theo thời gian của chuỗi Markov. Cụ thể, trace plot cho chuỗi \(m\) sẽ vẽ giá trị quan sát của chuỗi (trục y) theo số lần lặp tương ứng (trục x).

Bạn sẽ dựng trace plot của chuỗi \(m\) theo hai cách: dùng hàm có sẵn plot() áp dụng cho đối tượng mcmc.list là sleep_sim và, để kiểm soát đồ họa tốt hơn (và kiểm soát phân tích ở các chương sau), dùng ggplot() áp dụng cho đối tượng data.frame là sleep_chains. Cả sleep_sim và sleep_chains đều đã có trong workspace của bạn:

sleep_sim <- coda.samples(model = sleep_jags, variable.names = c("m", "s"), n.iter = 10000)
sleep_chains <- data.frame(sleep_sim[[1]], iter = 1:10000)

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

Áp dụng plot() cho sleep_sim với density = FALSE để dựng trace plot cho các chuỗi \(m\) và \(s\). LƯU Ý: 10.000 Iterations được ghi nhận bắt đầu sau giai đoạn "burn-in" khi các mẫu bị loại bỏ. Vì vậy, số đếm Iterations không bắt đầu từ 1!
Áp dụng ggplot() với lớp geom_line() cho sleep_chains để dựng lại trace plot của chuỗi $m`.
Phóng to: dựng một trace plot với ggplot() cho 100 lần lặp đầu tiên của chuỗi $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.

# Use plot() to construct trace plots of the m and s chains


# Use ggplot() to construct a trace plot of the m chain
ggplot(___, aes(x = ___, y = ___)) + 
    geom_line()

# Trace plot the first 100 iterations of the m chain
ggplot(___, aes(x = ___, y = ___)) + 
    geom_line()

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 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

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

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