全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 winbugs及其他软件专版
3038 5
2014-04-08
悬赏 5 个论坛币 已解决
新接触这个计量软件,想估计带跳的CKLS短期利率动态模型,自己攥的程序如下:
model {
k~dnorm(0,1)
mur~dnorm(2.3,0.25)
gama~dnorm(1,0.25)
a~dgamma(2.5,0.025)
sigmar2<-1/a
kesai~dnorm(u,tau1)
u~dnorm(2,0.25)
tau1<-1/sigmaj2
b~dgamma(2.5,0.025)
sigmaj2<-1/b
zt~dnorm(0,1)
J~dpois(lamda)
lamda~dnorm(0,1)
for (i in 1:N-1) { r[i+1]<-r+k*(mur-r)+sigmar2*pow(r,gama)*zt+J*kesai}
}

data
list(N=10, r=c(2.4155,2.4116,2.4093,2.4101,2.4041,2.3994,2.3984,2.3954,2.3830,2.3400  )
)

inits
list(k=0.01,mur=2.3,sigmar2=0.25,gama=1,u=1.5,sigmaj2=0.25,zt=0.01,lamda=0.01)



可是编译的时候,出现了multiple definition of node r[2]的问题,不知何解,各位大神,请不吝赐教!拜谢!!

最佳答案

TimeT 查看完整内容

程序改了下,见多处#号的解释。(声明:我未花时间研究你程序是否达到经济目的,只看程序是否能运行,改后能运行。) 你的最重要的问题就是你的那句For后面的语句,r
二维码

扫码加我 拉你入群

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

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

全部回复
2014-4-8 15:23:59
程序改了下,见多处#号的解释。(声明:我未花时间研究你程序是否达到经济目的,只看程序是否能运行,改后能运行。)
你的最重要的问题就是你的那句For后面的语句,r[i+1]<-...,这样是确定变量(不是随机变量)的赋值,你已用此句给r[i+1]都赋值了,在load data时r[...]又赋值,就重复赋值了。我是以为你的r[i+1]是随机变量,所以改成随机变量的“赋值”方式了。
水平有限,你试试。

model {
k~dnorm(0,1)
mur~dnorm(2.3,0.25)
gama~dnorm(1,0.25)
a~dgamma(2.5,0.025)
sigmar2<-1/a
kesai~dnorm(u,tau1)
u~dnorm(2,0.25)
tau1<-1/sigmaj2
b~dgamma(2.5,0.025)
sigmaj2<-1/b
# zt~dnorm(0,1) # deleted, because the changes made within the for loop below
J~dpois(lamda)
lamda~dnorm(0,1)I(0.001,) #in J~dpois(lamda), lamda must>0 so add I(0.001, ) to ensure that
for (i in 1:N-1) {
#r[i+1]<-r[i ]+k*(mur-r[i ])+sigmar2*pow(r,gama)*zt+J*kesai # deleted to change to the following 4 lines
rsTemp[i ]<- sigmar2*pow(r[i ],gama) # added
rTauTemp[i ]<- 1/rsTemp[i ]/rsTemp[i ] # added
rmuTemp[i ] <- k*(mur-r[i ])+J*kesai +r[i ] # added
r[i+1]~dnorm(rmuTemp[i ],rTauTemp[i ]) # added
        }
}


data
list(N=10, r=c(2.4155,2.4116,2.4093,2.4101,2.4041,2.3994,2.3984,2.3954,2.3830,2.3400  )
)

inits
list(k=0.01,mur=2.3,a=4,gama=1,u=1.5,b=4,lamda=0.01) #cannot set the following as inits, so delete: sigmar2=0.25,sigmaj2=0.25,zt=0.01
二维码

扫码加我 拉你入群

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

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

2014-4-10 10:53:17
程序发现了小错误,继续改。。。磨人啊。。
二维码

扫码加我 拉你入群

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

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

2014-4-12 15:17:04
TimeT 发表于 2014-4-12 13:45
程序改了下,见多处#号的解释。(声明:我未花时间研究你程序是否达到经济目的,只看程序是否能运行,改后能 ...
感谢!按照您的建议试了一下,问题解决了!谢谢!
二维码

扫码加我 拉你入群

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

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

2014-4-12 15:22:20
TimeT 发表于 2014-4-12 13:45
程序改了下,见多处#号的解释。(声明:我未花时间研究你程序是否达到经济目的,只看程序是否能运行,改后能 ...
第一次发帖,不知道那5个悬赏论坛币怎么给您啊?汗。。。。
二维码

扫码加我 拉你入群

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

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

2018-11-7 23:44:29
TimeT 发表于 2014-4-12 13:45
程序改了下,见多处#号的解释。(声明:我未花时间研究你程序是否达到经济目的,只看程序是否能运行,改后能 ...
大神我也遇到一样的问题了
1541605471409497.jpeg
二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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