pi<-seq(0,1,0.0001)
y_n<-60
n<-100

posterior<-((pi^y_n)*(1-pi)^(n-y_n))/beta(y_n+1,n-y_n+1)

の95％信頼区間を求めたいんですがその関数がよくわからないです誰かわかる人いらっしゃいますか？