Displaying 1 result from an estimated 1 matches for "survship".
2007 Mar 02
0
Matrix looping
...is my problem.
I have three ages: Nage<-c(1,2,3)
I have an weight matrix: Wt<-c( 0.04952867, 0.23808432, 0.34263880)
I have an age schedule of maturity: Mat<-c(0,1,1) where 0 is not mature,
and 1 is mature
I have a vulnerability schedule: Vul<-c(0,1,1)
I have an survivorship schedule: Survship<-c(1,0.4,0.16)
I also have leading parameters R0<-130.66; recK<-3.068; a<-5.48;
b<-0.0282; S<-0.4
I have annual catches for 100 years, ct<-runif(100,5,20)
Now I want a matrix of 100 years x 3 ages
yr<-c(1:100)
Nt<-matrix(0,nrow=length(yr)+1),ncol=length(Nage))
Now the f...