`panel5pt3.fn` <- function(){ ## Analysis of Microtus data using Model Mt -- a closed population model ## with time-specific probabilities of capture Y<-microtus.data[,1:5] nind<-nrow(Y) lik<-function(parms){ pvec<- expit(parms[1:5]) n0<- exp(parms[6]) N<-nind + n0 loglik<-rep(NA,nind) for(i in 1:nind){ yvec<-Y[i,] loglik[i]<- sum(yvec*log(pvec) + (1-yvec)*log(1-pvec)) } zerocell<- n0*sum(log(1-pvec)) fact<-lgamma(N+1)-lgamma(n0+1) -1*(fact + zerocell + sum(loglik)) } nlm(lik,rep(0,6),hessian=TRUE) }