绘制给定似然和先验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结果。任何帮助将不胜感激！

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")