นี่จะเป็นคำถามที่ @ Xi'an มากกว่าคำตอบ
ฉันจะตอบว่าความแปรปรวนหลัง
กับจำนวนของการทดลองที่จำนวนที่ประสบความสำเร็จและค่าสัมประสิทธิ์ของเบต้าก่อนที่เกินความแปรปรวนก่อน
เป็นไปได้ในรูปแบบทวินามตามตัวอย่างด้านล่างซึ่งมีโอกาสและ ก่อนหน้านี้อยู่ในทางตรงกันข้ามโดยสิ้นเชิงเพื่อให้หลังเป็น "ไกลเกินไปในระหว่าง" ดูเหมือนจะขัดแย้งกับคำพูดของ Gelman
V(θ|y)=α1β1(α1+β1)2(α1+β1+1)=(α0+k)(n−k+β0)(α0+n+β0)2(α0+n+β0+1),
nkα0,β0V(θ)=α0β0(α0+β0)2(α0+β0+1)
n <- 10
k <- 1
alpha0 <- 100
beta0 <- 20
theta <- seq(0.01,0.99,by=0.005)
likelihood <- theta^k*(1-theta)^(n-k)
prior <- function(theta,alpha0,beta0) return(dbeta(theta,alpha0,beta0))
posterior <- dbeta(theta,alpha0+k,beta0+n-k)
plot(theta,likelihood,type="l",ylab="density",col="lightblue",lwd=2)
likelihood_scaled <- dbeta(theta,k+1,n-k+1)
plot(theta,likelihood_scaled,type="l",ylim=c(0,max(c(likelihood_scaled,posterior,prior(theta,alpha0,beta0)))),ylab="density",col="lightblue",lwd=2)
lines(theta,prior(theta,alpha0,beta0),lty=2,col="gold",lwd=2)
lines(theta,posterior,lty=3,col="darkgreen",lwd=2)
legend("top",c("Likelihood","Prior","Posterior"),lty=c(1,2,3),lwd=2,col=c("lightblue","gold","darkgreen"))
> (postvariance <- (alpha0+k)*(n-k+beta0)/((alpha0+n+beta0)^2*(alpha0+n+beta0+1)))
[1] 0.001323005
> (priorvariance <- (alpha0*beta0)/((alpha0+beta0)^2*(alpha0+beta0+1)))
[1] 0.001147842
ดังนั้นตัวอย่างนี้แสดงให้เห็นถึงความแปรปรวนด้านหลังที่มีขนาดใหญ่กว่าในแบบจำลองทวินาม
แน่นอนว่านี่ไม่ใช่ความแปรปรวนด้านหลังที่คาดหวัง นั่นคือสิ่งที่ความแตกต่างอยู่?
ตัวเลขที่สอดคล้องกันคือ