[BACK]Return to biaspar.gp CVS log [TXT] [DIR] Up to [Local Repository] / imach / html / doc

Diff for /imach/html/doc/biaspar.gp between versions 1.3 and 1.4

version 1.3, 2006/02/09 15:20:33 version 1.4, 2015/08/04 07:14:39
Line 1 Line 0
   
 # Imach version 0.98b, January 2006, INED-EUROREVES  
 # biaspar.gp  
 set missing 'NaNq'  
 cd "C:/Program Files/IMaCh/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.245240  
 p2=0.092358  
 p3=-10.671890  
 p4=0.060969  
 p5=-2.645345  
 p6=-0.022325  
 p7=-4.773317  
 p8=0.007859  
   
 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*( 7.081e-001* 2.074e-003*cos(t)+ 7.061e-001* 1.411e-003*sin(t)),  2.328e-002 +2.000*(-7.061e-001* 2.074e-003*cos(t)+ 7.081e-001* 1.411e-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.375e-001* 2.333e-003*cos(t)+ 8.433e-001* 1.486e-003*sin(t)),  3.688e-002 +2.000*(-8.433e-001* 2.333e-003*cos(t)+ 5.375e-001* 1.486e-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.212e-001* 2.525e-003*cos(t)+ 9.070e-001* 1.399e-003*sin(t)),  5.839e-002 +2.000*(-9.070e-001* 2.525e-003*cos(t)+ 4.212e-001* 1.399e-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.981e-001* 2.739e-003*cos(t)+ 9.173e-001* 1.383e-003*sin(t)),  9.233e-002 +2.000*(-9.173e-001* 2.739e-003*cos(t)+ 3.981e-001* 1.383e-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.975e-001* 4.130e-003*cos(t)+ 9.176e-001* 2.278e-003*sin(t)),  1.457e-001 +2.000*(-9.176e-001* 4.130e-003*cos(t)+ 3.975e-001* 2.278e-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.928e-001* 8.973e-003*cos(t)+ 9.562e-001* 4.707e-003*sin(t)),  2.292e-001 +2.000*(-9.562e-001* 8.973e-003*cos(t)+ 2.928e-001* 4.707e-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.165e-001* 1.969e-002*cos(t)+ 9.763e-001* 8.913e-003*sin(t)),  3.591e-001 +2.000*(-9.763e-001* 1.969e-002*cos(t)+ 2.165e-001* 8.913e-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.139e-002*cos(t)+-1.878e-002* 1.672e-003*sin(t)),  2.328e-002 +2.000*( 1.878e-002* 3.139e-002*cos(t)+ 9.998e-001* 1.672e-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.995e-001* 2.153e-002*cos(t)+-3.254e-002* 2.005e-003*sin(t)),  3.688e-002 +2.000*( 3.254e-002* 2.153e-002*cos(t)+ 9.995e-001* 2.005e-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.984e-001* 1.374e-002*cos(t)+-5.600e-002* 2.239e-003*sin(t)),  5.839e-002 +2.000*( 5.600e-002* 1.374e-002*cos(t)+ 9.984e-001* 2.239e-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.950e-001* 8.279e-003*cos(t)+-9.943e-002* 2.449e-003*sin(t)),  9.233e-002 +2.000*( 9.943e-002* 8.279e-003*cos(t)+ 9.950e-001* 2.449e-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.670e-001* 6.440e-003*cos(t)+-2.549e-001* 3.655e-003*sin(t)),  1.457e-001 +2.000*( 2.549e-001* 6.440e-003*cos(t)+ 9.670e-001* 3.655e-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.855e-001* 9.526e-003*cos(t)+-8.106e-001* 6.807e-003*sin(t)),  2.292e-001 +2.000*( 8.106e-001* 9.526e-003*cos(t)+ 5.855e-001* 6.807e-003*sin(t)) not  
 # Age 95, p21 - p12  
 set label "95" at  9.952e-002, 3.591e-001 center  
 replot  9.952e-002+ 2.000*( 2.226e-001* 1.969e-002*cos(t)+-9.749e-001* 9.591e-003*sin(t)),  3.591e-001 +2.000*( 9.749e-001* 1.969e-002*cos(t)+ 2.226e-001* 9.591e-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.122e-002*cos(t)+-3.110e-002* 1.646e-003*sin(t)),  2.328e-002 +2.000*( 3.110e-002* 2.122e-002*cos(t)+ 9.995e-001* 1.646e-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.755e-002*cos(t)+-4.814e-002* 1.950e-003*sin(t)),  3.688e-002 +2.000*( 4.814e-002* 1.755e-002*cos(t)+ 9.988e-001* 1.950e-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.972e-001* 1.375e-002*cos(t)+-7.534e-002* 2.131e-003*sin(t)),  5.839e-002 +2.000*( 7.534e-002* 1.375e-002*cos(t)+ 9.972e-001* 2.131e-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.926e-001* 1.011e-002*cos(t)+-1.215e-001* 2.277e-003*sin(t)),  9.233e-002 +2.000*( 1.215e-001* 1.011e-002*cos(t)+ 9.926e-001* 2.277e-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.723e-001* 7.665e-003*cos(t)+-2.338e-001* 3.559e-003*sin(t)),  1.457e-001 +2.000*( 2.338e-001* 7.665e-003*cos(t)+ 9.723e-001* 3.559e-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.216e-001* 9.497e-003*cos(t)+-7.833e-001* 7.226e-003*sin(t)),  2.292e-001 +2.000*( 7.833e-001* 9.497e-003*cos(t)+ 6.216e-001* 7.226e-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.331e-001* 1.967e-002*cos(t)+-9.725e-001* 1.154e-002*sin(t)),  3.591e-001 +2.000*( 9.725e-001* 1.967e-002*cos(t)+ 2.331e-001* 1.154e-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.138e-002*cos(t)+-2.829e-003* 1.772e-003*sin(t)),  1.459e-002 +2.000*( 2.829e-003* 3.138e-002*cos(t)+ 1.000e+000* 1.772e-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.152e-002*cos(t)+-3.296e-003* 1.771e-003*sin(t)),  1.976e-002 +2.000*( 3.296e-003* 2.152e-002*cos(t)+ 1.000e+000* 1.771e-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.372e-002*cos(t)+-3.023e-003* 1.655e-003*sin(t)),  2.674e-002 +2.000*( 3.023e-003* 1.372e-002*cos(t)+ 1.000e+000* 1.655e-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.242e-003*cos(t)+-4.787e-003* 1.672e-003*sin(t)),  3.614e-002 +2.000*( 4.787e-003* 8.242e-003*cos(t)+ 1.000e+000* 1.672e-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.991e-001* 6.301e-003*cos(t)+-4.129e-002* 2.647e-003*sin(t)),  4.875e-002 +2.000*( 4.129e-002* 6.301e-003*cos(t)+ 9.991e-001* 2.647e-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.920e-001* 7.882e-003*cos(t)+-1.263e-001* 5.157e-003*sin(t)),  6.557e-002 +2.000*( 1.263e-001* 7.882e-003*cos(t)+ 9.920e-001* 5.157e-003*sin(t)) not  
 # Age 95, p21 - p13  
 set label "95" at  9.952e-002, 8.778e-002 center  
 replot  9.952e-002+ 2.000*( 8.713e-001* 1.061e-002*cos(t)+-4.908e-001* 9.377e-003*sin(t)),  8.778e-002 +2.000*( 4.908e-001* 1.061e-002*cos(t)+ 8.713e-001* 9.377e-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.992e-001* 2.122e-002*cos(t)+ 3.961e-002* 1.564e-003*sin(t)),  1.459e-002 +2.000*(-3.961e-002* 2.122e-002*cos(t)+ 9.992e-001* 1.564e-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.987e-001* 1.755e-002*cos(t)+ 5.131e-002* 1.529e-003*sin(t)),  1.976e-002 +2.000*(-5.131e-002* 1.755e-002*cos(t)+ 9.987e-001* 1.529e-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.977e-001* 1.374e-002*cos(t)+ 6.811e-002* 1.369e-003*sin(t)),  2.674e-002 +2.000*(-6.811e-002* 1.374e-002*cos(t)+ 9.977e-001* 1.369e-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.009e-002*cos(t)+ 1.004e-001* 1.338e-003*sin(t)),  3.614e-002 +2.000*(-1.004e-001* 1.009e-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.803e-001* 7.637e-003*cos(t)+ 1.977e-001* 2.231e-003*sin(t)),  4.875e-002 +2.000*(-1.977e-001* 7.637e-003*cos(t)+ 9.803e-001* 2.231e-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.247e-001* 8.660e-003*cos(t)+ 3.806e-001* 4.366e-003*sin(t)),  6.557e-002 +2.000*(-3.806e-001* 8.660e-003*cos(t)+ 9.247e-001* 4.366e-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.590e-001* 1.328e-002*cos(t)+ 5.121e-001* 8.035e-003*sin(t)),  8.778e-002 +2.000*(-5.121e-001* 1.328e-002*cos(t)+ 8.590e-001* 8.035e-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*( 5.483e-002* 3.141e-002*cos(t)+ 9.985e-001* 2.117e-002*sin(t)),  1.936e-001 +2.000*(-9.985e-001* 3.141e-002*cos(t)+ 5.483e-002* 2.117e-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*( 9.437e-002* 2.155e-002*cos(t)+ 9.955e-001* 1.749e-002*sin(t)),  1.734e-001 +2.000*(-9.955e-001* 2.155e-002*cos(t)+ 9.437e-002* 1.749e-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*( 7.002e-001* 1.395e-002*cos(t)+ 7.140e-001* 1.347e-002*sin(t)),  1.552e-001 +2.000*(-7.140e-001* 1.395e-002*cos(t)+ 7.002e-001* 1.347e-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.962e-001* 1.005e-002*cos(t)+ 8.694e-002* 8.227e-003*sin(t)),  1.389e-001 +2.000*(-8.694e-002* 1.005e-002*cos(t)+ 9.962e-001* 8.227e-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.869e-001* 7.529e-003*cos(t)+ 1.612e-001* 6.261e-003*sin(t)),  1.243e-001 +2.000*(-1.612e-001* 7.529e-003*cos(t)+ 9.869e-001* 6.261e-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.505e-001* 8.378e-003*cos(t)+ 5.260e-001* 7.632e-003*sin(t)),  1.112e-001 +2.000*(-5.260e-001* 8.378e-003*cos(t)+ 8.505e-001* 7.632e-003*sin(t)) not  
 # Age 95, p23 - p21  
 set label "95" at  2.085e-001, 9.952e-002 center  
 replot  2.085e-001+ 2.000*( 9.743e-001* 1.222e-002*cos(t)+ 2.251e-001* 1.022e-002*sin(t)),  9.952e-002 +2.000*(-2.251e-001* 1.222e-002*cos(t)+ 9.743e-001* 1.022e-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;  
Removed from v.1.3  
changed lines
  Added in v.1.4

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>