全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 winbugs及其他软件专版
1973 2
2016-03-13
求助啊啊
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
求助大神
二维码

扫码加我 拉你入群

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

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

全部回复
2016-3-13 22:16:19
我检查过数据的长度,加过xt1的定义,就从x1说重复,去掉x1就说从x2开始重复TAT
跪求给点思路,真的不知道错在哪里了TAT
二维码

扫码加我 拉你入群

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

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

2022-11-7 21:48:15
楼主,您好可以加个联系方式一起讨论讨论吗,我最近也在做这个研究
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群