r - 绘制 N 给定似然性和先验的后验分布

Question

假设似然是 Y|n~Binomial(n, theta) 并且先验是 n~Poisson(5)，我试图计算 Y=0 和 theta=0.2 的样本大小 n 的后验分布。我的代码如下：

Y <- 0
theta <- 0.2
n_grid <- seq(0,1,length=1000)
like <- dbinom(Y,n_grid,theta)
prior <- dpois(n_grid,5)
fy <- sum(like*prior)
post <- like*prior/fy
plot(n_grid,post,type="l")

在计算似然函数和先验时，我不断得到 NaN 结果。任何帮助，将不胜感激！

标签： rstatisticsbayesian

Answer 1

所以我意识到回答我自己的问题可能是非常规的，但我想出了我的解决方案，并认为我会发布答案来帮助其他人。

Y          <- 0
theta      <- 0.2
N          <- 0:0.01:100
like       <- dbinom(Y,N,theta)
prior      <- dpois(N,5)
fy         <- sum(like*prior)
post       <- like*prior/fy

plot(N,post,type="l")