结构突变贝叶斯分析关于GIBBS抽样openbugs编程问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

各位大神，最近想基于小样本做一个结构突变，研读东北财经大学硕士论文《我国GDP与碳排放关系的结构突变贝叶斯分析》时，利用文章的程序来做，但发现初值怎么赋值都有问题，实在不知道是什么地方出了问题，求各位大神帮帮忙啊！
model
{
for(i in 3:N){
y~dnorm(mu,tau[J])
mu<-alpha[J]+beta[J]*i+gamma*x1+lamda*x2
J<-1+step(i-x.change1)+step(i-x.change2)
x1<-y[i-1]
x2<-y[i-2]
}
for(j in 1:3){
alpha[j]~dnorm(0,0.001)
beta[j]~dnorm(0,0.001)
tau[j]~dgamma(1.001,0.001)
}
for(i in 1:N-2){
punif<-1/N-2
}
k[1]~dcat(punif[])
k[2]~dcat(punif[])
x.change1<-min(k[1],k[2])+1
x.change2<-max(k[1],k[2])+1
gamma~dnorm(0,0.001)
lamda~dnorm(0,0.001)
}
list(N=31,y=c(1439.85965,1506.94235,1593.38617,1724.48938,1765.57417,1878.66476,2010.69812,2148.31963,2183.13826,2177.70031,2295.66967,2375.6207,2498.76806,2681.7985,722.71694,2841.26338,3129.77766,3197.55622,3115.97441,3271.8091,3353.64009,3776.82765,4235.69684,4744.9877,5463.70378,5936.03632,6326.36456,6684.65068,7573.37999,7997.04373,8715.30724))#y值是1981-2011年中国二氧化碳排放量，数据来源：eia
list(gamma=0,lamda=0,alpha=c(0,0,0),beta=c(0,0,0),tau=c(0,0,0),k=c(1,28))
赋初值时显示错误error for node k[1] of type GraphCat,StdNode first argument must be a proportion

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：openbugs Gibbs抽样 Gibbs 贝叶斯分析 结构突变 结构 东北财经

帮你顶一下吧

大神问题最终解决了吗？

请问最后做出来了吗？

请问作出了吗 愿意有偿求教

结构突变点有人会吗 愿意有偿求教    qq 2587713066