openbugs做突变点贝叶斯推断，报错：multiple definitions of node xt[1]

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

求助啊啊
openbugs做ARMA（1，1）突变点贝叶斯推断，一直报错：multiple definitions of node xt[1]

代码
library(R2OpenBUGS)
N=1000
xt<-c(x,y)
e<-c(e1,e2)
data<-list("N","xt","e")
parameters<-c("k","alpha","sigma","beta")
inits<-function()
{
  list(k=500,alpha=c(0.2,0.6),beta=c(0.4,0.8),sigma=c(0.1,0.1))
}
model<-function(){
  xt[1]<-0.01
  k~dunif(1,N)
  for(i in 1:N){
    e[i]~dnorm(0,sigma[J[i]])
    xt[i]<-alpha[J[i]]*xt[i-1]+beta[J[i]]*e[i-1]+e[i]
    J[i]<-1+step(i-k)
  }
  for(j in 1:2){
    alpha[j]~dnorm(0,0.001)
    beta[j]~dnorm(0,0.001)
    sigma[j]~dnorm(0,0.001)
  }

}
write.model(model,"model.txt")
model.sim<-bugs(data,inits,parameters,"model.txt",n.chains = 3,n.iter = 10000,debug=TRUE)



虽然觉得做伸手党不好，但是已经卡了好几天了，要死了TAT
求助大神

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：DEFINITIONS definition openbugs Multiple multip multiple

我检查过数据的长度，加过xt1的定义，就从x1说重复，去掉x1就说从x2开始重复TAT
跪求给点思路，真的不知道错在哪里了TAT

楼主，您好可以加个联系方式一起讨论讨论吗，我最近也在做这个研究