楼主: 醉风尘
1960 4

[学科前沿] WinBUGS作GC-MSV参数估计时stat是灰色的 [推广有奖]

  • 0关注
  • 1粉丝

大专生

11%

还不是VIP/贵宾

-

威望
0
论坛币
623 个
通用积分
0.7500
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
70 点
帖子
1
精华
0
在线时间
88 小时
注册时间
2018-4-29
最后登录
2023-11-9

100论坛币
在用winbugs作GC-MSV的参数估计时,"chaeck model"和"load data"都没问题,程序也可以跑得起来,也set参数了,但是最后参数的"stat"选项的按钮是灰色的怎么回事呢?希望有高手帮忙解答,非常感谢!
程序代码如下:
model
{
for (i in 1:N) {
       
        Ysigma2[i,1,1] <- exp(th[i,2]);
        Ysigma2[i,2,2] <- exp(th[i,1])
        Ysigma2[i,1,2] <- rhoep*exp(0.5*th[i,1]+0.5*th[i,2])
        Ysigma2[i,2,1] <- Ysigma2[i,1,2];
        Y[i,1:2] ~ dmnorm(muy[],Ysigma2[i,,]);
        }
muy[1] <- 0;
muy[2] <- 0;
thmean[1,1] <- mu1;
thmean[1,2] <- mu2;
th[1,1] ~ dnorm(thmean[1,1],itaua2);
th[1,2] ~ dnorm(thmean[1,2],itaub2);
for (i in 2:N) {
        thmean[i,1] <- mu1 + phi1*(th[i-1,1]-mu1);
        thmean[i,2] <- mu2 + phi2*(th[i-1,2]-mu2)+phi21*(th[i-1,1]-mu1);
        th[i,1] ~ dnorm(thmean[i,1],itaua2);
        th[i,2] ~ dnorm(thmean[i,2],itaub2);
}
phi1star ~ dbeta(20,1.5)
phi1 <- 2*phi1star - 1;
phi2star ~ dbeta(20,1.5)
phi2 <- 2*phi2star - 1;
phi21 ~ dnorm(0,0.1);

itaua2 ~ dgamma(2.5,0.025);
taua <- sqrt(1/itaua2);
itaub2 ~ dgamma(2.5,0.025);
taub <- sqrt(1/itaub2);
mu1 ~ dnorm(0,0.04);
mu2 ~ dnorm(0,0.04);
rhoep ~ dunif(-1,1);
}

list(phi1star=0.99, phi2star=0.99, mu1=0, mu2=0, itaua2=100, itaub2=100, rhoep=0.7, phi21=0)
list(N = 5,Y=structure(.Data=c(0.001838 ,-0.000133 ,-0.000378 ,0.000278 ,-0.005169 ,0.000001 ,0.004320 ,-0.000086 ,0.003322 ,0.000086),
.Dim=c(5,2))
)
实际数据有两千多个,这里只是摘了一部分,是模型代码写的吗,还是数据或者初值定义的问题?

关键词:winbugs GC-MSV模型 代码 编程
沙发
15997351897 发表于 2020-3-7 16:43:23 |只看作者 |坛友微信交流群
请问程序代码是自己编的吗?

使用道具

藤椅
15997351897 发表于 2020-3-7 16:44:17 |只看作者 |坛友微信交流群
winbugs模型先验分布参数设置怎么写代码啊?

使用道具

板凳
oopp0099 发表于 2020-7-31 00:45:43 来自手机 |只看作者 |坛友微信交流群
醉风尘 发表于 2019-11-30 00:20
在用winbugs作GC-MSV的参数估计时,&quot;chaeck model&quot;和&quot;load data&quot;都没问题,程序也可以跑得起来,也set参数了 ...
qq留下联系

使用道具

报纸
JiaYuxuan 学生认证  发表于 2022-4-19 15:07:27 |只看作者 |坛友微信交流群
应该是迭代次数设置的太小了,增加迭代次数,stats统计信息就可以查看了

使用道具

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

本版微信群
加好友,备注jltj
拉您入交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-4-27 07:44