Internal estimation function for LSRA and LSRA2 functions
LSRA_opt(param, FF_a, Chist, M_a, Mat_age_a, Wt_age_a, sel_a, Recdevs_a, h_a,
Umax = 0.5, mode = 1)a numeric value representing log(R0)
numeric value, recent fishign mortality rate (apical F)
a vector of historical catch observations [nyears]
numeric value, natural mortality rate
a vector of maturity at age [nage]
a vector of weight at age [nage]
a vector of selectivity at age [nage]
a vector of recruitment deviations [nyears]
a numeric value of steepness values of the Bev-Holt Stock-Recruitment relationship
maximum harvest rate per year
1-5 see below
depends on mode but could be 1:objective function 2:trajectory of Fs 3: SSB depletion 4:log(R0) 5:diagnostic plots