全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
5888 9
2012-03-20
model
{
for (t in 1:n) {
yisigma2[t]<-1/exp(theta[t])
y[t] ~ dnorm(0,yisigma2[t])
                    }
mu~dnorm(0,0.1)
phistar~dbeta(20,1.5)
itau2~dgamma(2.5,0.025)
beta<-exp(mu/2)
phi<-2*phistar-1
tau<-sqrt(1/itau2)
theta0~dnorm(mu,itau2)
thmean[1]<-mu+phi*(theta0-mu)
theta[1]~dnorm(thmean[1],itau2)
for(t in 2:n) {
thmean[t]<-mu+phi*(theta[t-1]-mu)
theta[t]~dnorm(thmean[t],itau2)
                   }
}
#data
list( n=20,y=c(-0.0017,-0.0093,-0.0113,-0.0079,0.0016,0.0024,0.0100,0.0139,0.0025 ,
-0.0076,0.0105,0.0043,0.0022,0.0092,0.0038, 0.0025,0.0031,0.0024,-0.0103 ,0.0092))
#init
list(phistar=0.95,mu=0,itau2=50)


上面是我的程序,为什么我点击(load inits)左下角的状态显示(expected  variable name)
请大家帮帮忙看看,谢谢啊
二维码

扫码加我 拉你入群

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

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

全部回复
2012-3-20 21:02:47
有点看不懂。。。
二维码

扫码加我 拉你入群

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

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

2012-3-21 10:48:56
heilongwxl 发表于 2012-3-20 21:02
有点看不懂。。。
你是不懂这个软件,还是不懂我编的这个程序
二维码

扫码加我 拉你入群

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

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

2012-3-27 20:57:16
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[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)
}
}
使用这个程序能估计SV模型;但我看不懂中间的过程。
求高手帮助
二维码

扫码加我 拉你入群

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

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

2012-5-14 11:11:11
我也看不懂  程序看不懂  但我知道好像差不多
二维码

扫码加我 拉你入群

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

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

2012-10-20 16:12:02
我试过运行你的程序,是可以运行的。load inits时,会告知“this chain contains unitialized variables”,这不是什么错,你只要再点gen inits键即可。

我只是看你这个模型太累,能解释为何如此设分布和变量?
二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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