我如何在這個 r 函數中找到后驗分布？

Question

我試圖弄清楚如何從無信息先驗分布中找到后驗分布。

我的代碼是 -

like<-function(x) {dbinom(4. ,10, x)}
integrate(like,0,1)
pi<-seq(0,1, by=0.001)
like<-dbinom(4,10,pi)
post<-function(x) {dbinom(4,10,x)/0.2}
integrate(post,.6,.8)
plot(pi,post,type = "l", col = "red", ylab = "Posterior f(pi)")

我不斷收到錯誤消息 -

xy.coords(x, y, xlabel, ylabel, log) 中的錯誤：“x”和“y”長度不同

我將如何解決這個問題？

Answer 1

也許你正在尋找這個：

#Index
like<-function(x) {dbinom(4. ,10, x)}
integrate(like,0,1)
pi<-seq(0,1, by=0.001)
like<-dbinom(4,10,pi)
post<-function(x) {dbinom(4,10,x)/0.2}
integrate(post,.6,.8)
plot(post(pi),type = "l", col = "red", ylab = "Posterior f(pi)")