y=26
n=100
a=1
b=1
newa=a+y
newb=b+n-y
theta=seq(0.01,0.99,l=500)
plot(theta,dbeta(theta,newa,newb),type="l")
qbeta(0.5,newa,newb)
qbeta(0.9,newa,newb)
qbeta(0.1,newa,newb)
newa/(newa+newb)
(newa-1)/(newa+newb-2)
0.2647059*(1-0.2647059)/(newa+newb+1)
sqrt(0.001889677)
qbeta(0.25,newa,newb)
qbeta(0.75,newa,newb)
newa
newb
0.293313-0.2344279
qbeta(0.025,newa,newb)
qbeta(0.975,newa,newb)