1: cd "D:/cygwin/home/agnes/hid/hidinstit/merge/bruit/"
2: set xlabel "Age"
3: set ylabel "Probability"
4: set ter gif small size 400,300
5: plot [65:100] "vplrparaaa.txt" every :::0::0 u 1:2 "%lf %lf (%lf) %*lf (%*lf)" t"Stationary prevalence" w l 0,"vplrparaaa.txt" every :::0::0 u 1:($2+2*$3) "%lf %lf (%lf) %*lf (%*lf)" t"95% CI" w l 1,"vplrparaaa.txt" every :::0::0 u 1:($2-2*$3) "%lf %lf (%lf) %*lf (%*lf)" t"" w l 1,"prparaaa.txt" every :::0::0 u 1:($2) t"Observed prevalence " w l 2
6: set out "vparaaa11.gif"
7: replot
8:
9: set xlabel "Age"
10: set ylabel "Probability"
11: set ter gif small size 400,300
12: plot [65:100] "vplrparaaa.txt" every :::0::0 u 1:2 "%lf %*lf (%*lf) %lf (%lf)" t"Stationary prevalence" w l 0,"vplrparaaa.txt" every :::0::0 u 1:($2+2*$3) "%lf %*lf (%*lf) %lf (%lf)" t"95% CI" w l 1,"vplrparaaa.txt" every :::0::0 u 1:($2-2*$3) "%lf %*lf (%*lf) %lf (%lf)" t"" w l 1,"prparaaa.txt" every :::0::0 u 1:($6) t"Observed prevalence " w l 2
13: set out "vparaaa21.gif"
14: replot
15:
16: set ylabel "Years"
17: set ter gif small size 400,300
18: plot [65:100] "trparaaa.txt" every :::0::0 u 1:2 "%lf %lf (%lf) %*lf (%*lf) %*lf (%*lf)" t"TLE" w l ,"trparaaa.txt" every :::0::0 u 1:($2-$3*2) "%lf %lf (%lf) %*lf (%*lf) %*lf (%*lf)" t"" w l 0,"trparaaa.txt" every :::0::0 u 1:($2+$3*2) "%lf %lf (%lf) %*lf (%*lf) %*lf (%*lf)" t"" w l 0,"trparaaa.txt" every :::0::0 u 1:2 "%lf %*lf (%*lf) %lf (%lf) %*lf (%*lf)" t"LE in state (1)" w l ,"trparaaa.txt" every :::0::0 u 1:($2-$3*2) "%lf %*lf (%*lf) %lf (%lf) %*lf (%*lf)" t"" w l 0,"trparaaa.txt" every :::0::0 u 1:($2+$3*2) "%lf %*lf (%*lf) %lf (%lf) %*lf (%*lf)" t"" w l 0,"trparaaa.txt" every :::0::0 u 1:2 "%lf %*lf (%*lf) %*lf (%*lf) %lf (%lf)" t"LE in state (2)" w l ,"trparaaa.txt" every :::0::0 u 1:($2-$3*2) "%lf %*lf (%*lf) %*lf (%*lf) %lf (%lf)" t"" w l 0,"trparaaa.txt" every :::0::0 u 1:($2+$3*2) "%lf %*lf (%*lf) %*lf (%*lf) %lf (%lf)" t"" w l 0
19: set out "eparaaa1.gif"
20: replot
21:
22: set ter gif small size 400,300
23: plot [65:100] "erparaaa.txt" every :::0::0 u 1:2 t "e11" w l,"erparaaa.txt" every :::0::0 u 1:3 t "e12" w l
24: set out "expparaaa11.gif"
25: replot
26:
27: set ter gif small size 400,300
28: plot [65:100] "erparaaa.txt" every :::0::0 u 1:4 t "e21" w l,"erparaaa.txt" every :::0::0 u 1:5 t "e22" w l
29: set out "expparaaa21.gif"
30: replot
31:
32: set xlabel "Age"
33: set ylabel "Probability"
34: set ter gif small size 400,300
35: plot [65:100] "pijrparaaa.txt" u ($1==1 ? ($3):1/0):($5/($4+$5)) t"prev(1,2)" w l,"pijrparaaa.txt" u ($1==1 ? ($3):1/0):($8/($7+$8)) t"prev(2,2)" w l
36: set out "pparaaa11.gif"
37: replot
38:
39: p1=-6.490654
40: p2=0.045790
41: p3=-8.307042
42: p4=0.058647
43: p5=-3.638365
44: p6=-0.021300
45: p7=-5.549177
46: p8=0.034565
47:
48: set ter gif small size 400,300
49: set log y
50: plot [65:100] 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"
51: set out "peparaaa1.gif"
52: replot
53:
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>