Tôi đang đăng lại một "câu trả lời" cho một câu hỏi mà tôi đã đưa ra hai tuần trước tại đây: Tại sao Jeffreys trước lại hữu ích? Tuy nhiên, đây thực sự là một câu hỏi (và tôi cũng không có quyền đăng bình luận), vì vậy tôi hy vọng sẽ ổn khi làm điều này:
Trong liên kết ở trên, người ta đã thảo luận rằng tính năng thú vị của Jeffreys trước đó là, khi xác định lại mô hình, phân phối hậu quả mang lại xác suất sau tuân theo các hạn chế do biến đổi áp đặt. Say, như đã thảo luận ở đó, khi di chuyển từ xác suất thành công trong ví dụ Beta-Bernoulli để tỷ lệ cược , nó phải là trường hợp mà một đáp ứng sau .
Tôi muốn tạo ra một ví dụ số của bất biến của Jeffreys trước cho chuyển để tỷ lệ cược , và, thú vị hơn, thiếu của chúng thuộc priors khác (nói, Haldane, thống nhất, hoặc những người tùy ý).
Bây giờ, nếu hậu tố cho xác suất thành công là Beta (đối với bất kỳ Beta nào trước đó, không chỉ Jeffreys), thì hậu thế của tỷ lệ cược tuân theo phân phối Beta của loại thứ hai (xem Wikipedia) với cùng tham số . Sau đó, như được nhấn mạnh trong ví dụ bằng số bên dưới, không quá ngạc nhiên (ít nhất là với tôi) rằng có bất kỳ sự lựa chọn nào cho bất kỳ lựa chọn Beta nào trước đó (chơi xung quanh alpha0_U
và beta0_U
), không chỉ Jeffreys, nhận đầu ra của chương trình.
library(GB2)
# has the Beta density of the 2nd kind, the distribution of theta/(1-theta) if theta~Beta(alpha,beta)
theta_1 = 2/3 # a numerical example as in the above post
theta_2 = 1/3
odds_1 = theta_1/(1-theta_1) # the corresponding odds
odds_2 = theta_2/(1-theta_2)
n = 10 # some data
k = 4
alpha0_J = 1/2 # Jeffreys prior for the Beta-Bernoulli case
beta0_J = 1/2
alpha1_J = alpha0_J + k # the corresponding parameters of the posterior
beta1_J = beta0_J + n - k
alpha0_U = 0 # some other prior
beta0_U = 0
alpha1_U = alpha0_U + k # resulting posterior parameters for the other prior
beta1_U = beta0_U + n - k
# posterior probability that theta is between theta_1 and theta_2:
pbeta(theta_1,alpha1_J,beta1_J) - pbeta(theta_2,alpha1_J,beta1_J)
# the same for the corresponding odds, based on the beta distribution of the second kind
pgb2(odds_1, 1, 1,alpha1_J,beta1_J) - pgb2(odds_2, 1, 1,alpha1_J,beta1_J)
# same for the other prior and resulting posterior
pbeta(theta_1,alpha1_U,beta1_U) - pbeta(theta_2,alpha1_U,beta1_U)
pgb2(odds_1, 1, 1,alpha1_U,beta1_U) - pgb2(odds_2, 1, 1,alpha1_U,beta1_U)
Điều này mang lại cho tôi các câu hỏi sau đây:
- Tôi có phạm sai lầm không?
- Nếu không, có một kết quả như không thiếu sự bất biến trong các gia đình liên hợp, hoặc một cái gì đó như thế? (Kiểm tra nhanh khiến tôi nghi ngờ rằng tôi có thể chẳng hạn, cũng không tạo ra sự bất biến trong trường hợp bình thường.)
- Bạn có biết một (tốt nhất là đơn giản) ví dụ mà chúng ta làm được thiếu bất biến?