--- imach/html/doc/biaspar.gp	2004/06/16 21:44:30	1.1
+++ imach/html/doc/biaspar.gp	2005/10/25 11:53:26	1.2
@@ -1,245 +1,245 @@
-
-# Imach version 0.07a, May 2004, INED-EUROREVES 
-# biaspar.gp
-set missing 'NaNq'
-cd "D:/imachcvs/imach/html/doc/" 
-
-set out "biaspar/vbiaspar11.png" 
-
-#set out "vbiaspar11.png" 
-set xlabel "Age" 
-set ylabel "Probability" 
-set ter png small
-set size 0.65,0.65
-plot [70:95] "biaspar/vplrbiaspar.txt" every :::0::0 u 1:2 "%lf %lf (%lf) %*lf (%*lf)" t"Stable prevalence" w l 0,"biaspar/vplrbiaspar.txt" every :::0::0 u 1:($2+1.96*$3) "%lf %lf (%lf) %*lf (%*lf)" t"95% CI" w l 1,"biaspar/vplrbiaspar.txt" every :::0::0 u 1:($2-1.96*$3) "%lf %lf (%lf) %*lf (%*lf)" t"" w l 1,"biaspar/prbiaspar.txt" every :::0::0 u 1:($2) t"Observed prevalence " w l 2
-set out "biaspar/vbiaspar21.png" 
-
-#set out "vbiaspar21.png" 
-set xlabel "Age" 
-set ylabel "Probability" 
-set ter png small
-set size 0.65,0.65
-plot [70:95] "biaspar/vplrbiaspar.txt" every :::0::0 u 1:2 "%lf %*lf (%*lf) %lf (%lf)" t"Stable prevalence" w l 0,"biaspar/vplrbiaspar.txt" every :::0::0 u 1:($2+1.96*$3) "%lf %*lf (%*lf) %lf (%lf)" t"95% CI" w l 1,"biaspar/vplrbiaspar.txt" every :::0::0 u 1:($2-1.96*$3) "%lf %*lf (%*lf) %lf (%lf)" t"" w l 1,"biaspar/prbiaspar.txt" every :::0::0 u 1:($6) t"Observed prevalence " w l 2
-set out "biaspar/ebiaspar1.png" 
-set ylabel "Years" 
-set ter png small
-set size 0.65,0.65
-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 0,"biaspar/trbiaspar.txt" every :::0::0 u 1:($2+$3*2) "%lf %lf (%lf) %*lf (%*lf) %*lf (%*lf)" t"" w l 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 0,"biaspar/trbiaspar.txt" every :::0::0 u 1:($2+$3*2) "%lf %*lf (%*lf) %lf (%lf) %*lf (%*lf)" t"" w l 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 0,"biaspar/trbiaspar.txt" every :::0::0 u 1:($2+$3*2) "%lf %*lf (%*lf) %*lf (%*lf) %lf (%lf)" t"" w l 0
-set out "biaspar/expbiaspar11.png" 
-set ter png small
-set size 0.65,0.65
-plot [70:95] "biaspar/erbiaspar.txt" every :::0::0 u 1:2 t "e11" w l ,"biaspar/erbiaspar.txt" every :::0::0 u 1:4 t "e12" w l
-set out "biaspar/expbiaspar21.png" 
-set ter png small
-set size 0.65,0.65
-plot [70:95] "biaspar/erbiaspar.txt" every :::0::0 u 1:6 t "e21" w l ,"biaspar/erbiaspar.txt" every :::0::0 u 1:8 t "e22" w l
-set out "biaspar/pbiaspar11.png" 
-set xlabel "Age" 
-set ylabel "Probability" 
-set ter png small
-set size 0.65,0.65
-unset log y
-plot [70:95] "biaspar/pijrbiaspar.txt" u ($1==1 ? ($3):1/0):($5/($4+$5)) t"prev(1,2)" w l,"biaspar/pijrbiaspar.txt" u ($1==1 ? ($3):1/0):($8/($7+$8)) t"prev(2,2)" w l
-
-set out "biaspar/pbiaspar21.png" 
-set xlabel "Age" 
-set ylabel "Probability" 
-set ter png small
-set size 0.65,0.65
-unset log y
-plot [70:95] "biaspar/pijrbiaspar.txt" u ($1==1 ? ($3):1/0):($6/($4+$5)) t"prev(2,3)" w l,"biaspar/pijrbiaspar.txt" u ($1==1 ? ($3):1/0):($12/($10+$11)) t"prev(3,3)" w l
-p1=-12.245160 
-p2=0.092357 
-p3=-10.672078 
-p4=0.060971 
-p5=-2.645815 
-p6=-0.022320 
-p7=-4.773208 
-p8=0.007857 
-
-set out "biaspar/pebiaspar11.png" 
-
-set title "Probability"
-
-set ter png small
-set size 0.65,0.65
-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" 
-set out "biaspar/pebiaspar12.png" 
-
-set ylabel "Quasi-incidence per year"
-
-set ter png small
-set size 0.65,0.65
-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" 
-# Routine varprob
-set parametric;unset label
-set log y;set log x; set xlabel "p13 (year-1)";set ylabel "p12 (year-1)"
-set ter png small
-set size 0.65,0.65
-set out "biaspar/varpijgrbiaspar113-12.png"
-set label "65" at  1.459e-002, 2.328e-002 center
-# Age 65, p13 - p12
-plot [-pi:pi]  1.459e-002+ 2.000*( 6.866e-001* 2.057e-003*cos(t)+ 7.270e-001* 1.402e-003*sin(t)),  2.328e-002 +2.000*(-7.270e-001* 2.057e-003*cos(t)+ 6.866e-001* 1.402e-003*sin(t)) not
-# Age 70, p13 - p12
-set label "70" at  1.976e-002, 3.688e-002 center
-replot  1.976e-002+ 2.000*( 5.191e-001* 2.326e-003*cos(t)+ 8.547e-001* 1.471e-003*sin(t)),  3.688e-002 +2.000*(-8.547e-001* 2.326e-003*cos(t)+ 5.191e-001* 1.471e-003*sin(t)) not
-# Age 75, p13 - p12
-set label "75" at  2.674e-002, 5.839e-002 center
-replot  2.674e-002+ 2.000*( 4.119e-001* 2.525e-003*cos(t)+ 9.112e-001* 1.387e-003*sin(t)),  5.839e-002 +2.000*(-9.112e-001* 2.525e-003*cos(t)+ 4.119e-001* 1.387e-003*sin(t)) not
-# Age 80, p13 - p12
-set label "80" at  3.614e-002, 9.233e-002 center
-replot  3.614e-002+ 2.000*( 3.972e-001* 2.741e-003*cos(t)+ 9.177e-001* 1.382e-003*sin(t)),  9.233e-002 +2.000*(-9.177e-001* 2.741e-003*cos(t)+ 3.972e-001* 1.382e-003*sin(t)) not
-# Age 85, p13 - p12
-set label "85" at  4.875e-002, 1.457e-001 center
-replot  4.875e-002+ 2.000*( 3.964e-001* 4.133e-003*cos(t)+ 9.181e-001* 2.272e-003*sin(t)),  1.457e-001 +2.000*(-9.181e-001* 4.133e-003*cos(t)+ 3.964e-001* 2.272e-003*sin(t)) not
-# Age 90, p13 - p12
-set label "90" at  6.557e-002, 2.292e-001 center
-replot  6.557e-002+ 2.000*( 2.870e-001* 8.981e-003*cos(t)+ 9.579e-001* 4.671e-003*sin(t)),  2.292e-001 +2.000*(-9.579e-001* 8.981e-003*cos(t)+ 2.870e-001* 4.671e-003*sin(t)) not
-# Age 95, p13 - p12
-set label "95" at  8.778e-002, 3.591e-001 center
-replot  8.778e-002+ 2.000*( 2.106e-001* 1.972e-002*cos(t)+ 9.776e-001* 8.813e-003*sin(t)),  3.591e-001 +2.000*(-9.776e-001* 1.972e-002*cos(t)+ 2.106e-001* 8.813e-003*sin(t)) not
-set out "biaspar/varpijgrbiaspar113-12.png";replot;
-set parametric;unset label
-set log y;set log x; set xlabel "p21 (year-1)";set ylabel "p12 (year-1)"
-set ter png small
-set size 0.65,0.65
-set out "biaspar/varpijgrbiaspar121-12.png"
-set label "65" at  1.936e-001, 2.328e-002 center
-# Age 65, p21 - p12
-plot [-pi:pi]  1.936e-001+ 2.000*( 9.998e-001* 3.151e-002*cos(t)+-1.959e-002* 1.668e-003*sin(t)),  2.328e-002 +2.000*( 1.959e-002* 3.151e-002*cos(t)+ 9.998e-001* 1.668e-003*sin(t)) not
-# Age 70, p21 - p12
-set label "70" at  1.734e-001, 3.688e-002 center
-replot  1.734e-001+ 2.000*( 9.994e-001* 2.161e-002*cos(t)+-3.389e-002* 2.001e-003*sin(t)),  3.688e-002 +2.000*( 3.389e-002* 2.161e-002*cos(t)+ 9.994e-001* 2.001e-003*sin(t)) not
-# Age 75, p21 - p12
-set label "75" at  1.552e-001, 5.839e-002 center
-replot  1.552e-001+ 2.000*( 9.983e-001* 1.378e-002*cos(t)+-5.800e-002* 2.235e-003*sin(t)),  5.839e-002 +2.000*( 5.800e-002* 1.378e-002*cos(t)+ 9.983e-001* 2.235e-003*sin(t)) not
-# Age 80, p21 - p12
-set label "80" at  1.389e-001, 9.233e-002 center
-replot  1.389e-001+ 2.000*( 9.949e-001* 8.292e-003*cos(t)+-1.011e-001* 2.447e-003*sin(t)),  9.233e-002 +2.000*( 1.011e-001* 8.292e-003*cos(t)+ 9.949e-001* 2.447e-003*sin(t)) not
-# Age 85, p21 - p12
-set label "85" at  1.243e-001, 1.457e-001 center
-replot  1.243e-001+ 2.000*( 9.669e-001* 6.441e-003*cos(t)+-2.553e-001* 3.657e-003*sin(t)),  1.457e-001 +2.000*( 2.553e-001* 6.441e-003*cos(t)+ 9.669e-001* 3.657e-003*sin(t)) not
-# Age 90, p21 - p12
-set label "90" at  1.112e-001, 2.292e-001 center
-replot  1.112e-001+ 2.000*( 5.879e-001* 9.566e-003*cos(t)+-8.090e-001* 6.788e-003*sin(t)),  2.292e-001 +2.000*( 8.090e-001* 9.566e-003*cos(t)+ 5.879e-001* 6.788e-003*sin(t)) not
-# Age 95, p21 - p12
-set label "95" at  9.953e-002, 3.591e-001 center
-replot  9.953e-002+ 2.000*( 2.290e-001* 1.977e-002*cos(t)+-9.734e-001* 9.572e-003*sin(t)),  3.591e-001 +2.000*( 9.734e-001* 1.977e-002*cos(t)+ 2.290e-001* 9.572e-003*sin(t)) not
-set out "biaspar/varpijgrbiaspar121-12.png";replot;
-set parametric;unset label
-set log y;set log x; set xlabel "p23 (year-1)";set ylabel "p12 (year-1)"
-set ter png small
-set size 0.65,0.65
-set out "biaspar/varpijgrbiaspar123-12.png"
-set label "65" at  1.640e-001, 2.328e-002 center
-# Age 65, p23 - p12
-plot [-pi:pi]  1.640e-001+ 2.000*( 9.995e-001* 2.107e-002*cos(t)+-3.113e-002* 1.654e-003*sin(t)),  2.328e-002 +2.000*( 3.113e-002* 2.107e-002*cos(t)+ 9.995e-001* 1.654e-003*sin(t)) not
-# Age 70, p23 - p12
-set label "70" at  1.708e-001, 3.688e-002 center
-replot  1.708e-001+ 2.000*( 9.988e-001* 1.745e-002*cos(t)+-4.823e-002* 1.959e-003*sin(t)),  3.688e-002 +2.000*( 4.823e-002* 1.745e-002*cos(t)+ 9.988e-001* 1.959e-003*sin(t)) not
-# Age 75, p23 - p12
-set label "75" at  1.778e-001, 5.839e-002 center
-replot  1.778e-001+ 2.000*( 9.971e-001* 1.368e-002*cos(t)+-7.559e-002* 2.139e-003*sin(t)),  5.839e-002 +2.000*( 7.559e-002* 1.368e-002*cos(t)+ 9.971e-001* 2.139e-003*sin(t)) not
-# Age 80, p23 - p12
-set label "80" at  1.851e-001, 9.233e-002 center
-replot  1.851e-001+ 2.000*( 9.925e-001* 1.009e-002*cos(t)+-1.220e-001* 2.279e-003*sin(t)),  9.233e-002 +2.000*( 1.220e-001* 1.009e-002*cos(t)+ 9.925e-001* 2.279e-003*sin(t)) not
-# Age 85, p23 - p12
-set label "85" at  1.926e-001, 1.457e-001 center
-replot  1.926e-001+ 2.000*( 9.722e-001* 7.671e-003*cos(t)+-2.342e-001* 3.560e-003*sin(t)),  1.457e-001 +2.000*( 2.342e-001* 7.671e-003*cos(t)+ 9.722e-001* 3.560e-003*sin(t)) not
-# Age 90, p23 - p12
-set label "90" at  2.004e-001, 2.292e-001 center
-replot  2.004e-001+ 2.000*( 6.128e-001* 9.486e-003*cos(t)+-7.902e-001* 7.229e-003*sin(t)),  2.292e-001 +2.000*( 7.902e-001* 9.486e-003*cos(t)+ 6.128e-001* 7.229e-003*sin(t)) not
-# Age 95, p23 - p12
-set label "95" at  2.085e-001, 3.591e-001 center
-replot  2.085e-001+ 2.000*( 2.248e-001* 1.970e-002*cos(t)+-9.744e-001* 1.148e-002*sin(t)),  3.591e-001 +2.000*( 9.744e-001* 1.970e-002*cos(t)+ 2.248e-001* 1.148e-002*sin(t)) not
-set out "biaspar/varpijgrbiaspar123-12.png";replot;
-set parametric;unset label
-set log y;set log x; set xlabel "p21 (year-1)";set ylabel "p13 (year-1)"
-set ter png small
-set size 0.65,0.65
-set out "biaspar/varpijgrbiaspar121-13.png"
-set label "65" at  1.936e-001, 1.459e-002 center
-# Age 65, p21 - p13
-plot [-pi:pi]  1.936e-001+ 2.000*( 1.000e+000* 3.150e-002*cos(t)+-1.880e-003* 1.741e-003*sin(t)),  1.459e-002 +2.000*( 1.880e-003* 3.150e-002*cos(t)+ 1.000e+000* 1.741e-003*sin(t)) not
-# Age 70, p21 - p13
-set label "70" at  1.734e-001, 1.976e-002 center
-replot  1.734e-001+ 2.000*( 1.000e+000* 2.160e-002*cos(t)+-1.974e-003* 1.743e-003*sin(t)),  1.976e-002 +2.000*( 1.974e-003* 2.160e-002*cos(t)+ 1.000e+000* 1.743e-003*sin(t)) not
-# Age 75, p21 - p13
-set label "75" at  1.552e-001, 2.674e-002 center
-replot  1.552e-001+ 2.000*( 1.000e+000* 1.376e-002*cos(t)+-1.404e-003* 1.636e-003*sin(t)),  2.674e-002 +2.000*( 1.404e-003* 1.376e-002*cos(t)+ 1.000e+000* 1.636e-003*sin(t)) not
-# Age 80, p21 - p13
-set label "80" at  1.389e-001, 3.614e-002 center
-replot  1.389e-001+ 2.000*( 1.000e+000* 8.253e-003*cos(t)+-3.684e-003* 1.671e-003*sin(t)),  3.614e-002 +2.000*( 3.684e-003* 8.253e-003*cos(t)+ 1.000e+000* 1.671e-003*sin(t)) not
-# Age 85, p21 - p13
-set label "85" at  1.243e-001, 4.875e-002 center
-replot  1.243e-001+ 2.000*( 9.992e-001* 6.301e-003*cos(t)+-4.077e-002* 2.642e-003*sin(t)),  4.875e-002 +2.000*( 4.077e-002* 6.301e-003*cos(t)+ 9.992e-001* 2.642e-003*sin(t)) not
-# Age 90, p21 - p13
-set label "90" at  1.112e-001, 6.557e-002 center
-replot  1.112e-001+ 2.000*( 9.932e-001* 7.891e-003*cos(t)+-1.164e-001* 5.116e-003*sin(t)),  6.557e-002 +2.000*( 1.164e-001* 7.891e-003*cos(t)+ 9.932e-001* 5.116e-003*sin(t)) not
-# Age 95, p21 - p13
-set label "95" at  9.953e-002, 8.778e-002 center
-replot  9.953e-002+ 2.000*( 9.056e-001* 1.057e-002*cos(t)+-4.241e-001* 9.328e-003*sin(t)),  8.778e-002 +2.000*( 4.241e-001* 1.057e-002*cos(t)+ 9.056e-001* 9.328e-003*sin(t)) not
-set out "biaspar/varpijgrbiaspar121-13.png";replot;
-set parametric;unset label
-set log y;set log x; set xlabel "p23 (year-1)";set ylabel "p13 (year-1)"
-set ter png small
-set size 0.65,0.65
-set out "biaspar/varpijgrbiaspar123-13.png"
-set label "65" at  1.640e-001, 1.459e-002 center
-# Age 65, p23 - p13
-plot [-pi:pi]  1.640e-001+ 2.000*( 9.993e-001* 2.108e-002*cos(t)+ 3.790e-002* 1.549e-003*sin(t)),  1.459e-002 +2.000*(-3.790e-002* 2.108e-002*cos(t)+ 9.993e-001* 1.549e-003*sin(t)) not
-# Age 70, p23 - p13
-set label "70" at  1.708e-001, 1.976e-002 center
-replot  1.708e-001+ 2.000*( 9.988e-001* 1.745e-002*cos(t)+ 4.951e-002* 1.516e-003*sin(t)),  1.976e-002 +2.000*(-4.951e-002* 1.745e-002*cos(t)+ 9.988e-001* 1.516e-003*sin(t)) not
-# Age 75, p23 - p13
-set label "75" at  1.778e-001, 2.674e-002 center
-replot  1.778e-001+ 2.000*( 9.978e-001* 1.368e-002*cos(t)+ 6.667e-002* 1.362e-003*sin(t)),  2.674e-002 +2.000*(-6.667e-002* 1.368e-002*cos(t)+ 9.978e-001* 1.362e-003*sin(t)) not
-# Age 80, p23 - p13
-set label "80" at  1.851e-001, 3.614e-002 center
-replot  1.851e-001+ 2.000*( 9.950e-001* 1.006e-002*cos(t)+ 1.004e-001* 1.338e-003*sin(t)),  3.614e-002 +2.000*(-1.004e-001* 1.006e-002*cos(t)+ 9.950e-001* 1.338e-003*sin(t)) not
-# Age 85, p23 - p13
-set label "85" at  1.926e-001, 4.875e-002 center
-replot  1.926e-001+ 2.000*( 9.801e-001* 7.644e-003*cos(t)+ 1.985e-001* 2.219e-003*sin(t)),  4.875e-002 +2.000*(-1.985e-001* 7.644e-003*cos(t)+ 9.801e-001* 2.219e-003*sin(t)) not
-# Age 90, p23 - p13
-set label "90" at  2.004e-001, 6.557e-002 center
-replot  2.004e-001+ 2.000*( 9.273e-001* 8.612e-003*cos(t)+ 3.743e-001* 4.351e-003*sin(t)),  6.557e-002 +2.000*(-3.743e-001* 8.612e-003*cos(t)+ 9.273e-001* 4.351e-003*sin(t)) not
-# Age 95, p23 - p13
-set label "95" at  2.085e-001, 8.778e-002 center
-replot  2.085e-001+ 2.000*( 8.652e-001* 1.310e-002*cos(t)+ 5.014e-001* 8.033e-003*sin(t)),  8.778e-002 +2.000*(-5.014e-001* 1.310e-002*cos(t)+ 8.652e-001* 8.033e-003*sin(t)) not
-set out "biaspar/varpijgrbiaspar123-13.png";replot;
-set parametric;unset label
-set log y;set log x; set xlabel "p23 (year-1)";set ylabel "p21 (year-1)"
-set ter png small
-set size 0.65,0.65
-set out "biaspar/varpijgrbiaspar123-21.png"
-set label "65" at  1.640e-001, 1.936e-001 center
-# Age 65, p23 - p21
-plot [-pi:pi]  1.640e-001+ 2.000*( 3.032e-002* 3.151e-002*cos(t)+ 9.995e-001* 2.105e-002*sin(t)),  1.936e-001 +2.000*(-9.995e-001* 3.151e-002*cos(t)+ 3.032e-002* 2.105e-002*sin(t)) not
-# Age 70, p23 - p21
-set label "70" at  1.708e-001, 1.734e-001 center
-replot  1.708e-001+ 2.000*( 4.964e-002* 2.161e-002*cos(t)+ 9.988e-001* 1.741e-002*sin(t)),  1.734e-001 +2.000*(-9.988e-001* 2.161e-002*cos(t)+ 4.964e-002* 1.741e-002*sin(t)) not
-# Age 75, p23 - p21
-set label "75" at  1.778e-001, 1.552e-001 center
-replot  1.778e-001+ 2.000*( 5.539e-001* 1.385e-002*cos(t)+ 8.326e-001* 1.356e-002*sin(t)),  1.552e-001 +2.000*(-8.326e-001* 1.385e-002*cos(t)+ 5.539e-001* 1.356e-002*sin(t)) not
-# Age 80, p23 - p21
-set label "80" at  1.851e-001, 1.389e-001 center
-replot  1.851e-001+ 2.000*( 9.977e-001* 1.002e-002*cos(t)+ 6.827e-002* 8.244e-003*sin(t)),  1.389e-001 +2.000*(-6.827e-002* 1.002e-002*cos(t)+ 9.977e-001* 8.244e-003*sin(t)) not
-# Age 85, p23 - p21
-set label "85" at  1.926e-001, 1.243e-001 center
-replot  1.926e-001+ 2.000*( 9.865e-001* 7.536e-003*cos(t)+ 1.638e-001* 6.259e-003*sin(t)),  1.243e-001 +2.000*(-1.638e-001* 7.536e-003*cos(t)+ 9.865e-001* 6.259e-003*sin(t)) not
-# Age 90, p23 - p21
-set label "90" at  2.004e-001, 1.112e-001 center
-replot  2.004e-001+ 2.000*( 8.440e-001* 8.341e-003*cos(t)+ 5.363e-001* 7.658e-003*sin(t)),  1.112e-001 +2.000*(-5.363e-001* 8.341e-003*cos(t)+ 8.440e-001* 7.658e-003*sin(t)) not
-# Age 95, p23 - p21
-set label "95" at  2.085e-001, 9.953e-002 center
-replot  2.085e-001+ 2.000*( 9.795e-001* 1.210e-002*cos(t)+ 2.015e-001* 1.028e-002*sin(t)),  9.953e-002 +2.000*(-2.015e-001* 1.210e-002*cos(t)+ 9.795e-001* 1.028e-002*sin(t)) not
-set out "biaspar/varpijgrbiaspar123-21.png";replot;
-# Routine varevsij
-set noparametric;set nolabel; set ter png small;set size 0.65, 0.65
- set log y; set nolog x;set xlabel "Age"; set ylabel "Force of mortality (year-1)";
- plot "biaspar/prmorprev1-stablbased-rbiaspar.txt"  u 1:($3) not w l 1 
- replot "biaspar/prmorprev1-stablbased-rbiaspar.txt"  u 1:(($3+1.96*$4)) t "95% interval" w l 2 
- replot "biaspar/prmorprev1-stablbased-rbiaspar.txt"  u 1:(($3-1.96*$4)) not w l 2 
-set out "biaspar/varmuptjgr-stablbased-biaspar1.png";replot;
+
+# Imach version 0.98, September 2005, INED-EUROREVES 
+# biaspar.gp
+set missing 'NaNq'
+cd "/Users/brouard/bin/html/doc" 
+
+set out "biaspar/vbiaspar11.png" 
+
+#set out "vbiaspar11.png" 
+set xlabel "Age" 
+set ylabel "Probability" 
+set ter png small
+set size 0.65,0.65
+plot [70:95] "biaspar/vplrbiaspar.txt" every :::0::0 u 1:2 "%lf %lf (%lf) %*lf (%*lf)" t"Stable prevalence" w l 0,"biaspar/vplrbiaspar.txt" every :::0::0 u 1:($2+1.96*$3) "%lf %lf (%lf) %*lf (%*lf)" t"95% CI" w l 1,"biaspar/vplrbiaspar.txt" every :::0::0 u 1:($2-1.96*$3) "%lf %lf (%lf) %*lf (%*lf)" t"" w l 1,"biaspar/prbiaspar.txt" every :::0::0 u 1:($2) t"Observed prevalence " w l 2
+set out "biaspar/vbiaspar21.png" 
+
+#set out "vbiaspar21.png" 
+set xlabel "Age" 
+set ylabel "Probability" 
+set ter png small
+set size 0.65,0.65
+plot [70:95] "biaspar/vplrbiaspar.txt" every :::0::0 u 1:2 "%lf %*lf (%*lf) %lf (%lf)" t"Stable prevalence" w l 0,"biaspar/vplrbiaspar.txt" every :::0::0 u 1:($2+1.96*$3) "%lf %*lf (%*lf) %lf (%lf)" t"95% CI" w l 1,"biaspar/vplrbiaspar.txt" every :::0::0 u 1:($2-1.96*$3) "%lf %*lf (%*lf) %lf (%lf)" t"" w l 1,"biaspar/prbiaspar.txt" every :::0::0 u 1:($6) t"Observed prevalence " w l 2
+set out "biaspar/ebiaspar1.png" 
+set ylabel "Years" 
+set ter png small
+set size 0.65,0.65
+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 0,"biaspar/trbiaspar.txt" every :::0::0 u 1:($2+$3*2) "%lf %lf (%lf) %*lf (%*lf) %*lf (%*lf)" t"" w l 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 0,"biaspar/trbiaspar.txt" every :::0::0 u 1:($2+$3*2) "%lf %*lf (%*lf) %lf (%lf) %*lf (%*lf)" t"" w l 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 0,"biaspar/trbiaspar.txt" every :::0::0 u 1:($2+$3*2) "%lf %*lf (%*lf) %*lf (%*lf) %lf (%lf)" t"" w l 0
+set out "biaspar/expbiaspar11.png" 
+set ter png small
+set size 0.65,0.65
+plot [70:95] "biaspar/erbiaspar.txt" every :::0::0 u 1:2 t "e11" w l ,"biaspar/erbiaspar.txt" every :::0::0 u 1:4 t "e12" w l
+set out "biaspar/expbiaspar21.png" 
+set ter png small
+set size 0.65,0.65
+plot [70:95] "biaspar/erbiaspar.txt" every :::0::0 u 1:6 t "e21" w l ,"biaspar/erbiaspar.txt" every :::0::0 u 1:8 t "e22" w l
+set out "biaspar/pbiaspar11.png" 
+set xlabel "Age" 
+set ylabel "Probability" 
+set ter png small
+set size 0.65,0.65
+unset log y
+plot [70:95] "biaspar/pijrbiaspar.txt" u ($1==1 ? ($3):1/0):($5/($4+$5)) t"prev(1,2)" w l,"biaspar/pijrbiaspar.txt" u ($1==1 ? ($3):1/0):($8/($7+$8)) t"prev(2,2)" w l
+
+set out "biaspar/pbiaspar21.png" 
+set xlabel "Age" 
+set ylabel "Probability" 
+set ter png small
+set size 0.65,0.65
+unset log y
+plot [70:95] "biaspar/pijrbiaspar.txt" u ($1==1 ? ($3):1/0):($6/($4+$5)) t"prev(2,3)" w l,"biaspar/pijrbiaspar.txt" u ($1==1 ? ($3):1/0):($12/($10+$11)) t"prev(3,3)" w l
+p1=-12.245504 
+p2=0.092361 
+p3=-10.670953 
+p4=0.060958 
+p5=-2.645034 
+p6=-0.022329 
+p7=-4.775501 
+p8=0.007884 
+
+set out "biaspar/pebiaspar11.png" 
+
+set title "Probability"
+
+set ter png small
+set size 0.65,0.65
+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" 
+set out "biaspar/pebiaspar12.png" 
+
+set ylabel "Quasi-incidence per year"
+
+set ter png small
+set size 0.65,0.65
+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" 
+# Routine varprob
+set parametric;unset label
+set log y;set log x; set xlabel "p13 (year-1)";set ylabel "p12 (year-1)"
+set ter png small
+set size 0.65,0.65
+set out "biaspar/varpijgrbiaspar113-12.png"
+set label "65" at   1.460e-02,  2.328e-02 center
+# Age 65, p13 - p12
+plot [-pi:pi]   1.460e-02+ 2.000*(  7.019e-01*  2.099e-03*cos(t)+  7.123e-01*  1.386e-03*sin(t)),   2.328e-02 +2.000*( -7.123e-01*  2.099e-03*cos(t)+  7.019e-01*  1.386e-03*sin(t)) not
+# Age 70, p13 - p12
+set label "70" at   1.977e-02,  3.688e-02 center
+replot   1.977e-02+ 2.000*(  5.418e-01*  2.359e-03*cos(t)+  8.405e-01*  1.461e-03*sin(t)),   3.688e-02 +2.000*( -8.405e-01*  2.359e-03*cos(t)+  5.418e-01*  1.461e-03*sin(t)) not
+# Age 75, p13 - p12
+set label "75" at   2.675e-02,  5.839e-02 center
+replot   2.675e-02+ 2.000*(  4.274e-01*  2.544e-03*cos(t)+  9.041e-01*  1.383e-03*sin(t)),   5.839e-02 +2.000*( -9.041e-01*  2.544e-03*cos(t)+  4.274e-01*  1.383e-03*sin(t)) not
+# Age 80, p13 - p12
+set label "80" at   3.615e-02,  9.232e-02 center
+replot   3.615e-02+ 2.000*(  3.992e-01*  2.744e-03*cos(t)+  9.169e-01*  1.381e-03*sin(t)),   9.232e-02 +2.000*( -9.169e-01*  2.744e-03*cos(t)+  3.992e-01*  1.381e-03*sin(t)) not
+# Age 85, p13 - p12
+set label "85" at   4.875e-02,  1.457e-01 center
+replot   4.875e-02+ 2.000*(  3.981e-01*  4.136e-03*cos(t)+  9.174e-01*  2.268e-03*sin(t)),   1.457e-01 +2.000*( -9.174e-01*  4.136e-03*cos(t)+  3.981e-01*  2.268e-03*sin(t)) not
+# Age 90, p13 - p12
+set label "90" at   6.556e-02,  2.292e-01 center
+replot   6.556e-02+ 2.000*(  2.987e-01*  9.027e-03*cos(t)+  9.543e-01*  4.655e-03*sin(t)),   2.292e-01 +2.000*( -9.543e-01*  9.027e-03*cos(t)+  2.987e-01*  4.655e-03*sin(t)) not
+# Age 95, p13 - p12
+set label "95" at   8.777e-02,  3.591e-01 center
+replot   8.777e-02+ 2.000*(  2.247e-01*  1.984e-02*cos(t)+  9.744e-01*  8.795e-03*sin(t)),   3.591e-01 +2.000*( -9.744e-01*  1.984e-02*cos(t)+  2.247e-01*  8.795e-03*sin(t)) not
+set out "biaspar/varpijgrbiaspar113-12.png";replot;
+set parametric;unset label
+set log y;set log x; set xlabel "p21 (year-1)";set ylabel "p12 (year-1)"
+set ter png small
+set size 0.65,0.65
+set out "biaspar/varpijgrbiaspar121-12.png"
+set label "65" at   1.936e-01,  2.328e-02 center
+# Age 65, p21 - p12
+plot [-pi:pi]   1.936e-01+ 2.000*(  9.998e-01*  3.109e-02*cos(t)+ -1.822e-02*  1.692e-03*sin(t)),   2.328e-02 +2.000*(  1.822e-02*  3.109e-02*cos(t)+  9.998e-01*  1.692e-03*sin(t)) not
+# Age 70, p21 - p12
+set label "70" at   1.734e-01,  3.688e-02 center
+replot   1.734e-01+ 2.000*(  9.995e-01*  2.135e-02*cos(t)+ -3.169e-02*  2.026e-03*sin(t)),   3.688e-02 +2.000*(  3.169e-02*  2.135e-02*cos(t)+  9.995e-01*  2.026e-03*sin(t)) not
+# Age 75, p21 - p12
+set label "75" at   1.552e-01,  5.839e-02 center
+replot   1.552e-01+ 2.000*(  9.985e-01*  1.365e-02*cos(t)+ -5.500e-02*  2.256e-03*sin(t)),   5.839e-02 +2.000*(  5.500e-02*  1.365e-02*cos(t)+  9.985e-01*  2.256e-03*sin(t)) not
+# Age 80, p21 - p12
+set label "80" at   1.389e-01,  9.232e-02 center
+replot   1.389e-01+ 2.000*(  9.951e-01*  8.259e-03*cos(t)+ -9.930e-02*  2.453e-03*sin(t)),   9.232e-02 +2.000*(  9.930e-02*  8.259e-03*cos(t)+  9.951e-01*  2.453e-03*sin(t)) not
+# Age 85, p21 - p12
+set label "85" at   1.243e-01,  1.457e-01 center
+replot   1.243e-01+ 2.000*(  9.672e-01*  6.435e-03*cos(t)+ -2.542e-01*  3.661e-03*sin(t)),   1.457e-01 +2.000*(  2.542e-01*  6.435e-03*cos(t)+  9.672e-01*  3.661e-03*sin(t)) not
+# Age 90, p21 - p12
+set label "90" at   1.112e-01,  2.292e-01 center
+replot   1.112e-01+ 2.000*(  5.677e-01*  9.492e-03*cos(t)+ -8.232e-01*  6.841e-03*sin(t)),   2.292e-01 +2.000*(  8.232e-01*  9.492e-03*cos(t)+  5.677e-01*  6.841e-03*sin(t)) not
+# Age 95, p21 - p12
+set label "95" at   9.952e-02,  3.591e-01 center
+replot   9.952e-02+ 2.000*(  2.083e-01*  1.976e-02*cos(t)+ -9.781e-01*  9.575e-03*sin(t)),   3.591e-01 +2.000*(  9.781e-01*  1.976e-02*cos(t)+  2.083e-01*  9.575e-03*sin(t)) not
+set out "biaspar/varpijgrbiaspar121-12.png";replot;
+set parametric;unset label
+set log y;set log x; set xlabel "p23 (year-1)";set ylabel "p12 (year-1)"
+set ter png small
+set size 0.65,0.65
+set out "biaspar/varpijgrbiaspar123-12.png"
+set label "65" at   1.639e-01,  2.328e-02 center
+# Age 65, p23 - p12
+plot [-pi:pi]   1.639e-01+ 2.000*(  9.995e-01*  2.122e-02*cos(t)+ -3.216e-02*  1.649e-03*sin(t)),   2.328e-02 +2.000*(  3.216e-02*  2.122e-02*cos(t)+  9.995e-01*  1.649e-03*sin(t)) not
+# Age 70, p23 - p12
+set label "70" at   1.707e-01,  3.688e-02 center
+replot   1.707e-01+ 2.000*(  9.988e-01*  1.755e-02*cos(t)+ -4.952e-02*  1.952e-03*sin(t)),   3.688e-02 +2.000*(  4.952e-02*  1.755e-02*cos(t)+  9.988e-01*  1.952e-03*sin(t)) not
+# Age 75, p23 - p12
+set label "75" at   1.777e-01,  5.839e-02 center
+replot   1.777e-01+ 2.000*(  9.970e-01*  1.375e-02*cos(t)+ -7.687e-02*  2.132e-03*sin(t)),   5.839e-02 +2.000*(  7.687e-02*  1.375e-02*cos(t)+  9.970e-01*  2.132e-03*sin(t)) not
+# Age 80, p23 - p12
+set label "80" at   1.850e-01,  9.232e-02 center
+replot   1.850e-01+ 2.000*(  9.925e-01*  1.011e-02*cos(t)+ -1.224e-01*  2.275e-03*sin(t)),   9.232e-02 +2.000*(  1.224e-01*  1.011e-02*cos(t)+  9.925e-01*  2.275e-03*sin(t)) not
+# Age 85, p23 - p12
+set label "85" at   1.926e-01,  1.457e-01 center
+replot   1.926e-01+ 2.000*(  9.723e-01*  7.666e-03*cos(t)+ -2.339e-01*  3.562e-03*sin(t)),   1.457e-01 +2.000*(  2.339e-01*  7.666e-03*cos(t)+  9.723e-01*  3.562e-03*sin(t)) not
+# Age 90, p23 - p12
+set label "90" at   2.004e-01,  2.292e-01 center
+replot   2.004e-01+ 2.000*(  6.167e-01*  9.540e-03*cos(t)+ -7.872e-01*  7.206e-03*sin(t)),   2.292e-01 +2.000*(  7.872e-01*  9.540e-03*cos(t)+  6.167e-01*  7.206e-03*sin(t)) not
+# Age 95, p23 - p12
+set label "95" at   2.085e-01,  3.591e-01 center
+replot   2.085e-01+ 2.000*(  2.384e-01*  1.981e-02*cos(t)+ -9.712e-01*  1.149e-02*sin(t)),   3.591e-01 +2.000*(  9.712e-01*  1.981e-02*cos(t)+  2.384e-01*  1.149e-02*sin(t)) not
+set out "biaspar/varpijgrbiaspar123-12.png";replot;
+set parametric;unset label
+set log y;set log x; set xlabel "p21 (year-1)";set ylabel "p13 (year-1)"
+set ter png small
+set size 0.65,0.65
+set out "biaspar/varpijgrbiaspar121-13.png"
+set label "65" at   1.936e-01,  1.460e-02 center
+# Age 65, p21 - p13
+plot [-pi:pi]   1.936e-01+ 2.000*(  1.000e+00*  3.109e-02*cos(t)+ -1.676e-03*  1.773e-03*sin(t)),   1.460e-02 +2.000*(  1.676e-03*  3.109e-02*cos(t)+  1.000e+00*  1.773e-03*sin(t)) not
+# Age 70, p21 - p13
+set label "70" at   1.734e-01,  1.977e-02 center
+replot   1.734e-01+ 2.000*(  1.000e+00*  2.134e-02*cos(t)+ -1.717e-03*  1.772e-03*sin(t)),   1.977e-02 +2.000*(  1.717e-03*  2.134e-02*cos(t)+  1.000e+00*  1.772e-03*sin(t)) not
+# Age 75, p21 - p13
+set label "75" at   1.552e-01,  2.675e-02 center
+replot   1.552e-01+ 2.000*(  1.000e+00*  1.363e-02*cos(t)+ -1.187e-03*  1.657e-03*sin(t)),   2.675e-02 +2.000*(  1.187e-03*  1.363e-02*cos(t)+  1.000e+00*  1.657e-03*sin(t)) not
+# Age 80, p21 - p13
+set label "80" at   1.389e-01,  3.615e-02 center
+replot   1.389e-01+ 2.000*(  1.000e+00*  8.222e-03*cos(t)+ -3.844e-03*  1.674e-03*sin(t)),   3.615e-02 +2.000*(  3.844e-03*  8.222e-03*cos(t)+  1.000e+00*  1.674e-03*sin(t)) not
+# Age 85, p21 - p13
+set label "85" at   1.243e-01,  4.875e-02 center
+replot   1.243e-01+ 2.000*(  9.992e-01*  6.297e-03*cos(t)+ -4.035e-02*  2.643e-03*sin(t)),   4.875e-02 +2.000*(  4.035e-02*  6.297e-03*cos(t)+  9.992e-01*  2.643e-03*sin(t)) not
+# Age 90, p21 - p13
+set label "90" at   1.112e-01,  6.556e-02 center
+replot   1.112e-01+ 2.000*(  9.932e-01*  7.824e-03*cos(t)+ -1.163e-01*  5.151e-03*sin(t)),   6.556e-02 +2.000*(  1.163e-01*  7.824e-03*cos(t)+  9.932e-01*  5.151e-03*sin(t)) not
+# Age 95, p21 - p13
+set label "95" at   9.952e-02,  8.777e-02 center
+replot   9.952e-02+ 2.000*(  8.754e-01*  1.047e-02*cos(t)+ -4.835e-01*  9.400e-03*sin(t)),   8.777e-02 +2.000*(  4.835e-01*  1.047e-02*cos(t)+  8.754e-01*  9.400e-03*sin(t)) not
+set out "biaspar/varpijgrbiaspar121-13.png";replot;
+set parametric;unset label
+set log y;set log x; set xlabel "p23 (year-1)";set ylabel "p13 (year-1)"
+set ter png small
+set size 0.65,0.65
+set out "biaspar/varpijgrbiaspar123-13.png"
+set label "65" at   1.639e-01,  1.460e-02 center
+# Age 65, p23 - p13
+plot [-pi:pi]   1.639e-01+ 2.000*(  9.992e-01*  2.122e-02*cos(t)+  3.972e-02*  1.562e-03*sin(t)),   1.460e-02 +2.000*( -3.972e-02*  2.122e-02*cos(t)+  9.992e-01*  1.562e-03*sin(t)) not
+# Age 70, p23 - p13
+set label "70" at   1.707e-01,  1.977e-02 center
+replot   1.707e-01+ 2.000*(  9.987e-01*  1.756e-02*cos(t)+  5.147e-02*  1.527e-03*sin(t)),   1.977e-02 +2.000*( -5.147e-02*  1.756e-02*cos(t)+  9.987e-01*  1.527e-03*sin(t)) not
+# Age 75, p23 - p13
+set label "75" at   1.777e-01,  2.675e-02 center
+replot   1.777e-01+ 2.000*(  9.977e-01*  1.375e-02*cos(t)+  6.833e-02*  1.368e-03*sin(t)),   2.675e-02 +2.000*( -6.833e-02*  1.375e-02*cos(t)+  9.977e-01*  1.368e-03*sin(t)) not
+# Age 80, p23 - p13
+set label "80" at   1.850e-01,  3.615e-02 center
+replot   1.850e-01+ 2.000*(  9.949e-01*  1.009e-02*cos(t)+  1.005e-01*  1.339e-03*sin(t)),   3.615e-02 +2.000*( -1.005e-01*  1.009e-02*cos(t)+  9.949e-01*  1.339e-03*sin(t)) not
+# Age 85, p23 - p13
+set label "85" at   1.926e-01,  4.875e-02 center
+replot   1.926e-01+ 2.000*(  9.804e-01*  7.637e-03*cos(t)+  1.972e-01*  2.228e-03*sin(t)),   4.875e-02 +2.000*( -1.972e-01*  7.637e-03*cos(t)+  9.804e-01*  2.228e-03*sin(t)) not
+# Age 90, p23 - p13
+set label "90" at   2.004e-01,  6.556e-02 center
+replot   2.004e-01+ 2.000*(  9.254e-01*  8.649e-03*cos(t)+  3.789e-01*  4.358e-03*sin(t)),   6.556e-02 +2.000*( -3.789e-01*  8.649e-03*cos(t)+  9.254e-01*  4.358e-03*sin(t)) not
+# Age 95, p23 - p13
+set label "95" at   2.085e-01,  8.777e-02 center
+replot   2.085e-01+ 2.000*(  8.603e-01*  1.326e-02*cos(t)+  5.098e-01*  8.022e-03*sin(t)),   8.777e-02 +2.000*( -5.098e-01*  1.326e-02*cos(t)+  8.603e-01*  8.022e-03*sin(t)) not
+set out "biaspar/varpijgrbiaspar123-13.png";replot;
+set parametric;unset label
+set log y;set log x; set xlabel "p23 (year-1)";set ylabel "p21 (year-1)"
+set ter png small
+set size 0.65,0.65
+set out "biaspar/varpijgrbiaspar123-21.png"
+set label "65" at   1.639e-01,  1.936e-01 center
+# Age 65, p23 - p21
+plot [-pi:pi]   1.639e-01+ 2.000*(  3.376e-02*  3.110e-02*cos(t)+  9.994e-01*  2.119e-02*sin(t)),   1.936e-01 +2.000*( -9.994e-01*  3.110e-02*cos(t)+  3.376e-02*  2.119e-02*sin(t)) not
+# Age 70, p23 - p21
+set label "70" at   1.707e-01,  1.734e-01 center
+replot   1.707e-01+ 2.000*(  5.801e-02*  2.135e-02*cos(t)+  9.983e-01*  1.752e-02*sin(t)),   1.734e-01 +2.000*( -9.983e-01*  2.135e-02*cos(t)+  5.801e-02*  1.752e-02*sin(t)) not
+# Age 75, p23 - p21
+set label "75" at   1.777e-01,  1.552e-01 center
+replot   1.777e-01+ 2.000*(  8.046e-01*  1.382e-02*cos(t)+  5.938e-01*  1.352e-02*sin(t)),   1.552e-01 +2.000*( -5.938e-01*  1.382e-02*cos(t)+  8.046e-01*  1.352e-02*sin(t)) not
+# Age 80, p23 - p21
+set label "80" at   1.850e-01,  1.389e-01 center
+replot   1.850e-01+ 2.000*(  9.977e-01*  1.005e-02*cos(t)+  6.821e-02*  8.213e-03*sin(t)),   1.389e-01 +2.000*( -6.821e-02*  1.005e-02*cos(t)+  9.977e-01*  8.213e-03*sin(t)) not
+# Age 85, p23 - p21
+set label "85" at   1.926e-01,  1.243e-01 center
+replot   1.926e-01+ 2.000*(  9.869e-01*  7.531e-03*cos(t)+  1.615e-01*  6.256e-03*sin(t)),   1.243e-01 +2.000*( -1.615e-01*  7.531e-03*cos(t)+  9.869e-01*  6.256e-03*sin(t)) not
+# Age 90, p23 - p21
+set label "90" at   2.004e-01,  1.112e-01 center
+replot   2.004e-01+ 2.000*(  8.742e-01*  8.336e-03*cos(t)+  4.856e-01*  7.619e-03*sin(t)),   1.112e-01 +2.000*( -4.856e-01*  8.336e-03*cos(t)+  8.742e-01*  7.619e-03*sin(t)) not
+# Age 95, p23 - p21
+set label "95" at   2.085e-01,  9.952e-02 center
+replot   2.085e-01+ 2.000*(  9.838e-01*  1.218e-02*cos(t)+  1.795e-01*  1.016e-02*sin(t)),   9.952e-02 +2.000*( -1.795e-01*  1.218e-02*cos(t)+  9.838e-01*  1.016e-02*sin(t)) not
+set out "biaspar/varpijgrbiaspar123-21.png";replot;
+# Routine varevsij
+set noparametric;set nolabel; set ter png small;set size 0.65, 0.65
+ set log y; set nolog x;set xlabel "Age"; set ylabel "Force of mortality (year-1)";
+ plot "biaspar/prmorprev1-stablbased-rbiaspar.txt"  u 1:($3) not w l 1 
+ replot "biaspar/prmorprev1-stablbased-rbiaspar.txt"  u 1:(($3+1.96*$4)) t "95% interval" w l 2 
+ replot "biaspar/prmorprev1-stablbased-rbiaspar.txt"  u 1:(($3-1.96*$4)) not w l 2 
+set out "biaspar/varmuptjgr-stablbased-biaspar1.png";replot;