OpenBUGS中的数组索引不是整数错误

时间:2016-07-29 14:22:05

标签: bayesian markov-chains winbugs r2winbugs openbugs

我正在尝试在OpenBUGS中运行Markov模型来估计过渡概率。但我纠正错误“数组索引不是整数” 这是我的型号代码:

model{
p[1,1] ~ dnorm(0,1)
p[2,2] ~ dnorm(0,1)
p[3,3] ~ dnorm(0,1)
p[1,2] ~ dnorm(0,1)
p[2,1] ~ dnorm(0,1)
p[3,1] ~ dnorm(0,1)
p[1,3] <- 1-p[1,1]-p[1,2]
p[2,3] <- 1-p[2,1]-p[2,2]
p[3,2] <- 1-p[3,1]-p[3,3]

    for(i in 1:N){
      for(t in 2:T){
        y[i,t]~dcat(p[y[i,t-1],])
}
}
}

以下是我的R代码数据:

Data=structure(.Data=c(0.1, 0.2, 0.5, 0.2, 0.2, 0.8, 0.2, 0.4, 0.3, 0.2, 0.1, 0.6, 0.1, 0.5, 0.6), .Dim=c(5,3))
N=nrow(Data)
T=ncol(Data)
y=as.matrix(Data)

#Load the data
sp.data = list(y=y, N=N, T=T)

Sp.sim0<- bugs (sp.data, inits=NULL, 
       parameters=c('p'),                       
                   model.file= "model2.txt",n.chains=2,  DIC=TRUE,n.iter=1000, codaPkg=FALSE,debug=TRUE)

print(Sp.sim0)
attach.bugs(Sp.sim0)

有谁可以帮助我解决错误?

0 个答案:

没有答案