请winbugs高手帮忙看下我的数据录入错在哪里

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

各位好，在数据处理过程中需要用到winbugs，但是我本人对这个软件不是特别熟悉。请各位高手帮忙看下我的数据输入错在哪里了？不胜感激！数据如下：
model{
for(i in 1:N){
for(j in 1:T[i]){
mean_diff[i,j] ~dnorm(theta[i,j],prec[i,j])
prec[i,j]<-1/(sd_diff[i,j]* sd_diff[i,j])
theta[i,j]<-beta[i]*time_period[i,j]
}
beta[i] ~ dnorm(beta.p, tau)
}
beta.p ~ dnorm(0, 0.00001)
tau <-1/var
var <-pow(sd,2)
sd ~ dnorm(0, 0.01)I(0,)
pr.gr.zero <-step(beta.p)- equals(0, beta.p)
}
list(N = c(1, 2, 3, 3, 3 ,3 ,3 ,3, 3, 3, 3, 3, 3, 3, 3, 4, 5, 6),
T = c(1,1, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 1, 1, 1),
mean_diff = c(-0.50, -0.10, -0.15, -0.40, -0.28, -0.18, -0.12, -0.19, -0.06, 0.32, 0.27, 0.25, 0.15, -0.12, 0.18, 0.31, 1.02, -0.01),
sd_diff = c(2.59, 2.12, 1.97, 2.19, 2.41, 2.33, 2.22, 2.31, 2.29, 2.29, 2.21, 2.37, 2.36, 2.26, 2.36, 4.76, 5.30, 4.35),
time_period = c(4, 3, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 30, 24, 54))

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：winbugs WINBUG 数据录入 BUGS bug 录入

N错了吧，N=数据个数

601361436 发表于 2013-5-10 11:43
N错了吧，N=数据个数

谢谢你的回复，不过我把它改成N = 18后，依然是错的，错误的原因是expected key word structure.

是不是少了一个括号