s**u
发帖数: 383 | 1 经高手指导，compile 没有问题了， 但是initial 又出问题了， 说是NIL
dereference (read)， 查了manual, 说是 'NIL dereference (read)' can occur at
compilation in some circumstances when an inappropriate transformation is
made, for example an array into a scalar.用gen initial也出现同样的error.
请高人指导，谢谢。
model {
for (i in 1:k-1) {
for (j in i+1:k){
win1[i,j] ~ dbin(p[i,j], n)
logit(p[i,j]) <- d[i]-d[j]}
# p[i] <- phi(d[ind1[i]] - d[ind2[i]])
}
d[1:k] ~dmnorm(dmu[ ], prec.Sigma[,])
## Inverse Wishart Prior for Sigma
for (i in 1:k) {dmu[i] <-0}
prec.Sigma[1:k,1:k] ~dwish(Omega[,],k)
Sigma[1:k,1:k]<-inverse(prec.Sigma[,])
for (i in 1:k){ Omega[i,i] <- 0.1}
}
Data:
list(n= 4, k= 4, win1= structure(.Data = c(0, 4, 4, 4,
0, 0, 0, 1,
0, 4, 0, 3,
0, 3, 1, 0), .Dim = c(4, 4)))
inits
list(d=c(0.1,0.1,0.1,0.1), prec.Sigma= structure(.Data = c(1, 0, 0, 0,
0, 1,0, 0,
0, 0, 1, 0,
0, 0, 0, 1), .Dim = c(4, 4))) | q*****q
发帖数: 158 | 2 （k-1）
(i+1)
先把这些加上括号再看看。。。
at
【在 s**u 的大作中提到】
: 经高手指导，compile 没有问题了， 但是initial 又出问题了， 说是NIL
: dereference (read)， 查了manual, 说是 'NIL dereference (read)' can occur at
: compilation in some circumstances when an inappropriate transformation is
: made, for example an array into a scalar.用gen initial也出现同样的error.
: 请高人指导，谢谢。
: model {
: for (i in 1:k-1) {
: for (j in i+1:k){
: win1[i,j] ~ dbin(p[i,j], n)
: logit(p[i,j]) <- d[i]-d[j]}
|
|
