【50论坛币急求解答】请教winbugs问题

想用openbugs估计自己图中分布的参数吗，就是logistic回归基础上增加了一个lambda，请问我的model写的对吗？
model
   {
    for (i in 1:N) {
    y~dbin(p,1)
   p<-1/pow((exp(-beta0-beta1*x)+1),lambda)
       }
beta1~dnorm(0,1)
beta0~dnorm(0,1)
lambda~dgamma(2,1)
}

list(y=c(1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0),
x=c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23),
N=23)

#Inits

list(beta1=1,beta0=0,lambda=1)

顶一下

See this revision. It works on my computer.


model
   {
    for (i in 1:N) {
    y[i]~dbin(p[i],1)
   p[i]<-1/pow((exp(-beta0-beta1*x[i])+1),lambda)
       }
beta1~dnorm(0,1)
beta0~dnorm(0,1)
lambda~dgamma(2,1)
}

list(y=c(1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0),
x=c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23),
N=23)

#Inits

list(beta1=1,beta0=0,lambda=1)