• 综合主站

  • WinBUGS在统计分析中的应用（第一部分）

[未知用户] 数字格式不对。参见下面的回复。

To齐韬：你可能需要修改正文中的代码，我在OpenBUGS 3.0.7中运行你的代码会出错，错误信息和这里一样。我检查了一下，是数字格式的问题。
请问winbugs能计算递归吗？想用之解决一个随机变量自我n重卷积的问题。用递归来算，比如一个正态分布s~(mu,tau)，自我卷积n次求结果。该怎么弄呢？谢谢大虾们指教。
[未知用户] WinBUGS是个编程工具么？……不了解。你自己Google一下吧。
1 个月 后
[未知用户] 您好，请问进行MH算法需要另外编辑程序吗？
[未知用户] 不需要，要是MH都要编程，那WinBUGS就没有存在必要了……呃，确切地说，我说的是OpenBUGS……我从没用过WinBUGS，因为它已经确定没有继续的技术支持了，它的开发已经停止。OpenBUGS还在更新中。
6 天 后
请教一个问题：我运行第五步，load data的时候，软件的左下角显示expected variable name, 然后点击compile的时候，软件界面上没有任何变化，load inits和gen Init按钮是灰色的，不能点击运行。是怎么回事啊，请各位朋友帮帮忙，谢谢了。
[未知用户] 你要选中所有的数据，然后点load data，或者至少你要选中list这个词，这样BUGS才知道数据的起止。如果选择范围不恰当，BUGS就不能正确载入数据。

另一个可能性是你的模型中使用了数据没有提供的变量名，这种情况下，你要自己检查模型中的变量，如果有遗漏，就在数据的list()中加上。
5 天 后
model
{
for (i in 1 : I) {
for(j in 1 : J) {
for(k in 1:K) {
epsilon[i,j,k] ~ dnorm(0.0,tau)
y[i,j,k] <-alpha[i,k] + beta[i,k] * theta[j,k] + epsilon[i,j,k]
}
}
}
...
list(
I = 5,J = 5,K = 4,
y = structure(
.Data = c(
6.8,7.2,7.35,7.63,

.Dim = c(5,5,4)
)
)

编译出现multiple definitions of node y[1,1,1]
如果去掉数据部分的y定义就可以编译过去。
希望大家帮我解决下这个问题，万分感谢！
[未知用户] 去掉y即使模型语法正确，程序肯定也无法运行。

你再仔细看看你这里y的定义，c的回括号打到哪里去了？

还有，一共4个数字，请问该如何定义为一个5x5x4的数组？
model                                                               
{  
#distribution of Ys                                                             
###################                                                             
for (i in 1:N) {
    ysigmadet<-exp(th[i,1]+th[i,2])*(1-rhoep*rhoep);        
    Yisigma2[i,1,1] <- exp(th[i,2])/ysigmadet;
    Yisigma2[i,2,2] <- exp(th[i,1])/ysigmadet;
    Yisigma2[i,1,2] <- -rhoep*exp(0.5*th[i,1]+0.5*th[i,2])/ysigmadet;             Yisigma2[i,2,1] <- Yisigma2[i,1,2];
    Y[i,1:2]~ dmnorm(muy[],Yisigma2[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);
sig1[1]<-exp(0.5*th[1,1]);
sig2[1]<-exp(0.5*th[1,2]);
q[1]~dnorm(psi0,itau2);
rhoep[1]<-(exp(q[1])-1)/(exp(q[1])+1);
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);
		th[i,1]~dnorm(thmean[i,1],itaua2);
		th[i,2]~dnorm(thmean[i,2],itaub2);
		sig1<-exp(0.5*th[i,1]);
		sig2<-exp(0.5*th[i,2]);
		qmean<-psi0+psi*(q[i-1]-psi0);
		q[i]~dnorm(qmean[i],itau2);
		rhoep[i]<-(exp(q[i])-1)/(exp(q[i])+1);
   }                                                                            
#distribution of phi, mu, rhoep                                                    
###########################                                                     
phi1star ~ dbeta(20,1.5);   
phi1 <- 2*phi1star -1;  
phi2star ~ dbeta(20,1.5);   
phi2 <- 2*phi2star -1;
psistar ~ dbeta(20,1.5);   
psi <- 2*psistar -1;

itaua2 ~ dgamma(2.5,0.025);                                                      
taua <- sqrt(1/itaua2);
itaub2 ~ dgamma(2.5,0.025);                                                      
taub <- sqrt(1/itaub2);                                             
itau2 ~ dgamma(2.5,0.025);                                                      
tau <- sqrt(1/itau2);                                             
mu1 ~ dnorm(0,0.04);
mu2 ~ dnorm(0,0.04);
psi0~dnorm(0.7,0.1);
} 
list(phi1star=0.99, phi2star=0.99,mu1=0,mu2=0,itaua2=100,itaub2=100,psistar=0.99,psi0=1.9,itau2=100)
这是一个DC-MSV模型，大家能帮我看一下，到底是哪个初始值没有定义吗？
在编译完成之后，提示说有未定义的初始值，是怎么回事？
急急急，先谢过各位了！！！！！
[未知用户] 这是个体力活儿，自个儿把模型中所有的参数都列出来，然后对比一下是不是所有的参数都有初始值吧
20 天 后
[未知用户] 请问就是把您给的代码复制进WINBUGS,到compile时,出现educational version cannot do this model.什么原因阿?对这个软件不甚掌握,论文却急用请指教,谢谢
[未知用户] 用OpenBUGS吧，我从不用WinBUGS，后者早已经停止开发了
3 个月 后
每次启动winBUGS1.4.3都要弹出Licence Agreement，这是不是没安装好呢？
第七步时，它说edicational version cannot do this model，是不是必须注册了才行呀
[未知用户] 用OpenBUGS！！WinBUGS早已经停止开发了。OpenBUGS开源、免费、还在开发更新中。
16 天 后
大家好 OpenBUGS和WINbugs的程序一样吗
[未知用户] 这个是什么的程序呀 SO 长