全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
2860 4
2012-05-01
winbugs的code如下:
model
       {
      for ( i in 1:N )
         {
           r[i] ~ dbin(p[i], n[i])
           b[i] ~ dnorm(0, tau)
           logit(p[i]) <- alpha0 + alpha1*x1[i] + alpha2*x2[i] + alpha12*x1[i]*x2[i] + b[i]
         }
           alpha0 ~ dnorm(0.0, 1.0E-6)
           alpha1 ~ dnorm(0.0, 1.0E-6)
           alpha2 ~ dnorm(0.0, 1.0E-6)
           alpha12 ~ dnorm(0.0, 1.0E-6)
           tau ~ dgamma(0.001, 0.001)
           sigma <- 1/sqrt(tau)
       }
list(
    r=c(10,23,23,26,17,5,53,55,32,46,10,8,10,8,23,0,3,22,15,32,3)
    n=c(39,62,81,51,39,6,74,72,51,79,13,16,30,28,45,4,12,41,30,51,7)
    x1=c(0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1)
    x2=c(0,0,0,0,0,1,1,1,1,1,1,0,0,0,0,0,1,1,1,1,1)
    N=21
     )
list(
     alpha0=0,
     alpha1=0,
     alpha2=0,
     alpha12=0,
     tau=1
     )
在check model和load data的时候,都没问题,但是在complie的时候提示“variable N is not defined”,如何解决啊?谢谢。
二维码

扫码加我 拉你入群

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

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

全部回复
2012-5-1 12:55:21
叫高手来解决吧
二维码

扫码加我 拉你入群

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

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

2012-5-1 13:24:30
list(
     r=c(10,23,23,26,17,5,53,55,32,46,10,8,10,8,23,0,3,22,15,32,3),
     n=c(39,62,81,51,39,6,74,72,51,79,13,16,30,28,45,4,12,41,30,51,7),
     x1=c(0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1),
     x2=c(0,0,0,0,0,1,1,1,1,1,1,0,0,0,0,0,1,1,1,1,1),
     N=21)
二维码

扫码加我 拉你入群

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

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

2012-5-1 14:34:46
epoh 发表于 2012-5-1 13:24
list(
     r=c(10,23,23,26,17,5,53,55,32,46,10,8,10,8,23,0,3,22,15,32,3),
     n=c(39,62,81,51,39, ...
谢谢
二维码

扫码加我 拉你入群

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

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

2018-8-16 11:47:01
楼主怎么解决的?可以交流吗?
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群