İki hafta önce burada verdiğim bir soruya "cevap" yazıyorum: Jeffreys neden daha önce yararlı? Bu gerçekten bir soruydu (ve ben de o zaman yorum gönderme hakkım yoktu), bu yüzden umarım bunu yapmak için sorun yok:
Yukarıdaki bağlantıda, Jeffreys'in daha önce ilginç olan özelliğinin, modeli yeniden parametrelendirirken ortaya çıkan posterior dağılımın, dönüşümün getirdiği kısıtlamalara uyan posterior olasılıklar verdiği tartışılmaktadır. Orada anlatıldığı gibi başarı olasılığı taşırken ki, oran Beta-Bernoulli örnekte , bu durumda olması gereken bir arka tatmin olduğu .
Öncesinde transforme edilmesi için Jeffreys değişmezliği sayısal bir örnek oluşturmak için istediği oranları için ve daha ilginç bir şekilde, bunun, diğer önsel (örneğin, Haldane, tekdüze veya isteğe bağlı olanlar) yoksundur.
Şimdi, başarı olasılığı için posterior Beta ise (sadece Jeffreys için değil, önceki herhangi bir Beta için), olasılıkların posterioru aynı parametrelerle ikinci tür bir Beta dağılımını takip eder (Wikipedia'ya bakın) . Aşağıdaki sayısal örnekte vurgulanmıştır Sonra, çok değişmezliği olduğunu (en azından benim için) şaşırtıcı değildir öncesinde Beta herhangi bir seçim için (oynayabilir alpha0_U
ve beta0_U
), sadece Jeffreys, krş program çıktısı.
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)
Bu beni şu sorulara getiriyor:
- Bir hata yapabilir miyim?
- Hayır ise, eşlenik ailelerde değişmezlik olmaması gibi bir sonuç var mı, ya da bunun gibi bir şey mi? (Hızlı inceleme beni örneğin normal-normal durumda değişmezlik üretemediğimden şüpheleniyor.)
- Eğer biz hangi bir (tercihen basit) örnek biliyor musunuz yapmak değişmezliği eksikliği olsun?