RDSURVIV

Computes estimates of survival and capture probability using capture-recapture data with robust-design models and recoveries.

Kendall, W.L., J.D. Nichols, and J.E. Hines. 1997. Estimating Temporary Emigration using Capture-Recapture data with Pollocks Robust Design. Ecology 78:563-578.
Lindberg, M. S., W. L. Kendall, J. E. Hines, and M. G. Anderson. 2001. Combining band recovery data and Pollocks robust design to model temporary and permanent emigration. Biometrics 57:273-281.

Run RDSURVIV online

Hines, J.E. 1996. RDSURVIV User's Manual. Biological Resources Div., USGS, Laurel, MD

To convert capture-history data to RDSURVIV input, paste capture-histories into input box below, then click 'convert' button.

00 00 11 00 00 10 10
00 00 00 01 00 11 01
00 00 00 01 00 00 00
00 00 00 01 01 11 11
00 00 00 01 00 00 00
00 00 00 00 10 10 11
00 00 00 00 11 10 11
11 01 00 01 10 11 00
00 00 00 00 00 10 00
00 00 00 00 00 10 00
00 00 00 00 00 11 00
00 00 00 00 00 11 11
00 00 00 00 00 10 00
00 00 00 00 00 10 11
11 01 00 00 00 00 00
11 00 00 00 00 00 00
11 00 00 00 00 00 00
11 10 00 00 00 00 00
11 01 00 00 00 11 11
01 01 00 00 00 10 01
11 01 00 00 00 00 00
11 00 00 00 00 11 00
11 11 00 00 00 10 11
11 00 00 00 00 00 00
10 00 00 00 00 00 00
11 01 00 00 00 00 00
11 01 00 00 00 00 11
11 00 00 00 00 00 00
11 11 00 00 00 11 00
10 00 00 00 00 00 00
01 00 00 00 00 00 00
01 00 00 00 00 00 00
00 11 00 00 00 00 00
00 00 11 00 00 00 00
00 10 01 00 00 00 10
00 01 01 00 00 00 00
00 00 00 00 11 00 00
00 00 11 00 01 11 11
00 00 01 00 11 11 11
00 01 00 00 00 00 00
00 01 00 01 00 10 11
01 00 00 00 00 00 00

You can enter an ASCII text file in the box below. The file should look like a RDSURVIV input file (eg:

PROC TITLE ;
PROC MODEL NPAR=45 ADDCELL PRIMARY PERIODS=7 SECONDARY=2 2 2 2 2 2 2;
COHORT=20 /* u(1) */;
14:  /* 11 */;
2:  /* 10 */;
COHORT=5 /* u(2) */;
1:  /* 11 */;
1:  /* 10 */;
COHORT=10 /* m(1,2) */;
2:  /* 11 */;
1:  /* 10 */;
COHORT=4 /* u(3) */;
3:  /* 11 */;
0:  /* 10 */;
COHORT=0 /* m(1,3) */;
0:  /* 11 */;
0:  /* 10 */;
COHORT=2 /* m(2,3) */;
0:  /* 11 */;
0:  /* 10 */;
COHORT=4 /* u(4) */;
0:  /* 11 */;
0:  /* 10 */;
COHORT=0 /* m(1,4) */;
0:  /* 11 */;
0:  /* 10 */;
COHORT=2 /* m(2,4) */;
0:  /* 11 */;
0:  /* 10 */;
COHORT=0 /* m(3,4) */;
0:  /* 11 */;
0:  /* 10 */;
COHORT=3 /* u(5) */;
2:  /* 11 */;
1:  /* 10 */;
COHORT=0 /* m(1,5) */;
0:  /* 11 */;
0:  /* 10 */;
COHORT=0 /* m(2,5) */;
0:  /* 11 */;
0:  /* 10 */;
COHORT=2 /* m(3,5) */;
1:  /* 11 */;
0:  /* 10 */;
COHORT=2 /* m(4,5) */;
0:  /* 11 */;
1:  /* 10 */;
COHORT=6 /* u(6) */;
2:  /* 11 */;
4:  /* 10 */;
COHORT=1 /* m(1,6) */;
1:  /* 11 */;
0:  /* 10 */;
COHORT=4 /* m(2,6) */;
2:  /* 11 */;
2:  /* 10 */;
COHORT=1 /* m(3,6) */;
0:  /* 11 */;
1:  /* 10 */;
COHORT=2 /* m(4,6) */;
1:  /* 11 */;
1:  /* 10 */;
COHORT=6 /* m(5,6) */;
4:  /* 11 */;
2:  /* 10 */;
COHORT=0 /* u(7) */;
0:  /* 11 */;
0:  /* 10 */;
COHORT=0 /* m(1,7) */;
0:  /* 11 */;
0:  /* 10 */;
COHORT=1 /* m(2,7) */;
1:  /* 11 */;
0:  /* 10 */;
COHORT=1 /* m(3,7) */;
0:  /* 11 */;
1:  /* 10 */;
COHORT=0 /* m(4,7) */;
0:  /* 11 */;
0:  /* 10 */;
COHORT=0 /* m(5,7) */;
0:  /* 11 */;
0:  /* 10 */;
COHORT=13 /* m(6,7) */;
10:  /* 11 */;
1:  /* 10 */;
COHORT=20   /* R(1)  */;
10:  /*  m(1,2)  */;
0:  /*  m(1,3)  */;
0:  /*  m(1,4)  */;
0:  /*  m(1,5)  */;
1:  /*  m(1,6)  */;
0:  /*  m(1,7)  */;
COHORT=15   /* R(2)  */;
2:  /*  m(2,3)  */;
2:  /*  m(2,4)  */;
0:  /*  m(2,5)  */;
4:  /*  m(2,6)  */;
1:  /*  m(2,7)  */;
COHORT=6   /* R(3)  */;
0:  /*  m(3,4)  */;
2:  /*  m(3,5)  */;
1:  /*  m(3,6)  */;
1:  /*  m(3,7)  */;
COHORT=6   /* R(4)  */;
2:  /*  m(4,5)  */;
2:  /*  m(4,6)  */;
0:  /*  m(4,7)  */;
COHORT=7   /* R(5)  */;
6:  /*  m(5,6)  */;
0:  /*  m(5,7)  */;
COHORT=20   /* R(6)  */;
13:  /*  m(6,7)  */;
LABELS;
S(1)=p(1,1); S(2)=p(1,2); S(3)=c(1,2); 
S(4)=p(2,1); S(5)=p(2,2); S(6)=c(2,2);
S(7)=p(3,1); S(8)=p(3,2); S(9)=c(3,2);
S(10)=p(4,1); S(11)=p(4,2); S(12)=c(4,2);
S(13)=p(5,1); S(14)=p(5,2); S(15)=c(5,2); 
S(16)=p(6,1); S(17)=p(6,2); S(18)=c(6,2);
S(19)=p(7,1); S(20)=p(7,2); S(21)=c(7,2);
S(22)=PHI(1); S(23)=PHI(2); S(24)=PHI(3); 
S(25)=PHI(4); S(26)=PHI(5); S(27)=PHI(6);
S(28)=GAM'(2); S(29)=GAM'(3); S(30)=GAM'(4);
S(31)=GAM'(5); S(32)=GAM'(6); S(33)=GAM'(7);
S(34)=GAM(3); S(35)=GAM(4); S(36)=GAM(5);
S(37)=GAM(6); S(38)=GAM(7);
S(39)=THETA(1); S(40)=THETA(2); S(41)=THETA(3);
S(42)=THETA(4); S(43)=THETA(5); S(44)=THETA(6); S(45)=THETA(7);
PROC ESTIMATE MAXFN=32000 NOVAR NAME=M(0)Sng;
CONSTRAINTS;
S(2)=S(1) /* p(1,2)=p(1,1) */;
S(3)=S(1) /* c(1,2)=p(1,1) */;
S(5)=S(4) /* p(2,2)=p(2,1) */;
S(6)=S(4) /* c(2,2)=p(2,1) */;
S(8)=S(7) /* p(3,2)=p(3,1) */;
S(9)=S(7) /* c(3,2)=p(3,1) */;
S(11)=S(10) /* p(4,2)=p(4,1) */;
S(12)=S(10) /* c(4,2)=p(4,1) */;
S(14)=S(13) /* p(5,2)=p(5,1) */;
S(15)=S(13) /* c(5,2)=p(5,1) */;
S(17)=S(16) /* p(6,2)=p(6,1) */;
S(18)=S(16) /* c(6,2)=p(6,1) */;
S(20)=S(19) /* p(7,2)=p(7,1) */;
S(21)=S(19) /* c(7,2)=p(7,1) */;
S(23)=S(22) /* PHI(2)=PHI(1) */;
S(24)=S(22) /* PHI(3)=PHI(1) */;
S(25)=S(22) /* PHI(4)=PHI(1) */;
S(26)=S(22) /* PHI(5)=PHI(1) */;
S(27)=S(22) /* PHI(6)=PHI(1) */;
S(28)=0 /* GAM'(2)=0 -> no emigration */;
S(29)=0 /* GAM'(3)=0 -> no emigration */;
S(30)=0 /* GAM'(4)=0 -> no emigration */;
S(31)=0 /* GAM'(5)=0 -> no emigration */;
S(32)=0 /* GAM'(6)=0 -> no emigration */;
S(33)=0 /* GAM'(7)=0 -> no emigration */;
S(34)=0 /* GAM(3)=0 -> no emigration */;
S(35)=0 /* GAM(4)=0 -> no emigration */;
S(36)=0 /* GAM(5)=0 -> no emigration */;
S(37)=0 /* GAM(6)=0 -> no emigration */;
S(38)=0 /* GAM(7)=0 -> no emigration */;
S(39)=1 /* THETA(1) not used in this model */;
S(40)=1 /* THETA(2) not used in this model */;
S(41)=1 /* THETA(3) not used in this model */;
S(42)=1 /* THETA(4) not used in this model */;
S(43)=1 /* THETA(5) not used in this model */;
S(44)=1 /* THETA(6) not used in this model */;
S(45)=1 /* THETA(7) not used in this model */;
PROC ESTIMATE MAXFN=32000 NOVAR NAME=M(b)Sng;
INITIAL; RETAIN=A; CONSTRAINTS;
S(2)=S(1) /* p(1,2)=p(1,1) */;
S(5)=S(4) /* p(2,2)=p(2,1) */;
S(8)=S(7) /* p(3,2)=p(3,1) */;
S(11)=S(10) /* p(4,2)=p(4,1) */;
S(14)=S(13) /* p(5,2)=p(5,1) */;
S(17)=S(16) /* p(6,2)=p(6,1) */;
S(20)=S(19) /* p(7,2)=p(7,1) */;
S(28)=0 /* GAM'(2)=0 -> no emigration */;
S(29)=0 /* GAM'(3)=0 -> no emigration */;
S(30)=0 /* GAM'(4)=0 -> no emigration */;
S(31)=0 /* GAM'(5)=0 -> no emigration */;
S(32)=0 /* GAM'(6)=0 -> no emigration */;
S(33)=0 /* GAM'(7)=0 -> no emigration */;
S(34)=0 /* GAM(3)=0 -> no emigration */;
S(35)=0 /* GAM(4)=0 -> no emigration */;
S(36)=0 /* GAM(5)=0 -> no emigration */;
S(37)=0 /* GAM(6)=0 -> no emigration */;
S(38)=0 /* GAM(7)=0 -> no emigration */;
S(39)=1 /* THETA(1) not used in this model */;
S(40)=1 /* THETA(2) not used in this model */;
S(41)=1 /* THETA(3) not used in this model */;
S(42)=1 /* THETA(4) not used in this model */;
S(43)=1 /* THETA(5) not used in this model */;
S(44)=1 /* THETA(6) not used in this model */;
S(45)=1 /* THETA(7) not used in this model */;
PROC ESTIMATE MAXFN=32000 NOVAR NAME=M(t)Stng;
INITIAL; RETAIN=A; CONSTRAINTS;
S(3)=S(2) /* c(1,2)=c(1,2) */;
S(6)=S(5) /* c(2,2)=c(2,2) */;
S(9)=S(8) /* c(3,2)=c(3,2) */;
S(12)=S(11) /* c(4,2)=c(4,2) */;
S(15)=S(14) /* c(5,2)=c(5,2) */;
S(18)=S(17) /* c(6,2)=c(6,2) */;
S(21)=S(20) /* c(7,2)=c(7,2) */;
S(28)=0 /* GAM'(2)=0 -> no emigration */;
S(29)=0 /* GAM'(3)=0 -> no emigration */;
S(30)=0 /* GAM'(4)=0 -> no emigration */;
S(31)=0 /* GAM'(5)=0 -> no emigration */;
S(32)=0 /* GAM'(6)=0 -> no emigration */;
S(33)=0 /* GAM'(7)=0 -> no emigration */;
S(34)=0 /* GAM(3)=0 -> no emigration */;
S(35)=0 /* GAM(4)=0 -> no emigration */;
S(36)=0 /* GAM(5)=0 -> no emigration */;
S(37)=0 /* GAM(6)=0 -> no emigration */;
S(38)=0 /* GAM(7)=0 -> no emigration */;
S(39)=1 /* THETA(1) not used in this model */;
S(40)=1 /* THETA(2) not used in this model */;
S(41)=1 /* THETA(3) not used in this model */;
S(42)=1 /* THETA(4) not used in this model */;
S(43)=1 /* THETA(5) not used in this model */;
S(44)=1 /* THETA(6) not used in this model */;
S(45)=1 /* THETA(7) not used in this model */;
PROC ESTIMATE MAXFN=32000 NOVAR NAME=M(bt)Stng;
INITIAL; RETAIN=A;
S(45)=1; S(44)=1; S(43)=1; S(42)=1; S(41)=1; S(40)=1; S(39)=1;
CONSTRAINTS;
S(3)=-1 /* c(1,2)=-1 -> c(1,2)=p(1,2)**THETA(1) */;
S(6)=-1 /* c(2,2)=-1 -> c(2,2)=p(2,2)**THETA(2) */;
S(9)=-1 /* c(3,2)=-1 -> c(3,2)=p(3,2)**THETA(3) */;
S(12)=-1 /* c(4,2)=-1 -> c(4,2)=p(4,2)**THETA(4) */;
S(15)=-1 /* c(5,2)=-1 -> c(5,2)=p(5,2)**THETA(5) */;
S(18)=-1 /* c(6,2)=-1 -> c(6,2)=p(6,2)**THETA(6) */;
S(21)=-1 /* c(7,2)=-1 -> c(7,2)=p(7,2)**THETA(7) */;
S(28)=0 /* GAM'(2)=0 -> no emigration */;
S(29)=0 /* GAM'(3)=0 -> no emigration */;
S(30)=0 /* GAM'(4)=0 -> no emigration */;
S(31)=0 /* GAM'(5)=0 -> no emigration */;
S(32)=0 /* GAM'(6)=0 -> no emigration */;
S(33)=0 /* GAM'(7)=0 -> no emigration */;
S(34)=0 /* GAM(3)=0 -> no emigration */;
S(35)=0 /* GAM(4)=0 -> no emigration */;
S(36)=0 /* GAM(5)=0 -> no emigration */;
S(37)=0 /* GAM(6)=0 -> no emigration */;
S(38)=0 /* GAM(7)=0 -> no emigration */;
S(39)<20 /* THETA(1) */;
S(40)<20 /* THETA(2) */;
S(41)<20 /* THETA(3) */;
S(42)<20 /* THETA(4) */;
S(43)<20 /* THETA(5) */;
S(44)<20 /* THETA(6) */;
S(45)<20 /* THETA(7) */;
PROC TEST; PROC STOP;