1:
2: # Imach version 0.07a, May 2004, INED-EUROREVES
3: # biaspar.gp
4: set missing 'NaNq'
5: cd "D:/imachcvs/imach/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.245160
52: p2=0.092357
53: p3=-10.672078
54: p4=0.060971
55: p5=-2.645815
56: p6=-0.022320
57: p7=-4.773208
58: p8=0.007857
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.459e-002, 2.328e-002 center
83: # Age 65, p13 - p12
84: 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
85: # Age 70, p13 - p12
86: set label "70" at 1.976e-002, 3.688e-002 center
87: 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
88: # Age 75, p13 - p12
89: set label "75" at 2.674e-002, 5.839e-002 center
90: 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
91: # Age 80, p13 - p12
92: set label "80" at 3.614e-002, 9.233e-002 center
93: 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
94: # Age 85, p13 - p12
95: set label "85" at 4.875e-002, 1.457e-001 center
96: 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
97: # Age 90, p13 - p12
98: set label "90" at 6.557e-002, 2.292e-001 center
99: 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
100: # Age 95, p13 - p12
101: set label "95" at 8.778e-002, 3.591e-001 center
102: 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
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-001, 2.328e-002 center
110: # Age 65, p21 - p12
111: 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
112: # Age 70, p21 - p12
113: set label "70" at 1.734e-001, 3.688e-002 center
114: 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
115: # Age 75, p21 - p12
116: set label "75" at 1.552e-001, 5.839e-002 center
117: 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
118: # Age 80, p21 - p12
119: set label "80" at 1.389e-001, 9.233e-002 center
120: 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
121: # Age 85, p21 - p12
122: set label "85" at 1.243e-001, 1.457e-001 center
123: 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
124: # Age 90, p21 - p12
125: set label "90" at 1.112e-001, 2.292e-001 center
126: 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
127: # Age 95, p21 - p12
128: set label "95" at 9.953e-002, 3.591e-001 center
129: 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
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.640e-001, 2.328e-002 center
137: # Age 65, p23 - p12
138: 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
139: # Age 70, p23 - p12
140: set label "70" at 1.708e-001, 3.688e-002 center
141: 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
142: # Age 75, p23 - p12
143: set label "75" at 1.778e-001, 5.839e-002 center
144: 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
145: # Age 80, p23 - p12
146: set label "80" at 1.851e-001, 9.233e-002 center
147: 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
148: # Age 85, p23 - p12
149: set label "85" at 1.926e-001, 1.457e-001 center
150: 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
151: # Age 90, p23 - p12
152: set label "90" at 2.004e-001, 2.292e-001 center
153: 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
154: # Age 95, p23 - p12
155: set label "95" at 2.085e-001, 3.591e-001 center
156: 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
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-001, 1.459e-002 center
164: # Age 65, p21 - p13
165: 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
166: # Age 70, p21 - p13
167: set label "70" at 1.734e-001, 1.976e-002 center
168: 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
169: # Age 75, p21 - p13
170: set label "75" at 1.552e-001, 2.674e-002 center
171: 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
172: # Age 80, p21 - p13
173: set label "80" at 1.389e-001, 3.614e-002 center
174: 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
175: # Age 85, p21 - p13
176: set label "85" at 1.243e-001, 4.875e-002 center
177: 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
178: # Age 90, p21 - p13
179: set label "90" at 1.112e-001, 6.557e-002 center
180: 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
181: # Age 95, p21 - p13
182: set label "95" at 9.953e-002, 8.778e-002 center
183: 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
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.640e-001, 1.459e-002 center
191: # Age 65, p23 - p13
192: 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
193: # Age 70, p23 - p13
194: set label "70" at 1.708e-001, 1.976e-002 center
195: 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
196: # Age 75, p23 - p13
197: set label "75" at 1.778e-001, 2.674e-002 center
198: 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
199: # Age 80, p23 - p13
200: set label "80" at 1.851e-001, 3.614e-002 center
201: 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
202: # Age 85, p23 - p13
203: set label "85" at 1.926e-001, 4.875e-002 center
204: 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
205: # Age 90, p23 - p13
206: set label "90" at 2.004e-001, 6.557e-002 center
207: 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
208: # Age 95, p23 - p13
209: set label "95" at 2.085e-001, 8.778e-002 center
210: 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
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.640e-001, 1.936e-001 center
218: # Age 65, p23 - p21
219: 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
220: # Age 70, p23 - p21
221: set label "70" at 1.708e-001, 1.734e-001 center
222: 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
223: # Age 75, p23 - p21
224: set label "75" at 1.778e-001, 1.552e-001 center
225: 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
226: # Age 80, p23 - p21
227: set label "80" at 1.851e-001, 1.389e-001 center
228: 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
229: # Age 85, p23 - p21
230: set label "85" at 1.926e-001, 1.243e-001 center
231: 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
232: # Age 90, p23 - p21
233: set label "90" at 2.004e-001, 1.112e-001 center
234: 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
235: # Age 95, p23 - p21
236: set label "95" at 2.085e-001, 9.953e-002 center
237: 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
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;
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>