Как представить «Джеффриса до» «1 / N» в WinBUGS

Question

Джеффрис предложил подходящий априор для параметра 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    
}