- 阅读权限
- 255
- 威望
- 0 级
- 论坛币
- 623 个
- 通用积分
- 0.7500
- 学术水平
- 0 点
- 热心指数
- 0 点
- 信用等级
- 0 点
- 经验
- 70 点
- 帖子
- 1
- 精华
- 0
- 在线时间
- 88 小时
- 注册时间
- 2018-4-29
- 最后登录
- 2023-11-9
大专生
还不是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))
)
实际数据有两千多个,这里只是摘了一部分,是模型代码写的吗,还是数据或者初值定义的问题? |
|