`panel7pt1.fn` <-
function(ni=310,nb=10,nthin=1,nc=2,nz=500){
library("R2WinBUGS")

x<- sin( (impala[,5]/360)*(2*pi) ) *impala[,4]
x<-x/100

#x<- simdist.fn()
nind<-length(x)
y<-c(rep(1,nind),rep(0,nz))
x<-c(x,rep(NA,nz))

sink("model.txt")
cat("
model {
sigma~dunif(0,10)
sigma2<-sigma*sigma
psi~dunif(0,1)

for(i in 1:(nind+nz)){
 w[i]~dbern(psi)
 x[i]~dunif(0,4)
 logp[i]<-   -((x[i]*x[i])/sigma2)
  p[i]<-exp(logp[i])
 mu[i]<-w[i]*p[i] 
 y[i]~dbern(mu[i])
 
}
N<-sum(w[1:(nind+nz)])
D<- N/48
}
",fill=TRUE)
sink()

wst<-y
data <- list ("nind","x","nz","y")
inits <- function(){
  list (psi=runif(1),w=wst,sigma=runif(1,.4,2))

}
parameters <- c("N","sigma","D")


out <- bugs (data, inits, parameters, "model.txt", n.thin=nthin,n.chains=nc, n.burnin=nb,n.iter=ni,debug=FALSE,DIC=FALSE)







out


}