全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
4663 6
2010-09-04
最近在做论文,写了个SV模型如下,load data时我选了数据,但运行还是显示expected variable name,请各位高手帮忙,谢谢~~
# model{
for (i in 1:n)
{p[i]<-1/exp(theta[i])
y[i]~dnorm(0,p[i])
}
mu~dnorm(0,0.01)
phi1~dbeta(20,1.5)
itau2~dgamma(2.5,0.025)
phi<-2*phi1-1
tau<-sqrt(1/itau2)
theta0~dnorm(mu,itau2)
thmean[i]<-mu+phi*(theta0-mu)
theta[i]~dnorm(thmean[1],itau2)
for (j in 2:n)
{thmean[j]<-mu+phi*(theta[j-1]-mu)
theta[j]~dnorm(thmean[j],itau2)
}
}

list (y=c(5.917,6.0587,5.9334,5.8906,5.8189,6.0395,5.9027,6.0142,6.468,6.3844,6.4806,6.7442,6.9832,9.3499,9.1096,9.0589,9.1446,9.5877,8.1597),n=19)
二维码

扫码加我 拉你入群

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

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

全部回复
2010-9-4 16:55:59
model{
### likelihood: joint distribution of y
for (i in 1:n)
{p[i]<-1/exp(theta[i])
y[i]~dnorm(0,p[i])
}
### prior distributions
phi1~dbeta(20,1.5)
phi<-2*phi1-1
mu~dnorm(0,0.01)
itau2~dgamma(2.5,0.025)
tau<-sqrt(1/itau2)

theta0~dnorm(mu,itau2)
thmean[1]<-mu+phi*(theta0-mu)
theta[1]~dnorm(thmean[1],itau2)
for (j in 2:n)
{thmean[j]<-mu+phi*(theta[j-1]-mu)
theta[j]~dnorm(thmean[j],itau2)
}
}
#data
list(n=19,y=c(5.917,6.0587,5.9334,5.8906,5.8189,6.0395,5.9027,6.0142,6.468,6.3844,6.4806,6.7442,6.9832,9.3499,9.1096,9.0589,9.1446,9.5877,8.1597))
二维码

扫码加我 拉你入群

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

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

2010-9-4 17:51:50
谢谢~~这位大哥真是好人~~ 我刚试了,可以了
2# epoh
二维码

扫码加我 拉你入群

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

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

2010-9-15 15:57:21
我是初学者,请较各位高手,在基本SV模型中,怎样进行波动率预测?是theta还是exp(Theta/2)。。。还有我在winBUGS运行SV模型后点DIC 为什么是cannot calculate DIC for model,不知道是什么原因?各位高手帮帮忙~~
二维码

扫码加我 拉你入群

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

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

2011-1-4 10:20:19
都是高手啊,我要好好学习.
二维码

扫码加我 拉你入群

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

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

2011-8-10 19:58:19
y80820201 发表于 2010-9-15 15:57
我是初学者,请较各位高手,在基本SV模型中,怎样进行波动率预测?是theta还是exp(Theta/2)。。。还有我在 ...
不能计算DIC有两个原因,一是参数估计时不收敛,二是模型太复杂,或者分布很复杂,严重偏离正态分布,比如是混合正态分布,也不能计算DIC的值
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

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