1

我正在尝试在 OpenBUGS 中运行马尔可夫模型来估计过渡概率。但我收到一个错误“数组索引不是整数”这是我的模型代码：

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)

谁能帮我解决错误？

4

0 回答 0