R/project.leslie.R

In demogR: Analysis of Age-Structured Demographic Models

`project.leslie` <-
function(A,no,tmax,pop.sum=FALSE){
  if(length(no) != dim(A)[1])
    stop("Projection matrix and population vector have different number of states!")
  N <- matrix(0, nrow=length(no),ncol=tmax+1)
  N[,1] <- no
  pop <- no
  for(t in 1:tmax){
    pop <- A%*%pop
    N[,t+1] <- pop
  }
  if(pop.sum){
    N <- apply(N,2,sum)
  }
  return(N)
}