Jeffreys建议为参数N提供合适的先验,其中N = 0、1、2,...为p(N)∝1 / N。
下面是一个示例:
假设一个公平的硬币出现了y = 10次。有多少次 硬币被扔了吗?用N表示这个未知数量。
我们可以解析地计算后验分布,它是负二项分布,平均值为2y = 20。
但是,当我使用WinBUGS对其进行建模时,我不会如何表示先前的分布“ 1 / N”。
代码如下:
model{
y <- 10
y ~ dbin(0.5, N)
N ~ ??1/N?? #How to represent that p(N) is proportional to 1/N
}