全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
3614 3
2011-06-15
请教一个问题:
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这样定义。如果要实现上述代码的功能,该怎么办?
二维码

扫码加我 拉你入群

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

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

全部回复
2011-6-18 01:35:48
前面一个mu[i,j]的定义表明这是一个确定的logical的variable,如果你所有的数据都是给定的;
而后一个mu[i,j]的定义又说这是一个Random variable, 是stochastic的,那这当然是矛盾的。

你需要确定哪个是变量,哪个是数据,重新定义model
二维码

扫码加我 拉你入群

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

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

2011-6-20 08:31:23
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]这一部分都有给定的数据。
二维码

扫码加我 拉你入群

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

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

2011-6-26 11:23:57
你这样的逻辑是不对的,即使等号右侧已经是固定的值了,等号左边为什么会是Random的?除非你定义两个值mu1和mu2,然后建立二者之间的关系
二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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