Open zandjyo opened 2 years ago
Need to develop a script to output tuning values for annually varying selectivity parameters. Been using Grant's method for balancing variance of value with the mean of the variance of the parameter where the sum should be near 1.0 for example:
p1.1<-var(data.table(mods1[[i]]$parameters)[Label%like%"Size_DblN_peak_Fishery"&Pr_type=="dev"]$Value) p1.2<-mean(data.table(mods1[[i]]$parameters)[Label%like%"Size_DblN_peak_Fishery"&Pr_type=="dev"]$Parm_StDev)^2 p1=p1.1+p1.2
where p1 should be near 1.0
This is done iteratively for all the varying selectivity parameters... There has to be a better way...
Need to develop a script to output tuning values for annually varying selectivity parameters. Been using Grant's method for balancing variance of value with the mean of the variance of the parameter where the sum should be near 1.0 for example:
p1.1<-var(data.table(mods1[[i]]$parameters)[Label%like%"Size_DblN_peak_Fishery"&Pr_type=="dev"]$Value) p1.2<-mean(data.table(mods1[[i]]$parameters)[Label%like%"Size_DblN_peak_Fishery"&Pr_type=="dev"]$Parm_StDev)^2 p1=p1.1+p1.2
where p1 should be near 1.0
This is done iteratively for all the varying selectivity parameters... There has to be a better way...