底辺私立医大を卒業した医者って頭悪いよね？ Part16

N=950
y=c(rep(NA,N-1),30) # c(NA,NA,....,30)
thresh=cbind(rep(15,N),rep(25,N),rep(30,N)) # 15,25,30
ybin=c(rep(1,N-1),3) # 0:[,15), 1:[15,25), 2:[25,30), 3:[30,]
data=list(N=N,y=y,ybin=ybin,thresh=thresh)
inits=list(theta=0.5)
model.string='
model{
for (i in 1:N){
ybin[i] ~ dinterval( y[i] , thresh[i, ] )
y[i] ~ dbin(theta,50)
}
theta ~ dbeta(1,1)I(0.25,1)
}
'