WINBUGS软件对SV模型的参数估计

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

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）
请大家帮帮忙看看，谢谢啊

扫码加我 拉你入群

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

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

分享0 收藏4 回帖

关键词：Winbugs软件 winbugs WINBUG SV模型 BUGS 模型 软件

有点看不懂。。。

heilongwxl 发表于 2012-3-20 21:02
有点看不懂。。。

你是不懂这个软件，还是不懂我编的这个程序

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模型；但我看不懂中间的过程。
求高手帮助

我也看不懂  程序看不懂  但我知道好像差不多

我试过运行你的程序，是可以运行的。load inits时，会告知“this chain contains unitialized variables”，这不是什么错，你只要再点gen inits键即可。

我只是看你这个模型太累，能解释为何如此设分布和变量？

金融坦然 发表于 2012-3-27 20:57
model{
for(i in 1:n)
{p

我觉得看你这个模型太累，能解释为何如此设分布和变量？我有个简单一点的WINBUGS程序，也是对付SV的参数估计。不过似乎你的这个要复杂的多。

TaskShare 发表于 2012-10-20 16:46
我觉得看你这个模型太累，能解释为何如此设分布和变量？我有个简单一点的WINBUGS程序，也是对付SV的参数估 ...

你好 请问你可以把你的SV参数估计程序借我用一下吗  谢谢

TaskShare 发表于 2012-10-20 16:46
我觉得看你这个模型太累，能解释为何如此设分布和变量？我有个简单一点的WINBUGS程序，也是对付SV的参数估 ...

请问你winbugs做SV的程序还有吗？我可以参考一下吗？新学者，完全一头雾水

偏斜分布SN，在winbugs中怎么表示，有个参数是服从sn分布的，