求助wingbugs程序cheked model 和load data都可以，但是compile却提示N is not define

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

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.000001)
    alpha1~dnorm(0,0.000001)
    alpha2~dnorm(0,0.000001)
    alpha12~dnorm(0,0.000001)
    tau~dgamma(0.001,0.001)
    sigma <- 1/sqrt(tau)
    }
Data
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
)
Inits
list(
alpha0=0,
alpha1=0,
alpha2=0,
alpha12=0,
tau=1
)

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：wingbugs COMPILE DEFINE model Data 程序