Suy luận về lưu lượng theo ngày trong tuần

Đầu ra mô phỏng RJAGS 10.000 vòng lặp, rail_sim_1, có trong workspace của bạn cùng với một data frame chứa kết quả chuỗi Markov:

> head(rail_chains_1, 2)
         a b.1.       b.2.        s 
1 420.6966    0  -54.30783 118.2328
2 399.5823    0  -52.02570 119.9499

Các chuỗi này cung cấp 10.000 bộ giá trị cho a, lưu lượng điển hình trên đường mòn vào ngày cuối tuần, và b.2., phần chênh lệch giữa lưu lượng điển hình ngày thường so với cuối tuần. Ví dụ, bộ tham số đầu tiên cho thấy thường có 420.6966 người đi vào ngày cuối tuần và ít hơn 54.30783 người vào ngày thường. Vậy nên thường có 420.6966 - 54.30783 = 366.3888 người đi vào ngày thường. Bạn sẽ dùng dữ liệu mô phỏng này để suy luận về lưu lượng trên đường mòn vào ngày thường.

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

Kết hợp các giá trị chuỗi a và b.2. để tạo một chuỗi 10.000 giá trị cho lưu lượng điển hình vào ngày thường. Lưu thành weekday_mean trong rail_chains_1.
Dùng ggplot() để vẽ biểu đồ mật độ cho các giá trị chuỗi weekday_mean.
Tạo khoảng tin cậy 95% cho lưu lượng điển hình vào ngày thường.

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.

# Construct a chain of values for the typical weekday volume
rail_chains_1 <- rail_chains_1 %>% 
    mutate(weekday_mean = ___)

# Construct a density plot of the weekday chain
ggplot(___, aes(x = ___)) + 
    geom_density()

# 95% credible interval for typical weekday volume

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

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

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