楼主: flashcrazy
1576 3

winbugs模型设定问题 [推广有奖]

  • 0关注
  • 9粉丝

VIP

已卖:560份资源

讲师

75%

还不是VIP/贵宾

-

威望
0
论坛币
3697 个
通用积分
2.4841
学术水平
19 点
热心指数
21 点
信用等级
17 点
经验
2690 点
帖子
111
精华
0
在线时间
1079 小时
注册时间
2007-7-1
最后登录
2021-12-20

楼主
flashcrazy 发表于 2010-12-31 17:38:52 |AI写论文
5论坛币
初接触bugs,遇到麻烦希望大家帮忙解决。
weight和g[],wd[],tq[]都是设定的数据,需要估出theta[j], bg[i,1], bwd[i,1], btq[i,1]
我是想设bg[i,1]服从正态,bg[i,2]不需要估计,设定其为= 0
但compiling时,提示multiple definitions of node bg[1,1],水平太差,请达人指教!

model{
  for (j in 1:nperson){
     
theta[j] ~ dnorm(0, 1);


  
for (i in 1:nitem){
      
Y[j, i] ~ dbern(prob[j, i]);


ppg[j,i] <- exp(weightg[i]*theta[j]- bg[i,g[i]]);

ppwd[j,i] <- exp(weightwd[i]*theta[j]- bwd[i,wd[i]]);

pptq[j,i] <- exp(weighttq[i]*theta[j]- btq[i,tq[i]]);

num[j,i] <- (pw[i] + pg[i]*ppg[j,i] + pwd[i]*ppwd[j,i] + ptq[i]*pptq[j,i])

den[j,i] <- (1+weightg[i]*ppg[j,i] + weightwd[i]*ppwd[j,i] + weighttq[i]*pptq[j,i]);
      
prob[j, i] <- num[j,i]/den[j,i];



for(k in 1:nt){
   
bg[i,k] ~ dnorm(0, .5);


bwd[i,k] ~ dnorm(0, .5);


btq[i,k] ~ dnorm(0, .5);
   
}

bg[i,2] <- 0;

bwd[i,2] <- 0;

btq[i,2] <-0;

}

}
}

关键词:winbugs WINBUG BUGS bug Win 模型 winbugs

本帖被以下文库推荐

沙发
flashcrazy 发表于 2010-12-31 20:26:12
有什么建议没

藤椅
flashcrazy 发表于 2011-1-1 13:19:36
问题已解决

板凳
liuqi99 发表于 2013-9-1 03:00:03
问题怎么解决的?

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

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2025-12-9 08:15