Hello,
I am trying to Apply multinomial logit model in joint diseases mapping in openbugs.
What does the error message "Sorry, something went wrong in procedure Stack.Value in module GraphStack" mean? It's given when I run the following model:
(the model runs with small values "osse" node, for example osse=structure(.Data=c(69, 56, 199,131, 179, 360,257, 327, 535,28, 73, 26,23, 44, 64,279, 363, 391,186, 365, 276,65, 41, 100,177, 188, 320,77, 65, 153,91, 90, 182,29, 40, 47,17, 14, 31,92, 212, 213,125, 153, 140,26, 23, 36,278, 293, 751,19, 21, 61,229, 294, 330,18, 29, 81,57, 81, 193,68, 108, 172,91, 98, 171,187, 152, 427,30, 30, 56,49, 34, 65,37, 33, 71,1189, 1195, 1403,67, 96, 92,25, 43, 106), .Dim=c(30,3)).)
model
{
# shared clustering terms
us[1:I] ~ car.normal(mapa[],weia[],numa[],tau.us)
# specific clustering terms
up[1,1:I] ~ car.normal(mapa[],weia[],numa[],tau.u[1])
up[2,1:I] ~ car.normal(mapa[],weia[],numa[],tau.u[2])
# heterogeneity terms
for (k in 1:2){
for (i in 1:I){
v[k,i]~ dnorm(0.00,tau.v[k])
}
}
# i area, k disease
for (i in 1:I) {
osse[i,1:K] ~ dmulti(p[i,1:K],n[i] )
n[i] <- sum(osse[i,])
for (k in 1:K) {
p[i,k] <- phi[i,k] / sum(phi[i,])
log(phi[i,k]) <-alpha[k]+u[k,i]+v[k,i]
devi.each[i,k]<- (osse[i,k]+0.05)*log((osse[i,k]+0.05)/(n[i]*p[i,k]+0.05))
}
}
# deviance
dev <- 2*(sum(devi.each[,]))
# clustering terms (shared + specific)
for (i in 1:I){
u[1,i]<-us[i]*delta+up[1,i]
u[2,i]<-us[i]/delta+up[2,i]
}
logdelta ~ dnorm(0,5.9)
delta <- exp(logdelta)
# disease number 3 (lung) is the reference
alpha[3]<-0
for (i in 1:I){
u[3,i]<-0
v[3,i]<-0
}
for (i in 1:I)
{
for (k in 1:2)
{ exp.up[k,i]<-exp(up[k,i])}
exp.us[i]<-exp(us[i])
}
# Weia
for(j in 1:sumNumNeig)
{weia[j]<-1}
for (k in 1:2){
alpha[k] ~ dflat()
tau.v[k] ~ dgamma(0.5,0.0005)
tau.u[k] ~ dgamma(0.5,0.0005)
}
tau.us ~ dgamma(0.5,0.0005)
}
list( I=30, numa = c(3, 3, 3, 4, 4, 8, 6, 6, 4, 3,
6, 4, 3, 5, 4, 4, 4, 3, 5, 5,
4, 7, 6, 5, 4, 8, 3, 5, 6, 7),
mapa= c(30, 9, 3,
30, 21, 3,
30, 2, 1,
20, 19, 12, 7,
22, 20, 19, 6,
29, 26, 25, 23, 22, 20, 7, 5,
29, 20, 14, 12, 6, 4,
30, 28, 24, 23, 11, 9,
30, 24, 8, 1,
26, 24, 18,
30, 23, 22, 21, 15, 8,
27, 14, 7, 4,
22, 19, 15,
29, 27, 16, 12, 7,
22, 21, 13, 11,
29, 27, 17, 14,
29, 26, 18, 16,
26, 17, 10,
22, 20, 13, 5, 4,
19, 7, 6, 5, 4,
30, 15, 11, 2,
23, 19, 15, 13, 11, 6, 5,
28, 25, 22, 11, 8, 6,
28, 26, 10, 9, 8,
28, 26, 23, 6,
29, 28, 25, 24, 18, 17, 10, 6,
16, 14, 12,
26, 25, 24, 23, 8,
26, 17, 16, 14, 7, 6,
21, 11, 9, 8, 3, 2, 1),
sumNumNeig = 142,
I=30,
K=3,
osse=structure(.Data=c(996,995,2806,645,755,1613,298,256,1094,1661,1518,1260,84,51,160,177,199,101,5387,3603,4674,
224,235,329,59,130,118,1820,1380,3243,93,94,218,944,905,921,198,202,514,218,218,348,160,166,214,1039,1364,1216,285,223,475,280,203,376,314,375,1001,427,619,586,532,511,690,100,122,265,495,391,803,1138,901,1512,344,412,742,
1043,674,1965,299,344,409,178,160,154,364,489,670,327,297,275), .Dim=c(30,3))
)
#intial
list(tau.us=0.5,
alpha=c(0.5,0.5,NA),
tau.v=c(0.5,0.5),
tau.u=c(0.5, 0.5),
logdelta=0,
v=structure(.Data=c(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA), .Dim=c(3,30)))
)
Thanks for any help!
-------------------------------------------------------------------
This list is for discussion of modelling issues and the BUGS software.
For help with crashes and error messages, first mail [log in to unmask]
To mail the BUGS list, mail to [log in to unmask]
Before mailing, please check the archive at www.jiscmail.ac.uk/lists/bugs.html
Please do not mail attachments to the list.
To leave the BUGS list, send LEAVE BUGS to [log in to unmask]
If this fails, mail [log in to unmask], NOT the whole list
|