winbugs 中的multiple definition of node问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

请教一个问题：
winbugs中
for (i in 1:m){
for (j in 1:k){
modis[i,j]~dnorm(mu[i,j],tau1)
amsre[i,j]~dnorm(mu[i,j],tau2)
mu[i,j]<-sst.mu[i,j]+a*mean1[i,j]+b*mean2[i,j]+enta[i,j]
mu[i,j]~dnorm(27.31,5.08)
}
}
mean1和mean2是两种尺度的观测数据，mu[i,j]是两种尺度观测数据的关联方程
又已知mu[i,j]~dnorm(27.31,5.08)
这样写代码总是提示multiple definition of mu，winbugs中只允许constant node这样定义。如果要实现上述代码的功能，该怎么办？

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：definition Multiple winbugs multip WINBUG winbugs

前面一个mu[i,j]的定义表明这是一个确定的logical的variable，如果你所有的数据都是给定的；
而后一个mu[i,j]的定义又说这是一个Random variable, 是stochastic的，那这当然是矛盾的。

你需要确定哪个是变量，哪个是数据，重新定义model

mu[i,j]<-sst.mu[i,j]+a*mean1[i,j]+b*mean2[i,j]+enta[i,j]

如果通过上述公式计算出来的mu，遵循mu[i,j]~dnorm(27.31,5.08)，在winbugs中怎么实现？sst.mu[i,j]+a*mean1[i,j]+b*mean2[i,j]+enta[i,j]这一部分都有给定的数据。

你这样的逻辑是不对的，即使等号右侧已经是固定的值了，等号左边为什么会是Random的？除非你定义两个值mu1和mu2，然后建立二者之间的关系