我如何在这个 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)")