全部版块 我的主页
› 论坛 › 数据科学与人工智能 › 数据分析与数据科学 › R语言论坛
winBUGS程序使用数据输入问题
KarlShanks
2035 6
2018-04-03
各位大佬，我在复现一篇paper 的时候，不知道怎么输入合适的数据，才用winbugs,不熟悉，求指教
代码是paper给的原代码，没有问题，但是数据有点奇怪，按照paper需要的数据，我读不进去。点load data一直没响应
下面是moedel 代码：
model{
       for (g in 1:ngames) {
           y1[g] ~ dpois(theta[g,1])
           y2[g] ~ dpois(theta[g,2])
          ynew[g,1] ~ dpois(theta[g,1])
          ynew[g,2] ~ dpois(theta[g,2])
          log(theta[g,1]) <- home+att[hometeam[g]]+def[awayteam[g]]
          log(theta[g,2]) <- att[awayteam[g]]+def[hometeam[g]]
          }
    }

home ~  dnorm(0,0.0001)

for(t in 1:nteams){
   att.star[t] ~ dnorm(mu.att,tau.att)
   def.star[t] ~dnorm(mu.def,tau.def)
   att[t] <- att.star[t]-mean(att.star[])
   def[t] <- def.star[t]-mean(def.star[])
}

mu.att ~ dnorm(0,0.0001)
mu.def ~ dnorm(0,0.0001)
tau.att ~ dgamma(0.01,0.01)
tau.def ~ dgamma(0.01,0.01)

for( t in 1:nteams){
   grp.att[t] ~ dcat(p.att[t,])
   grp.def[t] ~ dcat(p.def[t,])
   att[t] ~ dt(mu.att[grp.att[t]],tau.att[grp.att[t]],4)
   def[t] ~ dt(mu.def[grp.def[t]],tau.def[grp.att[t]],4)
   att.star[t]<- att[t]-mean(att[])
   def.star[t]<-def[t]-mean(def[])
   p.att[t,1:3] ~ ddirch(prior.att[t,])
   p.def[t,1:3] ~ ddirch(prior.def[t,])
}

mu.att[1] ~ dj1.dnorm.trunc(0,0.001,-3,0)
mu.def[1] ~ dj1.dnorm.trunc(0,0.0001,0,3)
tau.att[1] ~ dgamma(0.01,0.01)
tau.def[1] ~ dgamma(0.01,0.01)
mu.att[2] <-0
mu.def[2]<-0
tau.att[2] ~ dgamma(0.01,0.01)
tau.def[2] ~ dgamma(0.01,0.01)
mu.att[3] ~ dj1.dnorm.trunc(0,0.001,0,3)
mu.def[3] ~ dj1.dnorm.trunc(0,0.0001,-3,0)
tau.att[3] ~ dgamma(0.01,0.01)
tau.def[3] ~ dgamma(0.01,0.01)
}

y1y2y3y4

1

1

1

1

0

0

2

11

1

1

3

14

0

0

4

10



不知道怎么输入进去，求助

二维码

扫码加我 拉你入群

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

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

全部回复
13644444
2018-4-3 01:11:34
哇哇哇哇哇哇哇哇哇哇哇哇哇哇哇哇哇哇哇
二维码

扫码加我 拉你入群

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

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

KarlShanks
2018-4-3 01:14:17
需要输入的数据有四个，y1 y2是主队客队进球数据，y3 y4 是球队编号。一共380场比赛记录
二维码

扫码加我 拉你入群

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

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

KarlShanks
2018-4-3 01:19:32
list(ngames=380,
y1=c(1,0,1,0,0,1,0,1,4,0,3,2,0,0,0,2,1,2,3,0,0,3,3,0,2,8,2,1,0,1,0,2,3,1,1,0,1,0,0,1,4,1,3,3,0,2,3,4,4,3,4,2,2,3,1,0,1,1,2,2,0,0,2,2,1,1,6,2,2,3,1,4,3,1,2,1,2,2,2,1,0,1,1,2,3,1,1,1,1,3,0,3,3,2,3,0,0,3,1,3,3,0,0,1,2,2,5,1,3,2,3,1,2,0,0,2,3,1,2,1,0,3,1,2,3,1,1,1,0,0,4,5,0,1,3,0,0,2,1,1,1,1,4,4,2,1,2,2,2,1,1,2,0,1,1,0,1,0,1,1,2,1,1,0,3,2,2,0,1,0,1,1,5,0,1,0,1,0,1,3,1,1,1,2,1,2,1,1,0,0,1,2,1,1,3,1,1,1,3,2,1,3,1,0,3,1,1,1,3,0,3,1,5,0,3,1,2,2,1,3,1,2,1,3,0,2,0,0,1,3,7,3,2,1,0,1,1,3,2,0,3,1,2,2,2,1,2,5,0,1,3,3,2,0,4,0,5,1,2,1,1,2,1,1,1,0,5,1,1,1,2,1,1,1,0,2,1,1,2,3,0,1,1,0,2,0,2,3,1,3,2,0,1,2,1,3,0,1,1,2,2,2,3,2,2,2,2,3,0,0,0,2,1,2,2,0,3,1,2,1,4,1,2,1,2,4,3,1,0,1,0,3,1,4,1,0,0,2,1,2,3,1,0,4,0,1,4,0,1,2,4,0,4,6,2,1,0,1,1,0,3,1,2,2,2,0,1,0,0,4,2,3,3,2,2,0,1,2,2,3),
y2=c(2,0,1,0,4,1,0,2,0,2,1,1,1,1,0,3,1,0,0,0,1,1,1,0,0,2,1,5,1,0,5,2,0,0,2,2,1,1,0,1,3,2,1,0,3,2,1,0,0,0,1,1,0,1,1,0,2,1,1,0,4,2,0,2,2,5,0,0,1,1,1,1,1,1,0,3,1,2,0,2,2,1,0,2,1,2,3,6,0,5,1,1,3,2,1,2,2,1,3,0,2,1,0,0,1,3,0,3,1,1,1,2,3,0,1,1,3,2,0,1,2,0,1,1,1,2,3,1,0,1,2,1,3,1,0,4,1,1,0,0,2,0,1,2,0,2,1,1,1,2,1,0,0,1,2,2,0,2,0,2,2,2,0,5,0,3,3,0,1,2,1,1,0,0,1,0,1,2,1,1,0,1,3,3,1,2,1,0,1,2,2,1,2,2,0,0,4,2,0,0,1,1,0,0,0,1,2,0,2,1,1,1,2,0,1,2,0,3,2,2,0,0,1,1,3,2,2,0,1,0,1,0,0,2,1,2,1,3,1,0,2,2,0,1,1,2,3,0,1,5,0,0,2,1,0,0,2,2,0,1,2,0,1,0,0,2,0,1,3,0,1,0,0,0,2,0,1,0,1,0,3,1,0,5,0,1,1,2,1,0,1,0,2,1,1,1,2,3,0,4,0,1,3,1,0,3,0,1,2,2,3,1,1,2,1,0,0,0,0,1,0,1,0,2,3,0,0,0,3,6,0,0,0,0,2,0,0,0,1,1,0,0,1,4,2,2,0,3,1,2,4,0,0,1,0,0,1,1,4,2,3,1,2,1,0,2,0,0,1,1,1,1,2,0,2,1,0,0,3,2),
hometeam=c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,11,1,12,3,15,6,19,4,20,8,10,16,13,9,7,14,18,2,17,5,11,1,16,13,15,18,2,19,14,20,10,12,3,9,4,7,8,6,5,17,11,1,13,19,14,18,16,2,15,20,12,3,9,17,5,7,6,10,4,8,11,16,3,4,18,10,1,2,19,5,12,13,9,17,14,15,8,6,20,7,10,11,1,3,19,4,5,16,2,18,13,9,17,7,14,15,8,6,12,20,10,16,12,19,17,7,14,8,3,15,11,1,9,4,5,20,6,13,18,2,10,16,3,19,17,15,8,7,14,12,1,13,2,4,6,18,11,9,5,20,1,18,11,13,2,9,4,5,6,20,16,12,3,19,7,14,8,10,17,15,3,10,16,12,19,17,7,15,14,8,11,1,2,5,18,9,20,6,13,4,20,11,1,12,3,19,20,8,4,15,6,16,13,2,17,5,7,14,18,10,9,13,19,15,20,18,11,1,16,2,14,10,9,17,5,7,8,6,12,4,3,1,16,13,2,19,14,15,20,11,18,12,9,4,5,8,6,10,17,7,1,3,9,5,7,8,6,2,4,20,11,16,12,13,14,18,19,17,15,10,3,2,6,4,18,1,13,9,5,20,10,16,12,3,17,7,14,15,8,19,11,13,2,9,5,6,18,4,20,1,15,16,12,3,17,7,14,8,10,19,11,13,2,4,20,1,9,5,6,18,17,14,15,8,19,10,10,11,1,16,2,4,5,3,19,18,11,13,17,7,14,15,8,6,12,20,9,3,7,16,12,10,11,16,2,19,4,5,18,1,3,12,13,9,17,7,14,15,20,8,6),
awayteam=c(18,11,14,10,16,17,12,19,15,3,1,8,5,4,6,9,7,2,20,18,13,17,16,14,5,10,2,9,7,15,19,11,6,3,12,20,1,8,4,4,10,17,6,8,5,9,12,7,3,16,15,18,13,1,19,2,20,11,14,6,9,3,17,8,4,12,5,7,10,13,19,11,15,1,2,16,14,20,18,8,14,17,6,15,7,20,13,9,12,10,19,18,1,11,16,3,2,5,4,8,17,12,15,14,13,9,7,20,6,18,4,10,5,2,19,16,1,3,11,2,13,18,4,5,1,6,20,9,11,19,15,17,12,8,16,10,7,14,3,13,11,5,18,4,2,6,20,1,9,8,17,16,15,12,7,3,10,19,14,16,17,10,15,19,7,8,14,3,12,4,2,1,6,11,13,9,18,20,5,4,5,18,11,1,2,6,20,9,13,15,7,10,17,12,3,8,14,16,19,13,13,2,14,7,16,18,17,5,10,9,3,1,4,12,6,8,15,11,19,20,9,7,12,6,3,5,4,10,8,17,1,2,16,18,14,15,13,19,11,20,5,6,12,7,3,10,17,4,9,8,16,1,18,2,14,11,20,19,15,11,10,16,13,17,12,15,18,14,19,2,5,7,20,3,1,8,6,9,4,13,15,8,17,19,14,10,12,3,7,11,1,20,6,18,9,5,13,4,2,12,8,17,14,10,7,16,3,15,19,4,2,6,11,13,18,20,1,9,5,7,14,12,16,17,3,8,15,19,10,9,18,1,5,11,6,12,14,17,15,6,7,20,8,13,9,16,2,3,10,16,18,11,4,5,1,19,2,13,20,4,17,20,8,14,15,9,7,13,6,12,1,4,5,11,16,19,3,2,10,18)
)

输进去没反应。
我不知道ynew[ ]这个该不该输入进去，因为没有ynew只是在做预测，不需要训练
二维码

扫码加我 拉你入群

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

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

KarlShanks
2018-4-3 01:20:29
(ngames=380,
y1=c(1,0,1,0,0,1,0,1,4,0,3,2,0,0,0,2,1,2,3,0,0,3,3,0,2,8,2,1,0,1,0,2,3,1,1,0,1,0,0,1,4,1,3,3,0,2,3,4,4,3,4,2,2,3,1,0,1,1,2,2,0,0,2,2,1,1,6,2,2,3,1,4,3,1,2,1,2,2,2,1,0,1,1,2,3,1,1,1,1,3,0,3,3,2,3,0,0,3,1,3,3,0,0,1,2,2,5,1,3,2,3,1,2,0,0,2,3,1,2,1,0,3,1,2,3,1,1,1,0,0,4,5,0,1,3,0,0,2,1,1,1,1,4,4,2,1,2,2,2,1,1,2,0,1,1,0,1,0,1,1,2,1,1,0,3,2,2,0,1,0,1,1,5,0,1,0,1,0,1,3,1,1,1,2,1,2,1,1,0,0,1,2,1,1,3,1,1,1,3,2,1,3,1,0,3,1,1,1,3,0,3,1,5,0,3,1,2,2,1,3,1,2,1,3,0,2,0,0,1,3,7,3,2,1,0,1,1,3,2,0,3,1,2,2,2,1,2,5,0,1,3,3,2,0,4,0,5,1,2,1,1,2,1,1,1,0,5,1,1,1,2,1,1,1,0,2,1,1,2,3,0,1,1,0,2,0,2,3,1,3,2,0,1,2,1,3,0,1,1,2,2,2,3,2,2,2,2,3,0,0,0,2,1,2,2,0,3,1,2,1,4,1,2,1,2,4,3,1,0,1,0,3,1,4,1,0,0,2,1,2,3,1,0,4,0,1,4,0,1,2,4,0,4,6,2,1,0,1,1,0,3,1,2,2,2,0,1,0,0,4,2,3,3,2,2,0,1,2,2,3),
y2=c(2,0,1,0,4,1,0,2,0,2,1,1,1,1,0,3,1,0,0,0,1,1,1,0,0,2,1,5,1,0,5,2,0,0,2,2,1,1,0,1,3,2,1,0,3,2,1,0,0,0,1,1,0,1,1,0,2,1,1,0,4,2,0,2,2,5,0,0,1,1,1,1,1,1,0,3,1,2,0,2,2,1,0,2,1,2,3,6,0,5,1,1,3,2,1,2,2,1,3,0,2,1,0,0,1,3,0,3,1,1,1,2,3,0,1,1,3,2,0,1,2,0,1,1,1,2,3,1,0,1,2,1,3,1,0,4,1,1,0,0,2,0,1,2,0,2,1,1,1,2,1,0,0,1,2,2,0,2,0,2,2,2,0,5,0,3,3,0,1,2,1,1,0,0,1,0,1,2,1,1,0,1,3,3,1,2,1,0,1,2,2,1,2,2,0,0,4,2,0,0,1,1,0,0,0,1,2,0,2,1,1,1,2,0,1,2,0,3,2,2,0,0,1,1,3,2,2,0,1,0,1,0,0,2,1,2,1,3,1,0,2,2,0,1,1,2,3,0,1,5,0,0,2,1,0,0,2,2,0,1,2,0,1,0,0,2,0,1,3,0,1,0,0,0,2,0,1,0,1,0,3,1,0,5,0,1,1,2,1,0,1,0,2,1,1,1,2,3,0,4,0,1,3,1,0,3,0,1,2,2,3,1,1,2,1,0,0,0,0,1,0,1,0,2,3,0,0,0,3,6,0,0,0,0,2,0,0,0,1,1,0,0,1,4,2,2,0,3,1,2,4,0,0,1,0,0,1,1,4,2,3,1,2,1,0,2,0,0,1,1,1,1,2,0,2,1,0,0,3,2),
hometeam=c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,11,1,12,3,15,6,19,4,20,8,10,16,13,9,7,14,18,2,17,5,11,1,16,13,15,18,2,19,14,20,10,12,3,9,4,7,8,6,5,17,11,1,13,19,14,18,16,2,15,20,12,3,9,17,5,7,6,10,4,8,11,16,3,4,18,10,1,2,19,5,12,13,9,17,14,15,8,6,20,7,10,11,1,3,19,4,5,16,2,18,13,9,17,7,14,15,8,6,12,20,10,16,12,19,17,7,14,8,3,15,11,1,9,4,5,20,6,13,18,2,10,16,3,19,17,15,8,7,14,12,1,13,2,4,6,18,11,9,5,20,1,18,11,13,2,9,4,5,6,20,16,12,3,19,7,14,8,10,17,15,3,10,16,12,19,17,7,15,14,8,11,1,2,5,18,9,20,6,13,4,20,11,1,12,3,19,20,8,4,15,6,16,13,2,17,5,7,14,18,10,9,13,19,15,20,18,11,1,16,2,14,10,9,17,5,7,8,6,12,4,3,1,16,13,2,19,14,15,20,11,18,12,9,4,5,8,6,10,17,7,1,3,9,5,7,8,6,2,4,20,11,16,12,13,14,18,19,17,15,10,3,2,6,4,18,1,13,9,5,20,10,16,12,3,17,7,14,15,8,19,11,13,2,9,5,6,18,4,20,1,15,16,12,3,17,7,14,8,10,19,11,13,2,4,20,1,9,5,6,18,17,14,15,8,19,10,10,11,1,16,2,4,5,3,19,18,11,13,17,7,14,15,8,6,12,20,9,3,7,16,12,10,11,16,2,19,4,5,18,1,3,12,13,9,17,7,14,15,20,8,6),
awayteam=c(18,11,14,10,16,17,12,19,15,3,1,8,5,4,6,9,7,2,20,18,13,17,16,14,5,10,2,9,7,15,19,11,6,3,12,20,1,8,4,4,10,17,6,8,5,9,12,7,3,16,15,18,13,1,19,2,20,11,14,6,9,3,17,8,4,12,5,7,10,13,19,11,15,1,2,16,14,20,18,8,14,17,6,15,7,20,13,9,12,10,19,18,1,11,16,3,2,5,4,8,17,12,15,14,13,9,7,20,6,18,4,10,5,2,19,16,1,3,11,2,13,18,4,5,1,6,20,9,11,19,15,17,12,8,16,10,7,14,3,13,11,5,18,4,2,6,20,1,9,8,17,16,15,12,7,3,10,19,14,16,17,10,15,19,7,8,14,3,12,4,2,1,6,11,13,9,18,20,5,4,5,18,11,1,2,6,20,9,13,15,7,10,17,12,3,8,14,16,19,13,13,2,14,7,16,18,17,5,10,9,3,1,4,12,6,8,15,11,19,20,9,7,12,6,3,5,4,10,8,17,1,2,16,18,14,15,13,19,11,20,5,6,12,7,3,10,17,4,9,8,16,1,18,2,14,11,20,19,15,11,10,16,13,17,12,15,18,14,19,2,5,7,20,3,1,8,6,9,4,13,15,8,17,19,14,10,12,3,7,11,1,20,6,18,9,5,13,4,2,12,8,17,14,10,7,16,3,15,19,4,2,6,11,13,18,20,1,9,5,7,14,12,16,17,3,8,15,19,10,9,18,1,5,11,6,12,14,17,15,6,7,20,8,13,9,16,2,3,10,16,18,11,4,5,1,19,2,13,20,4,17,20,8,14,15,9,7,13,6,12,1,4,5,11,16,19,3,2,10,18)
)
ynew 这个是预测，我不知道该怎么输入进去
二维码

扫码加我 拉你入群

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

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

KarlShanks
2018-4-3 01:22:02
KarlShanks 发表于 2018-4-3 01:14
需要输入的数据有四个，y1 y2是主队客队进球数据，y3 y4 是球队编号。一共380场比赛记录
ynew那个数据是预测，我不知道需不需要输入，因为不是文章模型里面提过的数据
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

微信QQ空间QQ微博
扫码加好友，拉您进群
各岗位、行业、专业交流群