Return to biaspar.gp CVS log

Up to [Local Repository] / imach / html / doc

Annotation of imach/html/doc/biaspar.gp, revision 1.2

1.2     ! brouard     1: 
        !             2: # Imach version 0.98, September 2005, INED-EUROREVES 
        !             3: # biaspar.gp
        !             4: set missing 'NaNq'
        !             5: cd "/Users/brouard/bin/html/doc" 
        !             6: 
        !             7: set out "biaspar/vbiaspar11.png" 
        !             8: 
        !             9: #set out "vbiaspar11.png" 
        !            10: set xlabel "Age" 
        !            11: set ylabel "Probability" 
        !            12: set ter png small
        !            13: set size 0.65,0.65
        !            14: 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
        !            15: set out "biaspar/vbiaspar21.png" 
        !            16: 
        !            17: #set out "vbiaspar21.png" 
        !            18: set xlabel "Age" 
        !            19: set ylabel "Probability" 
        !            20: set ter png small
        !            21: set size 0.65,0.65
        !            22: 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
        !            23: set out "biaspar/ebiaspar1.png" 
        !            24: set ylabel "Years" 
        !            25: set ter png small
        !            26: set size 0.65,0.65
        !            27: 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
        !            28: set out "biaspar/expbiaspar11.png" 
        !            29: set ter png small
        !            30: set size 0.65,0.65
        !            31: 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
        !            32: set out "biaspar/expbiaspar21.png" 
        !            33: set ter png small
        !            34: set size 0.65,0.65
        !            35: 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
        !            36: set out "biaspar/pbiaspar11.png" 
        !            37: set xlabel "Age" 
        !            38: set ylabel "Probability" 
        !            39: set ter png small
        !            40: set size 0.65,0.65
        !            41: unset log y
        !            42: 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
        !            43: 
        !            44: set out "biaspar/pbiaspar21.png" 
        !            45: set xlabel "Age" 
        !            46: set ylabel "Probability" 
        !            47: set ter png small
        !            48: set size 0.65,0.65
        !            49: unset log y
        !            50: 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
        !            51: p1=-12.245504 
        !            52: p2=0.092361 
        !            53: p3=-10.670953 
        !            54: p4=0.060958 
        !            55: p5=-2.645034 
        !            56: p6=-0.022329 
        !            57: p7=-4.775501 
        !            58: p8=0.007884 
        !            59: 
        !            60: set out "biaspar/pebiaspar11.png" 
        !            61: 
        !            62: set title "Probability"
        !            63: 
        !            64: set ter png small
        !            65: set size 0.65,0.65
        !            66: set log y
        !            67: 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" 
        !            68: set out "biaspar/pebiaspar12.png" 
        !            69: 
        !            70: set ylabel "Quasi-incidence per year"
        !            71: 
        !            72: set ter png small
        !            73: set size 0.65,0.65
        !            74: set log y
        !            75: 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" 
        !            76: # Routine varprob
        !            77: set parametric;unset label
        !            78: set log y;set log x; set xlabel "p13 (year-1)";set ylabel "p12 (year-1)"
        !            79: set ter png small
        !            80: set size 0.65,0.65
        !            81: set out "biaspar/varpijgrbiaspar113-12.png"
        !            82: set label "65" at   1.460e-02,  2.328e-02 center
        !            83: # Age 65, p13 - p12
        !            84: 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
        !            85: # Age 70, p13 - p12
        !            86: set label "70" at   1.977e-02,  3.688e-02 center
        !            87: 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
        !            88: # Age 75, p13 - p12
        !            89: set label "75" at   2.675e-02,  5.839e-02 center
        !            90: 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
        !            91: # Age 80, p13 - p12
        !            92: set label "80" at   3.615e-02,  9.232e-02 center
        !            93: 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
        !            94: # Age 85, p13 - p12
        !            95: set label "85" at   4.875e-02,  1.457e-01 center
        !            96: 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
        !            97: # Age 90, p13 - p12
        !            98: set label "90" at   6.556e-02,  2.292e-01 center
        !            99: 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
        !           100: # Age 95, p13 - p12
        !           101: set label "95" at   8.777e-02,  3.591e-01 center
        !           102: 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
        !           103: set out "biaspar/varpijgrbiaspar113-12.png";replot;
        !           104: set parametric;unset label
        !           105: set log y;set log x; set xlabel "p21 (year-1)";set ylabel "p12 (year-1)"
        !           106: set ter png small
        !           107: set size 0.65,0.65
        !           108: set out "biaspar/varpijgrbiaspar121-12.png"
        !           109: set label "65" at   1.936e-01,  2.328e-02 center
        !           110: # Age 65, p21 - p12
        !           111: 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
        !           112: # Age 70, p21 - p12
        !           113: set label "70" at   1.734e-01,  3.688e-02 center
        !           114: 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
        !           115: # Age 75, p21 - p12
        !           116: set label "75" at   1.552e-01,  5.839e-02 center
        !           117: 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
        !           118: # Age 80, p21 - p12
        !           119: set label "80" at   1.389e-01,  9.232e-02 center
        !           120: 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
        !           121: # Age 85, p21 - p12
        !           122: set label "85" at   1.243e-01,  1.457e-01 center
        !           123: 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
        !           124: # Age 90, p21 - p12
        !           125: set label "90" at   1.112e-01,  2.292e-01 center
        !           126: 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
        !           127: # Age 95, p21 - p12
        !           128: set label "95" at   9.952e-02,  3.591e-01 center
        !           129: 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
        !           130: set out "biaspar/varpijgrbiaspar121-12.png";replot;
        !           131: set parametric;unset label
        !           132: set log y;set log x; set xlabel "p23 (year-1)";set ylabel "p12 (year-1)"
        !           133: set ter png small
        !           134: set size 0.65,0.65
        !           135: set out "biaspar/varpijgrbiaspar123-12.png"
        !           136: set label "65" at   1.639e-01,  2.328e-02 center
        !           137: # Age 65, p23 - p12
        !           138: 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
        !           139: # Age 70, p23 - p12
        !           140: set label "70" at   1.707e-01,  3.688e-02 center
        !           141: 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
        !           142: # Age 75, p23 - p12
        !           143: set label "75" at   1.777e-01,  5.839e-02 center
        !           144: 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
        !           145: # Age 80, p23 - p12
        !           146: set label "80" at   1.850e-01,  9.232e-02 center
        !           147: 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
        !           148: # Age 85, p23 - p12
        !           149: set label "85" at   1.926e-01,  1.457e-01 center
        !           150: 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
        !           151: # Age 90, p23 - p12
        !           152: set label "90" at   2.004e-01,  2.292e-01 center
        !           153: 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
        !           154: # Age 95, p23 - p12
        !           155: set label "95" at   2.085e-01,  3.591e-01 center
        !           156: 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
        !           157: set out "biaspar/varpijgrbiaspar123-12.png";replot;
        !           158: set parametric;unset label
        !           159: set log y;set log x; set xlabel "p21 (year-1)";set ylabel "p13 (year-1)"
        !           160: set ter png small
        !           161: set size 0.65,0.65
        !           162: set out "biaspar/varpijgrbiaspar121-13.png"
        !           163: set label "65" at   1.936e-01,  1.460e-02 center
        !           164: # Age 65, p21 - p13
        !           165: 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
        !           166: # Age 70, p21 - p13
        !           167: set label "70" at   1.734e-01,  1.977e-02 center
        !           168: 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
        !           169: # Age 75, p21 - p13
        !           170: set label "75" at   1.552e-01,  2.675e-02 center
        !           171: 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
        !           172: # Age 80, p21 - p13
        !           173: set label "80" at   1.389e-01,  3.615e-02 center
        !           174: 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
        !           175: # Age 85, p21 - p13
        !           176: set label "85" at   1.243e-01,  4.875e-02 center
        !           177: 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
        !           178: # Age 90, p21 - p13
        !           179: set label "90" at   1.112e-01,  6.556e-02 center
        !           180: 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
        !           181: # Age 95, p21 - p13
        !           182: set label "95" at   9.952e-02,  8.777e-02 center
        !           183: 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
        !           184: set out "biaspar/varpijgrbiaspar121-13.png";replot;
        !           185: set parametric;unset label
        !           186: set log y;set log x; set xlabel "p23 (year-1)";set ylabel "p13 (year-1)"
        !           187: set ter png small
        !           188: set size 0.65,0.65
        !           189: set out "biaspar/varpijgrbiaspar123-13.png"
        !           190: set label "65" at   1.639e-01,  1.460e-02 center
        !           191: # Age 65, p23 - p13
        !           192: 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
        !           193: # Age 70, p23 - p13
        !           194: set label "70" at   1.707e-01,  1.977e-02 center
        !           195: 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
        !           196: # Age 75, p23 - p13
        !           197: set label "75" at   1.777e-01,  2.675e-02 center
        !           198: 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
        !           199: # Age 80, p23 - p13
        !           200: set label "80" at   1.850e-01,  3.615e-02 center
        !           201: 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
        !           202: # Age 85, p23 - p13
        !           203: set label "85" at   1.926e-01,  4.875e-02 center
        !           204: 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
        !           205: # Age 90, p23 - p13
        !           206: set label "90" at   2.004e-01,  6.556e-02 center
        !           207: 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
        !           208: # Age 95, p23 - p13
        !           209: set label "95" at   2.085e-01,  8.777e-02 center
        !           210: 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
        !           211: set out "biaspar/varpijgrbiaspar123-13.png";replot;
        !           212: set parametric;unset label
        !           213: set log y;set log x; set xlabel "p23 (year-1)";set ylabel "p21 (year-1)"
        !           214: set ter png small
        !           215: set size 0.65,0.65
        !           216: set out "biaspar/varpijgrbiaspar123-21.png"
        !           217: set label "65" at   1.639e-01,  1.936e-01 center
        !           218: # Age 65, p23 - p21
        !           219: 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
        !           220: # Age 70, p23 - p21
        !           221: set label "70" at   1.707e-01,  1.734e-01 center
        !           222: 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
        !           223: # Age 75, p23 - p21
        !           224: set label "75" at   1.777e-01,  1.552e-01 center
        !           225: 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
        !           226: # Age 80, p23 - p21
        !           227: set label "80" at   1.850e-01,  1.389e-01 center
        !           228: 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
        !           229: # Age 85, p23 - p21
        !           230: set label "85" at   1.926e-01,  1.243e-01 center
        !           231: 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
        !           232: # Age 90, p23 - p21
        !           233: set label "90" at   2.004e-01,  1.112e-01 center
        !           234: 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
        !           235: # Age 95, p23 - p21
        !           236: set label "95" at   2.085e-01,  9.952e-02 center
        !           237: 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
        !           238: set out "biaspar/varpijgrbiaspar123-21.png";replot;
        !           239: # Routine varevsij
        !           240: set noparametric;set nolabel; set ter png small;set size 0.65, 0.65
        !           241:  set log y; set nolog x;set xlabel "Age"; set ylabel "Force of mortality (year-1)";
        !           242:  plot "biaspar/prmorprev1-stablbased-rbiaspar.txt"  u 1:($3) not w l 1 
        !           243:  replot "biaspar/prmorprev1-stablbased-rbiaspar.txt"  u 1:(($3+1.96*$4)) t "95% interval" w l 2 
        !           244:  replot "biaspar/prmorprev1-stablbased-rbiaspar.txt"  u 1:(($3-1.96*$4)) not w l 2 
        !           245: set out "biaspar/varmuptjgr-stablbased-biaspar1.png";replot;