如何在WinBUGS中表示“杰弗里斯的先验”“ 1 / N”

时间:2019-06-29 14:20:57

标签: bayesian winbugs

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    
}

0 个答案:

没有答案