Vấn đề
Tôi đang viết một hàm R thực hiện phân tích Bayes để ước tính mật độ sau được cung cấp trước và dữ liệu được thông báo. Tôi muốn chức năng gửi cảnh báo nếu người dùng cần xem xét lại trước.
Trong câu hỏi này, tôi quan tâm đến việc học cách đánh giá trước. Các câu hỏi trước đây đã đề cập đến các cơ chế nêu rõ các linh mục được thông báo ( ở đây và ở đây .)
Các trường hợp sau đây có thể yêu cầu đánh giá lại trước:
- dữ liệu đại diện cho một trường hợp cực đoan không được tính khi nêu trước
- lỗi trong dữ liệu (ví dụ: nếu dữ liệu tính theo đơn vị g khi trước đó tính bằng kg)
- sai trước được chọn từ một nhóm các linh mục có sẵn vì một lỗi trong mã
Trong trường hợp đầu tiên, các linh mục thường không đủ sức khuếch tán để dữ liệu nói chung sẽ áp đảo chúng trừ khi các giá trị dữ liệu nằm trong một phạm vi không được hỗ trợ (ví dụ <0 cho logN hoặc Gamma). Các trường hợp khác là lỗi hoặc lỗi.
Câu hỏi
- Có bất kỳ vấn đề nào liên quan đến tính hợp lệ của việc sử dụng dữ liệu để đánh giá trước không?
- bất kỳ thử nghiệm cụ thể phù hợp nhất cho vấn đề này?
Ví dụ
Dưới đây là hai bộ dữ liệu được kết hợp kém với trước vì chúng đến từ các quần thể có (đỏ) hoặc (xanh dương).
Dữ liệu màu xanh có thể là kết hợp dữ liệu + trước hợp lệ trong khi dữ liệu màu đỏ sẽ yêu cầu phân phối trước được hỗ trợ cho các giá trị âm.
set.seed(1)
x<- seq(0.01,15,by=0.1)
plot(x, dlnorm(x), type = 'l', xlim = c(-15,15),xlab='',ylab='')
points(rnorm(50,0,5),jitter(rep(0,50),factor =0.2), cex = 0.3, col = 'red')
points(rnorm(50,8,0.5),jitter(rep(0,50),factor =0.4), cex = 0.3, col = 'blue')