R/estN.LSC.R

In CARE1: Statistical package for population size estimation in capture-recapture models.

estN.LSC <-
function(z){
t <- as.integer(log(length(z) + 1) / log(2))
Mat <- matrix(0, ncol = t, nrow = 2 ^ t - 1)
for(i in 1:(2 ^ t - 1)) Mat[i,] <- rev(as.integer(intToBits(i))[1:t])
M <- sum(z)
D <- M - mean(sapply(1:t, function(j) Sub.S(z, t, Mat, j)))
Chat <- 1 - mean(sapply(1:t, function(j) Sub.S(z, t, Mat, j)) / sapply(1:t, function(j) Sub.n(z, t, Mat, j)))

MatA <- matrix(0, t, t)
MatB <- matrix(0, t, t)
MatC <- matrix(0, t, t)

for(i in 1:(t - 1)){
for(j in (i + 1):t){
MatA[i, j] = Sub.A(z, t, Mat, i, j) * Sub.r(z, t, Mat, i, j)
}
}

Nhat <- D / Chat + 1 / (t * Chat) * sum(MatA)
names(Nhat) <- "Nhat"
return(Nhat)
}