#Primero creamos una matriz de proyeccion y un vector de clases de estados con sus abundancias#
#Esta es la matriz#
A <- matrix(c(0, 30, 60, 10, 0.25, 0, 0, 0, 0, 0.5, 0, 0, 0, 0, 0.35, 0), nr = 4, byrow = TRUE)
A
## [,1] [,2] [,3] [,4]
## [1,] 0.00 30.0 60.00 10
## [2,] 0.25 0.0 0.00 0
## [3,] 0.00 0.5 0.00 0
## [4,] 0.00 0.0 0.35 0
#Este es el vector con las abundancias#
Nt <- matrix(c(500, 100, 300, 50), ncol = 1)
Nt
## [,1]
## [1,] 500
## [2,] 100
## [3,] 300
## [4,] 50
#Ahora realizamos la multiplicación de matrices entre la matriz de proyección y Nt para saber el número de individuos en el siguiente año para cada estado#
Nt1= A %*% Nt
Nt1
## [,1]
## [1,] 21500
## [2,] 125
## [3,] 50
## [4,] 105
Nt2=A%*%Nt1
Nt2
## [,1]
## [1,] 7800.0
## [2,] 5375.0
## [3,] 62.5
## [4,] 17.5