Estimating survival and exploitation rates from a tagging study on juvenile Atlantic halibut

Table 8. Estimates of survival and tag-recovery from the Hoenig et al (1988) model denoted by M(*) F(t) assuming
Instantaneous natural survival rate (M)
	Estimate	SE
All years	0.084	0.026
Instantaneous Tag-recovery rates ()
Year	Estimate	SE
1995	0.036	0.008
1996	0.042	0.009
1997	0.048	0.009
1998	0.033	0.005
1999	0.036	0.004
2000	0.033	0.004
2001	0.026	0.004
2002	0.022	0.004
2003	0.024	0.004
2004	0.013	0.003
2005	0.006	0.002
2006	0.007	0.003
2007	0.004	0.002

See Hoenig et al (1988) CJFAS 55, 1466-1476

f(i) - instantaneous fishing mortality rate

this rate.

It is assumed that tag shedding is zero as no double tags were applied to estimate

this rate.

This code can be used directly with the online surviv program at Patuxent.

Visit: http://www.mbr-pwrc.usgs.gov/software/surviv.html

Recoveries took place from 1999 to 2007 (13 years).

Parms 1 to 13 are the F

parm 14 is M */

------ Cut and paste code below into the code box at ---------

http://www.mbr-pwrc.usgs.gov/software/surviv.html
-----------------------------------------------------------------------

proc title 'Halibut model - Hoenig et al (1998) model M(*) F(t) .1';

inline prd=1;

inline sump=0;

inline do j=i,13 /* loop for each recovery year */;

inline tmp=f(j)+s(14);

inline srv=exp(-tmp);

inline p(j-i+1,i)=prd*(1-srv)*f(j)/tmp;

inline sump=sump+p(j-i+1,i);

inline prd=prd*srv;

inline end do;

inline p(14-i+1,i)=1-sump /* compute last cell (not recovered) */;

inline end do;

cohort=648 ; 22 :; 24 :; 29 :; 21 :; 20 :; 9 :; 7 :; 1 :;

3 :; 0 :; 0 :; 1 :; 0 :;

cohort=34 ; 0 :; 0 :; 0 :; 2 :; 0 :; 0 :; 2 :; 0 :; 0 :; 1 :; 1 :; 0 :;

cohort=200 ; 4 :; 3 :; 7 :; 6 :; 3 :; 2 :; 0 :; 1 :; 0 :; 0 :; 0 :;

cohort=887 ; 24 :; 36 :; 35 :; 15 :; 10 :; 7 :; 2 :; 4 :; 3 :; 1 :;

cohort=892 ; 12 :; 23 :; 27 :; 11 :; 12 :; 9 :; 3 :; 2 :; 0 :;

cohort=349 ; 1 :; 9 :; 9 :; 8 :; 2 :; 0 :; 0 :; 0 :;

cohort=799 ; 11 :; 19 :; 22 :; 12 :; 3 :; 5 :; 5 :;

s(1)=f(1);

s(2)=f(2);

s(3)=f(3);

s(4)=f(4);

s(5)=f(5);

s(6)=f(6);

s(7)=f(7);

s(8)=f(8);

s(9)=f(9);

s(10)=f(10);

s(11)=f(11);

s(12)=f(12);

s(13)=f(13);

s(14)=m;
/* Now to estimate the model */

proc estimate novar maxfn=32001 name=base;