imach/html/doc/biaspar.gp - view

File: [Local Repository] / imach / html / doc / biaspar.gp
Revision 1.2: download - view: text, annotated - select for diffs
Tue Oct 25 11:53:26 2005 UTC (18 years, 8 months ago) by brouard
Branches: MAIN
CVS tags: HEAD

Version 0.98 MacOS/X Panther

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;