2

我正在尝试编写保留模型,但遇到了一些问题:

  1. 当我按“加载初始化”时,我收到消息

    该链包含单元化变量

    尽管我已经初始化了所有变量。

  2. 当我尝试更新我的模型时,忽略第一个问题并按“gen inits”,我收到消息

    gamma full 条件的形状参数太小

    尽管我使用的伽马分布参数与 Ntzoufras 和 Dellaportas (2002) 在“包含索赔计数不确定性的突出责任的贝叶斯模型”中建议的参数相同。

这是 WinBUGS 代码:

#MODEL
model {

for (i in 1:n) {

for (j in 1:n){

x[i,j] ~ dlnorm(m[i,j], t)

m[i,j] <- c + a[i] + b[j]

} 

a[1] <- -sum(a[2:n])

b[1] <- -sum(b[2:n])

c ~ dnorm(0,0.001)

for (i in 2:n) {

a[i] ~ dnorm(0, 0.01)

b[i] ~ dnorm(0,0.01) }

t ~ dgamma(0.001, 0.001)

sigma2 <- 1/t

R[1] <- 0

for (i in 2:n){

R[i] <- sum(x[i,(n+2-i):n])

}

Total <- sum(R[])

}
#DATA
list(n=10)

x[,1] x[,2] [,3] x[,4] x[,5] x[,6] x[,7] x[,8] x[,9] x[,10]

357848  766940  610542  482940 527326 574398 146342 139950 227229 67948

352118  884021  933894 1183289 445745 320996 527804 266172 425046 NA

290507 1001799  926219 1016654 750816 146923 495992 280405 NA NA

310608 1108250  776189 1562400 272482 352053 206286 NA NA NA

443160  693190  991983  769488 504851 470639 NA NA NA NA

396132  937085  847498  805037 705960 NA NA NA NA NA

440832  847631 1131398 1063269 NA NA NA NA NA NA

359480 1061648 1443370 NA NA NA NA NA NA NA

376686  986608 NA NA NA NA NA NA NA NA

344014 NA NA NA NA NA NA NA NA NA

END
#INITS



list(a=c(NA,2,2,2,2,2,2,2,2,2),
b=c(NA,0,0,0,0,0,0,0,0,0), c=10, t=10)
4

0 回答 0