-# Imach version 0.98q5, August 2015,INED-EUROREVES-Institut de longevite-Japan Society for the Promotion of Science (Grant-in-Aid for Scientific Research 25293121), Intel Software 2015
+# IMaCh-0.99r19
# biaspar.gp
set datafile missing 'NaNq'
-cd "/Users/nbrouard/Documents/imach/cvs/imach/build/osx/html/doc"
+cd "/Users/nbrouard/Documents/imach/cvs/imach/html/doc"
-# 1st: Period (stable) prevalence with CI: 'vpl' files
+#Diagram of the model
-set out "biaspar/vbiaspar1_1.png"
+delta=0.03;delta2=0.07;unset arrow;
+yoff=(2 > 2? 0:1);
-#set out "vbiaspar1_1.png"
+#Peripheral arrows
+set for [i=1:2] for [j=1:2] arrow i*10+j from cos(pi*((1-(2/2)*2./2)/2+(i-1)*2./2))-(i!=j?(i-j)/abs(i-j)*delta:0), yoff +sin(pi*((1-(2/2)*2./2)/2+(i-1)*2./2)) + (i!=j?(i-j)/abs(i-j)*delta:0) rto -0.95*(cos(pi*((1-(2/2)*2./2)/2+(i-1)*2./2))+(i!=j?(i-j)/abs(i-j)*delta:0) - cos(pi*((1-(2/2)*2./2)/2+(j-1)*2./2)) + (i!=j?(i-j)/abs(i-j)*delta2:0)), -0.95*(sin(pi*((1-(2/2)*2./2)/2+(i-1)*2./2)) + (i!=j?(i-j)/abs(i-j)*delta:0) - sin(pi*((1-(2/2)*2./2)/2+(j-1)*2./2))+( i!=j?(i-j)/abs(i-j)*delta2:0)) ls (i < j? 1:2)
+
+#Centripete arrows (turning in other direction (1-i) instead of (i-1))
+set for [i=1:2] arrow (2+1)*10+i from cos(pi*((1-(2/2)*2./2)/2+(1-i)*2./2))-(i!=j?(i-j)/abs(i-j)*delta:0), yoff +sin(pi*((1-(2/2)*2./2)/2+(1-i)*2./2)) + (i!=j?(i-j)/abs(i-j)*delta:0) rto -0.80*(cos(pi*((1-(2/2)*2./2)/2+(1-i)*2./2))+(i!=j?(i-j)/abs(i-j)*delta:0) ), -0.80*(sin(pi*((1-(2/2)*2./2)/2+(1-i)*2./2)) + (i!=j?(i-j)/abs(i-j)*delta:0) + yoff ) ls 4
+
+#show arrow
+unset label
+
+#States labels, starting from 2 (2-i) instead of (1-i), was (i-1)
+set for [i=1:2] label i sprintf("State %d",i) center at cos(pi*((1-(2/2)*2./2)/2+(2-i)*2./2)), yoff+sin(pi*((1-(2/2)*2./2)/2+(2-i)*2./2)) font "helvetica, 16" tc rgbcolor "blue"
+
+set label 2+1 sprintf("State %d",2+1) center at 0.,0. font "helvetica, 16" tc rgbcolor "red"
+
+#show label
+unset border;unset xtics; unset ytics;
+
+
+set ter svg size 640, 480;set out "biaspar/D_biaspar_.svg"
+unset log y; plot [-1.2:1.2][yoff-1.2:1.2] 1/0 not; set out;reset;
+
+# Contributions to the Likelihood, mle >=1. For mle=4 no interpolation, pure matrix products.
+#
+
+ set log y; unset log x;set xlabel "Age"; set ylabel "Likelihood (-2Log(L))";
+set ter pngcairo size 640, 480
+set out "biaspar/ILK_biaspar-dest.png";
+set log y;plot "biaspar/ILK_biaspar.txt" u 2:(-$13):6 t "All sample, transitions colored by destination" with dots lc variable; set out;
+
+set out "biaspar/ILK_biaspar-ori.png";
+set log y;plot "biaspar/ILK_biaspar.txt" u 2:(-$13):5 t "All sample, transitions colored by origin" with dots lc variable; set out;
+
+
+set out "biaspar/ILK_biaspar-p1j.png";set ylabel "Probability for each individual/wave";unset log;
+# plot weighted, mean weight should have point size of 0.5
+ plot "biaspar/ILK_biaspar.txt" u 2:($5 == 1 && $6==1 ? $10 : 1/0):($12/4.):6 t "p11" with points pointtype 7 ps variable lc variable \
+,\
+ "" u 2:($5 == 1 && $6==2 ? $10 : 1/0):($12/4.):6 t "p12" with points pointtype 7 ps variable lc variable ,\
+ "" u 2:($5 == 1 && $6==3 ? $10 : 1/0):($12/4.):6 t "p13" with points pointtype 7 ps variable lc variable ;
+set out; unset ylabel;
+
+set out "biaspar/ILK_biaspar-p2j.png";set ylabel "Probability for each individual/wave";unset log;
+# plot weighted, mean weight should have point size of 0.5
+ plot "biaspar/ILK_biaspar.txt" u 2:($5 == 2 && $6==1 ? $10 : 1/0):($12/4.):6 t "p21" with points pointtype 7 ps variable lc variable \
+,\
+ "" u 2:($5 == 2 && $6==2 ? $10 : 1/0):($12/4.):6 t "p22" with points pointtype 7 ps variable lc variable ,\
+ "" u 2:($5 == 2 && $6==3 ? $10 : 1/0):($12/4.):6 t "p23" with points pointtype 7 ps variable lc variable ;
+set out; unset ylabel;
+
+set out;unset log
+
+# 1st: Forward (stable period) prevalence with CI: 'VPL_' files and live state =1
+#
+
+set out "biaspar/V_biaspar_1-1-1.svg"
+
+#set out "V_biaspar_1-1-1.svg"
+set title "Alive state 1 ()" font "Helvetica,12"
set xlabel "Age"
set ylabel "Probability"
-set ter png small size 320, 240
-plot [70:95] "biaspar/vplrbiaspar.txt" every :::0::0 u 1:2 "%lf %lf (%lf) %*lf (%*lf)" t"Period (stable) prevalence" w l lt 0,"biaspar/vplrbiaspar.txt" every :::0::0 u 1:($2+1.96*$3) "%lf %lf (%lf) %*lf (%*lf)" t"95% CI" w l lt 1,"biaspar/vplrbiaspar.txt" every :::0::0 u 1:($2-1.96*$3) "%lf %lf (%lf) %*lf (%*lf)" t"" w l lt 1,"biaspar/prbiaspar.txt" every :::0::0 u 1:($2) t"Observed prevalence " w l lt 2
-set out "biaspar/vbiaspar2_1.png"
+set ter svg size 640, 480
+plot [70:95] "biaspar/VPL_biaspar.txt" every :::0::0 u 1:($2==1 ? $3:1/0) "%lf %lf %lf (%lf) %*lf (%*lf)" t"Forward prevalence" w l lt 0,"biaspar/VPL_biaspar.txt" every :::0::0 u 1:($2==1 ? $3+1.96*$4 : 1/0) "%lf %lf %lf (%lf) %*lf (%*lf)" t"95% CI" w l lt 1,"biaspar/VPL_biaspar.txt" every :::0::0 u 1:($2==1 ? $3-1.96*$4 : 1/0) "%lf %lf %lf (%lf) %*lf (%*lf)" t"" w l lt 1,"biaspar/P_biaspar.txt" u 1:(($2)) t 'Observed prevalence in state 1' with line lt 3
+set out ;unset title;
+
+# 1st: Forward (stable period) prevalence with CI: 'VPL_' files and live state =2
+#
-#set out "vbiaspar2_1.png"
+set out "biaspar/V_biaspar_2-1-1.svg"
+
+#set out "V_biaspar_2-1-1.svg"
+set title "Alive state 2 ()" font "Helvetica,12"
set xlabel "Age"
set ylabel "Probability"
-set ter png small size 320, 240
-plot [70:95] "biaspar/vplrbiaspar.txt" every :::0::0 u 1:2 "%lf %*lf (%*lf) %lf (%lf)" t"Period (stable) prevalence" w l lt 0,"biaspar/vplrbiaspar.txt" every :::0::0 u 1:($2+1.96*$3) "%lf %*lf (%*lf) %lf (%lf)" t"95% CI" w l lt 1,"biaspar/vplrbiaspar.txt" every :::0::0 u 1:($2-1.96*$3) "%lf %*lf (%*lf) %lf (%lf)" t"" w l lt 1,"biaspar/prbiaspar.txt" every :::0::0 u 1:($6) t"Observed prevalence " w l lt 2
-# 2nd: Total life expectancy with CI: 't' files
+set ter svg size 640, 480
+plot [70:95] "biaspar/VPL_biaspar.txt" every :::0::0 u 1:($2==1 ? $3:1/0) "%lf %lf %*lf (%*lf) %lf (%lf)" t"Forward prevalence" w l lt 0,"biaspar/VPL_biaspar.txt" every :::0::0 u 1:($2==1 ? $3+1.96*$4 : 1/0) "%lf %lf %*lf (%*lf) %lf (%lf)" t"95% CI" w l lt 1,"biaspar/VPL_biaspar.txt" every :::0::0 u 1:($2==1 ? $3-1.96*$4 : 1/0) "%lf %lf %*lf (%*lf) %lf (%lf)" t"" w l lt 1,"biaspar/P_biaspar.txt" u 1:(($5)) t 'Observed prevalence in state 2' with line lt 3
+set out ;unset title;
+
+# 2nd: Total life expectancy with CI: 't' files
+#
+
+set out "biaspar/E_biaspar_1-1.svg"
-set out "biaspar/ebiaspar1.png"
+set label "popbased 0 ()" at graph 0.98,0.5 center rotate font "Helvetica,12"
set ylabel "Years"
-set ter png small size 320, 240
-plot [70:95] "biaspar/trbiaspar.txt" every :::0::0 u 1:2 "%lf %lf (%lf) %*lf (%*lf) %*lf (%*lf)" t"TLE" w l ,"biaspar/trbiaspar.txt" every :::0::0 u 1:($2-$3*2) "%lf %lf (%lf) %*lf (%*lf) %*lf (%*lf)" t"" w l lt 0,"biaspar/trbiaspar.txt" every :::0::0 u 1:($2+$3*2) "%lf %lf (%lf) %*lf (%*lf) %*lf (%*lf)" t"" w l lt 0,"biaspar/trbiaspar.txt" every :::0::0 u 1:2 "%lf %*lf (%*lf) %lf (%lf) %*lf (%*lf)" t"LE in state (1)" w l ,"biaspar/trbiaspar.txt" every :::0::0 u 1:($2-$3*2) "%lf %*lf (%*lf) %lf (%lf) %*lf (%*lf)" t"" w l lt 0,"biaspar/trbiaspar.txt" every :::0::0 u 1:($2+$3*2) "%lf %*lf (%*lf) %lf (%lf) %*lf (%*lf)" t"" w l lt 0,"biaspar/trbiaspar.txt" every :::0::0 u 1:2 "%lf %*lf (%*lf) %*lf (%*lf) %lf (%lf)" t"LE in state (2)" w l ,"biaspar/trbiaspar.txt" every :::0::0 u 1:($2-$3*2) "%lf %*lf (%*lf) %*lf (%*lf) %lf (%lf)" t"" w l lt 0,"biaspar/trbiaspar.txt" every :::0::0 u 1:($2+$3*2) "%lf %*lf (%*lf) %*lf (%*lf) %lf (%lf)" t"" w l lt 0
-set out "biaspar/expbiaspar11.png"
-set ter png small size 320, 240
-plot [70:95] "biaspar/erbiaspar.txt" every :::0::0 u 1:2 t "e11" w l ,"biaspar/erbiaspar.txt" every :::0::0 u 1:3 t "e12" w l ,"biaspar/erbiaspar.txt" every :::0::0 u 1:4 t "e1." w l
-set out "biaspar/expbiaspar21.png"
-set ter png small size 320, 240
-plot [70:95] "biaspar/erbiaspar.txt" every :::0::0 u 1:5 t "e21" w l ,"biaspar/erbiaspar.txt" every :::0::0 u 1:6 t "e22" w l ,"biaspar/erbiaspar.txt" every :::0::0 u 1:7 t "e2." w l
-#
-#
-#CV preval stable (period): 'pij' files, cov=1 state=1
-set out "biaspar/pbiaspar1_1.png"
+set ter svg size 640, 480
+plot [70:95] "biaspar/T_biaspar.txt" every :::0::0 u 1:($2==0 && $4!=0 ?$4 : 1/0) "%lf %lf %lf %lf (%lf) %*lf (%*lf) %*lf (%*lf)" t"TLE" w l lt 1, \
+"biaspar/T_biaspar.txt" every :::0::0 u 1:($2==0 && $4!=0 ? $4-$5*2 : 1/0) "%lf %lf %lf %lf (%lf) %*lf (%*lf) %*lf (%*lf)" t"" w l lt 0,"biaspar/T_biaspar.txt" every :::0::0 u 1:($2==0 && $4!=0 ? $4+$5*2 : 1/0) "%lf %lf %lf %lf (%lf) %*lf (%*lf) %*lf (%*lf)" t"" w l lt 0,\
+"biaspar/T_biaspar.txt" every :::0::0 u 1:($2==0 && $4!=0 ?$4 : 1/0) "%lf %lf %lf %*lf (%*lf) %lf (%lf) %*lf (%*lf)" t"LE in state (1)" w l lt 3, \
+"biaspar/T_biaspar.txt" every :::0::0 u 1:($2==0 && $4!=0 ? $4-$5*2 : 1/0) "%lf %lf %lf %*lf (%*lf) %lf (%lf) %*lf (%*lf)" t"" w l lt 0,"biaspar/T_biaspar.txt" every :::0::0 u 1:($2==0 && $4!=0 ? $4+$5*2 : 1/0) "%lf %lf %lf %*lf (%*lf) %lf (%lf) %*lf (%*lf)" t"" w l lt 0,\
+"biaspar/T_biaspar.txt" every :::0::0 u 1:($2==0 && $4!=0 ?$4 : 1/0) "%lf %lf %lf %*lf (%*lf) %*lf (%*lf) %lf (%lf)" t"LE in state (2)" w l lt 4, \
+"biaspar/T_biaspar.txt" every :::0::0 u 1:($2==0 && $4!=0 ? $4-$5*2 : 1/0) "%lf %lf %lf %*lf (%*lf) %*lf (%*lf) %lf (%lf)" t"" w l lt 0,"biaspar/T_biaspar.txt" every :::0::0 u 1:($2==0 && $4!=0 ? $4+$5*2 : 1/0) "%lf %lf %lf %*lf (%*lf) %*lf (%*lf) %lf (%lf)" t"" w l lt 0
+set label "popbased 1 ()" at graph 0.98,0.5 center rotate font "Helvetica,12"
+
+replot "biaspar/T_biaspar.txt" every :::0::0 u 1:($2==1 && $4!=0 ?$4 : 1/0) "%lf %lf %lf %lf (%lf) %*lf (%*lf) %*lf (%*lf)" t"TLE" w l lt 1, \
+"biaspar/T_biaspar.txt" every :::0::0 u 1:($2==1 && $4!=0 ? $4-$5*2 : 1/0) "%lf %lf %lf %lf (%lf) %*lf (%*lf) %*lf (%*lf)" t"" w l lt 0,"biaspar/T_biaspar.txt" every :::0::0 u 1:($2==1 && $4!=0 ? $4+$5*2 : 1/0) "%lf %lf %lf %lf (%lf) %*lf (%*lf) %*lf (%*lf)" t"" w l lt 0,\
+"biaspar/T_biaspar.txt" every :::0::0 u 1:($2==1 && $4!=0 ?$4 : 1/0) "%lf %lf %lf %*lf (%*lf) %lf (%lf) %*lf (%*lf)" t"LE in state (1)" w l lt 3, \
+"biaspar/T_biaspar.txt" every :::0::0 u 1:($2==1 && $4!=0 ? $4-$5*2 : 1/0) "%lf %lf %lf %*lf (%*lf) %lf (%lf) %*lf (%*lf)" t"" w l lt 0,"biaspar/T_biaspar.txt" every :::0::0 u 1:($2==1 && $4!=0 ? $4+$5*2 : 1/0) "%lf %lf %lf %*lf (%*lf) %lf (%lf) %*lf (%*lf)" t"" w l lt 0,\
+"biaspar/T_biaspar.txt" every :::0::0 u 1:($2==1 && $4!=0 ?$4 : 1/0) "%lf %lf %lf %*lf (%*lf) %*lf (%*lf) %lf (%lf)" t"LE in state (2)" w l lt 4, \
+"biaspar/T_biaspar.txt" every :::0::0 u 1:($2==1 && $4!=0 ? $4-$5*2 : 1/0) "%lf %lf %lf %*lf (%*lf) %*lf (%*lf) %lf (%lf)" t"" w l lt 0,"biaspar/T_biaspar.txt" every :::0::0 u 1:($2==1 && $4!=0 ? $4+$5*2 : 1/0) "%lf %lf %lf %*lf (%*lf) %*lf (%*lf) %lf (%lf)" t"" w l lt 0
+set out;set out "biaspar/E_biaspar_1-1.svg"; replot; set out; unset label;
+
+
+# 3d: Life expectancy with EXP_ files: combination=1 state=1
+#
+
+set out "biaspar/EXP_biaspar_1-1-1.svg"
+set label "()" at graph 0.98,0.5 center rotate font "Helvetica,12"
+set ter svg size 640, 480
+plot [70:95] "biaspar/E_biaspar.txt" every :::0::0 u 1:2 t "e11" w l ,"biaspar/E_biaspar.txt" every :::0::0 u 1:3 t "e12" w l ,"biaspar/E_biaspar.txt" every :::0::0 u 1:4 t "e1." w l
+
+# 3d: Life expectancy with EXP_ files: combination=1 state=2
+#
+
+set out "biaspar/EXP_biaspar_2-1-1.svg"
+set label "()" at graph 0.98,0.5 center rotate font "Helvetica,12"
+set ter svg size 640, 480
+plot [70:95] "biaspar/E_biaspar.txt" every :::0::0 u 1:5 t "e21" w l ,"biaspar/E_biaspar.txt" every :::0::0 u 1:6 t "e22" w l ,"biaspar/E_biaspar.txt" every :::0::0 u 1:7 t "e2." w l
+unset label;
+
+#
+#
+# Survival functions in state j : 'LIJ_' files, cov=1 state=1
+#
+
+set out "biaspar/LIJ_biaspar_1-1-1.svg"
+set label "Alive state 1 ()" at graph 0.98,0.5 center rotate font "Helvetica,12"
+set xlabel "Age"
+set ylabel "Probability to be alive"
+set ter svg size 640, 480
+unset log y
+plot [70:95] "biaspar/PIJ_biaspar.txt" u ($1==1 ? ($3):1/0):($4/($4+$5+$6)) t "l(1,1)" w l, '' u ($1==1 ? ($3):1/0):($7/($7+$8+$9)) t "l(2,1)" w l
+set out; unset label;
+
+#
+#
+# Survival functions in state j : 'LIJ_' files, cov=1 state=2
+#
+
+set out "biaspar/LIJ_biaspar_2-1-1.svg"
+set label "Alive state 2 ()" at graph 0.98,0.5 center rotate font "Helvetica,12"
+set xlabel "Age"
+set ylabel "Probability to be alive"
+set ter svg size 640, 480
+unset log y
+plot [70:95] "biaspar/PIJ_biaspar.txt" u ($1==1 ? ($3):1/0):($5/($4+$5+$6)) t "l(1,2)" w l, '' u ($1==1 ? ($3):1/0):($8/($7+$8+$9)) t "l(2,2)" w l
+set out; unset label;
+
+#
+#
+# Survival functions in state j and all livestates from state i by final state j: 'lij' files, cov=1 state=1
+#
+
+set out "biaspar/LIJT_biaspar_1-1-1.svg"
+set label "Alive state 1 ()" at graph 0.98,0.5 center rotate font "Helvetica,12"
+set xlabel "Age"
+set ylabel "Probability to be alive"
+set ter svg size 640, 480
+unset log y
+plot [70:95] "biaspar/PIJ_biaspar.txt" u (($1==1 && (floor($2)%5 == 0)) ? ($3):1/0):($4) t "l(1,1)" w l, '' u (($1==1 && (floor($2)%5 == 0)) ? ($3):1/0):($5) t "l(1,2)" w l, '' u (($1==1 && (floor($2)%5 == 0)) ? ($3):1/0):($4 +$5) t"l(1,.)" w l
+set out; unset label;
+
+#
+#
+# Survival functions in state j and all livestates from state i by final state j: 'lij' files, cov=1 state=2
+#
+
+set out "biaspar/LIJT_biaspar_2-1-1.svg"
+set label "Alive state 2 ()" at graph 0.98,0.5 center rotate font "Helvetica,12"
+set xlabel "Age"
+set ylabel "Probability to be alive"
+set ter svg size 640, 480
+unset log y
+plot [70:95] "biaspar/PIJ_biaspar.txt" u (($1==1 && (floor($2)%5 == 0)) ? ($3):1/0):($7) t "l(2,1)" w l, '' u (($1==1 && (floor($2)%5 == 0)) ? ($3):1/0):($8) t "l(2,2)" w l, '' u (($1==1 && (floor($2)%5 == 0)) ? ($3):1/0):($7 +$8) t"l(2,.)" w l
+set out; unset label;
+
+#
+#
+#CV preval stable (forward): 'pij' files, covariatecombination#=1 state=1
+#
+
+set out "biaspar/P_biaspar_1-1-1.svg"
+set label "Alive state 1 ()" at graph 0.98,0.5 center rotate font "Helvetica,12"
set xlabel "Age"
set ylabel "Probability"
-set ter png small size 320, 240
+set ter svg size 640, 480
unset log y
-plot [70:95] "biaspar/pijrbiaspar.txt" u ($1==1 ? ($3):1/0):($4/($4+$5)) t "prev(1,1)" w l, '' u ($1==1 ? ($3):1/0):($7/($7+$8)) t "prev(2,1)" w l
+plot [70:95] "biaspar/PIJ_biaspar.txt" u ($1==1 ? ($3):1/0):($4/($4+$5)) t "prev(1,1)" w l, '' u ($1==1 ? ($3):1/0):($7/($7+$8)) t "prev(2,1)" w l
+set out; unset label;
#
#
-#CV preval stable (period): 'pij' files, cov=1 state=2
-set out "biaspar/pbiaspar2_1.png"
+#CV preval stable (forward): 'pij' files, covariatecombination#=1 state=2
+#
+
+set out "biaspar/P_biaspar_2-1-1.svg"
+set label "Alive state 2 ()" at graph 0.98,0.5 center rotate font "Helvetica,12"
set xlabel "Age"
set ylabel "Probability"
-set ter png small size 320, 240
+set ter svg size 640, 480
+unset log y
+plot [70:95] "biaspar/PIJ_biaspar.txt" u ($1==1 ? ($3):1/0):($5/($4+$5)) t "prev(1,2)" w l, '' u ($1==1 ? ($3):1/0):($8/($7+$8)) t "prev(2,2)" w l
+set out; unset label;
+
+#
+#
+#Projection of prevalence to forward stable prevalence (period): 'PROJ_' files, covariatecombination#=1 state=1
+#
+# hpijx=probability over h years, hp.jx is weighted by observed prev
+
+set out "biaspar/PROJ_biaspar_1-1-1.svg"
+set label "Alive state 1 ()" at graph 0.98,0.5 center rotate font "Helvetica,12"
+set xlabel "Age"
+set ylabel "Prevalence"
+set ter svg size 640, 480
+unset log y
+plot [70:95] "biaspar/F_biaspar.txt" u 2:( $5/(1.-$11)):1 t 'pw.1' with line lc variable ,\
+ '' u 2:( (($1-$2) == 1919 ) ? $5/(1.-$11) : 1/0):1 with labels center not
+set out; unset label;
+
+#
+#
+#Projection of prevalence to forward stable prevalence (period): 'PROJ_' files, covariatecombination#=1 state=2
+#
+# hpijx=probability over h years, hp.jx is weighted by observed prev
+
+set out "biaspar/PROJ_biaspar_2-1-1.svg"
+set label "Alive state 2 ()" at graph 0.98,0.5 center rotate font "Helvetica,12"
+set xlabel "Age"
+set ylabel "Prevalence"
+set ter svg size 640, 480
unset log y
-plot [70:95] "biaspar/pijrbiaspar.txt" u ($1==1 ? ($3):1/0):($5/($4+$5)) t "prev(1,2)" w l, '' u ($1==1 ? ($3):1/0):($8/($7+$8)) t "prev(2,2)" w l
+plot [70:95] "biaspar/F_biaspar.txt" u 2:( $8/(1.-$11)):1 t 'pw.2' with line lc variable ,\
+ '' u 2:( (($1-$2) == 1919 ) ? $8/(1.-$11) : 1/0):1 with labels center not
+set out; unset label;
##############
-#MLE estimated parameters
+#9eme MLE estimated parameters
#############
# initial state 1
# current state 2
-p1=-12.202711; p2=0.091764;
+p1=-12.317578; p2=0.093112;
# current state 3
-p3=-10.299238; p4=0.056626;
+p3=-10.163170; p4=0.054838;
# initial state 2
# current state 1
-p5=-1.545619; p6=-0.035395;
+p5=-1.668569; p6=-0.033906;
# current state 3
-p7=-5.565963; p8=0.016951;
+p7=-5.686920; p8=0.018135;
##############
#
##############
-#Graphics of of probabilities or incidences
+#10eme Graphics of probabilities or incidences
#############
# logi(p12/p11)=a12+b12*age+c12age*age+d12*V1+e12*V1*age
# logi(p12/p11)=p1 +p2*age +p3*age*age+ p4*V1+ p5*V1*age
# +exp(a13+b13*age+c13age*age+d13*V1+e13*V1*age))
# +exp(a14+b14*age+c14age*age+d14*V1+e14*V1*age)+...)
#
-# ng=1
-# jk=1 to 2^0=1
-# jk=1
+#Number of graphics: first is logit, 2nd is probabilities, third is incidences per year
+#model=
+# Type of graphic ng=1
+# k1=1 to 2^0=1
-set out "biaspar/pebiaspar1_1.png"
-set title "Probability"
+# Combination of dummy k1=1 which is
+#
+
+set out "biaspar/PE_biaspar_1-1-1.svg"
+set key outside
+set title "()" font "Helvetica,12"
+
+set ter svg size 640, 480
+set ylabel "Value of the logit of the model"
+
+unset log y
+plot [70:95] p1+p2*x w l lw 2 lt (3*1+2)%3+1 dt 1 t "logit(p12)" , p3+p4*x w l lw 2 lt (3*1+3)%3+1 dt 1 t "logit(p13)" , p5+p6*x w l lw 2 lt (3*2+1)%3+1 dt 2 t "logit(p21)" , p7+p8*x w l lw 2 lt (3*2+3)%3+1 dt 2 t "logit(p23)"
+ set out; unset title;set key default;
+#Number of graphics: first is logit, 2nd is probabilities, third is incidences per year
+#model=
+# Type of graphic ng=2
+# k1=1 to 2^0=1
+
+
+# Combination of dummy k1=1 which is
+#
+
+set out "biaspar/PE_biaspar_1-2-1.svg"
+set key outside
+set title "()" font "Helvetica,12"
+
+set ter svg size 640, 480
+set ylabel "Probability"
-set ter png small size 320, 240
set log y
-plot [70:95] exp(p1+p2*x)/(1+exp(p1+p2*x)+exp(p3+p4*x)) t "p12" , exp(p3+p4*x)/(1+exp(p1+p2*x)+exp(p3+p4*x)) t "p13" , exp(p5+p6*x)/(1+exp(p5+p6*x)+exp(p7+p8*x)) t "p21" , exp(p7+p8*x)/(1+exp(p5+p6*x)+exp(p7+p8*x)) t "p23" # ng=2
-# jk=1 to 2^0=1
-# jk=1
+plot [70:95] (1.)/(1+exp(p1+p2*x)+exp(p3+p4*x)) w l lw 2 lt (3*1+1)%3+1 dt 1 t "p11" , exp(p1+p2*x)/(1+exp(p1+p2*x)+exp(p3+p4*x)) w l lw 2 lt (3*1+2)%3+1 dt 1 t "p12" , exp(p3+p4*x)/(1+exp(p1+p2*x)+exp(p3+p4*x)) w l lw 2 lt (3*1+3)%3+1 dt 1 t "p13" , exp(p5+p6*x)/(1+exp(p5+p6*x)+exp(p7+p8*x)) w l lw 2 lt (3*2+1)%3+1 dt 2 t "p21" , (1.)/(1+exp(p5+p6*x)+exp(p7+p8*x)) w l lw 2 lt (3*2+2)%3+1 dt 2 t "p22" , exp(p7+p8*x)/(1+exp(p5+p6*x)+exp(p7+p8*x)) w l lw 2 lt (3*2+3)%3+1 dt 2 t "p23"
+ set out; unset title;set key default;
+#Number of graphics: first is logit, 2nd is probabilities, third is incidences per year
+#model=
+# Type of graphic ng=3
+# k1=1 to 2^0=1
+
-set out "biaspar/pebiaspar1_2.png"
+# Combination of dummy k1=1 which is
+#
+
+set out "biaspar/PE_biaspar_1-3-1.svg"
+set key outside
+set title "()" font "Helvetica,12"
+set ter svg size 640, 480
set ylabel "Quasi-incidence per year"
-set ter png small size 320, 240
set log y
-plot [70:95] 12.000000*exp(p1+p2*x)/(1+exp(p1+p2*x)+exp(p3+p4*x)) t "p12" , 12.000000*exp(p3+p4*x)/(1+exp(p1+p2*x)+exp(p3+p4*x)) t "p13" , 12.000000*exp(p5+p6*x)/(1+exp(p5+p6*x)+exp(p7+p8*x)) t "p21" , 12.000000*exp(p7+p8*x)/(1+exp(p5+p6*x)+exp(p7+p8*x)) t "p23"
+plot [70:95] (1.)/(1+exp(p1+p2*x)+exp(p3+p4*x)) w l lw 2 lt (3*1+1)%3+1 dt 1 t "i11" , 12.000000*exp(p1+p2*x)/(1+exp(p1+p2*x)+exp(p3+p4*x)) w l lw 2 lt (3*1+2)%3+1 dt 1 t "i12" , 12.000000*exp(p3+p4*x)/(1+exp(p1+p2*x)+exp(p3+p4*x)) w l lw 2 lt (3*1+3)%3+1 dt 1 t "i13" , 12.000000*exp(p5+p6*x)/(1+exp(p5+p6*x)+exp(p7+p8*x)) w l lw 2 lt (3*2+1)%3+1 dt 2 t "i21" , (1.)/(1+exp(p5+p6*x)+exp(p7+p8*x)) w l lw 2 lt (3*2+2)%3+1 dt 2 t "i22" , 12.000000*exp(p7+p8*x)/(1+exp(p5+p6*x)+exp(p7+p8*x)) w l lw 2 lt (3*2+3)%3+1 dt 2 t "i23"
+ set out; unset title;set key default;
+
# Routine varprob
+# Ellipsoids of confidence
+#
+
set parametric;unset label
set log y;set log x; set xlabel "p13 (year-1)";set ylabel "p12 (year-1)"
-set ter png small size 320, 240
-set out "biaspar/varpijgrbiaspar113-12.png"
-set label "65" at 1.597e-02, 2.337e-02 center
-# Age 65, p13 - p12
-plot [-pi:pi] 1.597e-02+ 2.000*( 7.325e-01* 2.084e-03*cos(t)+ 6.807e-01* 1.465e-03*sin(t)), 2.337e-02 +2.000*( -6.807e-01* 2.084e-03*cos(t)+ 7.325e-01* 1.465e-03*sin(t)) not
+set ter svg size 640, 480
+set out "biaspar/VARPIJGR_biaspar_113-12.svg"
+set label "70" at 2.140e-02, 3.617e-02 center
# Age 70, p13 - p12
-set label "70" at 2.117e-02, 3.691e-02 center
-replot 2.117e-02+ 2.000*( 5.240e-01* 2.315e-03*cos(t)+ 8.517e-01* 1.521e-03*sin(t)), 3.691e-02 +2.000*( -8.517e-01* 2.315e-03*cos(t)+ 5.240e-01* 1.521e-03*sin(t)) not
+plot [-pi:pi] 2.140e-02+ 2.000*( 5.514e-01* 2.264e-03*cos(t)+ 8.342e-01* 1.523e-03*sin(t)), 3.617e-02 +2.000*( -8.342e-01* 2.264e-03*cos(t)+ 5.514e-01* 1.523e-03*sin(t)) not
# Age 75, p13 - p12
-set label "75" at 2.803e-02, 5.827e-02 center
-replot 2.803e-02+ 2.000*( 3.929e-01* 2.501e-03*cos(t)+ 9.196e-01* 1.403e-03*sin(t)), 5.827e-02 +2.000*( -9.196e-01* 2.501e-03*cos(t)+ 3.929e-01* 1.403e-03*sin(t)) not
+set label "75" at 2.808e-02, 5.747e-02 center
+replot 2.808e-02+ 2.000*( 3.943e-01* 2.448e-03*cos(t)+ 9.190e-01* 1.398e-03*sin(t)), 5.747e-02 +2.000*( -9.190e-01* 2.448e-03*cos(t)+ 3.943e-01* 1.398e-03*sin(t)) not
# Age 80, p13 - p12
-set label "80" at 3.707e-02, 9.186e-02 center
-replot 3.707e-02+ 2.000*( 3.828e-01* 2.718e-03*cos(t)+ 9.238e-01* 1.387e-03*sin(t)), 9.186e-02 +2.000*( -9.238e-01* 2.718e-03*cos(t)+ 3.828e-01* 1.387e-03*sin(t)) not
+set label "80" at 3.681e-02, 9.123e-02 center
+replot 3.681e-02+ 2.000*( 3.741e-01* 2.676e-03*cos(t)+ 9.274e-01* 1.373e-03*sin(t)), 9.123e-02 +2.000*( -9.274e-01* 2.676e-03*cos(t)+ 3.741e-01* 1.373e-03*sin(t)) not
# Age 85, p13 - p12
-set label "85" at 4.893e-02, 1.445e-01 center
-replot 4.893e-02+ 2.000*( 4.026e-01* 4.082e-03*cos(t)+ 9.154e-01* 2.265e-03*sin(t)), 1.445e-01 +2.000*( -9.154e-01* 4.082e-03*cos(t)+ 4.026e-01* 2.265e-03*sin(t)) not
+set label "85" at 4.816e-02, 1.445e-01 center
+replot 4.816e-02+ 2.000*( 3.906e-01* 4.057e-03*cos(t)+ 9.206e-01* 2.238e-03*sin(t)), 1.445e-01 +2.000*( -9.206e-01* 4.057e-03*cos(t)+ 3.906e-01* 2.238e-03*sin(t)) not
# Age 90, p13 - p12
-set label "90" at 6.441e-02, 2.268e-01 center
-replot 6.441e-02+ 2.000*( 2.874e-01* 8.787e-03*cos(t)+ 9.578e-01* 4.578e-03*sin(t)), 2.268e-01 +2.000*( -9.578e-01* 8.787e-03*cos(t)+ 2.874e-01* 4.578e-03*sin(t)) not
+set label "90" at 6.283e-02, 2.283e-01 center
+replot 6.283e-02+ 2.000*( 2.754e-01* 8.795e-03*cos(t)+ 9.613e-01* 4.488e-03*sin(t)), 2.283e-01 +2.000*( -9.613e-01* 8.795e-03*cos(t)+ 2.754e-01* 4.488e-03*sin(t)) not
# Age 95, p13 - p12
-set label "95" at 8.441e-02, 3.543e-01 center
-replot 8.441e-02+ 2.000*( 2.019e-01* 1.921e-02*cos(t)+ 9.794e-01* 8.453e-03*sin(t)), 3.543e-01 +2.000*( -9.794e-01* 1.921e-02*cos(t)+ 2.019e-01* 8.453e-03*sin(t)) not
-set out "biaspar/varpijgrbiaspar113-12.png";replot;
+set label "95" at 8.159e-02, 3.590e-01 center
+replot 8.159e-02+ 2.000*( 1.909e-01* 1.934e-02*cos(t)+ 9.816e-01* 8.213e-03*sin(t)), 3.590e-01 +2.000*( -9.816e-01* 1.934e-02*cos(t)+ 1.909e-01* 8.213e-03*sin(t)) not
+set out;
+set out "biaspar/VARPIJGR_biaspar_113-12.svg";replot;set out;
+# Ellipsoids of confidence
+#
+
set parametric;unset label
set log y;set log x; set xlabel "p21 (year-1)";set ylabel "p12 (year-1)"
-set ter png small size 320, 240
-set out "biaspar/varpijgrbiaspar121-12.png"
-set label "65" at 2.481e-01, 2.337e-02 center
-# Age 65, p21 - p12
-plot [-pi:pi] 2.481e-01+ 2.000*( 9.999e-01* 3.884e-02*cos(t)+ -1.528e-02* 1.677e-03*sin(t)), 2.337e-02 +2.000*( 1.528e-02* 3.884e-02*cos(t)+ 9.999e-01* 1.677e-03*sin(t)) not
+set ter svg size 640, 480
+set out "biaspar/VARPIJGR_biaspar_121-12.svg"
+set label "70" at 2.047e-01, 3.617e-02 center
# Age 70, p21 - p12
-set label "70" at 2.084e-01, 3.691e-02 center
-replot 2.084e-01+ 2.000*( 9.996e-01* 2.493e-02*cos(t)+ -2.848e-02* 2.005e-03*sin(t)), 3.691e-02 +2.000*( 2.848e-02* 2.493e-02*cos(t)+ 9.996e-01* 2.005e-03*sin(t)) not
+plot [-pi:pi] 2.047e-01+ 2.000*( 9.996e-01* 2.493e-02*cos(t)+ -2.854e-02* 1.942e-03*sin(t)), 3.617e-02 +2.000*( 2.854e-02* 2.493e-02*cos(t)+ 9.996e-01* 1.942e-03*sin(t)) not
# Age 75, p21 - p12
-set label "75" at 1.749e-01, 5.827e-02 center
-replot 1.749e-01+ 2.000*( 9.986e-01* 1.487e-02*cos(t)+ -5.306e-02* 2.233e-03*sin(t)), 5.827e-02 +2.000*( 5.306e-02* 1.487e-02*cos(t)+ 9.986e-01* 2.233e-03*sin(t)) not
+set label "75" at 1.730e-01, 5.747e-02 center
+replot 1.730e-01+ 2.000*( 9.986e-01* 1.494e-02*cos(t)+ -5.266e-02* 2.181e-03*sin(t)), 5.747e-02 +2.000*( 5.266e-02* 1.494e-02*cos(t)+ 9.986e-01* 2.181e-03*sin(t)) not
# Age 80, p21 - p12
-set label "80" at 1.467e-01, 9.186e-02 center
-replot 1.467e-01+ 2.000*( 9.949e-01* 8.452e-03*cos(t)+ -1.006e-01* 2.434e-03*sin(t)), 9.186e-02 +2.000*( 1.006e-01* 8.452e-03*cos(t)+ 9.949e-01* 2.434e-03*sin(t)) not
+set label "80" at 1.462e-01, 9.123e-02 center
+replot 1.462e-01+ 2.000*( 9.951e-01* 8.486e-03*cos(t)+ -9.892e-02* 2.403e-03*sin(t)), 9.123e-02 +2.000*( 9.892e-02* 8.486e-03*cos(t)+ 9.951e-01* 2.403e-03*sin(t)) not
# Age 85, p21 - p12
-set label "85" at 1.230e-01, 1.445e-01 center
-replot 1.230e-01+ 2.000*( 9.691e-01* 6.397e-03*cos(t)+ -2.468e-01* 3.619e-03*sin(t)), 1.445e-01 +2.000*( 2.468e-01* 6.397e-03*cos(t)+ 9.691e-01* 3.619e-03*sin(t)) not
+set label "85" at 1.235e-01, 1.445e-01 center
+replot 1.235e-01+ 2.000*( 9.684e-01* 6.423e-03*cos(t)+ -2.493e-01* 3.599e-03*sin(t)), 1.445e-01 +2.000*( 2.493e-01* 6.423e-03*cos(t)+ 9.684e-01* 3.599e-03*sin(t)) not
# Age 90, p21 - p12
-set label "90" at 1.031e-01, 2.268e-01 center
-replot 1.031e-01+ 2.000*( 5.313e-01* 9.166e-03*cos(t)+ -8.472e-01* 6.591e-03*sin(t)), 2.268e-01 +2.000*( 8.472e-01* 9.166e-03*cos(t)+ 5.313e-01* 6.591e-03*sin(t)) not
+set label "90" at 1.042e-01, 2.283e-01 center
+replot 1.042e-01+ 2.000*( 5.573e-01* 9.289e-03*cos(t)+ -8.303e-01* 6.600e-03*sin(t)), 2.283e-01 +2.000*( 8.303e-01* 9.289e-03*cos(t)+ 5.573e-01* 6.600e-03*sin(t)) not
# Age 95, p21 - p12
-set label "95" at 8.634e-02, 3.543e-01 center
-replot 8.634e-02+ 2.000*( 1.792e-01* 1.914e-02*cos(t)+ -9.838e-01* 8.537e-03*sin(t)), 3.543e-01 +2.000*( 9.838e-01* 1.914e-02*cos(t)+ 1.792e-01* 8.537e-03*sin(t)) not
-set out "biaspar/varpijgrbiaspar121-12.png";replot;
+set label "95" at 8.796e-02, 3.590e-01 center
+replot 8.796e-02+ 2.000*( 1.959e-01* 1.934e-02*cos(t)+ -9.806e-01* 8.711e-03*sin(t)), 3.590e-01 +2.000*( 9.806e-01* 1.934e-02*cos(t)+ 1.959e-01* 8.711e-03*sin(t)) not
+set out;
+set out "biaspar/VARPIJGR_biaspar_121-12.svg";replot;set out;
+# Ellipsoids of confidence
+#
+
set parametric;unset label
set log y;set log x; set xlabel "p23 (year-1)";set ylabel "p12 (year-1)"
-set ter png small size 320, 240
-set out "biaspar/varpijgrbiaspar123-12.png"
-set label "65" at 1.338e-01, 2.337e-02 center
-# Age 65, p23 - p12
-plot [-pi:pi] 1.338e-01+ 2.000*( 9.993e-01* 1.768e-02*cos(t)+ -3.696e-02* 1.655e-03*sin(t)), 2.337e-02 +2.000*( 3.696e-02* 1.768e-02*cos(t)+ 9.993e-01* 1.655e-03*sin(t)) not
+set ter svg size 640, 480
+set out "biaspar/VARPIJGR_biaspar_123-12.svg"
+set label "70" at 1.406e-01, 3.617e-02 center
# Age 70, p23 - p12
-set label "70" at 1.460e-01, 3.691e-02 center
-replot 1.460e-01+ 2.000*( 9.985e-01* 1.541e-02*cos(t)+ -5.441e-02* 1.957e-03*sin(t)), 3.691e-02 +2.000*( 5.441e-02* 1.541e-02*cos(t)+ 9.985e-01* 1.957e-03*sin(t)) not
+plot [-pi:pi] 1.406e-01+ 2.000*( 9.987e-01* 1.493e-02*cos(t)+ -5.178e-02* 1.920e-03*sin(t)), 3.617e-02 +2.000*( 5.178e-02* 1.493e-02*cos(t)+ 9.987e-01* 1.920e-03*sin(t)) not
# Age 75, p23 - p12
-set label "75" at 1.591e-01, 5.827e-02 center
-replot 1.591e-01+ 2.000*( 9.967e-01* 1.275e-02*cos(t)+ -8.076e-02* 2.137e-03*sin(t)), 5.827e-02 +2.000*( 8.076e-02* 1.275e-02*cos(t)+ 9.967e-01* 2.137e-03*sin(t)) not
+set label "75" at 1.542e-01, 5.747e-02 center
+replot 1.542e-01+ 2.000*( 9.970e-01* 1.241e-02*cos(t)+ -7.775e-02* 2.112e-03*sin(t)), 5.747e-02 +2.000*( 7.775e-02* 1.241e-02*cos(t)+ 9.970e-01* 2.112e-03*sin(t)) not
# Age 80, p23 - p12
-set label "80" at 1.734e-01, 9.186e-02 center
-replot 1.734e-01+ 2.000*( 9.924e-01* 9.884e-03*cos(t)+ -1.231e-01* 2.277e-03*sin(t)), 9.186e-02 +2.000*( 1.231e-01* 9.884e-03*cos(t)+ 9.924e-01* 2.277e-03*sin(t)) not
+set label "80" at 1.690e-01, 9.123e-02 center
+replot 1.690e-01+ 2.000*( 9.927e-01* 9.663e-03*cos(t)+ -1.210e-01* 2.265e-03*sin(t)), 9.123e-02 +2.000*( 1.210e-01* 9.663e-03*cos(t)+ 9.927e-01* 2.265e-03*sin(t)) not
# Age 85, p23 - p12
-set label "85" at 1.889e-01, 1.445e-01 center
-replot 1.889e-01+ 2.000*( 9.741e-01* 7.696e-03*cos(t)+ -2.260e-01* 3.522e-03*sin(t)), 1.445e-01 +2.000*( 2.260e-01* 7.696e-03*cos(t)+ 9.741e-01* 3.522e-03*sin(t)) not
+set label "85" at 1.851e-01, 1.445e-01 center
+replot 1.851e-01+ 2.000*( 9.733e-01* 7.552e-03*cos(t)+ -2.295e-01* 3.515e-03*sin(t)), 1.445e-01 +2.000*( 2.295e-01* 7.552e-03*cos(t)+ 9.733e-01* 3.515e-03*sin(t)) not
# Age 90, p23 - p12
-set label "90" at 2.057e-01, 2.268e-01 center
-replot 2.057e-01+ 2.000*( 6.276e-01* 9.309e-03*cos(t)+ -7.786e-01* 7.132e-03*sin(t)), 2.268e-01 +2.000*( 7.786e-01* 9.309e-03*cos(t)+ 6.276e-01* 7.132e-03*sin(t)) not
+set label "90" at 2.028e-01, 2.283e-01 center
+replot 2.028e-01+ 2.000*( 6.011e-01* 9.259e-03*cos(t)+ -7.992e-01* 7.110e-03*sin(t)), 2.283e-01 +2.000*( 7.992e-01* 9.259e-03*cos(t)+ 6.011e-01* 7.110e-03*sin(t)) not
# Age 95, p23 - p12
-set label "95" at 2.238e-01, 3.543e-01 center
-replot 2.238e-01+ 2.000*( 2.443e-01* 1.925e-02*cos(t)+ -9.697e-01* 1.172e-02*sin(t)), 3.543e-01 +2.000*( 9.697e-01* 1.925e-02*cos(t)+ 2.443e-01* 1.172e-02*sin(t)) not
-set out "biaspar/varpijgrbiaspar123-12.png";replot;
+set label "95" at 2.219e-01, 3.590e-01 center
+replot 2.219e-01+ 2.000*( 2.273e-01* 1.937e-02*cos(t)+ -9.738e-01* 1.173e-02*sin(t)), 3.590e-01 +2.000*( 9.738e-01* 1.937e-02*cos(t)+ 2.273e-01* 1.173e-02*sin(t)) not
+set out;
+set out "biaspar/VARPIJGR_biaspar_123-12.svg";replot;set out;
+# Ellipsoids of confidence
+#
+
set parametric;unset label
set log y;set log x; set xlabel "p21 (year-1)";set ylabel "p13 (year-1)"
-set ter png small size 320, 240
-set out "biaspar/varpijgrbiaspar121-13.png"
-set label "65" at 2.481e-01, 1.597e-02 center
-# Age 65, p21 - p13
-plot [-pi:pi] 2.481e-01+ 2.000*( 1.000e+00* 3.883e-02*cos(t)+ -2.210e-03* 1.821e-03*sin(t)), 1.597e-02 +2.000*( 2.210e-03* 3.883e-02*cos(t)+ 1.000e+00* 1.821e-03*sin(t)) not
+set ter svg size 640, 480
+set out "biaspar/VARPIJGR_biaspar_121-13.svg"
+set label "70" at 2.047e-01, 2.140e-02 center
# Age 70, p21 - p13
-set label "70" at 2.084e-01, 2.117e-02 center
-replot 2.084e-01+ 2.000*( 1.000e+00* 2.492e-02*cos(t)+ -2.683e-03* 1.773e-03*sin(t)), 2.117e-02 +2.000*( 2.683e-03* 2.492e-02*cos(t)+ 1.000e+00* 1.773e-03*sin(t)) not
+plot [-pi:pi] 2.047e-01+ 2.000*( 1.000e+00* 2.492e-02*cos(t)+ -2.081e-03* 1.781e-03*sin(t)), 2.140e-02 +2.000*( 2.081e-03* 2.492e-02*cos(t)+ 1.000e+00* 1.781e-03*sin(t)) not
# Age 75, p21 - p13
-set label "75" at 1.749e-01, 2.803e-02 center
-replot 1.749e-01+ 2.000*( 1.000e+00* 1.485e-02*cos(t)+ -2.605e-03* 1.622e-03*sin(t)), 2.803e-02 +2.000*( 2.605e-03* 1.485e-02*cos(t)+ 1.000e+00* 1.622e-03*sin(t)) not
+set label "75" at 1.730e-01, 2.808e-02 center
+replot 1.730e-01+ 2.000*( 1.000e+00* 1.492e-02*cos(t)+ -2.295e-03* 1.606e-03*sin(t)), 2.808e-02 +2.000*( 2.295e-03* 1.492e-02*cos(t)+ 1.000e+00* 1.606e-03*sin(t)) not
# Age 80, p21 - p13
-set label "80" at 1.467e-01, 3.707e-02 center
-replot 1.467e-01+ 2.000*( 1.000e+00* 8.412e-03*cos(t)+ -4.929e-03* 1.650e-03*sin(t)), 3.707e-02 +2.000*( 4.929e-03* 8.412e-03*cos(t)+ 1.000e+00* 1.650e-03*sin(t)) not
+set label "80" at 1.462e-01, 3.681e-02 center
+replot 1.462e-01+ 2.000*( 1.000e+00* 8.448e-03*cos(t)+ -5.527e-03* 1.619e-03*sin(t)), 3.681e-02 +2.000*( 5.527e-03* 8.448e-03*cos(t)+ 1.000e+00* 1.619e-03*sin(t)) not
# Age 85, p21 - p13
-set label "85" at 1.230e-01, 4.893e-02 center
-replot 1.230e-01+ 2.000*( 9.992e-01* 6.267e-03*cos(t)+ -3.981e-02* 2.636e-03*sin(t)), 4.893e-02 +2.000*( 3.981e-02* 6.267e-03*cos(t)+ 9.992e-01* 2.636e-03*sin(t)) not
+set label "85" at 1.235e-01, 4.816e-02 center
+replot 1.235e-01+ 2.000*( 9.994e-01* 6.288e-03*cos(t)+ -3.508e-02* 2.591e-03*sin(t)), 4.816e-02 +2.000*( 3.508e-02* 6.288e-03*cos(t)+ 9.994e-01* 2.591e-03*sin(t)) not
# Age 90, p21 - p13
-set label "90" at 1.031e-01, 6.441e-02 center
-replot 1.031e-01+ 2.000*( 9.921e-01* 7.441e-03*cos(t)+ -1.255e-01* 5.013e-03*sin(t)), 6.441e-02 +2.000*( 1.255e-01* 7.441e-03*cos(t)+ 9.921e-01* 5.013e-03*sin(t)) not
+set label "90" at 1.042e-01, 6.283e-02 center
+replot 1.042e-01+ 2.000*( 9.956e-01* 7.558e-03*cos(t)+ -9.337e-02* 4.919e-03*sin(t)), 6.283e-02 +2.000*( 9.337e-02* 7.558e-03*cos(t)+ 9.956e-01* 4.919e-03*sin(t)) not
# Age 95, p21 - p13
-set label "95" at 8.634e-02, 8.441e-02 center
-replot 8.634e-02+ 2.000*( 6.785e-01* 9.541e-03*cos(t)+ -7.346e-01* 8.651e-03*sin(t)), 8.441e-02 +2.000*( 7.346e-01* 9.541e-03*cos(t)+ 6.785e-01* 8.651e-03*sin(t)) not
-set out "biaspar/varpijgrbiaspar121-13.png";replot;
+set label "95" at 8.796e-02, 8.159e-02 center
+replot 8.796e-02+ 2.000*( 8.858e-01* 9.520e-03*cos(t)+ -4.641e-01* 8.679e-03*sin(t)), 8.159e-02 +2.000*( 4.641e-01* 9.520e-03*cos(t)+ 8.858e-01* 8.679e-03*sin(t)) not
+set out;
+set out "biaspar/VARPIJGR_biaspar_121-13.svg";replot;set out;
+# Ellipsoids of confidence
+#
+
set parametric;unset label
set log y;set log x; set xlabel "p23 (year-1)";set ylabel "p13 (year-1)"
-set ter png small size 320, 240
-set out "biaspar/varpijgrbiaspar123-13.png"
-set label "65" at 1.338e-01, 1.597e-02 center
-# Age 65, p23 - p13
-plot [-pi:pi] 1.338e-01+ 2.000*( 9.991e-01* 1.769e-02*cos(t)+ 4.256e-02* 1.662e-03*sin(t)), 1.597e-02 +2.000*( -4.256e-02* 1.769e-02*cos(t)+ 9.991e-01* 1.662e-03*sin(t)) not
+set ter svg size 640, 480
+set out "biaspar/VARPIJGR_biaspar_123-13.svg"
+set label "70" at 1.406e-01, 2.140e-02 center
# Age 70, p23 - p13
-set label "70" at 1.460e-01, 2.117e-02 center
-replot 1.460e-01+ 2.000*( 9.986e-01* 1.541e-02*cos(t)+ 5.267e-02* 1.580e-03*sin(t)), 2.117e-02 +2.000*( -5.267e-02* 1.541e-02*cos(t)+ 9.986e-01* 1.580e-03*sin(t)) not
+plot [-pi:pi] 1.406e-01+ 2.000*( 9.985e-01* 1.494e-02*cos(t)+ 5.556e-02* 1.579e-03*sin(t)), 2.140e-02 +2.000*( -5.556e-02* 1.494e-02*cos(t)+ 9.985e-01* 1.579e-03*sin(t)) not
# Age 75, p23 - p13
-set label "75" at 1.591e-01, 2.803e-02 center
-replot 1.591e-01+ 2.000*( 9.977e-01* 1.273e-02*cos(t)+ 6.800e-02* 1.375e-03*sin(t)), 2.803e-02 +2.000*( -6.800e-02* 1.273e-02*cos(t)+ 9.977e-01* 1.375e-03*sin(t)) not
+set label "75" at 1.542e-01, 2.808e-02 center
+replot 1.542e-01+ 2.000*( 9.976e-01* 1.240e-02*cos(t)+ 6.965e-02* 1.358e-03*sin(t)), 2.808e-02 +2.000*( -6.965e-02* 1.240e-02*cos(t)+ 9.976e-01* 1.358e-03*sin(t)) not
# Age 80, p23 - p13
-set label "80" at 1.734e-01, 3.707e-02 center
-replot 1.734e-01+ 2.000*( 9.950e-01* 9.861e-03*cos(t)+ 9.978e-02* 1.332e-03*sin(t)), 3.707e-02 +2.000*( -9.978e-02* 9.861e-03*cos(t)+ 9.950e-01* 1.332e-03*sin(t)) not
+set label "80" at 1.690e-01, 3.681e-02 center
+replot 1.690e-01+ 2.000*( 9.951e-01* 9.643e-03*cos(t)+ 9.913e-02* 1.314e-03*sin(t)), 3.681e-02 +2.000*( -9.913e-02* 9.643e-03*cos(t)+ 9.951e-01* 1.314e-03*sin(t)) not
# Age 85, p23 - p13
-set label "85" at 1.889e-01, 4.893e-02 center
-replot 1.889e-01+ 2.000*( 9.810e-01* 7.672e-03*cos(t)+ 1.938e-01* 2.231e-03*sin(t)), 4.893e-02 +2.000*( -1.938e-01* 7.672e-03*cos(t)+ 9.810e-01* 2.231e-03*sin(t)) not
+set label "85" at 1.851e-01, 4.816e-02 center
+replot 1.851e-01+ 2.000*( 9.815e-01* 7.522e-03*cos(t)+ 1.914e-01* 2.205e-03*sin(t)), 4.816e-02 +2.000*( -1.914e-01* 7.522e-03*cos(t)+ 9.815e-01* 2.205e-03*sin(t)) not
# Age 90, p23 - p13
-set label "90" at 2.057e-01, 6.441e-02 center
-replot 2.057e-01+ 2.000*( 9.336e-01* 8.471e-03*cos(t)+ 3.582e-01* 4.337e-03*sin(t)), 6.441e-02 +2.000*( -3.582e-01* 8.471e-03*cos(t)+ 9.336e-01* 4.337e-03*sin(t)) not
+set label "90" at 2.028e-01, 6.283e-02 center
+replot 2.028e-01+ 2.000*( 9.340e-01* 8.361e-03*cos(t)+ 3.573e-01* 4.223e-03*sin(t)), 6.283e-02 +2.000*( -3.573e-01* 8.361e-03*cos(t)+ 9.340e-01* 4.223e-03*sin(t)) not
# Age 95, p23 - p13
-set label "95" at 2.238e-01, 8.441e-02 center
-replot 2.238e-01+ 2.000*( 9.049e-01* 1.305e-02*cos(t)+ 4.257e-01* 8.023e-03*sin(t)), 8.441e-02 +2.000*( -4.257e-01* 1.305e-02*cos(t)+ 9.049e-01* 8.023e-03*sin(t)) not
-set out "biaspar/varpijgrbiaspar123-13.png";replot;
+set label "95" at 2.219e-01, 8.159e-02 center
+replot 2.219e-01+ 2.000*( 9.094e-01* 1.298e-02*cos(t)+ 4.159e-01* 7.733e-03*sin(t)), 8.159e-02 +2.000*( -4.159e-01* 1.298e-02*cos(t)+ 9.094e-01* 7.733e-03*sin(t)) not
+set out;
+set out "biaspar/VARPIJGR_biaspar_123-13.svg";replot;set out;
+# Ellipsoids of confidence
+#
+
set parametric;unset label
set log y;set log x; set xlabel "p23 (year-1)";set ylabel "p21 (year-1)"
-set ter png small size 320, 240
-set out "biaspar/varpijgrbiaspar123-21.png"
-set label "65" at 1.338e-01, 2.481e-01 center
-# Age 65, p23 - p21
-plot [-pi:pi] 1.338e-01+ 2.000*( 8.647e-03* 3.884e-02*cos(t)+ 1.000e+00* 1.767e-02*sin(t)), 2.481e-01 +2.000*( -1.000e+00* 3.884e-02*cos(t)+ 8.647e-03* 1.767e-02*sin(t)) not
+set ter svg size 640, 480
+set out "biaspar/VARPIJGR_biaspar_123-21.svg"
+set label "70" at 1.406e-01, 2.047e-01 center
# Age 70, p23 - p21
-set label "70" at 1.460e-01, 2.084e-01 center
-replot 1.460e-01+ 2.000*( 1.255e-02* 2.493e-02*cos(t)+ 9.999e-01* 1.539e-02*sin(t)), 2.084e-01 +2.000*( -9.999e-01* 2.493e-02*cos(t)+ 1.255e-02* 1.539e-02*sin(t)) not
+plot [-pi:pi] 1.406e-01+ 2.000*( 3.099e-02* 2.492e-02*cos(t)+ 9.995e-01* 1.490e-02*sin(t)), 2.047e-01 +2.000*( -9.995e-01* 2.492e-02*cos(t)+ 3.099e-02* 1.490e-02*sin(t)) not
# Age 75, p23 - p21
-set label "75" at 1.591e-01, 1.749e-01 center
-replot 1.591e-01+ 2.000*( 4.580e-02* 1.486e-02*cos(t)+ 9.990e-01* 1.270e-02*sin(t)), 1.749e-01 +2.000*( -9.990e-01* 1.486e-02*cos(t)+ 4.580e-02* 1.270e-02*sin(t)) not
+set label "75" at 1.542e-01, 1.730e-01 center
+replot 1.542e-01+ 2.000*( 8.468e-02* 1.494e-02*cos(t)+ 9.964e-01* 1.235e-02*sin(t)), 1.730e-01 +2.000*( -9.964e-01* 1.494e-02*cos(t)+ 8.468e-02* 1.235e-02*sin(t)) not
# Age 80, p23 - p21
-set label "80" at 1.734e-01, 1.467e-01 center
-replot 1.734e-01+ 2.000*( 9.957e-01* 9.824e-03*cos(t)+ 9.255e-02* 8.399e-03*sin(t)), 1.467e-01 +2.000*( -9.255e-02* 9.824e-03*cos(t)+ 9.957e-01* 8.399e-03*sin(t)) not
+set label "80" at 1.690e-01, 1.462e-01 center
+replot 1.690e-01+ 2.000*( 9.896e-01* 9.619e-03*cos(t)+ 1.436e-01* 8.421e-03*sin(t)), 1.462e-01 +2.000*( -1.436e-01* 9.619e-03*cos(t)+ 9.896e-01* 8.421e-03*sin(t)) not
# Age 85, p23 - p21
-set label "85" at 1.889e-01, 1.230e-01 center
-replot 1.889e-01+ 2.000*( 9.875e-01* 7.569e-03*cos(t)+ 1.578e-01* 6.226e-03*sin(t)), 1.230e-01 +2.000*( -1.578e-01* 7.569e-03*cos(t)+ 9.875e-01* 6.226e-03*sin(t)) not
+set label "85" at 1.851e-01, 1.235e-01 center
+replot 1.851e-01+ 2.000*( 9.870e-01* 7.423e-03*cos(t)+ 1.606e-01* 6.252e-03*sin(t)), 1.235e-01 +2.000*( -1.606e-01* 7.423e-03*cos(t)+ 9.870e-01* 6.252e-03*sin(t)) not
# Age 90, p23 - p21
-set label "90" at 2.057e-01, 1.031e-01 center
-replot 2.057e-01+ 2.000*( 9.508e-01* 8.134e-03*cos(t)+ 3.099e-01* 7.328e-03*sin(t)), 1.031e-01 +2.000*( -3.099e-01* 8.134e-03*cos(t)+ 9.508e-01* 7.328e-03*sin(t)) not
+set label "90" at 2.028e-01, 1.042e-01 center
+replot 2.028e-01+ 2.000*( 9.136e-01* 8.053e-03*cos(t)+ 4.065e-01* 7.432e-03*sin(t)), 1.042e-01 +2.000*( -4.065e-01* 8.053e-03*cos(t)+ 9.136e-01* 7.432e-03*sin(t)) not
# Age 95, p23 - p21
-set label "95" at 2.238e-01, 8.634e-02 center
-replot 2.238e-01+ 2.000*( 9.979e-01* 1.231e-02*cos(t)+ 6.495e-02* 9.056e-03*sin(t)), 8.634e-02 +2.000*( -6.495e-02* 1.231e-02*cos(t)+ 9.979e-01* 9.056e-03*sin(t)) not
-set out "biaspar/varpijgrbiaspar123-21.png";replot;
+set label "95" at 2.219e-01, 8.796e-02 center
+replot 2.219e-01+ 2.000*( 9.960e-01* 1.226e-02*cos(t)+ 8.959e-02* 9.318e-03*sin(t)), 8.796e-02 +2.000*( -8.959e-02* 1.226e-02*cos(t)+ 9.960e-01* 9.318e-03*sin(t)) not
+set out;
+set out "biaspar/VARPIJGR_biaspar_123-21.svg";replot;set out;
+# Routine varevsij
+unset title
+
+unset parametric;unset label; set ter svg size 640, 480
+ set log y; unset log x;set xlabel "Age"; set ylabel "Force of mortality (year-1)";
+set out "biaspar/VARMUPTJGR--STABLBASED_biaspar1.svg";
+ plot "biaspar/PRMORPREV-1-STABLBASED_biaspar.txt" u 1:($3) not w l lt 1
+ replot "biaspar/PRMORPREV-1-STABLBASED_biaspar.txt" u 1:(($3+1.96*$4)) t "95% interval" w l lt 2
+ replot "biaspar/PRMORPREV-1-STABLBASED_biaspar.txt" u 1:(($3-1.96*$4)) not w l lt 2
+set out;
+set out "biaspar/VARMUPTJGR--STABLBASED_biaspar1.svg";replot;set out;
+
# Routine varevsij
-unset parametric;unset label; set ter png small size 320, 240
+unset title
+
+unset parametric;unset label; set ter svg size 640, 480
set log y; unset log x;set xlabel "Age"; set ylabel "Force of mortality (year-1)";
- plot "biaspar/prmorprev1-stablbased-rbiaspar.txt" u 1:($3) not w l lt 1
- replot "biaspar/prmorprev1-stablbased-rbiaspar.txt" u 1:(($3+1.96*$4)) t "95% interval" w l lt 2
- replot "biaspar/prmorprev1-stablbased-rbiaspar.txt" u 1:(($3-1.96*$4)) not w l lt 2
-set out "biaspar/varmuptjgr-stablbased-biaspar1.png";replot;
+set out "biaspar/VARMUPTJGR--POPULBASED-NOMOBIL_biaspar1.svg";
+ plot "biaspar/PRMORPREV-1-POPULBASED-NOMOBIL_biaspar.txt" u 1:($3) not w l lt 1
+ replot "biaspar/PRMORPREV-1-POPULBASED-NOMOBIL_biaspar.txt" u 1:(($3+1.96*$4)) t "95% interval" w l lt 2
+ replot "biaspar/PRMORPREV-1-POPULBASED-NOMOBIL_biaspar.txt" u 1:(($3-1.96*$4)) not w l lt 2
+set out;
+set out "biaspar/VARMUPTJGR--POPULBASED-NOMOBIL_biaspar1.svg";replot;set out;
Log filename:biaspar.log
-
-Imach version 0.98q5, August 2015,INED-EUROREVES-Institut de longevite-Japan Society for the Promotion of Science (Grant-in-Aid for Scientific Research 25293121), Intel Software 2015
-$Revision$ $Date$
+Version 0.99r19 $Revision$ $Date$
Enter the parameter file name:
-pathimach=/Users/nbrouard/Documents/imach/cvs/imach/build/osx/
-pathtot=html/doc/biaspar.imach
- path=html/doc/
+pathimach=/Users/nbrouard/Documents/imach/cvs/imach/html/doc/
+pathtot=biaspar.imach
+ path=/Users/nbrouard/Documents/imach/cvs/imach/html/doc/
optionfile=biaspar.imach
optionfilext=imach
optionfilefiname='biaspar'
-Local time (at start): Tue Aug 18 18:20:44 2015
-# Imach version 0.98q2, April 2015, INED-EUROREVES \r
+Compiled with: Clang/LLVM GNU GCC/G++ for Mac OS 64-bit using GNU C version 40201.
+Running on: Darwin tugault.local 18.2.0 Darwin Kernel Version 18.2.0: Thu Dec 20 20:46:53 PST 2018; root:xnu-4903.241.1~1/RELEASE_X86_64 x86_64
+ Local time (at start): Wed May 22 22:32:36 2019
+#0.99r10\r
+#Number of iterations & function calls = 33 & 6007, -2 Log likelihood = 46617.105944307616\r
title=1st_example datafile=data1.txt lastobs=8600 firstpass=1 lastpass=4
-ftol=1.000000e-08 stepm=1 ncovcol=2 nlstate=2 ndeath=1 maxwav=4 mle=1 weight=0
-model=1+age+.
-#q Parameters\r
-12 -12.245240 0.092358
-13 -10.671890 0.060969
-21 -2.645345 -0.022325
-23 -4.773317 0.007859
-# Scales\r
-12 0.000000e+00 0.000000e+00
-13 0.000000e+00 0.000000e+00
-21 0.000000e+00 0.000000e+00
-23 0.000000e+00 0.000000e+00
-#covariance matrix#\r
-121 0.00000e+00
-122 0.00000e+00 0.00000e+00
-131 0.00000e+00 0.00000e+00 0.00000e+00
-132 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
-211 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
-212 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
-231 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
-232 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
+ftol=1.000000e-08 stepm=1 ncovcol=2 nqv=0 ntv=0 nqtv=0 nlstate=2 ndeath=1 maxwav=4 mle=1 weight=0
+model=1+age+
+# Parameters nlstate*nlstate*ncov a12*1 + b12 * age + ...\r
+12 -12.317575 0.093112
+13 -10.163167 0.054838
+21 -1.668568 -0.033906
+23 -5.686920 0.018135
+# Scales (for hessian or gradient estimation)\r
+12 2.000000e-04 2.000000e-06
+13 3.000000e-04 4.000000e-06
+21 1.000000e-04 1.000000e-06
+23 3.000000e-04 3.000000e-06
+# Covariance matrix \r
+# 121 Var(a12)\r
+# 122 Cov(b12,a12) Var(b12)\r
+# ...\r
+# 232 Cov(b23,a12) Cov(b23,b12) ... Var (b23)\r
+#121 Var(a12)\r
+#122 Cov(b12,a12) Var(b12)\r
+#131 Cov(a13,a12) Cov(a13,b12) Var(a13)\r
+#132 Cov(b13,a12) Cov(b13,b12) Cov(b13,a13) Var(b13)\r
+#211 Cov(a21,a12) Cov(a21,b12) Cov(a21,a13) Cov(a21,b13) Var(a21)\r
+#212 Cov(b21,a12) Cov(b21,b12) Cov(b21,a13) Cov(b21,b13) Cov(b21,a21) Var(b21)\r
+#231 Cov(a23,a12) Cov(a23,b12) Cov(a23,a13) Cov(a23,b13) Cov(a23,a21) Cov(a23,b21) Var(a23)\r
+#232 Cov(b23,a12) Cov(b23,b12) Cov(b23,a13) Cov(b23,b13) Cov(b23,a21) Cov(b23,b21) Cov(b23,a23) Var(b23)\r
+121 1.08463e-01
+122 -1.30062e-03 1.56897e-05
+131 -5.39098e-02 6.59806e-04 3.01977e-01
+132 6.43147e-04 -7.94760e-06 -3.73264e-03 4.64351e-05
+211 7.75092e-02 -9.24613e-04 1.91004e-02 -2.45003e-04 4.71113e-01
+212 -9.24119e-04 1.10769e-05 -2.40929e-04 3.09854e-06 -5.60404e-03 6.70191e-05
+231 4.86258e-02 -5.50028e-04 -9.42536e-02 1.06700e-03 -2.57577e-03 2.94283e-05 2.50103e-01
+232 -5.52529e-04 6.29410e-06 1.11646e-03 -1.27952e-05 3.55428e-05 -4.16204e-07 -2.83674e-03 3.23581e-05
+Covariate type in the data: V1, DummyV(V1)=0, FixedV(V1)=0
+Covariate type in the data: V2, DummyV(V2)=0, FixedV(V2)=0
+Comment line
+# Data. Two fixed dummy covariates (ncovcol=2), zero fixed quantitative covariate (nqv=0), zero time varying covariate (ntv=0), zero quantitative time varying covariate (nqtv=0)
+Comment line
+#id V1 V2 W DBirth DDeath D_1st S_1 D_2nd S_2 D_3rd S_3 D_4 S_4
+Model=
+Typevar: 0 for simple covariate (dummy, quantitative, fixed or varying), 1 for age product, 2 for product
+Fixed[k] 0=fixed (product or simple), 1 varying, 2 fixed with age product, 3 varying with age product
+Dummy[k] 0=dummy (0 1), 1 quantitative (single or product without age), 2 dummy with age product, 3 quant with age product
+ncoveff=0, nqfveff=0, ntveff=0, nqtveff=0, cptcovn=0
+ncovf=0, ncovv=0, ncova=0, nsd=0, nsq=0
Total number of individuals= 8270, Agemin = 70.00, Agemax= 104.17
-Delay (in months) between two waves Min=1 (for indiviudal 6017) Max=74 (5168) Mean=24.040080
+Information! Unknown status for individual 1 line=1 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 40 line=40 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(40)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 55 line=55 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(55)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 84 line=84 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(84)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 92 line=92 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(92)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 103 line=103 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(103)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 186 line=186 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(186)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 202 line=202 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(202)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 207 line=207 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(207)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 239 line=239 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(239)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 260 line=260 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(260)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 269 line=269 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(269)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 302 line=302 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(302)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 329 line=329 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(329)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 335 line=335 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(335)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 346 line=346 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(346)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 347 line=347 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(347)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 378 line=378 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(378)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 414 line=414 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(414)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 431 line=431 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(431)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 443 line=443 occurred at last wave 4 at known date 5/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(443)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 482 line=482 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(482)/wave(4) contributes to the likelihood as 1-p-13 .
+Error! Death for individual 503 line=503 occurred at 10/1990 after last wave 4 interviewed at 7/1990. Potential bias if other individuals are still alive at this date but ignored. This case (503)/wave (4) is skipped, no contribution to likelihood. Please add a new fictive wave at the date of last vital status scan, with a dead status or alive but unknown state status (-1). See documentation
- Age 70 1.=665 loss[1]=5.1% 2.=0 loss[2]=NaNQ% 1.=631 prev[1]=100.0% 2.=0 prev[2]=0.0% 1-1=34 11=563 12=33 13=35 33=6
-Age 71 1.=663 loss[1]=5.7% 2.=2 loss[2]=0.0% 1.=625 prev[1]=99.7% 2.=2 prev[2]=0.3% 1-1=38 11=530 12=47 13=48 21=1 22=1 3-1=2 33=31
-Age 72 1.=1181 loss[1]=5.6% 2.=35 loss[2]=5.7% 1.=1115 prev[1]=97.1% 2.=33 prev[2]=2.9% 1-1=66 11=969 12=70 13=76 2-1=2 21=9 22=18 23=6 3-1=2 33=87
-Age 73 1.=1119 loss[1]=5.1% 2.=48 loss[2]=8.3% 1.=1062 prev[1]=96.0% 2.=44 prev[2]=4.0% 1-1=57 11=912 12=71 13=79 2-1=4 21=12 22=26 23=6 3-1=12 33=94
-Age 74 1.=1563 loss[1]=4.2% 2.=90 loss[2]=4.4% 1.=1497 prev[1]=94.6% 2.=86 prev[2]=5.4% 1-1=66 11=1277 12=104 13=116 2-1=4 21=30 22=41 23=15 3-1=4 33=125
-Age 75 1.=1463 loss[1]=6.2% 2.=98 loss[2]=7.1% 1.=1372 prev[1]=93.8% 2.=91 prev[2]=6.2% 1-1=91 11=1155 12=98 13=119 2-1=7 21=15 22=55 23=21 3-1=4 33=124
-Age 76 1.=1390 loss[1]=4.5% 2.=112 loss[2]=3.6% 1.=1328 prev[1]=92.5% 2.=108 prev[2]=7.5% 1-1=62 11=1090 12=111 13=127 2-1=4 21=25 22=54 23=29 3-1=4 33=142
-Age 77 1.=1271 loss[1]=5.0% 2.=101 loss[2]=5.9% 1.=1208 prev[1]=92.7% 2.=95 prev[2]=7.3% 1-1=63 11=962 12=131 13=115 2-1=6 21=24 22=52 23=19 3-1=7 33=145
-Age 78 1.=1166 loss[1]=4.9% 2.=122 loss[2]=4.1% 1.=1109 prev[1]=90.5% 2.=117 prev[2]=9.5% 1-1=57 11=894 12=111 13=104 2-1=5 21=21 22=75 23=21 3-1=4 33=146
-Age 79 1.=1045 loss[1]=5.7% 2.=133 loss[2]=5.3% 1.=985 prev[1]=88.7% 2.=126 prev[2]=11.3% 1-1=60 11=773 12=109 13=103 2-1=7 21=31 22=68 23=27 3-1=8 33=136
-Age 80 1.=986 loss[1]=4.1% 2.=128 loss[2]=3.1% 1.=946 prev[1]=88.4% 2.=124 prev[2]=11.6% 1-1=40 11=718 12=140 13=88 2-1=4 21=21 22=71 23=32 3-1=6 33=121
-Age 81 1.=872 loss[1]=4.9% 2.=130 loss[2]=3.8% 1.=829 prev[1]=86.9% 2.=125 prev[2]=13.1% 1-1=43 11=629 12=102 13=98 2-1=5 21=24 22=74 23=27 3-1=5 33=126
-Age 82 1.=803 loss[1]=4.4% 2.=146 loss[2]=4.1% 1.=768 prev[1]=84.6% 2.=140 prev[2]=15.4% 1-1=35 11=546 12=131 13=91 2-1=6 21=15 22=77 23=48 3-1=11 33=100
-Age 83 1.=726 loss[1]=6.5% 2.=116 loss[2]=0.9% 1.=679 prev[1]=85.5% 2.=115 prev[2]=14.5% 1-1=47 11=456 12=124 13=99 2-1=1 21=19 22=71 23=25 3-1=8 33=133
-Age 84 1.=622 loss[1]=5.6% 2.=147 loss[2]=4.1% 1.=587 prev[1]=80.6% 2.=141 prev[2]=19.4% 1-1=35 11=394 12=112 13=81 2-1=6 21=18 22=82 23=41 3-1=5 33=114
-Age 85 1.=526 loss[1]=4.4% 2.=127 loss[2]=6.3% 1.=503 prev[1]=80.9% 2.=119 prev[2]=19.1% 1-1=23 11=315 12=122 13=66 2-1=8 21=13 22=67 23=39 3-1=4 33=137
-Age 86 1.=464 loss[1]=5.2% 2.=130 loss[2]=6.2% 1.=440 prev[1]=78.3% 2.=122 prev[2]=21.7% 1-1=24 11=271 12=103 13=66 2-1=8 21=18 22=70 23=34 3-1=3 33=112
-Age 87 1.=372 loss[1]=4.8% 2.=143 loss[2]=2.8% 1.=354 prev[1]=71.8% 2.=139 prev[2]=28.2% 1-1=18 11=219 12=79 13=56 2-1=4 21=15 22=71 23=53 3-1=1 33=81
-Age 88 1.=283 loss[1]=3.9% 2.=164 loss[2]=6.7% 1.=272 prev[1]=64.0% 2.=153 prev[2]=36.0% 1-1=11 11=172 12=63 13=37 2-1=11 21=25 22=84 23=44 3-1=8 33=102
-Age 89 1.=227 loss[1]=3.5% 2.=152 loss[2]=4.6% 1.=219 prev[1]=60.2% 2.=145 prev[2]=39.8% 1-1=8 11=135 12=45 13=39 2-1=7 21=21 22=73 23=51 3-1=2 33=77
-Age 90 1.=181 loss[1]=3.9% 2.=143 loss[2]=4.9% 1.=174 prev[1]=56.1% 2.=136 prev[2]=43.9% 1-1=7 11=104 12=41 13=29 2-1=7 21=18 22=82 23=36 3-1=1 33=75
-Age 91 1.=146 loss[1]=4.1% 2.=109 loss[2]=8.3% 1.=140 prev[1]=58.3% 2.=100 prev[2]=41.7% 1-1=6 11=78 12=36 13=26 2-1=9 21=11 22=51 23=38 3-1=3 33=77
-Age 92 1.=114 loss[1]=6.1% 2.=101 loss[2]=2.0% 1.=107 prev[1]=51.9% 2.=99 prev[2]=48.1% 1-1=7 11=47 12=43 13=17 2-1=2 21=11 22=52 23=36 3-1=2 33=61
-Age 93 1.=78 loss[1]=3.8% 2.=89 loss[2]=12.4% 1.=75 prev[1]=49.0% 2.=78 prev[2]=51.0% 1-1=3 11=36 12=25 13=14 2-1=11 21=10 22=46 23=22 3-1=2 33=62
-Age 94 1.=51 loss[1]=3.9% 2.=71 loss[2]=1.4% 1.=49 prev[1]=41.2% 2.=70 prev[2]=58.8% 1-1=2 11=24 12=8 13=17 2-1=1 21=4 22=39 23=27 3-1=1 33=37
-Age 95 1.=38 loss[1]=5.3% 2.=57 loss[2]=5.3% 1.=36 prev[1]=40.0% 2.=54 prev[2]=60.0% 1-1=2 11=16 12=13 13=7 2-1=3 21=4 22=31 23=19 3-1=1 33=39
-Age 96 1.=25 loss[1]=4.0% 2.=43 loss[2]=4.7% 1.=24 prev[1]=36.9% 2.=41 prev[2]=63.1% 1-1=1 11=11 12=9 13=4 2-1=2 21=2 22=22 23=17 3-1=2 33=25
-Age 97 1.=11 loss[1]=0.0% 2.=36 loss[2]=8.3% 1.=11 prev[1]=25.0% 2.=33 prev[2]=75.0% 11=4 12=7 2-1=3 21=2 22=22 23=9 33=23
-Age 98 1.=5 loss[1]=0.0% 2.=21 loss[2]=0.0% 1.=5 prev[1]=19.2% 2.=21 prev[2]=80.8% 11=1 12=4 21=1 22=13 23=7 3-1=1 33=17
-Age 99 1.=6 loss[1]=0.0% 2.=14 loss[2]=0.0% 1.=6 prev[1]=30.0% 2.=14 prev[2]=70.0% 11=3 12=1 13=2 21=2 22=2 23=10 33=6
-Age 100 1.=2 loss[1]=0.0% 2.=9 loss[2]=0.0% 1.=2 prev[1]=18.2% 2.=9 prev[2]=81.8% 11=1 12=1 21=2 22=5 23=2 33=2
-Age 101 1.=2 loss[1]=50.0% 2.=0 loss[2]=NaNQ% 1.=1 prev[1]=100.0% 2.=0 prev[2]=0.0% 1-1=1 12=1 3-1=1 33=3
-Age 102 1.=2 loss[1]=0.0% 2.=1 loss[2]=0.0% 1.=2 prev[1]=66.7% 2.=1 prev[2]=33.3% 12=2 22=1 33=1
-Age 103 1.=0 loss[1]=NaNQ% 2.=0 loss[2]=NaNQ% 1.=0 prev[1]=NaNQ% 2.=0 prev[2]=NaNQ%
-Age 104 1.=0 loss[1]=NaNQ% 2.=0 loss[2]=NaNQ% 1.=0 prev[1]=NaNQ% 2.=0 prev[2]=NaNQ%
-Age 105 1.=0 loss[1]=NaNQ% 2.=0 loss[2]=NaNQ% 1.=0 prev[1]=NaNQ% 2.=0 prev[2]=NaNQ%
-Age 106 1.=0 loss[1]=NaNQ% 2.=0 loss[2]=NaNQ% 1.=0 prev[1]=NaNQ% 2.=0 prev[2]=NaNQ% -1-1=69 -11=718 -12=195 -13=261
-Total 1.=18068 loss[1]=5.0% 2.=2818 loss[2]=4.9% 1.=17161 prev[1]=86.5% 2.=2681 prev[2]=13.5% -1-1=69 -11=718 -12=195 -13=261 1-1=907 11=13305 12=2097 13=1759 2-1=137 21=424 22=1496 23=761 3-1=114 33=2567
-Powell
+Information! Unknown status for individual 516 line=516 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(516)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 537 line=537 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(537)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 544 line=544 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(544)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 569 line=569 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(569)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 581 line=581 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(581)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 595 line=595 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(595)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 620 line=620 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(620)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 630 line=630 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(630)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 696 line=696 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(696)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 699 line=699 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(699)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 704 line=704 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(704)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 707 line=707 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(707)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 712 line=712 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(712)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 713 line=713 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(713)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 762 line=762 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(762)/wave(4) contributes to the likelihood as 1-p-13 .
+Warning! Death for individual 766 line=766 occurred at 1/1989 before last wave 4 interviewed at 7/1990 and should have been coded as death instead of '-1'. This case (766)/wave (4) is contributing to likelihood.
+Information! Unknown status for individual 770 line=770 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(770)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 798 line=798 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(798)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 853 line=853 occurred at last wave 4 at known date 9/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(853)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 855 line=855 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(855)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 918 line=918 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(918)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 942 line=942 occurred at last wave 4 at known date 9/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(942)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 946 line=946 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(946)/wave(4) contributes to the likelihood as 1-p-13 .
+Warning! Death for individual 951 line=951 occurred at 4/1989 before last wave 4 interviewed at 7/1990 and should have been coded as death instead of '-1'. This case (951)/wave (4) is contributing to likelihood.
+Information! Unknown status for individual 958 line=958 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(958)/wave(4) contributes to the likelihood as 1-p-13 .
+Warning! Death for individual 966 line=966 occurred at 3/1990 before last wave 4 interviewed at 7/1990 and should have been coded as death instead of '-1'. This case (966)/wave (4) is contributing to likelihood.
+Information! Unknown status for individual 970 line=970 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(970)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 982 line=982 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(982)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1004 line=1004 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1004)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1039 line=1039 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1039)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1081 line=1081 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1081)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1091 line=1091 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1091)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1099 line=1099 occurred at last wave 4 at known date 9/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1099)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1157 line=1157 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1157)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1177 line=1177 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1177)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1179 line=1179 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1179)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1197 line=1197 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1197)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1201 line=1201 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1201)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1215 line=1215 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1215)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1242 line=1242 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1242)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1268 line=1268 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1268)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1270 line=1270 occurred at last wave 4 at known date 5/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1270)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1284 line=1284 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1284)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1295 line=1295 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1295)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1297 line=1297 occurred at last wave 4 at known date 5/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1297)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1363 line=1363 occurred at last wave 4 at known date 5/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1363)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1406 line=1406 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1406)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1477 line=1477 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1477)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1543 line=1543 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1543)/wave(4) contributes to the likelihood as 1-p-13 .
+Warning! Death for individual 1571 line=1571 occurred at 11/1989 before last wave 4 interviewed at 8/1990 and should have been coded as death instead of '-1'. This case (1571)/wave (4) is contributing to likelihood.
+Information! Unknown status for individual 1580 line=1580 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1580)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1581 line=1581 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1581)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1593 line=1593 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1593)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1620 line=1620 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1620)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1672 line=1672 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1672)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1680 line=1680 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1680)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1737 line=1737 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1737)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1764 line=1764 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1764)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1777 line=1777 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1777)/wave(4) contributes to the likelihood as 1-p-13 .
+Warning! Death for individual 1781 line=1781 occurred at 4/1990 before last wave 4 interviewed at 7/1990 and should have been coded as death instead of '-1'. This case (1781)/wave (4) is contributing to likelihood.
+Information! Unknown status for individual 1803 line=1803 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1803)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1849 line=1849 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1849)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1881 line=1881 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1881)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1903 line=1903 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1903)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1937 line=1937 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1937)/wave(4) contributes to the likelihood as 1-p-13 .
+Warning! Death for individual 1948 line=1948 occurred at 9/1989 before last wave 4 interviewed at 7/1990 and should have been coded as death instead of '-1'. This case (1948)/wave (4) is contributing to likelihood.
+Information! Unknown status for individual 1949 line=1949 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1949)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1968 line=1968 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1968)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 1984 line=1984 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(1984)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 2009 line=2009 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2009)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 2019 line=2019 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2019)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 2025 line=2025 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2025)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 2027 line=2027 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2027)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 2032 line=2032 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2032)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 2036 line=2036 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2036)/wave(4) contributes to the likelihood as 1-p-13 .
+Warning! Death for individual 2062 line=2062 occurred at 1/1990 before last wave 4 interviewed at 7/1990 and should have been coded as death instead of '-1'. This case (2062)/wave (4) is contributing to likelihood.
+Information! Unknown status for individual 2069 line=2069 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2069)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 2096 line=2096 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2096)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 2189 line=2189 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2189)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 2198 line=2198 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2198)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 2311 line=2311 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2311)/wave(4) contributes to the likelihood as 1-p-13 .
+Error! Death for individual 2327 line=2327 occurred at 6/1991 after last wave 4 interviewed at 7/1990. Potential bias if other individuals are still alive at this date but ignored. This case (2327)/wave (4) is skipped, no contribution to likelihood. Please add a new fictive wave at the date of last vital status scan, with a dead status or alive but unknown state status (-1). See documentation
-Powell iter=1 -2*LL=46542.394740741831 0 sec. 1 sec. 1 -12.245240000000 2 0.092358000000 3 -10.671890000000 4 0.060969000000 5 -2.645345000000 6 -0.022325000000 7 -4.773317000000 8 0.007859000000
+Information! Unknown status for individual 2354 line=2354 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2354)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 2377 line=2377 occurred at last wave 4 at known date 4/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2377)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 2388 line=2388 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2388)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 2392 line=2392 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2392)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 2432 line=2432 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2432)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 2446 line=2446 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2446)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 2513 line=2513 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2513)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 2574 line=2574 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2574)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 2583 line=2583 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2583)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 2638 line=2638 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2638)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 2643 line=2643 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2643)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 2655 line=2655 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2655)/wave(4) contributes to the likelihood as 1-p-13 .
+Warning! Death for individual 2662 line=2662 occurred at 7/1989 before last wave 4 interviewed at 7/1990 and should have been coded as death instead of '-1'. This case (2662)/wave (4) is contributing to likelihood.
+Information! Unknown status for individual 2681 line=2681 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2681)/wave(4) contributes to the likelihood as 1-p-13 .
+Error! Death for individual 2743 line=2743 occurred at 6/1991 after last wave 4 interviewed at 7/1990. Potential bias if other individuals are still alive at this date but ignored. This case (2743)/wave (4) is skipped, no contribution to likelihood. Please add a new fictive wave at the date of last vital status scan, with a dead status or alive but unknown state status (-1). See documentation
-Considering the time needed for this last iteration #1: 0 seconds,
- - if your program needs 10 iterations to converge, convergence will be
- reached in 0 day(s) 0 hour(s) 0 minute(s) 0 second(s) i.e.
- on Tue Aug 18 18:20:45 2015 (current time is Tue Aug 18 18:20:45 2015);
- - if your program needs 20 iterations to converge, convergence will be
- reached in 0 day(s) 0 hour(s) 0 minute(s) 0 second(s) i.e.
- on Tue Aug 18 18:20:45 2015 (current time is Tue Aug 18 18:20:45 2015);
- - if your program needs 30 iterations to converge, convergence will be
- reached in 0 day(s) 0 hour(s) 0 minute(s) 0 second(s) i.e.
- on Tue Aug 18 18:20:45 2015 (current time is Tue Aug 18 18:20:45 2015);
-1.............2.....................3...............4...................5..........6.............7...............8..................
-Powell iter=2 -2*LL=46542.366177077667 23 sec. 24 sec. 1 -12.245474667968 2 0.092360078052 3 -10.671718090206 4 0.060967648859 5 -2.641226372487 6 -0.022378298324 7 -4.773347352998 8 0.007897501084
+Warning! Death for individual 2782 line=2782 occurred at 9/1989 before last wave 4 interviewed at 7/1990 and should have been coded as death instead of '-1'. This case (2782)/wave (4) is contributing to likelihood.
+Information! Unknown status for individual 2841 line=2841 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2841)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 2888 line=2888 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2888)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 2919 line=2919 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2919)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 2970 line=2970 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2970)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 2980 line=2980 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(2980)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3012 line=3012 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3012)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3130 line=3130 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3130)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3140 line=3140 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3140)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3187 line=3187 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3187)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3234 line=3234 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3234)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3274 line=3274 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3274)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3331 line=3331 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3331)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3394 line=3394 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3394)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3396 line=3396 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3396)/wave(4) contributes to the likelihood as 1-p-13 .
+Warning! Death for individual 3414 line=3414 occurred at 9/1989 before last wave 4 interviewed at 7/1990 and should have been coded as death instead of '-1'. This case (3414)/wave (4) is contributing to likelihood.
+Information! Unknown status for individual 3417 line=3417 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3417)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3430 line=3430 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3430)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3451 line=3451 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3451)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3453 line=3453 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3453)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3472 line=3472 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3472)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3501 line=3501 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3501)/wave(4) contributes to the likelihood as 1-p-13 .
+Warning! Death for individual 3517 line=3517 occurred at 6/1990 before last wave 4 interviewed at 7/1990 and should have been coded as death instead of '-1'. This case (3517)/wave (4) is contributing to likelihood.
+Information! Unknown status for individual 3539 line=3539 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3539)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3557 line=3557 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3557)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3564 line=3564 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3564)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3565 line=3565 occurred at last wave 4 at known date 10/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3565)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3591 line=3591 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3591)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3595 line=3595 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3595)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3652 line=3652 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3652)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3684 line=3684 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3684)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3777 line=3777 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3777)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3787 line=3787 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3787)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3790 line=3790 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3790)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3805 line=3805 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3805)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3841 line=3841 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3841)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3848 line=3848 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3848)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3863 line=3863 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3863)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3949 line=3949 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3949)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3960 line=3960 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3960)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 3981 line=3981 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(3981)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 4176 line=4176 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(4176)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 4186 line=4186 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(4186)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 4190 line=4190 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(4190)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 4288 line=4288 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(4288)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 4295 line=4295 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(4295)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 4320 line=4320 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(4320)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 4360 line=4360 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(4360)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 4381 line=4381 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(4381)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 4403 line=4403 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(4403)/wave(4) contributes to the likelihood as 1-p-13 .
+Error! Death for individual 4417 line=4417 occurred at 8/1990 after last wave 4 interviewed at 7/1990. Potential bias if other individuals are still alive at this date but ignored. This case (4417)/wave (4) is skipped, no contribution to likelihood. Please add a new fictive wave at the date of last vital status scan, with a dead status or alive but unknown state status (-1). See documentation
-Considering the time needed for this last iteration #2: 23 seconds,
- - if your program needs 10 iterations to converge, convergence will be
- reached in 0 day(s) 0 hour(s) 3 minute(s) 4 second(s) i.e.
- on Tue Aug 18 18:24:12 2015 (current time is Tue Aug 18 18:21:08 2015);
- - if your program needs 20 iterations to converge, convergence will be
- reached in 0 day(s) 0 hour(s) 6 minute(s) 54 second(s) i.e.
- on Tue Aug 18 18:28:02 2015 (current time is Tue Aug 18 18:21:08 2015);
- - if your program needs 30 iterations to converge, convergence will be
- reached in 0 day(s) 0 hour(s) 10 minute(s) 44 second(s) i.e.
- on Tue Aug 18 18:31:52 2015 (current time is Tue Aug 18 18:21:08 2015);
-1............2.................3............4...................5..........6.............7..........8.........................Gaining to use new average direction of P0 P8 instead of biggest increase direction 8 :
+Information! Unknown status for individual 4467 line=4467 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(4467)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 4476 line=4476 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(4476)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 4528 line=4528 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(4528)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 4558 line=4558 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(4558)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 4611 line=4611 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(4611)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 4614 line=4614 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(4614)/wave(4) contributes to the likelihood as 1-p-13 .
+Warning! Death for individual 4747 line=4747 occurred at 3/1990 before last wave 4 interviewed at 7/1990 and should have been coded as death instead of '-1'. This case (4747)/wave (4) is contributing to likelihood.
+Information! Unknown status for individual 4751 line=4751 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(4751)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 4758 line=4758 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(4758)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 4816 line=4816 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(4816)/wave(4) contributes to the likelihood as 1-p-13 .
+Error! Death for individual 4834 line=4834 occurred at 10/1990 after last wave 4 interviewed at 6/1990. Potential bias if other individuals are still alive at this date but ignored. This case (4834)/wave (4) is skipped, no contribution to likelihood. Please add a new fictive wave at the date of last vital status scan, with a dead status or alive but unknown state status (-1). See documentation
-Powell iter=3 -2*LL=46542.180101184662 22 sec. 46 sec. 1 -12.243554549912 2 0.092379481639 3 -10.687240979880 4 0.060948425212 5 -2.601681584093 6 -0.022889573696 7 -4.797716971796 8 0.008263023171
+Information! Unknown status for individual 4836 line=4836 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(4836)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 4864 line=4864 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(4864)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 4906 line=4906 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(4906)/wave(4) contributes to the likelihood as 1-p-13 .
+Error! Death for individual 4928 line=4928 occurred at 11/1990 after last wave 4 interviewed at 8/1990. Potential bias if other individuals are still alive at this date but ignored. This case (4928)/wave (4) is skipped, no contribution to likelihood. Please add a new fictive wave at the date of last vital status scan, with a dead status or alive but unknown state status (-1). See documentation
-Considering the time needed for this last iteration #3: 22 seconds,
- - if your program needs 10 iterations to converge, convergence will be
- reached in 0 day(s) 0 hour(s) 2 minute(s) 34 second(s) i.e.
- on Tue Aug 18 18:24:04 2015 (current time is Tue Aug 18 18:21:30 2015);
- - if your program needs 20 iterations to converge, convergence will be
- reached in 0 day(s) 0 hour(s) 6 minute(s) 14 second(s) i.e.
- on Tue Aug 18 18:27:44 2015 (current time is Tue Aug 18 18:21:30 2015);
- - if your program needs 30 iterations to converge, convergence will be
- reached in 0 day(s) 0 hour(s) 9 minute(s) 54 second(s) i.e.
- on Tue Aug 18 18:31:24 2015 (current time is Tue Aug 18 18:21:30 2015);
-1.................2.................3..........4....................5..........6.............7..........8.......
-Powell iter=4 -2*LL=46541.944882756965 20 sec. 66 sec. 1 -12.242484543366 2 0.092391161272 3 -10.684956992257 4 0.060933369989 5 -2.571767476006 6 -0.023216409131 7 -4.817698692685 8 0.008460497128
-1.........2..................3..........4....................5..........6.............7..........8..............Gaining to use new average direction of P0 P8 instead of biggest increase direction 3 :
+Information! Unknown status for individual 4946 line=4946 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(4946)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 4962 line=4962 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(4962)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 4982 line=4982 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(4982)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 4997 line=4997 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(4997)/wave(4) contributes to the likelihood as 1-p-13 .
+Warning! Death for individual 5009 line=5009 occurred at 6/1990 before last wave 4 interviewed at 6/1990 and should have been coded as death instead of '-1'. This case (5009)/wave (4) is contributing to likelihood.
+Warning! Death for individual 5011 line=5011 occurred at 2/1989 before last wave 4 interviewed at 8/1990 and should have been coded as death instead of '-1'. This case (5011)/wave (4) is contributing to likelihood.
+Information! Unknown status for individual 5063 line=5063 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5063)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 5089 line=5089 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5089)/wave(4) contributes to the likelihood as 1-p-13 .
+Warning! Death for individual 5114 line=5114 occurred at 3/1990 before last wave 4 interviewed at 8/1990 and should have been coded as death instead of '-1'. This case (5114)/wave (4) is contributing to likelihood.
+Information! Unknown status for individual 5122 line=5122 occurred at last wave 4 at known date 9/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5122)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 5135 line=5135 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5135)/wave(4) contributes to the likelihood as 1-p-13 .
+Warning! Death for individual 5142 line=5142 occurred at 8/1989 before last wave 4 interviewed at 8/1990 and should have been coded as death instead of '-1'. This case (5142)/wave (4) is contributing to likelihood.
+Information! Unknown status for individual 5162 line=5162 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5162)/wave(4) contributes to the likelihood as 1-p-13 .
+Warning! Death for individual 5189 line=5189 occurred at 12/1989 before last wave 4 interviewed at 6/1990 and should have been coded as death instead of '-1'. This case (5189)/wave (4) is contributing to likelihood.
+Information! Unknown status for individual 5216 line=5216 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5216)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 5223 line=5223 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5223)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 5239 line=5239 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5239)/wave(4) contributes to the likelihood as 1-p-13 .
+Warning! Death for individual 5251 line=5251 occurred at 2/1990 before last wave 4 interviewed at 7/1990 and should have been coded as death instead of '-1'. This case (5251)/wave (4) is contributing to likelihood.
+Information! Unknown status for individual 5265 line=5265 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5265)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 5307 line=5307 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5307)/wave(4) contributes to the likelihood as 1-p-13 .
+Warning! Death for individual 5355 line=5355 occurred at 3/1990 before last wave 4 interviewed at 6/1990 and should have been coded as death instead of '-1'. This case (5355)/wave (4) is contributing to likelihood.
+Error! Death for individual 5378 line=5378 occurred at 9/1990 after last wave 4 interviewed at 7/1990. Potential bias if other individuals are still alive at this date but ignored. This case (5378)/wave (4) is skipped, no contribution to likelihood. Please add a new fictive wave at the date of last vital status scan, with a dead status or alive but unknown state status (-1). See documentation
+
+Information! Unknown status for individual 5388 line=5388 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5388)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 5407 line=5407 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5407)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 5503 line=5503 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5503)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 5508 line=5508 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5508)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 5521 line=5521 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5521)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 5533 line=5533 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5533)/wave(4) contributes to the likelihood as 1-p-13 .
+Warning! Death for individual 5541 line=5541 occurred at 3/1989 before last wave 4 interviewed at 7/1990 and should have been coded as death instead of '-1'. This case (5541)/wave (4) is contributing to likelihood.
+Information! Unknown status for individual 5558 line=5558 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5558)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 5623 line=5623 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5623)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 5632 line=5632 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5632)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 5653 line=5653 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5653)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 5659 line=5659 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5659)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 5670 line=5670 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5670)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 5676 line=5676 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5676)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 5704 line=5704 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5704)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 5759 line=5759 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5759)/wave(4) contributes to the likelihood as 1-p-13 .
+Warning! Death for individual 5769 line=5769 occurred at 8/1989 before last wave 4 interviewed at 7/1990 and should have been coded as death instead of '-1'. This case (5769)/wave (4) is contributing to likelihood.
+Information! Unknown status for individual 5785 line=5785 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5785)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 5834 line=5834 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5834)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 5887 line=5887 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5887)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 5913 line=5913 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5913)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 5914 line=5914 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5914)/wave(4) contributes to the likelihood as 1-p-13 .
+Warning! Death for individual 5921 line=5921 occurred at 5/1989 before last wave 4 interviewed at 9/1990 and should have been coded as death instead of '-1'. This case (5921)/wave (4) is contributing to likelihood.
+Information! Unknown status for individual 5927 line=5927 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5927)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 5983 line=5983 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(5983)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 6050 line=6050 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(6050)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 6051 line=6051 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(6051)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 6091 line=6091 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(6091)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 6111 line=6111 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(6111)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 6180 line=6180 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(6180)/wave(4) contributes to the likelihood as 1-p-13 .
+Error! Death for individual 6294 line=6294 occurred at 11/1990 after last wave 4 interviewed at 7/1990. Potential bias if other individuals are still alive at this date but ignored. This case (6294)/wave (4) is skipped, no contribution to likelihood. Please add a new fictive wave at the date of last vital status scan, with a dead status or alive but unknown state status (-1). See documentation
-Powell iter=5 -2*LL=46538.155159115864 21 sec. 87 sec. 1 -12.277818683456 2 0.092782784461 3 -10.606927524558 4 0.060433498790 5 -1.636618302422 6 -0.034173907116 7 -5.415159199358 8 0.015209892785
-1.........2..............3........4...................5.........6..............7..........8........directest= 3.559482495482, t= -1.587615025675, f1= 46538.155159115864,f2= 46537.965437599931,f3= 46541.605702231980, del= 0.135251826250
-f1-2f2+f3= 3.829986147983, f1-f2-del= 0.054469689683, f1-f3= -3.450543116116
+Information! Unknown status for individual 6299 line=6299 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(6299)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 6316 line=6316 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(6316)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 6326 line=6326 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(6326)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 6335 line=6335 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(6335)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 6346 line=6346 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(6346)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 6349 line=6349 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(6349)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 6361 line=6361 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(6361)/wave(4) contributes to the likelihood as 1-p-13 .
+Warning! Death for individual 6392 line=6392 occurred at 6/1990 before last wave 4 interviewed at 7/1990 and should have been coded as death instead of '-1'. This case (6392)/wave (4) is contributing to likelihood.
+Information! Unknown status for individual 6398 line=6398 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(6398)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 6424 line=6424 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(6424)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 6459 line=6459 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(6459)/wave(4) contributes to the likelihood as 1-p-13 .
+Error! Death for individual 6563 line=6563 occurred at 6/1991 after last wave 4 interviewed at 7/1990. Potential bias if other individuals are still alive at this date but ignored. This case (6563)/wave (4) is skipped, no contribution to likelihood. Please add a new fictive wave at the date of last vital status scan, with a dead status or alive but unknown state status (-1). See documentation
-Powell iter=6 -2*LL=46537.965437599931 17 sec. 104 sec. 1 -12.286824634318 2 0.092796946661 3 -10.605118146544 4 0.060392081785 5 -1.601001695625 6 -0.034696856573 7 -5.445686958604 8 0.015541650236
-1............2...............3..........4...........................5.........6.............7.................8........
-Powell iter=7 -2*LL=46537.953896904983 21 sec. 125 sec. 1 -12.288663475348 2 0.092794258423 3 -10.602520523521 4 0.060372798332 5 -1.596765949187 6 -0.034763087057 7 -5.450260084532 8 0.015582787665
-1...........2...............3........4...................5..........6..............7............8.............Gaining to use new average direction of P0 P8 instead of biggest increase direction 4 :
+Information! Unknown status for individual 6603 line=6603 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(6603)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 6634 line=6634 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(6634)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 6672 line=6672 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(6672)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 6674 line=6674 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(6674)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 6730 line=6730 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(6730)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 6739 line=6739 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(6739)/wave(4) contributes to the likelihood as 1-p-13 .
+Warning! Death for individual 6744 line=6744 occurred at 12/1989 before last wave 4 interviewed at 7/1990 and should have been coded as death instead of '-1'. This case (6744)/wave (4) is contributing to likelihood.
+Information! Unknown status for individual 6860 line=6860 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(6860)/wave(4) contributes to the likelihood as 1-p-13 .
+Error! Death for individual 6925 line=6925 occurred at 10/1990 after last wave 4 interviewed at 7/1990. Potential bias if other individuals are still alive at this date but ignored. This case (6925)/wave (4) is skipped, no contribution to likelihood. Please add a new fictive wave at the date of last vital status scan, with a dead status or alive but unknown state status (-1). See documentation
-Powell iter=8 -2*LL=46537.933555419499 20 sec. 145 sec. 1 -12.288754155703 2 0.092753643773 3 -10.588920461186 4 0.060230465648 5 -1.588742425574 6 -0.034884297842 7 -5.460441275653 8 0.015679980832
-1..........2................3..........4.........5...........6..............7...........8............Gaining to use new average direction of P0 P8 instead of biggest increase direction 1 :
+Information! Unknown status for individual 6961 line=6961 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(6961)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 6967 line=6967 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(6967)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 6999 line=6999 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(6999)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 7032 line=7032 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(7032)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 7087 line=7087 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(7087)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 7097 line=7097 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(7097)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 7181 line=7181 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(7181)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 7209 line=7209 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(7209)/wave(4) contributes to the likelihood as 1-p-13 .
+Warning! Death for individual 7229 line=7229 occurred at 12/1989 before last wave 4 interviewed at 6/1990 and should have been coded as death instead of '-1'. This case (7229)/wave (4) is contributing to likelihood.
+Information! Unknown status for individual 7233 line=7233 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(7233)/wave(4) contributes to the likelihood as 1-p-13 .
+Warning! Death for individual 7273 line=7273 occurred at 1/1989 before last wave 4 interviewed at 7/1990 and should have been coded as death instead of '-1'. This case (7273)/wave (4) is contributing to likelihood.
+Information! Unknown status for individual 7290 line=7290 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(7290)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 7294 line=7294 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(7294)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 7319 line=7319 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(7319)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 7326 line=7326 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(7326)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 7422 line=7422 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(7422)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 7427 line=7427 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(7427)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 7480 line=7480 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(7480)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 7493 line=7493 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(7493)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 7508 line=7508 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(7508)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 7515 line=7515 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(7515)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 7552 line=7552 occurred at last wave 4 at known date 6/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(7552)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 7598 line=7598 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(7598)/wave(4) contributes to the likelihood as 1-p-13 .
+Error! Death for individual 7636 line=7636 occurred at 6/1991 after last wave 4 interviewed at 7/1990. Potential bias if other individuals are still alive at this date but ignored. This case (7636)/wave (4) is skipped, no contribution to likelihood. Please add a new fictive wave at the date of last vital status scan, with a dead status or alive but unknown state status (-1). See documentation
-Powell iter=9 -2*LL=46537.891454595432 19 sec. 164 sec. 1 -12.283266439012 2 0.092658255416 3 -10.562891684762 4 0.059948818330 5 -1.581316871784 6 -0.034987633432 7 -5.472744984969 8 0.015801808619
-1.....2..............3........4............5............6..............7...........8...........Gaining to use new average direction of P0 P8 instead of biggest increase direction 8 :
+Information! Unknown status for individual 7639 line=7639 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(7639)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 7724 line=7724 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(7724)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 7763 line=7763 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(7763)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 7823 line=7823 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(7823)/wave(4) contributes to the likelihood as 1-p-13 .
+Warning! Death for individual 7827 line=7827 occurred at 7/1988 before last wave 4 interviewed at 7/1990 and should have been coded as death instead of '-1'. This case (7827)/wave (4) is contributing to likelihood.
+Error! Death for individual 7832 line=7832 occurred at 9/1990 after last wave 4 interviewed at 7/1990. Potential bias if other individuals are still alive at this date but ignored. This case (7832)/wave (4) is skipped, no contribution to likelihood. Please add a new fictive wave at the date of last vital status scan, with a dead status or alive but unknown state status (-1). See documentation
-Powell iter=10 -2*LL=46537.765677731819 18 sec. 182 sec. 1 -12.262794789308 2 0.092367916069 3 -10.470813624311 4 0.058882333788 5 -1.562236459514 6 -0.035239826239 7 -5.510520845667 8 0.016209203343
-1......2.......................3.......4........5............6..............7............8............Gaining to use new average direction of P0 P8 instead of biggest increase direction 8 :
+Information! Unknown status for individual 7863 line=7863 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(7863)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 7873 line=7873 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(7873)/wave(4) contributes to the likelihood as 1-p-13 .
+Warning! Death for individual 7905 line=7905 occurred at 3/1990 before last wave 4 interviewed at 7/1990 and should have been coded as death instead of '-1'. This case (7905)/wave (4) is contributing to likelihood.
+Warning! Death for individual 7976 line=7976 occurred at 9/1989 before last wave 4 interviewed at 7/1990 and should have been coded as death instead of '-1'. This case (7976)/wave (4) is contributing to likelihood.
+Information! Unknown status for individual 7985 line=7985 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(7985)/wave(4) contributes to the likelihood as 1-p-13 .
+Warning! Death for individual 8043 line=8043 occurred at 10/1989 before last wave 4 interviewed at 6/1990 and should have been coded as death instead of '-1'. This case (8043)/wave (4) is contributing to likelihood.
+Information! Unknown status for individual 8069 line=8069 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(8069)/wave(4) contributes to the likelihood as 1-p-13 .
+Error! Death for individual 8130 line=8130 occurred at 8/1990 after last wave 4 interviewed at 7/1990. Potential bias if other individuals are still alive at this date but ignored. This case (8130)/wave (4) is skipped, no contribution to likelihood. Please add a new fictive wave at the date of last vital status scan, with a dead status or alive but unknown state status (-1). See documentation
-Powell iter=11 -2*LL=46537.498821119196 19 sec. 201 sec. 1 -12.214937649099 2 0.091832303757 3 -10.310758734665 4 0.056891524738 5 -1.540234661352 6 -0.035490430955 7 -5.566745000929 8 0.016892295483
-1.....2..............3.......4........5.................6...............7..............8.......
-Powell iter=12 -2*LL=46537.411404609535 18 sec. 219 sec. 1 -12.203395121992 2 0.091768319083 3 -10.300995242801 4 0.056661013905 5 -1.545167229526 6 -0.035400936285 7 -5.565324661507 8 0.016938176004
-1............2..................3..........4............5..................6..............7.............8.....................Gaining to use new average direction of P0 P8 instead of biggest increase direction 3 :
+Information! Unknown status for individual 8145 line=8145 occurred at last wave 4 at known date 7/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(8145)/wave(4) contributes to the likelihood as 1-p-13 .
+Information! Unknown status for individual 8178 line=8178 occurred at last wave 4 at known date 8/1990. Please, check if your unknown date of death 99/9999 means a live state -1 at wave 4. This case(8178)/wave(4) contributes to the likelihood as 1-p-13 .
+Error! Death for individual 8257 line=8257 occurred at 6/1991 after last wave 4 interviewed at 7/1990. Potential bias if other individuals are still alive at this date but ignored. This case (8257)/wave (4) is skipped, no contribution to likelihood. Please add a new fictive wave at the date of last vital status scan, with a dead status or alive but unknown state status (-1). See documentation
-Powell iter=13 -2*LL=46537.410703925241 22 sec. 241 sec. 1 -12.202513270154 2 0.091761198901 3 -10.299418789414 4 0.056628413435 5 -1.545343667372 6 -0.035398382663 7 -5.565921175501 8 0.016950736226
-1.........2................3............4.............5..........6....................7..............8........#Number of iterations & function calls = 13 & 1694, -2 Log likelihood = 46537.410262065809
+Delay (in months) between two waves Min=1 (for indiviudal 6017) Max=74 (5168) Mean=24.058442
+
+ Age 70 1.=665 loss[1]=5.1% 2.=0 loss[2]=NaNQ% 1.=631 prev[1]=100.0% 2.=0 prev[2]=0.0% 1-1=34 11=563 12=33 13=35
+Age 71 1.=663 loss[1]=5.7% 2.=2 loss[2]=0.0% 1.=625 prev[1]=99.7% 2.=2 prev[2]=0.3% 1-1=38 11=530 12=47 13=48 21=1 22=1
+Age 72 1.=1181 loss[1]=5.6% 2.=35 loss[2]=5.7% 1.=1115 prev[1]=97.1% 2.=33 prev[2]=2.9% 1-1=66 11=969 12=70 13=76 2-1=2 21=9 22=18 23=6
+Age 73 1.=1119 loss[1]=5.1% 2.=47 loss[2]=6.4% 1.=1062 prev[1]=96.0% 2.=44 prev[2]=4.0% 1-1=57 11=912 12=71 13=79 2-1=3 21=12 22=26 23=6
+Age 74 1.=1562 loss[1]=4.2% 2.=90 loss[2]=4.4% 1.=1497 prev[1]=94.6% 2.=86 prev[2]=5.4% 1-1=65 11=1277 12=104 13=116 2-1=4 21=30 22=41 23=15
+Age 75 1.=1461 loss[1]=6.1% 2.=97 loss[2]=6.2% 1.=1372 prev[1]=93.8% 2.=91 prev[2]=6.2% 1-1=89 11=1155 12=98 13=119 2-1=6 21=15 22=55 23=21
+Age 76 1.=1388 loss[1]=4.3% 2.=111 loss[2]=2.7% 1.=1328 prev[1]=92.5% 2.=108 prev[2]=7.5% 1-1=60 11=1090 12=111 13=127 2-1=3 21=25 22=54 23=29
+Age 77 1.=1270 loss[1]=4.9% 2.=100 loss[2]=5.0% 1.=1208 prev[1]=92.7% 2.=95 prev[2]=7.3% 1-1=62 11=962 12=131 13=115 2-1=5 21=24 22=52 23=19
+Age 78 1.=1163 loss[1]=4.6% 2.=122 loss[2]=4.1% 1.=1109 prev[1]=90.5% 2.=117 prev[2]=9.5% 1-1=54 11=894 12=111 13=104 2-1=5 21=21 22=75 23=21
+Age 79 1.=1045 loss[1]=5.7% 2.=131 loss[2]=3.8% 1.=985 prev[1]=88.7% 2.=126 prev[2]=11.3% 1-1=60 11=773 12=109 13=103 2-1=5 21=31 22=68 23=27
+Age 80 1.=986 loss[1]=4.1% 2.=128 loss[2]=3.1% 1.=946 prev[1]=88.4% 2.=124 prev[2]=11.6% 1-1=40 11=718 12=140 13=88 2-1=4 21=21 22=71 23=32
+Age 81 1.=870 loss[1]=4.7% 2.=130 loss[2]=3.8% 1.=829 prev[1]=86.9% 2.=125 prev[2]=13.1% 1-1=41 11=629 12=102 13=98 2-1=5 21=24 22=74 23=27
+Age 82 1.=799 loss[1]=3.9% 2.=145 loss[2]=3.4% 1.=768 prev[1]=84.6% 2.=140 prev[2]=15.4% 1-1=31 11=546 12=131 13=91 2-1=5 21=15 22=77 23=48
+Age 83 1.=722 loss[1]=6.0% 2.=115 loss[2]=0.0% 1.=679 prev[1]=85.5% 2.=115 prev[2]=14.5% 1-1=43 11=456 12=124 13=99 21=19 22=71 23=25
+Age 84 1.=620 loss[1]=5.3% 2.=147 loss[2]=4.1% 1.=587 prev[1]=80.6% 2.=141 prev[2]=19.4% 1-1=33 11=394 12=112 13=81 2-1=6 21=18 22=82 23=41
+Age 85 1.=526 loss[1]=4.4% 2.=126 loss[2]=5.6% 1.=503 prev[1]=80.9% 2.=119 prev[2]=19.1% 1-1=23 11=315 12=122 13=66 2-1=7 21=13 22=67 23=39
+Age 86 1.=464 loss[1]=5.2% 2.=128 loss[2]=4.7% 1.=440 prev[1]=78.3% 2.=122 prev[2]=21.7% 1-1=24 11=271 12=103 13=66 2-1=6 21=18 22=70 23=34
+Age 87 1.=370 loss[1]=4.3% 2.=143 loss[2]=2.8% 1.=354 prev[1]=71.8% 2.=139 prev[2]=28.2% 1-1=16 11=219 12=79 13=56 2-1=4 21=15 22=71 23=53
+Age 88 1.=283 loss[1]=3.9% 2.=163 loss[2]=6.1% 1.=272 prev[1]=64.0% 2.=153 prev[2]=36.0% 1-1=11 11=172 12=63 13=37 2-1=10 21=25 22=84 23=44
+Age 89 1.=226 loss[1]=3.1% 2.=151 loss[2]=4.0% 1.=219 prev[1]=60.2% 2.=145 prev[2]=39.8% 1-1=7 11=135 12=45 13=39 2-1=6 21=21 22=73 23=51
+Age 90 1.=181 loss[1]=3.9% 2.=143 loss[2]=4.9% 1.=174 prev[1]=56.1% 2.=136 prev[2]=43.9% 1-1=7 11=104 12=41 13=29 2-1=7 21=18 22=82 23=36
+Age 91 1.=146 loss[1]=4.1% 2.=107 loss[2]=6.5% 1.=140 prev[1]=58.3% 2.=100 prev[2]=41.7% 1-1=6 11=78 12=36 13=26 2-1=7 21=11 22=51 23=38
+Age 92 1.=114 loss[1]=6.1% 2.=101 loss[2]=2.0% 1.=107 prev[1]=51.9% 2.=99 prev[2]=48.1% 1-1=7 11=47 12=43 13=17 2-1=2 21=11 22=52 23=36
+Age 93 1.=78 loss[1]=3.8% 2.=88 loss[2]=11.4% 1.=75 prev[1]=49.0% 2.=78 prev[2]=51.0% 1-1=3 11=36 12=25 13=14 2-1=10 21=10 22=46 23=22
+Age 94 1.=50 loss[1]=2.0% 2.=70 loss[2]=0.0% 1.=49 prev[1]=41.2% 2.=70 prev[2]=58.8% 1-1=1 11=24 12=8 13=17 21=4 22=39 23=27
+Age 95 1.=38 loss[1]=5.3% 2.=56 loss[2]=3.6% 1.=36 prev[1]=40.0% 2.=54 prev[2]=60.0% 1-1=2 11=16 12=13 13=7 2-1=2 21=4 22=31 23=19
+Age 96 1.=25 loss[1]=4.0% 2.=42 loss[2]=2.4% 1.=24 prev[1]=36.9% 2.=41 prev[2]=63.1% 1-1=1 11=11 12=9 13=4 2-1=1 21=2 22=22 23=17
+Age 97 1.=11 loss[1]=0.0% 2.=36 loss[2]=8.3% 1.=11 prev[1]=25.0% 2.=33 prev[2]=75.0% 11=4 12=7 2-1=3 21=2 22=22 23=9
+Age 98 1.=5 loss[1]=0.0% 2.=21 loss[2]=0.0% 1.=5 prev[1]=19.2% 2.=21 prev[2]=80.8% 11=1 12=4 21=1 22=13 23=7
+Age 99 1.=6 loss[1]=0.0% 2.=14 loss[2]=0.0% 1.=6 prev[1]=30.0% 2.=14 prev[2]=70.0% 11=3 12=1 13=2 21=2 22=2 23=10
+Age 100 1.=2 loss[1]=0.0% 2.=9 loss[2]=0.0% 1.=2 prev[1]=18.2% 2.=9 prev[2]=81.8% 11=1 12=1 21=2 22=5 23=2
+Age 101 1.=1 loss[1]=0.0% 2.=0 loss[2]=NaNQ% 1.=1 prev[1]=100.0% 2.=0 prev[2]=0.0% 12=1
+Age 102 1.=2 loss[1]=0.0% 2.=1 loss[2]=0.0% 1.=2 prev[1]=66.7% 2.=1 prev[2]=33.3% 12=2 22=1
+Age 103 1.=0 loss[1]=NaNQ% 2.=0 loss[2]=NaNQ% 1.=0 prev[1]=NaNQ% 2.=0 prev[2]=NaNQ%
+Age 104 1.=0 loss[1]=NaNQ% 2.=0 loss[2]=NaNQ% 1.=0 prev[1]=NaNQ% 2.=0 prev[2]=NaNQ%
+1 1.=0 loss[1]=NaNQ% 2.=0 loss[2]=NaNQ% 1.=0 prev[1]=NaNQ% 2.=0 prev[2]=NaNQ%
+0 1.=0 loss[1]=NaNQ% 2.=0 loss[2]=NaNQ% 1.=0 prev[1]=NaNQ% 2.=0 prev[2]=NaNQ%
+Total 1.=18042 loss[1]=4.9% 2.=2799 loss[2]=4.2% 1.=17161 prev[1]=86.5% 2.=2681 prev[2]=13.5% 1-1=881 11=13305 12=2097 13=1759 2-1=118 21=424 22=1496 23=761
+
+
+12 -12.3175753 ln(2097/13305)= -1.8476322
+13 -10.1631670 ln(1759/13305)= -2.0233944
+21 -1.6685685 ln(424/1496)= -1.2608167
+23 -5.6869201 ln(761/1496)= -0.6759168
+
+ -1-1=0 -10=0 -11=0 -12=0 -13=0 0-1=0 00=0 01=0 02=0 03=0 1-1=881 10=0 11=13305 12=2097 13=1759 2-1=118 20=0 21=424 22=1496 23=761 3-1=0 30=0 31=0 32=0 33=0
+Powell
+
+Powell iter=1 -2*LL=46617.105944340874 0 sec. 10 sec.
+#model= 1 + age
+12 -12.3175753 0.0931122
+13 -10.1631670 0.0548384
+21 -1.6685685 -0.0339056
+23 -5.6869201 0.0181350
+1................2..................3.................4.......................5.........................6...................7................8......................
+Powell iter=2 -2*LL=46617.105944308234 13 sec. 23 sec.
+#model= 1 + age
+12 -12.3175784 0.0931122
+13 -10.1631705 0.0548384
+21 -1.6685685 -0.0339056
+23 -5.6869199 0.0181350
+
+Considering the time needed for this last iteration #2: 13 seconds,
+ - if your program needs 10 iterations to converge, convergence will be
+ reached in 0 day(s) 0 hour(s) 1 minute(s) 44 second(s) i.e.
+ on Wed May 22 22:34:43 2019 (current time is Wed May 22 22:32:59 2019);
+ - if your program needs 20 iterations to converge, convergence will be
+ reached in 0 day(s) 0 hour(s) 3 minute(s) 54 second(s) i.e.
+ on Wed May 22 22:36:53 2019 (current time is Wed May 22 22:32:59 2019);
+ - if your program needs 30 iterations to converge, convergence will be
+ reached in 0 day(s) 0 hour(s) 6 minute(s) 4 second(s) i.e.
+ on Wed May 22 22:39:03 2019 (current time is Wed May 22 22:32:59 2019);
+1.................................2......................3.........................4..............................5...............................6....................7.................................8.......................
+#Number of iterations & function calls = 2 & 424, -2 Log likelihood = 46617.105944307536
# Parameters nlstate*nlstate*ncov a12*1 + b12 * age + ...
-12 -12.2027105 0.0917637
-13 -10.2992383 0.0566262
-21 -1.5456185 -0.0353950
-23 -5.5659628 0.0169512
+12 -12.3175784 0.0931122
+13 -10.1631705 0.0548384
+21 -1.6685685 -0.0339056
+23 -5.6869199 0.0181350
Calculation of the hessian matrix. Wait...
-12345678.12.13.14.15.16.17.18.23.24.25.26.27.28.34.35.36.37.38.45.46.47.48.56.57.58.67.68.78
+1-2-3-4-5-6-7-8-.1-2.1-3.1-4.1-5.1-6.1-7.1-8.2-3.2-4.2-5.2-6.2-7.2-8.3-4.3-5.3-6.3-7.3-8.4-5.4-6.4-7.4-8.5-6.5-7.5-8.6-7.6-8.7-8
Inverting the hessian to get the covariance matrix. Wait...
#Hessian matrix#
-2.301e+03 1.898e+05 3.980e+02 3.280e+04 -3.628e+02 -3.007e+04 -3.111e+02 -2.563e+04
-1.898e+05 1.574e+07 3.252e+04 2.694e+06 -3.003e+04 -2.503e+06 -2.554e+04 -2.114e+06
-3.980e+02 3.252e+04 8.499e+02 6.839e+04 -8.040e+01 -6.688e+03 3.773e+02 3.120e+04
-3.280e+04 2.694e+06 6.839e+04 5.532e+06 -6.734e+03 -5.639e+05 3.117e+04 2.590e+06
--3.628e+02 -3.003e+04 -8.040e+01 -6.734e+03 4.570e+02 3.818e+04 4.994e+01 4.126e+03
--3.007e+04 -2.503e+06 -6.688e+03 -5.639e+05 3.818e+04 3.208e+06 4.126e+03 3.421e+05
--3.111e+02 -2.554e+04 3.773e+02 3.117e+04 4.994e+01 4.126e+03 1.024e+03 8.867e+04
--2.563e+04 -2.114e+06 3.120e+04 2.590e+06 4.126e+03 3.421e+05 8.867e+04 7.720e+06
-Parameters, T and confidence intervals
-12 -12.2027105 T= -36.632 CI=[ -12.8689514 ; -11.5364696] 0.0917637 T= 22.911 CI=[ 0.0837533 ; 0.0997740]
-13 -10.2992383 T= -18.738 CI=[ -11.3985038 ; -9.1999728] 0.0566262 T= 8.316 CI=[ 0.0430080 ; 0.0702444]
-21 -1.5456185 T= -2.259 CI=[ -2.9142844 ; -0.1769526] -0.0353950 T= -4.335 CI=[ -0.0517251 ; -0.0190649]
-23 -5.5659628 T= -11.196 CI=[ -6.5602821 ; -4.5716434] 0.0169512 T= 2.998 CI=[ 0.0056422 ; 0.0282603]
+2.296930e+03 1.895497e+05 3.935378e+02 3.244270e+04 -3.627633e+02 -3.009693e+04 -3.043646e+02 -2.506988e+04
+1.895497e+05 1.572740e+07 3.217180e+04 2.665640e+06 -3.004624e+04 -2.507644e+06 -2.497200e+04 -2.067052e+06
+3.935378e+02 3.217180e+04 8.583265e+02 6.900254e+04 -7.950011e+01 -6.614852e+03 3.761395e+02 3.110935e+04
+3.244270e+04 2.665640e+06 6.900254e+04 5.575042e+06 -6.660064e+03 -5.573884e+05 3.108386e+04 2.584342e+06
+-3.627633e+02 -3.004624e+04 -7.950011e+01 -6.660064e+03 4.568873e+02 3.819631e+04 4.751955e+01 3.915572e+03
+-3.009693e+04 -2.507644e+06 -6.614852e+03 -5.573884e+05 3.819631e+04 3.210369e+06 3.920001e+03 3.247002e+05
+-3.043646e+02 -2.497200e+04 3.761395e+02 3.108386e+04 4.751955e+01 3.920001e+03 1.019065e+03 8.831012e+04
+-2.506988e+04 -2.067052e+06 3.110935e+04 2.584342e+06 3.915572e+03 3.247002e+05 8.831012e+04 7.695216e+06
+Parameters, Wald tests and Wald-based confidence intervals
+ W is simply the result of the division of the parameter by the square root of covariance of the parameter.
+ And Wald-based confidence intervals plus and minus 1.96 * W
+ It might be better to visualize the covariance matrix. See the page 'Matrix of variance-covariance of one-step probabilities' and its graphs.
+12 -12.3175784 W= -37.237 CI=[ -12.9659254 ; -11.6692314] 0.0931122 W= 23.405 CI=[ 0.0853146 ; 0.1009098]
+13 -10.1631705 W= -18.488 CI=[ -11.2406101 ; -9.0857309] 0.0548384 W= 8.045 CI=[ 0.0414779 ; 0.0681989]
+21 -1.6685685 W= -2.395 CI=[ -3.0342962 ; -0.3028408] -0.0339056 W= -4.080 CI=[ -0.0501948 ; -0.0176164]
+23 -5.6869199 W= -11.373 CI=[ -6.6669549 ; -4.7068849] 0.0181350 W= 3.189 CI=[ 0.0069875 ; 0.0292825]
# Scales (for hessian or gradient estimation)
-12 1.00000e-04 1.00000e-06
-13 1.00000e-04 1.00000e-06
+12 2.00000e-04 2.00000e-06
+13 3.00000e-04 4.00000e-06
21 1.00000e-04 1.00000e-06
-23 3.00000e-04 1.00000e-06
+23 3.00000e-04 3.00000e-06
# Covariance matrix
# 121 Var(a12)
# 122 Cov(b12,a12) Var(b12)
#212 Cov(b21,a12) Cov(b21,b12) Cov(b21,a13) Cov(b21,b13) Cov(b21,a21) Var(b21)
#231 Cov(a23,a12) Cov(a23,b12) Cov(a23,a13) Cov(a23,b13) Cov(a23,a21) Cov(a23,b21) Var(a23)
#232 Cov(b23,a12) Cov(b23,b12) Cov(b23,a13) Cov(b23,b13) Cov(b23,a21) Cov(b23,b21) Cov(b23,a23) Var(b23)
-121 1.10969e-01
-122 -1.33029e-03 1.60415e-05
-131 -5.53768e-02 6.77306e-04 3.02096e-01
-132 6.58106e-04 -8.12775e-06 -3.73026e-03 4.63642e-05
-211 7.59223e-02 -9.02778e-04 2.53226e-02 -3.27358e-04 4.68312e-01
-212 -9.05016e-04 1.08138e-05 -3.15929e-04 4.09141e-06 -5.57262e-03 6.66682e-05
-231 5.24440e-02 -5.95206e-04 -8.93872e-02 1.00357e-03 6.06408e-03 -7.39139e-05 2.47168e-01
-232 -5.95756e-04 6.80630e-06 1.05928e-03 -1.20524e-05 -6.07266e-05 7.34461e-07 -2.80315e-03 3.19737e-05
+121 1.0942155e-01
+122 -1.31211e-03 1.5827559e-05
+131 -5.42260e-02 6.63610e-04 3.0218559e-01
+132 6.46897e-04 -7.99273e-06 -3.73518e-03 4.6465922e-05
+211 8.05586e-02 -9.61087e-04 1.88983e-02 -2.43008e-04 4.8553004e-01
+212 -9.60482e-04 1.15119e-05 -2.38550e-04 3.07513e-06 -5.77598e-03 6.9069819e-05
+231 4.89528e-02 -5.53986e-04 -9.43826e-02 1.06862e-03 -1.83095e-03 2.05848e-05 2.5001784e-01
+232 -5.56252e-04 6.33917e-06 1.11795e-03 -1.28139e-05 2.71734e-05 -3.16845e-07 -2.83579e-03 3.2347405e-05
+# agemin agemax for life expectancy, bage fage (if mle==0 ie no data nor Max likelihood).\r
begin-prev-date=1/1/1984 end-prev-date=1/6/1988 mov_average=0
-pop_based=0
-prevforecast=0 starting-proj-date=1/1/2000 final-proj-date=1/1/2000 mobil_average=0
-Computing period (stable) prevalence: result on file 'plrbiaspar.txt'
+pop_based=1
+prevforecast=1 starting-proj-date=1/1/1989 final-proj-date=1/1/2000 mobil_average=0
+#backcast=1 starting-back-date=1/1/1990 final-back-date=1/1/1970 mobil_average=1\r
+# Results model=1+age.\r
+result: .\r
+Computing forward period (stable) prevalence: result on file 'PL_biaspar.txt'
#************
-Computing pij: result on file 'pijrbiaspar.txt'
-Computing standard deviation of one-step probabilities: result on file 'probrbiaspar.txt'
-Computing matrix of variance covariance of one-step probabilities: result on file 'probcovrbiaspar.txt'
-and correlation matrix of one-step probabilities: result on file 'probcorrbiaspar.txt'
-65 13-12 mu 1.5972e-02 2.3367e-02 Var 3.3236e-06 3.1628e-06 cor -0.338 cov -1.0950e-06 Eig 4.341e-06 2.145e-06 1stv 0.733 -0.681 tan -0.929
-70 13-12 mu 2.1166e-02 3.6913e-02 Var 3.1497e-06 4.5213e-06 cor -0.360 cov -1.3581e-06 Eig 5.357e-06 2.314e-06 1stv 0.524 -0.852 tan -1.625
-75 13-12 mu 2.8027e-02 5.8265e-02 Var 2.6312e-06 5.5954e-06 cor -0.404 cov -1.5492e-06 Eig 6.257e-06 1.969e-06 1stv 0.393 -0.920 tan -2.341
-80 13-12 mu 3.7067e-02 9.1858e-02 Var 2.7234e-06 6.5867e-06 cor -0.456 cov -1.9326e-06 Eig 7.387e-06 1.923e-06 1stv 0.383 -0.924 tan -2.413
-85 13-12 mu 4.8930e-02 1.4455e-01 Var 7.0008e-06 1.4795e-05 cor -0.418 cov -4.2498e-06 Eig 1.666e-05 5.132e-06 1stv 0.403 -0.915 tan -2.274
-90 13-12 mu 6.4406e-02 2.2681e-01 Var 2.5605e-05 7.2563e-05 cor -0.359 cov -1.5484e-05 Eig 7.721e-05 2.096e-05 1stv 0.287 -0.958 tan -3.333
-95 13-12 mu 8.4408e-02 3.5434e-01 Var 8.3580e-05 3.5679e-04 cor -0.341 cov -5.8808e-05 Eig 3.689e-04 7.146e-05 1stv 0.202 -0.979 tan -4.852
-65 21-12 mu 2.4815e-01 2.3367e-02 Var 1.5081e-03 3.1628e-06 cor 0.333 cov 2.3009e-05 Eig 1.509e-03 2.811e-06 1stv 1.000 0.015 tan 0.015
-70 21-12 mu 2.0839e-01 3.6913e-02 Var 6.2123e-04 4.5213e-06 cor 0.332 cov 1.7587e-05 Eig 6.217e-04 4.020e-06 1stv 1.000 0.028 tan 0.028
-75 21-12 mu 1.7490e-01 5.8265e-02 Var 2.2058e-04 5.5954e-06 cor 0.326 cov 1.1456e-05 Eig 2.212e-04 4.987e-06 1stv 0.999 0.053 tan 0.053
-80 21-12 mu 1.4670e-01 9.1858e-02 Var 7.0767e-05 6.5867e-06 cor 0.304 cov 6.5593e-06 Eig 7.143e-05 5.923e-06 1stv 0.995 0.101 tan 0.101
-85 21-12 mu 1.2299e-01 1.4455e-01 Var 3.9227e-05 1.4795e-05 cor 0.276 cov 6.6547e-06 Eig 4.092e-05 1.310e-05 1stv 0.969 0.247 tan 0.255
-90 21-12 mu 1.0307e-01 2.2681e-01 Var 5.4893e-05 7.2563e-05 cor 0.289 cov 1.8261e-05 Eig 8.401e-05 4.344e-05 1stv 0.531 0.847 tan 1.595
-95 21-12 mu 8.6343e-02 3.5434e-01 Var 8.2297e-05 3.5679e-04 cor 0.302 cov 5.1723e-05 Eig 3.662e-04 7.287e-05 1stv 0.179 0.984 tan 5.489
-65 23-12 mu 1.3378e-01 2.3367e-02 Var 3.1221e-04 3.1628e-06 cor 0.364 cov 1.1447e-05 Eig 3.126e-04 2.739e-06 1stv 0.999 0.037 tan 0.037
-70 23-12 mu 1.4596e-01 3.6913e-02 Var 2.3686e-04 4.5213e-06 cor 0.388 cov 1.2697e-05 Eig 2.376e-04 3.829e-06 1stv 0.999 0.054 tan 0.054
-75 23-12 mu 1.5914e-01 5.8265e-02 Var 1.6142e-04 5.5954e-06 cor 0.423 cov 1.2709e-05 Eig 1.625e-04 4.566e-06 1stv 0.997 0.081 tan 0.081
-80 23-12 mu 1.7343e-01 9.1858e-02 Var 9.6292e-05 6.5867e-06 cor 0.449 cov 1.1306e-05 Eig 9.769e-05 5.184e-06 1stv 0.992 0.123 tan 0.124
-85 23-12 mu 1.8890e-01 1.4455e-01 Var 5.6832e-05 1.4795e-05 cor 0.356 cov 1.0309e-05 Eig 5.922e-05 1.240e-05 1stv 0.974 0.226 tan 0.232
-90 23-12 mu 2.0566e-01 2.2681e-01 Var 6.4965e-05 7.2563e-05 cor 0.255 cov 1.7485e-05 Eig 8.666e-05 5.087e-05 1stv 0.628 0.779 tan 1.241
-95 23-12 mu 2.2383e-01 3.5434e-01 Var 1.5120e-04 3.5679e-04 cor 0.238 cov 5.5313e-05 Eig 3.707e-04 1.373e-04 1stv 0.244 0.970 tan 3.969
-65 21-13 mu 2.4815e-01 1.5972e-02 Var 1.5081e-03 3.3236e-06 cor 0.047 cov 3.3259e-06 Eig 1.508e-03 3.316e-06 1stv 1.000 0.002 tan 0.002
-70 21-13 mu 2.0839e-01 2.1166e-02 Var 6.2123e-04 3.1497e-06 cor 0.037 cov 1.6580e-06 Eig 6.212e-04 3.145e-06 1stv 1.000 0.003 tan 0.003
-75 21-13 mu 1.7490e-01 2.8027e-02 Var 2.2058e-04 2.6312e-06 cor 0.024 cov 5.6775e-07 Eig 2.206e-04 2.630e-06 1stv 1.000 0.003 tan 0.003
-80 21-13 mu 1.4670e-01 3.7067e-02 Var 7.0767e-05 2.7234e-06 cor 0.024 cov 3.3537e-07 Eig 7.077e-05 2.722e-06 1stv 1.000 0.005 tan 0.005
-85 21-13 mu 1.2299e-01 4.8930e-02 Var 3.9227e-05 7.0008e-06 cor 0.078 cov 1.2860e-06 Eig 3.928e-05 6.950e-06 1stv 0.999 0.040 tan 0.040
-90 21-13 mu 1.0307e-01 6.4406e-02 Var 5.4893e-05 2.5605e-05 cor 0.100 cov 3.7654e-06 Eig 5.537e-05 2.513e-05 1stv 0.992 0.126 tan 0.127
-95 21-13 mu 8.6343e-02 8.4408e-02 Var 8.2297e-05 8.3580e-05 cor 0.097 cov 8.0750e-06 Eig 9.104e-05 7.484e-05 1stv 0.679 0.735 tan 1.083
-65 23-13 mu 1.3378e-01 1.5972e-02 Var 3.1221e-04 3.3236e-06 cor -0.409 cov -1.3182e-05 Eig 3.128e-04 2.762e-06 1stv 0.999 -0.043 tan -0.043
-70 23-13 mu 1.4596e-01 2.1166e-02 Var 2.3686e-04 3.1497e-06 cor -0.453 cov -1.2362e-05 Eig 2.375e-04 2.498e-06 1stv 0.999 -0.053 tan -0.053
-75 23-13 mu 1.5914e-01 2.8027e-02 Var 1.6142e-04 2.6312e-06 cor -0.528 cov -1.0874e-05 Eig 1.622e-04 1.890e-06 1stv 0.998 -0.068 tan -0.068
-80 23-13 mu 1.7343e-01 3.7067e-02 Var 9.6292e-05 2.7234e-06 cor -0.585 cov -9.4784e-06 Eig 9.724e-05 1.773e-06 1stv 0.995 -0.100 tan -0.100
-85 23-13 mu 1.8890e-01 4.8930e-02 Var 5.6832e-05 7.0008e-06 cor -0.513 cov -1.0241e-05 Eig 5.885e-05 4.978e-06 1stv 0.981 -0.194 tan -0.198
-90 23-13 mu 2.0566e-01 6.4406e-02 Var 6.4965e-05 2.5605e-05 cor -0.434 cov -1.7707e-05 Eig 7.176e-05 1.881e-05 1stv 0.934 -0.358 tan -0.384
-95 23-13 mu 2.2383e-01 8.4408e-02 Var 1.5120e-04 8.3580e-05 cor -0.363 cov -4.0852e-05 Eig 1.704e-04 6.436e-05 1stv 0.905 -0.426 tan -0.470
-65 23-21 mu 1.3378e-01 2.4815e-01 Var 3.1221e-04 1.5081e-03 cor -0.015 cov -1.0342e-05 Eig 1.508e-03 3.121e-04 1stv 0.009 -1.000 tan -115.645
-70 23-21 mu 1.4596e-01 2.0839e-01 Var 2.3686e-04 6.2123e-04 cor -0.013 cov -4.8257e-06 Eig 6.213e-04 2.368e-04 1stv 0.013 -1.000 tan -79.662
-75 23-21 mu 1.5914e-01 1.7490e-01 Var 1.6142e-04 2.2058e-04 cor -0.014 cov -2.7177e-06 Eig 2.207e-04 1.613e-04 1stv 0.046 -0.999 tan -21.813
-80 23-21 mu 1.7343e-01 1.4670e-01 Var 9.6292e-05 7.0767e-05 cor -0.029 cov -2.3931e-06 Eig 9.651e-05 7.054e-05 1stv 0.996 -0.093 tan -0.093
-85 23-21 mu 1.8890e-01 1.2299e-01 Var 5.6832e-05 3.9227e-05 cor -0.061 cov -2.8869e-06 Eig 5.729e-05 3.877e-05 1stv 0.987 -0.158 tan -0.160
-90 23-21 mu 2.0566e-01 1.0307e-01 Var 6.4965e-05 5.4893e-05 cor -0.062 cov -3.6726e-06 Eig 6.616e-05 5.370e-05 1stv 0.951 -0.310 tan -0.326
-95 23-21 mu 2.2383e-01 8.6343e-02 Var 1.5120e-04 8.2297e-05 cor -0.040 cov -4.5038e-06 Eig 1.515e-04 8.200e-05 1stv 0.998 -0.065 tan -0.065
-Computing Health Expectancies: result on file 'erbiaspar.txt'
-65|66|67|68|69|70|71|72|73|74|75|76|77|78|79|80|81|82|83|84|85|86|87|88|89|90|91|92|93|94|95|
-Computing Total Life expectancies with their standard errors: file 'trbiaspar.txt'
-Computing Health Expectancies and standard errors: result on file 'stderbiaspar.txt'
-Computing Covar. of Health Expectancies: result on file 'cverbiaspar.txt'
-Computing Variance-covariance of DFLEs: file 'vrbiaspar.txt'
-65|66|67|68|69|70|71|72|73|74|75|76|77|78|79|80|81|82|83|84|85|86|87|88|89|90|91|92|93|94|95|
-Computing total mortality p.j=w1*p1j+w2*p2j+..: result on file 'prmorprev1-stablbased-rbiaspar.txt'
-End of Imach
-Local time at start Tue Aug 18 18:20:44 2015
+Computing pij: result on file 'PIJ_biaspar.txt'
+Computing standard deviation of one-step probabilities: result on file 'PROB_rbiaspar.txt'
+Computing matrix of variance covariance of one-step probabilities: result on file 'PROBCOV_biaspar.txt'
+and correlation matrix of one-step probabilities: result on file 'PROBCOR_biaspar.txt'
+70 13-12 mu 2.1400e-02 3.6166e-02 Var 3.1732e-06 4.2735e-06 cor -0.351 cov -1.2918e-06 Eig 5.127e-06 2.319e-06 1stv 0.551 -0.834 tan -1.513
+75 13-12 mu 2.8084e-02 5.7475e-02 Var 2.5817e-06 5.3630e-06 cor -0.393 cov -1.4626e-06 Eig 5.991e-06 1.954e-06 1stv 0.394 -0.919 tan -2.331
+80 13-12 mu 3.6812e-02 9.1225e-02 Var 2.6231e-06 6.4219e-06 cor -0.446 cov -1.8300e-06 Eig 7.160e-06 1.885e-06 1stv 0.374 -0.927 tan -2.479
+85 13-12 mu 4.8162e-02 1.4452e-01 Var 6.7559e-06 1.4710e-05 cor -0.413 cov -4.1154e-06 Eig 1.646e-05 5.010e-06 1stv 0.391 -0.921 tan -2.357
+90 13-12 mu 6.2827e-02 2.2829e-01 Var 2.4481e-05 7.3017e-05 cor -0.358 cov -1.5149e-05 Eig 7.736e-05 2.014e-05 1stv 0.275 -0.961 tan -3.490
+95 13-12 mu 8.1592e-02 3.5900e-01 Var 7.8620e-05 3.6276e-04 cor -0.340 cov -5.7431e-05 Eig 3.739e-04 6.745e-05 1stv 0.191 -0.982 tan -5.142
+70 21-12 mu 2.0469e-01 3.6166e-02 Var 6.2086e-04 4.2735e-06 cor 0.342 cov 1.7621e-05 Eig 6.214e-04 3.770e-06 1stv 1.000 0.029 tan 0.029
+75 21-12 mu 1.7304e-01 5.7475e-02 Var 2.2257e-04 5.3630e-06 cor 0.332 cov 1.1485e-05 Eig 2.232e-04 4.757e-06 1stv 0.999 0.053 tan 0.053
+80 21-12 mu 1.4621e-01 9.1225e-02 Var 7.1363e-05 6.4219e-06 cor 0.305 cov 6.5203e-06 Eig 7.201e-05 5.774e-06 1stv 0.995 0.099 tan 0.099
+85 21-12 mu 1.2348e-01 1.4452e-01 Var 3.9495e-05 1.4710e-05 cor 0.283 cov 6.8331e-06 Eig 4.125e-05 1.295e-05 1stv 0.968 0.249 tan 0.257
+90 21-12 mu 1.0424e-01 2.2829e-01 Var 5.6829e-05 7.3017e-05 cor 0.307 cov 1.9770e-05 Eig 8.629e-05 4.356e-05 1stv 0.557 0.830 tan 1.490
+95 21-12 mu 8.7961e-02 3.5900e-01 Var 8.7335e-05 3.6276e-04 cor 0.322 cov 5.7307e-05 Eig 3.742e-04 7.589e-05 1stv 0.196 0.981 tan 5.006
+70 23-12 mu 1.4061e-01 3.6166e-02 Var 2.2238e-04 4.2735e-06 cor 0.368 cov 1.1339e-05 Eig 2.230e-04 3.686e-06 1stv 0.999 0.052 tan 0.052
+75 23-12 mu 1.5420e-01 5.7475e-02 Var 1.5306e-04 5.3630e-06 cor 0.405 cov 1.1590e-05 Eig 1.540e-04 4.459e-06 1stv 0.997 0.078 tan 0.078
+80 23-12 mu 1.6900e-01 9.1225e-02 Var 9.2084e-05 6.4219e-06 cor 0.436 cov 1.0600e-05 Eig 9.338e-05 5.130e-06 1stv 0.993 0.121 tan 0.122
+85 23-12 mu 1.8515e-01 1.4452e-01 Var 5.4682e-05 1.4710e-05 cor 0.352 cov 9.9793e-06 Eig 5.704e-05 1.236e-05 1stv 0.973 0.229 tan 0.236
+90 23-12 mu 2.0275e-01 2.2829e-01 Var 6.3261e-05 7.3017e-05 cor 0.249 cov 1.6893e-05 Eig 8.572e-05 5.056e-05 1stv 0.601 0.799 tan 1.330
+95 23-12 mu 2.2194e-01 3.5900e-01 Var 1.4978e-04 3.6276e-04 cor 0.226 cov 5.2587e-05 Eig 3.750e-04 1.375e-04 1stv 0.227 0.974 tan 4.283
+70 21-13 mu 2.0469e-01 2.1400e-02 Var 6.2086e-04 3.1732e-06 cor 0.029 cov 1.2856e-06 Eig 6.209e-04 3.170e-06 1stv 1.000 0.002 tan 0.002
+75 21-13 mu 1.7304e-01 2.8084e-02 Var 2.2257e-04 2.5817e-06 cor 0.021 cov 5.0493e-07 Eig 2.226e-04 2.581e-06 1stv 1.000 0.002 tan 0.002
+80 21-13 mu 1.4621e-01 3.6812e-02 Var 7.1363e-05 2.6231e-06 cor 0.028 cov 3.7993e-07 Eig 7.137e-05 2.621e-06 1stv 1.000 0.006 tan 0.006
+85 21-13 mu 1.2348e-01 4.8162e-02 Var 3.9495e-05 6.7559e-06 cor 0.070 cov 1.1505e-06 Eig 3.954e-05 6.715e-06 1stv 0.999 0.035 tan 0.035
+90 21-13 mu 1.0424e-01 6.2827e-02 Var 5.6829e-05 2.4481e-05 cor 0.082 cov 3.0607e-06 Eig 5.712e-05 2.419e-05 1stv 0.996 0.093 tan 0.094
+95 21-13 mu 8.7961e-02 8.1592e-02 Var 8.7335e-05 7.8620e-05 cor 0.076 cov 6.2933e-06 Eig 9.063e-05 7.532e-05 1stv 0.886 0.464 tan 0.524
+70 23-13 mu 1.4061e-01 2.1400e-02 Var 2.2238e-04 3.1732e-06 cor -0.461 cov -1.2236e-05 Eig 2.231e-04 2.492e-06 1stv 0.998 -0.056 tan -0.056
+75 23-13 mu 1.5420e-01 2.8084e-02 Var 1.5306e-04 2.5817e-06 cor -0.531 cov -1.0558e-05 Eig 1.538e-04 1.845e-06 1stv 0.998 -0.070 tan -0.070
+80 23-13 mu 1.6900e-01 3.6812e-02 Var 9.2084e-05 2.6231e-06 cor -0.579 cov -9.0014e-06 Eig 9.298e-05 1.726e-06 1stv 0.995 -0.099 tan -0.100
+85 23-13 mu 1.8515e-01 4.8162e-02 Var 5.4682e-05 6.7559e-06 cor -0.506 cov -9.7162e-06 Eig 5.658e-05 4.861e-06 1stv 0.982 -0.191 tan -0.195
+90 23-13 mu 2.0275e-01 6.2827e-02 Var 6.3261e-05 2.4481e-05 cor -0.442 cov -1.7380e-05 Eig 6.991e-05 1.783e-05 1stv 0.934 -0.357 tan -0.383
+95 23-13 mu 2.2194e-01 8.1592e-02 Var 1.4978e-04 7.8620e-05 cor -0.379 cov -4.1150e-05 Eig 1.686e-04 5.980e-05 1stv 0.909 -0.416 tan -0.457
+70 23-21 mu 1.4061e-01 2.0469e-01 Var 2.2238e-04 6.2086e-04 cor -0.033 cov -1.2366e-05 Eig 6.212e-04 2.220e-04 1stv 0.031 -1.000 tan -32.254
+75 23-21 mu 1.5420e-01 1.7304e-01 Var 1.5306e-04 2.2257e-04 cor -0.032 cov -5.9497e-06 Eig 2.231e-04 1.526e-04 1stv 0.085 -0.996 tan -11.767
+80 23-21 mu 1.6900e-01 1.4621e-01 Var 9.2084e-05 7.1363e-05 cor -0.038 cov -3.0714e-06 Eig 9.253e-05 7.092e-05 1stv 0.990 -0.144 tan -0.145
+85 23-21 mu 1.8515e-01 1.2348e-01 Var 5.4682e-05 3.9495e-05 cor -0.055 cov -2.5387e-06 Eig 5.510e-05 3.908e-05 1stv 0.987 -0.161 tan -0.163
+90 23-21 mu 2.0275e-01 1.0424e-01 Var 6.3261e-05 5.6829e-05 cor -0.060 cov -3.5686e-06 Eig 6.485e-05 5.524e-05 1stv 0.914 -0.407 tan -0.445
+95 23-21 mu 2.2194e-01 8.7961e-02 Var 1.4978e-04 8.7335e-05 cor -0.050 cov -5.6628e-06 Eig 1.503e-04 8.683e-05 1stv 0.996 -0.090 tan -0.090
+Prevalence combination of varying and fixed dummies 1
+
+Computing forecasting: result on file 'F_biaspar.txt', please wait...
+End of Computing forecasting
+Computing Health Expectancies: result on file 'E_biaspar.txt' ...70|71|72|73|74|75|76|77|78|79|80|81|82|83|84|85|86|87|88|89|90|91|92|93|94|95|
+done evsij
+Computing Total Life expectancies with their standard errors: file 'T_biaspar.txt' ...
+ Computing State-specific Expectancies and standard errors: result on file 'STDE_biaspar.txt'
+ Computing Covar. of State-specific Expectancies: result on file 'CVE_biaspar.txt'
+ Computing Variance-covariance of State-specific Expectancies: file 'V_biaspar.txt' ...
+#****** Result for:******
+ cvevsij 70|71|72|73|74|75|76|77|78|79|80|81|82|83|84|85|86|87|88|89|90|91|92|93|94|95|
+ end cvevsij
+ varevsij vpopbased=0
+Computing total mortality p.j=w1*p1j+w2*p2j+..: result on file 'PRMORPREV-1-STABLBASED_biaspar.txt'
+Computing age specific forward period (stable) prevalences in each health state
+varevsij vpopbased=1
+Computing total mortality p.j=w1*p1j+w2*p2j+..: result on file 'PRMORPREV-1-POPULBASED-NOMOBIL_biaspar.txt'
+Computing age specific forward period (stable) prevalences in each health state
+done selection
+done State-specific expectancies
+Computing Variance-covariance of forward period (stable) prevalence: file 'VPL_biaspar.txt' ...
+#****** ******
+done variance-covariance of forward period prevalence
+End of Imach with 14 errors and/or warnings 31. Please look at the log file for details.
+Local time at start Wed May 22 22:32:36 2019
-Local time at end Tue Aug 18 18:26:33 2015
+Local time at end Wed May 22 22:34:09 2019
-Total time used 0 day(s) 0 hour(s) 5 minute(s) 49 second(s)
-Total time was 349 Sec.
+Total time used 0 day(s) 0 hour(s) 1 minute(s) 33 second(s)
+Total time was 93 Sec.
git://henry.ined.fr / .git / commitdiff