楼主: jiawei5990
2851 3

求助!!关于用MCMC方法估计 MSSV的问题 [推广有奖]

  • 0关注
  • 1粉丝

本科生

82%

还不是VIP/贵宾

-

威望
0
论坛币
485 个
通用积分
25.2811
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
588 点
帖子
40
精华
0
在线时间
173 小时
注册时间
2013-1-30
最后登录
2025-5-4

楼主
jiawei5990 发表于 2014-11-29 23:25:31 |AI写论文
100论坛币
用openbugs做MSSV的参数估计,但迭代到后来发现参数的轨迹是这样的情况 QQ截图20141129232155.png 还有几个参数完全轨迹是平的  而且自相关系数一直很高 100阶了还自相关。。。。。不知道怎么回事,求高人解答啊 附上代码 自己照着SV模型的写的 不知道有没问题。。。。。
# prior distributions
alpha0~dnorm(0,0.1);
alpha1~dnorm(0,0.1);
p00~dbeta(20,1.5);
p11~dbeta(20,1.5);
phistar~dbeta(20,1.5);
insigma~dgamma(2.5,0.025);
phi<-2*phistar-1;
sigma<-sqrt(1/insigma);
Q[1,1]<-p00;
Q[1,2]<-1-p00;
Q[2,1]<-1-p11;
Q[2,2]<-p11;
s0~dbern(0.5);
s[1]~dbern(0.5);
x[1]<-s0+1;
V[1,1]<-0.5;
V[1,2]<-0.5;
M[1]<-0.5;
for (t in 2:2218)
{x[t]<-s[t-1]+1;
V[t,1]<-Q[x[t],1];
V[t,2]<-Q[x[t],2];
M[t]<-V[t,2];
s[t]~dbern(M[t]);
}


h0~dnorm(alpha0,insigma);
hmean[1]<-(1-phi)*alpha0+(s[1]-phi*s0)*alpha1+phi*h0;
h[1]~dnorm(hmean[1],insigma);
for (t in 2:2218) {hmean[t]<-(1-phi)*alpha0+(s[t]-phi*s[t-1])*alpha1+phi*h[t-1];
h[t]~dnorm(hmean[t],insigma);
}

关键词:mcmc MSS CMC distribution openbugs 而且 模型

沙发
jiawei5990 发表于 2014-11-29 23:26:34
不知道有没谁懂这个软件语言的  我也不知道 自己赋值的时候有没问题。。。。。。

藤椅
fucuo124651 发表于 2021-2-20 16:07:56 来自手机
你好,请问你有没有MSSV的open bugs代码?

板凳
tigerbang 发表于 2021-4-7 12:28:14
你好,请问你后面解决了吗?我想用winbugs写MSSV,代码有些和你的一样,但是还有些没写出来,可以请教下吗?QQ:1565891884  万分感谢!

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注jltj
拉您入交流群
GMT+8, 2025-12-31 09:33