`panel5pt6.fn` <- function(ni=5000,nb=2000,nthin=2,nc=3,nz=125){ # This R script fits a behavioral response model to the Microtus # data in WinBUGS, using data augmentation (Panel 5.6). The model is # a hierarchical version of the Yang and Chao (2005; Biometrics) model. # See Royle (2008 -- Analysis of capture-recapture models with # individual covariates using data augmentation, Biometrics). library("R2WinBUGS") Y<-microtus.data[,1:5] nind<-nrow(Y) J<-ncol(Y) wt<-microtus.data[,6] first<-rep(NA,nind) ## this block of code creates a dummy variable for previous capture prevcap<-matrix(0,nrow=nind,ncol=J) for(i in 1:nind){ first[i]<- min((1:J)[Y[i,]==1]) if(first[i]