本人新手,遇到一个问题请教大家。对模型进行compile时程序说“index out of range
SpatialCARNormal.Node.Set [00000D9CH]
.args GraphStochastic.Args Fields
.beg INTEGER 6
.car SpatialCARNormal.Node [016E4B20H]
.i INTEGER 0
.index INTEGER 4
.node SpatialCARNormal.Node [016E4B60H]
.numIslands INTEGER 1684805647
.numNeigh INTEGER 2
.off INTEGER -1
.p GraphNodes.Node [016E2C10H]
.res INTEGER 0
.size INTEGER 6
.start0 INTEGER”
其中,模型为:
#model3: Spatial + random effect
Model
{
#Likelihood
for (i in 1:6) {
Severe ~ dpois(mu)
log(mu) <- beta[1] + beta[2] * Inter_densitiy + beta[3] * lslseglen + beta[4] * hslseglen + beta[5] * HBWP + beta[6] * NHBWP + theta + phi
theta[i] ~ dnorm(0.0,tau.h)
xi[i] <- theta[i] + phi[i]
}
# Prior
for (k in 1:6) {beta[k] ~ dnorm(0.0,1.0E-5)}
phi[1:6] ~ car.normal (adj[], weights[], num[], tau.c)
tau.h ~ dgamma(0.001, 0.001) #precision: random effect
tau.c ~ dgamma(0.1,0.1) # precision: spatial
sd.h <- sd(theta[])
sd.c <- sd(phi[])
alpha <- sd.c/(sd.h + sd.c)
}
数据为:
list(
adj=
c(4, 5, 2, 1, 3, 2, 0, 5, 0, 4),
num=
c(2, 1, 2, 1, 2, 2),
weights=
c(0.21, 0.20, 0.14, 0.16, 0.14, 0.16, 0.23, 0.21, 0.23, 0.20))
请大家帮我看看是怎么回事,十分感谢