求助：我的winbugs程序有问题，哪位大侠帮我看看

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

小女子在作一个研究，具体x~ dmnorm(mu[], T[,]），mu是一个6维的列向量omega~ 逆wish(R[,], 7)，R是一个6维对角矩阵，并且其对角矩阵上的6个元素都服从dgamma(0.001, 0.001)，现在我想估计出omega及R，我编的程序如下：
model
{
  for(i in 1:N){x[i,1:6] ~ dmnorm(mu[], T[,])}
T[1:6, 1:6]~ dwish(R[,], 7)
R[1,1]~dgamma(0.001, 0.001)
R[2,2]~dgamma(0.001, 0.001)
R[3,3]~dgamma(0.001, 0.001)
R[4,4]~dgamma(0.001, 0.001)
R[5,5]~dgamma(0.001, 0.001)
R[6,6]~dgamma(0.001, 0.001)
omega[1:6,1:6]<-inverse(T[,])
}
list(N=15,  mu = c(0, 0, 0, 0, 0, 0), x=structure(
              .Data=c(-1.1108, -0.84162, 0.84162, -0.43073, 0.43073, 2.3263,
                      0.62293,  -0.43073, -0.083652, -0.25335, -1.1108, 0.84162,
                      -0.43073, 0.083652, -0.62293, 0.84162, -1.5011, 0.84162,
                      0.25335, 2.3263, 3.0902, -0.43073, -0.25335, 0.84162,
                      0.083652, 0.25335, -0.62293, -0.43073, -0.84162, 1.5011,
                      -0.84162, -0.083652, 1.1108, -0.43073, 0.62293, 0.84162,
                       0.43073, -0.62293, 0.083652, 0.43073, -0.43073, 0.84162,
                       0.84162, 1.5011,  0.43073, 1.5011, -0.62293, 0.84162,
                       -0.083652, 0.43073, -0.43073, 0.083652, 0.083652, 0.84162,
                       -0.62293, 1.1108, 0.25335, 1.1108, 2.3263, 0.84162,
                        2.3263, 0.84162, 1.5011, 3.0902, 1.5011, 0.84162,
                        -1.5011, 0.62293, -0.25335, -0.43073, 0.25335, 0.84162,
                         1.1108, -1.5011, -0.62293, -0.083652, 1.1108, 0.84162,
                         1.5011, -0.25335, -0.62293, 0.62293, -0.083652, 1.1108,
                         -0.25335, -1.1108, 0.62293, 0.25335, 0.84162, 0.84162),
                           .Dim=c(15,6)),            
             R(1,2:6)=c(0,0,0,0,0), R(2,1)=0, R(2,3:6)=c(0,0,0,0), R(3,1:2)=c(0,0), R(3,4:6)=c(0,0,0),
             R(4,1:3)=c(0,0,0), R(4,5:6)=(0,0), R(5,1:4)=C(0,0,0,0), R(5,6)=0, R(6,1:5)=C(0,0,0,0,0))

list(R[1,1] = 0.01, R[2,2] = 0.01, R[3,3] = 0.01, R[4,4] = 0.01, R[5,5] = 0.01, R[6,6] = 0.01,T = structure(
.Data = c(0.1, 0.3, 0.1, 0.2, 0.1, 0.1,

0.1, 0.1, 0.7, 0.1, 0.1, 0.3,

0.1, 0.3, 0.1, 0.1, 0.2, 0.1,

0.4, 0.1, 0.1, 0.5, 0.1, 0.1,

0.1, 0.6, 0.1, 0.1, 0.1, 0.7,

0.3, 0.1, 0.1, 0.1, 0.2, 0.1),

.Dim=c(6, 6)))

黑色字体没有问题，但在load初始值时，显示：expected an equals sign,我也不知道是哪里错了。请高人帮我指正。

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：winbugs WINBUG BUGS Win bug 求助 程序 大侠 winbugs

太深奥了，不懂，

N为什么是15？