Closed ejardim closed 8 years ago
Fixed. It uses the matching yearclass with, for example in bevholt
if (Rmodel == 1) { // beverton holt
for(int y=minYear+minAge; y<=maxYear; ++y){
predLogR = ra(y) + log(ssb(y-minAge)) - log(mfexp(rb(y)) + ssb(y-minAge));
varLogR = log(pow(srCV,2)+1);
nll += nldnorm(r(y), predLogR, varLogR);
}
}
Why is the fit being done on the previous year ssb ?