load inits无法高亮，请帮忙看一下程序，谢谢了？

大家好，我采用附件中的程序，可是无法加载原始数据，load inits无法高亮，不知道哪里出问题了，请帮忙看一下，谢谢了。


model
{
beta[1] <-0
for (tt in 2:nTx){
beta[tt]~dnorm(0,1.0E-6)
}
for(ss in 1:nStudies){
alpha[ss] ~ dnorm(0,1.0E-6)
}
for(ii in 1:LnObs ){
Lmu[ii] < - alpha[Lstudy[ii]]*multi[ii] + beta[Ltx[ii]]-beta[Lbase[ii]]
Lprec[ii] < - 1/pow(Lse[ii],2)
Lmean[ii] ~ dnorm(Lmu[ii],Lprec[ii])
}
for(ss in 1:BnObs){
logCumHaz[ss] < - alpha[Bstudy[ss]] + beta[Btx[ss]] - beta[Bbase[ss]]
cumFail[ss] < - 1-exp(-1*exp(logCumHaz[ss]))
Br[ss] ~ dbin(cumFail[ss], Bn[ss])
}
for (hh in 1:nTx) {
hr[hh] < -exp(beta[hh])
}
# Ranking plot
for (ll in 1:nTx) {
for (mm in 1:nTx) {
rk[ll,mm] < - equals(ranked(beta[],mm),beta[ll])
}
}
}
# Data
# Data set descriptors
list(LnObs = 5, BnObs = 8, nTx = 4, nStudies = 5)
# Log hazard ratio and log hazard data
Lstudy[] Ltx[] Lbase[] Lmean[] Lse[] multi[]
1 1 1 0 0.066 1
1 2 1 0.055 0.063 1
1 3 1 -0.154 0.070 1
1 4 1 -0.209 0.072 1
2 2 1 -0.276 0.203 0
END
# Binary data
Bstudy[] Btx[] Bbase[] Br[] Bn[]
3 3 1 1 229
3 1 1 1 227
4 2 1 4 374
4 3 1 3 372
4 4 1 2 358
4 1 1 7 361
5 3 1 1 554
5 1 1 2 270
END
# Initial values
list(alpha = c(-0.50,-0.50,-0.50,-0.50,-0.50), beta =
c(NA,-0.5,-0.5,-0.5))
list(alpha = c(0.50,0.50,0.50,0.50,0.50), beta = c
(NA,0.5,0.5,0.5))

请大家帮忙看一下吧，谢谢了啊

自己顶自己一下

自己再顶一下

……你这个可以跑的，第一次用不知道操作步骤吧……
要一步步来，先check model，然后选中data以后load data，然后compile，然后才load inits.
可以参照user manual里tutorial的说明进行操作。

生/ty灵火 发表于 2015-5-2 16:19
……你这个可以跑的，第一次用不知道操作步骤吧……
要一步步来，先check model，然后选中data以后load da ...

生/ty灵火 发表于 2015-5-2 16:19
……你这个可以跑的，第一次用不知道操作步骤吧……
要一步步来，先check model，然后选中data以后load da ...

我跑的步骤，请帮忙看一下，哪里出问题了，我就是按照手册来的，谢谢了啊

还是没有操作出来

这两行也要load

数据要加载三次：用红字标出分别加载一次。
# Data
# Data set descriptors
list(LnObs = 5, BnObs = 8, nTx = 4, nStudies = 5)
# Log hazard ratio and log hazard data
Lstudy[] Ltx[] Lbase[] Lmean[] Lse[] multi[]
1 1 1 0 0.066 1
1 2 1 0.055 0.063 1
1 3 1 -0.154 0.070 1
1 4 1 -0.209 0.072 1
2 2 1 -0.276 0.203 0
END
# Binary data
Bstudy[] Btx[] Bbase[] Br[] Bn[]
3 3 1 1 229
3 1 1 1 227
4 2 1 4 374
4 3 1 3 372
4 4 1 2 358
4 1 1 7 361
5 3 1 1 554
5 1 1 2 270
END
# Initial values
list(alpha = c(-0.50,-0.50,-0.50,-0.50,-0.50), beta =
c(NA,-0.5,-0.5,-0.5))
list(alpha = c(0.50,0.50,0.50,0.50,0.50), beta = c
(NA,0.5,0.5,0.5))