এটি উত্তরটির চেয়ে @ সিয়ানকে আরও প্রশ্ন করতে চলেছে।
আমি উত্তর দিতে যাচ্ছিলাম যে একটি উত্তর বৈকল্পিক
V(θ|y)=α1β1(α1+β1)2(α1+β1+1)=(α0+k)(n−k+β0)(α0+n+β0)2(α0+n+β0+1),
সঙ্গে
n পরীক্ষার সংখ্যা,
k সাফল্যের সংখ্যা এবং
α0,β0 পূর্বে বিটার সহগগুলি, পূর্বের বৈকল্পিককে অতিক্রম করে
V(θ)=α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
অতএব, এই উদাহরণটি দ্বিপদী মডেলের আরও বৃহত্তর পশ্চাত্পদটি প্রস্তাব করে।
অবশ্যই এটি প্রত্যাশিত উত্তর বৈকল্পিকতা নয়। এখানেই কি তাত্পর্য রয়েছে?
সংশ্লিষ্ট চিত্রটি হ'ল