m**c
发帖数: 88 | 1 model{
for (t in 1:n){
y[t,1:5]~ dmnorm( mu[], T[,] )
}
index ~ dcat(Prob[])
for (i in 1:H) {
x_spike[i] ~ dnorm(1,1.0E10)
x_slap[i] ~ dnorm(0.98,100)
lambda[i] <- Index[index,i]
x[i,1] <- (1- lambda[i])* x_spike[i] + lambda[i]* x_slap[i]
}
for (j in 1:P){
mu[j] <- inprod(A[j,],x[,])+ A0[j,1]
}
T[1:P,1:P] <- inverse(V[,])
}
DAtA
.................
.................
.................
END
这个WinBUGS程序我可以在Sample Monnitor Tool 的node框里选择保存 mu，x，x_
spike，x_slap，index等节点的采样样本，当时无法保存y的采样结果？怎么回事啊？
难道y不是节点（node）？y～N（mu，T）。mu可以保存，y怎么不行啊？哪位给看看啊? |
|
