@sgaichas @MaC-Perez I had been puzzled at the large nLL for the recruitment penalties.
Turns out I was populating a vector of length Nrecdevs * Nspp with the SUM of the nLL for each recruitment deviation, and then summing those. e.g. effectively upweighting the recdevs by a factor of 400. Oops.
This is due to the version of dnorm() that gets called with a dvar_vector resid and a dvar_vector std.dev. (ie it returns the sum not the individual nLLs)
Fixed in commit 745a725. Sorry about this!
@sgaichas @MaC-Perez I had been puzzled at the large nLL for the recruitment penalties.
Turns out I was populating a vector of length
Nrecdevs * Nsppwith the SUM of the nLL for each recruitment deviation, and then summing those. e.g. effectively upweighting the recdevs by a factor of 400. Oops.This is due to the version of
dnorm()that gets called with advar_vector residand advar_vector std.dev. (ie it returns the sum not the individual nLLs)Fixed in commit 745a725. Sorry about this!
