全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 MATLAB等数学软件专版
1548 3
2010-12-31
悬赏 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;

}

}
}
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2010-12-31 20:26:12
有什么建议没
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2011-1-1 13:19:36
问题已解决
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-9-1 03:00:03
问题怎么解决的?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群