1: cd "../mytry"
2: set xlabel "Age"
3: set ylabel "Probability"
4: set ter gif small size 400,300
5: plot [70:100] "vplrmypar.txt" u 1:2 "%lf %lf (%lf) %lf (%lf)" t"Stationary prevalence" w l 0,"vplrmypar.txt" u 1:($2+2*$3) "%lf %lf (%lf) %lf (%lf)" t"95% CI" w l 1,"vplrmypar.txt" u 1:($2-2*$3) "%lf %lf (%lf) %lf (%lf)" t"" w l 1,"prmypar.txt" u 1:($2) t"Observed prevalence " w l 2
6: set out "vmypar1.gif"
7: replot
8:
9: set xlabel "Age"
10: set ylabel "Probability"
11: set ter gif small size 400,300
12: plot [70:100] "vplrmypar.txt" u 1:4 "%lf %lf (%lf) %lf (%lf)" t"Stationary prevalence" w l 0,"vplrmypar.txt" u 1:($4+2*$5) "%lf %lf (%lf) %lf (%lf)" t"95% CI" w l 1,"vplrmypar.txt" u 1:($4-2*$5) "%lf %lf (%lf) %lf (%lf)" t"" w l 1,"prmypar.txt" u 1:($6) t"Observed prevalence " w l 2
13: set out "vmypar2.gif"
14: replot
15:
16: set ylabel "Years"
17: set ter gif small size 400,300
18: plot [70:100] "trmypar.txt" u 1:2 "%lf %lf (%lf) %lf (%lf) %lf (%lf)" t"TLE" w l ,"trmypar.txt" u 1:($2-2*$3) "%lf %lf (%lf) %lf (%lf) %lf (%lf)" t"" w l 0,"trmypar.txt" u 1:($2+2*$3) "%lf %lf (%lf) %lf (%lf) %lf (%lf)" t"" w l 0,"trmypar.txt" u 1:4 "%lf %lf (%lf) %lf (%lf) %lf (%lf)" t"LE in state (1)" w l ,"trmypar.txt" u 1:($4-2*$5) "%lf %lf (%lf) %lf (%lf) %lf (%lf)" t"" w l 0,"trmypar.txt" u 1:($4+2*$5) "%lf %lf (%lf) %lf (%lf) %lf (%lf)" t"" w l 0,"trmypar.txt" u 1:6 "%lf %lf (%lf) %lf (%lf) %lf (%lf)" t"LE in state (2)" w l ,"trmypar.txt" u 1:($6-2*$7) "%lf %lf (%lf) %lf (%lf) %lf (%lf)" t"" w l 0,"trmypar.txt" u 1:($6+2*$7) "%lf %lf (%lf) %lf (%lf) %lf (%lf)" t"" w l 0
19: set out "emypar.gif"
20: replot
21:
22: set ter gif small size 400,300
23: plot [70:100] "ermypar.txt" u 1:2 t "e11" w l,"ermypar.txt" u 1:3 t "e12" w l
24: set out "exmypar1.gif"
25: replot
26:
27: set ter gif small size 400,300
28: plot [70:100] "ermypar.txt" u 1:4 t "e21" w l,"ermypar.txt" u 1:5 t "e22" w l
29: set out "exmypar2.gif"
30: replot
31:
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>