Differences

This shows you the differences between two versions of the page.

--- admb_real:uncertain_catch [2010/09/14 15:30] – lennartfr
+++ admb_real:uncertain_catch [2010/09/27 12:06] (current) – lennartfr
@@ Line 1: / Line 1: @@
+====== NOTE: This section is largely obsolete! ======
+FIXME This model may well not be the one we will end up using, so the results in this section should be taken with a pinch of salt. //--Lennart, 20.09.2010.//
 ===Models with uncertainties in the catch data===
-It is unrealistic to assume that the catch data avaible from the ICES AFWG reports are without error. For this reason we make the following enhacement to the model:
+It is unrealistic to assume that the catch data avaible from the ICES AFWG reports are without error. For this reason we make the following enhancement to the model:
 Replace the recurrence relation:
-^N<sub>i,t</sub>      ^= ^(N<sub>i,t-1</sub>-C<sub>i,t-1</sub>)e<sup>-M</sup>  ^
+^N<sub>i,t</sub>      ^= ^(N<sub>i,t-1</sub> - C<sub>i,t-1</sub>)e<sup>-M</sup>  ^
 with:
-^N<sub>i,t</sub>      ^= ^(N<sub>i,t-1</sub>-C<sub>i,t-1</sub>)e<sup>-M</sup>  ^
+^N<sub>i,t</sub>      ^= ^(N<sub>i,t-1</sub>-C<sub>i,t-1</sub> exp[us<sub>u</sub>-0.5s<sub>u</sub><sup>2</sup>] )e<sup>-M</sup>  ^
-where u is a latent variable with expectation 0 , and standard deviation s<sub>u</sub>.
+where u is a latent variable vector with one component for each component of the catch data. It has expectation 0 , and standard deviation 1. The non-latent parameter s<sub>u</sub> controls the standard deviation of the catch itself.
+:!: When this model is added to the previously defined models they have too many degrees of freedom. Therefore the catchability should be controlled somehow. We try the model with two separate catchabilities, and lock the catchability of 0-2 year olds to be 0.1. A few preliminary experimental results are given below:
+|{{:admb_real:preliminary_two_q_vpa_85_86.png?200|}} |{{:admb_real:preliminary_two_q_vpa_87_88.png?200|}} |{{:admb_real:preliminary_two_q_vpa_89_90.png?200|}} |
+Although the population trajectories are not too different from the [[admb_real:vpa_estimates|VPA trajectories]] for the first and third plot, the parameter s<sub>u</sub> is at its lower bound. For the second plot the mortality is at its lower bound. Only two cohorts are present per run, since otherwise the ADMB fails to converge. ADMB also fails to converge for the cohorts born in 1991/1992, and 1993/1994, suggesting that further control of the catchability must be enforced somehow.