Return to biaspar.gp CVS log

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

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

1.5     ! brouard     1: 
        !             2: # 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
        !             3: # biaspar.gp
        !             4: set datafile missing 'NaNq'
        !             5: cd "/Users/nbrouard/Documents/imach/cvs/imach/build/osx/html/doc" 
        !             6: 
        !             7: # 1st: Period (stable) prevalence with CI: 'vpl' files
        !             8: 
        !             9: set out "biaspar/vbiaspar1_1.png" 
        !            10: 
        !            11: #set out "vbiaspar1_1.png" 
        !            12: set xlabel "Age" 
        !            13: set ylabel "Probability" 
        !            14: set ter png small size 320, 240
        !            15: 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
        !            16: set out "biaspar/vbiaspar2_1.png" 
        !            17: 
        !            18: #set out "vbiaspar2_1.png" 
        !            19: set xlabel "Age" 
        !            20: set ylabel "Probability" 
        !            21: set ter png small size 320, 240
        !            22: 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
        !            23: # 2nd: Total life expectancy with CI: 't' files
        !            24: 
        !            25: set out "biaspar/ebiaspar1.png" 
        !            26: set ylabel "Years" 
        !            27: set ter png small size 320, 240
        !            28: 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
        !            29: set out "biaspar/expbiaspar11.png" 
        !            30: set ter png small size 320, 240
        !            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:3 t "e12" w l ,"biaspar/erbiaspar.txt" every :::0::0 u 1:4 t "e1." w l
        !            32: set out "biaspar/expbiaspar21.png" 
        !            33: set ter png small size 320, 240
        !            34: 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
        !            35: #
        !            36: #
        !            37: #CV preval stable (period): 'pij' files, cov=1 state=1
        !            38: set out "biaspar/pbiaspar1_1.png" 
        !            39: set xlabel "Age" 
        !            40: set ylabel "Probability" 
        !            41: set ter png small size 320, 240
        !            42: unset log y
        !            43: 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
        !            44: 
        !            45: #
        !            46: #
        !            47: #CV preval stable (period): 'pij' files, cov=1 state=2
        !            48: set out "biaspar/pbiaspar2_1.png" 
        !            49: set xlabel "Age" 
        !            50: set ylabel "Probability" 
        !            51: set ter png small size 320, 240
        !            52: unset log y
        !            53: 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
        !            54: 
        !            55: ##############
        !            56: #MLE estimated parameters
        !            57: #############
        !            58: # initial state 1
        !            59: #   current state 2
        !            60: p1=-12.202711; p2=0.091764; 
        !            61: #   current state 3
        !            62: p3=-10.299238; p4=0.056626; 
        !            63: # initial state 2
        !            64: #   current state 1
        !            65: p5=-1.545619; p6=-0.035395; 
        !            66: #   current state 3
        !            67: p7=-5.565963; p8=0.016951; 
        !            68: ##############
        !            69: #
        !            70: 
        !            71: ##############
        !            72: #Graphics of of probabilities or incidences
        !            73: #############
        !            74: # logi(p12/p11)=a12+b12*age+c12age*age+d12*V1+e12*V1*age
        !            75: # logi(p12/p11)=p1 +p2*age +p3*age*age+ p4*V1+ p5*V1*age
        !            76: # logi(p13/p11)=a13+b13*age+c13age*age+d13*V1+e13*V1*age
        !            77: # logi(p13/p11)=p6 +p7*age +p8*age*age+ p9*V1+ p10*V1*age
        !            78: # p12+p13+p14+p11=1=p11(1+exp(a12+b12*age+c12age*age+d12*V1+e12*V1*age)
        !            79: #                      +exp(a13+b13*age+c13age*age+d13*V1+e13*V1*age)+...)
        !            80: # p11=1/(1+exp(a12+b12*age+c12age*age+d12*V1+e12*V1*age)
        !            81: #                      +exp(a13+b13*age+c13age*age+d13*V1+e13*V1*age)+...)
        !            82: # p12=exp(a12+b12*age+c12age*age+d12*V1+e12*V1*age)/
        !            83: #     (1+exp(a12+b12*age+c12age*age+d12*V1+e12*V1*age)
        !            84: #       +exp(a13+b13*age+c13age*age+d13*V1+e13*V1*age))
        !            85: #       +exp(a14+b14*age+c14age*age+d14*V1+e14*V1*age)+...)
        !            86: #
        !            87: # ng=1
        !            88: #   jk=1 to 2^0=1
        !            89: #    jk=1
        !            90: 
        !            91: set out "biaspar/pebiaspar1_1.png" 
        !            92: 
        !            93: set title "Probability"
        !            94: 
        !            95: set ter png small size 320, 240
        !            96: set log y
        !            97: 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
        !            98: #   jk=1 to 2^0=1
        !            99: #    jk=1
        !           100: 
        !           101: set out "biaspar/pebiaspar1_2.png" 
        !           102: 
        !           103: set ylabel "Quasi-incidence per year"
        !           104: 
        !           105: set ter png small size 320, 240
        !           106: set log y
        !           107: 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" 
        !           108: # Routine varprob
        !           109: set parametric;unset label
        !           110: set log y;set log x; set xlabel "p13 (year-1)";set ylabel "p12 (year-1)"
        !           111: set ter png small size 320, 240
        !           112: set out "biaspar/varpijgrbiaspar113-12.png"
        !           113: set label "65" at   1.597e-02,  2.337e-02 center
        !           114: # Age 65, p13 - p12
        !           115: 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
        !           116: # Age 70, p13 - p12
        !           117: set label "70" at   2.117e-02,  3.691e-02 center
        !           118: 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
        !           119: # Age 75, p13 - p12
        !           120: set label "75" at   2.803e-02,  5.827e-02 center
        !           121: 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
        !           122: # Age 80, p13 - p12
        !           123: set label "80" at   3.707e-02,  9.186e-02 center
        !           124: 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
        !           125: # Age 85, p13 - p12
        !           126: set label "85" at   4.893e-02,  1.445e-01 center
        !           127: 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
        !           128: # Age 90, p13 - p12
        !           129: set label "90" at   6.441e-02,  2.268e-01 center
        !           130: 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
        !           131: # Age 95, p13 - p12
        !           132: set label "95" at   8.441e-02,  3.543e-01 center
        !           133: 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
        !           134: set out "biaspar/varpijgrbiaspar113-12.png";replot;
        !           135: set parametric;unset label
        !           136: set log y;set log x; set xlabel "p21 (year-1)";set ylabel "p12 (year-1)"
        !           137: set ter png small size 320, 240
        !           138: set out "biaspar/varpijgrbiaspar121-12.png"
        !           139: set label "65" at   2.481e-01,  2.337e-02 center
        !           140: # Age 65, p21 - p12
        !           141: 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
        !           142: # Age 70, p21 - p12
        !           143: set label "70" at   2.084e-01,  3.691e-02 center
        !           144: 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
        !           145: # Age 75, p21 - p12
        !           146: set label "75" at   1.749e-01,  5.827e-02 center
        !           147: 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
        !           148: # Age 80, p21 - p12
        !           149: set label "80" at   1.467e-01,  9.186e-02 center
        !           150: 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
        !           151: # Age 85, p21 - p12
        !           152: set label "85" at   1.230e-01,  1.445e-01 center
        !           153: 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
        !           154: # Age 90, p21 - p12
        !           155: set label "90" at   1.031e-01,  2.268e-01 center
        !           156: 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
        !           157: # Age 95, p21 - p12
        !           158: set label "95" at   8.634e-02,  3.543e-01 center
        !           159: 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
        !           160: set out "biaspar/varpijgrbiaspar121-12.png";replot;
        !           161: set parametric;unset label
        !           162: set log y;set log x; set xlabel "p23 (year-1)";set ylabel "p12 (year-1)"
        !           163: set ter png small size 320, 240
        !           164: set out "biaspar/varpijgrbiaspar123-12.png"
        !           165: set label "65" at   1.338e-01,  2.337e-02 center
        !           166: # Age 65, p23 - p12
        !           167: 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
        !           168: # Age 70, p23 - p12
        !           169: set label "70" at   1.460e-01,  3.691e-02 center
        !           170: 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
        !           171: # Age 75, p23 - p12
        !           172: set label "75" at   1.591e-01,  5.827e-02 center
        !           173: 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
        !           174: # Age 80, p23 - p12
        !           175: set label "80" at   1.734e-01,  9.186e-02 center
        !           176: 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
        !           177: # Age 85, p23 - p12
        !           178: set label "85" at   1.889e-01,  1.445e-01 center
        !           179: 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
        !           180: # Age 90, p23 - p12
        !           181: set label "90" at   2.057e-01,  2.268e-01 center
        !           182: 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
        !           183: # Age 95, p23 - p12
        !           184: set label "95" at   2.238e-01,  3.543e-01 center
        !           185: 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
        !           186: set out "biaspar/varpijgrbiaspar123-12.png";replot;
        !           187: set parametric;unset label
        !           188: set log y;set log x; set xlabel "p21 (year-1)";set ylabel "p13 (year-1)"
        !           189: set ter png small size 320, 240
        !           190: set out "biaspar/varpijgrbiaspar121-13.png"
        !           191: set label "65" at   2.481e-01,  1.597e-02 center
        !           192: # Age 65, p21 - p13
        !           193: 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
        !           194: # Age 70, p21 - p13
        !           195: set label "70" at   2.084e-01,  2.117e-02 center
        !           196: 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
        !           197: # Age 75, p21 - p13
        !           198: set label "75" at   1.749e-01,  2.803e-02 center
        !           199: 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
        !           200: # Age 80, p21 - p13
        !           201: set label "80" at   1.467e-01,  3.707e-02 center
        !           202: 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
        !           203: # Age 85, p21 - p13
        !           204: set label "85" at   1.230e-01,  4.893e-02 center
        !           205: 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
        !           206: # Age 90, p21 - p13
        !           207: set label "90" at   1.031e-01,  6.441e-02 center
        !           208: 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
        !           209: # Age 95, p21 - p13
        !           210: set label "95" at   8.634e-02,  8.441e-02 center
        !           211: 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
        !           212: set out "biaspar/varpijgrbiaspar121-13.png";replot;
        !           213: set parametric;unset label
        !           214: set log y;set log x; set xlabel "p23 (year-1)";set ylabel "p13 (year-1)"
        !           215: set ter png small size 320, 240
        !           216: set out "biaspar/varpijgrbiaspar123-13.png"
        !           217: set label "65" at   1.338e-01,  1.597e-02 center
        !           218: # Age 65, p23 - p13
        !           219: 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
        !           220: # Age 70, p23 - p13
        !           221: set label "70" at   1.460e-01,  2.117e-02 center
        !           222: 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
        !           223: # Age 75, p23 - p13
        !           224: set label "75" at   1.591e-01,  2.803e-02 center
        !           225: 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
        !           226: # Age 80, p23 - p13
        !           227: set label "80" at   1.734e-01,  3.707e-02 center
        !           228: 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
        !           229: # Age 85, p23 - p13
        !           230: set label "85" at   1.889e-01,  4.893e-02 center
        !           231: 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
        !           232: # Age 90, p23 - p13
        !           233: set label "90" at   2.057e-01,  6.441e-02 center
        !           234: 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
        !           235: # Age 95, p23 - p13
        !           236: set label "95" at   2.238e-01,  8.441e-02 center
        !           237: 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
        !           238: set out "biaspar/varpijgrbiaspar123-13.png";replot;
        !           239: set parametric;unset label
        !           240: set log y;set log x; set xlabel "p23 (year-1)";set ylabel "p21 (year-1)"
        !           241: set ter png small size 320, 240
        !           242: set out "biaspar/varpijgrbiaspar123-21.png"
        !           243: set label "65" at   1.338e-01,  2.481e-01 center
        !           244: # Age 65, p23 - p21
        !           245: 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
        !           246: # Age 70, p23 - p21
        !           247: set label "70" at   1.460e-01,  2.084e-01 center
        !           248: 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
        !           249: # Age 75, p23 - p21
        !           250: set label "75" at   1.591e-01,  1.749e-01 center
        !           251: 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
        !           252: # Age 80, p23 - p21
        !           253: set label "80" at   1.734e-01,  1.467e-01 center
        !           254: 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
        !           255: # Age 85, p23 - p21
        !           256: set label "85" at   1.889e-01,  1.230e-01 center
        !           257: 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
        !           258: # Age 90, p23 - p21
        !           259: set label "90" at   2.057e-01,  1.031e-01 center
        !           260: 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
        !           261: # Age 95, p23 - p21
        !           262: set label "95" at   2.238e-01,  8.634e-02 center
        !           263: 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
        !           264: set out "biaspar/varpijgrbiaspar123-21.png";replot;
        !           265: # Routine varevsij
        !           266: unset parametric;unset label; set ter png small size 320, 240
        !           267:  set log y; unset log x;set xlabel "Age"; set ylabel "Force of mortality (year-1)";
        !           268:  plot "biaspar/prmorprev1-stablbased-rbiaspar.txt"  u 1:($3) not w l lt 1 
        !           269:  replot "biaspar/prmorprev1-stablbased-rbiaspar.txt"  u 1:(($3+1.96*$4)) t "95% interval" w l lt 2 
        !           270:  replot "biaspar/prmorprev1-stablbased-rbiaspar.txt"  u 1:(($3-1.96*$4)) not w l lt 2 
        !           271: set out "biaspar/varmuptjgr-stablbased-biaspar1.png";replot;