`panel5pt2.fn` <-
function(){
## Analysis of the Microtus data using Model M0 -- a simple closed
## population model with constant probability of detection

Y<-microtus.data[,1:5]
nind<-nrow(Y)
nvec<-table(apply(Y,1,sum))

lik0<-function(parms){
  p<-  expit(parms[1])
  n0<- exp(parms[2])
  N<-nind + n0
  cp<-dbinom(0:5,5,p)
  -1*(lgamma(N+1) - lgamma(n0+1) + sum(c(n0,nvec)*log(cp) ))
}

 nlm(lik0,c(0,1),hessian=TRUE)

}