Return to biaspar.plt CVS log

Up to [Local Repository] / imach096d / doc

Annotation of imach096d/doc/biaspar.plt, revision 1.2

1.2     ! brouard     1: cd "E:/usr/local/src/imach/Copie4doc/" 
1.1       brouard     2: set xlabel "Age" 
                      3: set ylabel "Probability" 
                      4: set ter gif small size 400,300
                      5: plot [70:95] "vplrbiaspar.txt" every :::0::0 u 1:2 "%lf %lf (%lf) %*lf (%*lf)" t"Stationary prevalence" w l 0,"vplrbiaspar.txt" every :::0::0 u 1:($2+2*$3) "%lf %lf (%lf) %*lf (%*lf)" t"95% CI" w l 1,"vplrbiaspar.txt" every :::0::0 u 1:($2-2*$3) "%lf %lf (%lf) %*lf (%*lf)" t"" w l 1,"prbiaspar.txt" every :::0::0 u 1:($2) t"Observed prevalence " w l 2
                      6: set out "vbiaspar11.gif" 
                      7: replot
                      8: 
                      9: set xlabel "Age" 
                     10: set ylabel "Probability" 
                     11: set ter gif small size 400,300
                     12: plot [70:95] "vplrbiaspar.txt" every :::0::0 u 1:2 "%lf %*lf (%*lf) %lf (%lf)" t"Stationary prevalence" w l 0,"vplrbiaspar.txt" every :::0::0 u 1:($2+2*$3) "%lf %*lf (%*lf) %lf (%lf)" t"95% CI" w l 1,"vplrbiaspar.txt" every :::0::0 u 1:($2-2*$3) "%lf %*lf (%*lf) %lf (%lf)" t"" w l 1,"prbiaspar.txt" every :::0::0 u 1:($6) t"Observed prevalence " w l 2
                     13: set out "vbiaspar21.gif" 
                     14: replot
                     15: 
                     16: set ylabel "Years" 
                     17: set ter gif small size 400,300
                     18: plot [70:95] "trbiaspar.txt" every :::0::0 u 1:2 "%lf %lf (%lf) %*lf (%*lf) %*lf (%*lf)" t"TLE" w l ,"trbiaspar.txt" every :::0::0 u 1:($2-$3*2) "%lf %lf (%lf) %*lf (%*lf) %*lf (%*lf)" t"" w l 0,"trbiaspar.txt" every :::0::0 u 1:($2+$3*2) "%lf %lf (%lf) %*lf (%*lf) %*lf (%*lf)" t"" w l 0,"trbiaspar.txt" every :::0::0 u 1:2 "%lf %*lf (%*lf) %lf (%lf) %*lf (%*lf)" t"LE in state (1)" w l ,"trbiaspar.txt" every :::0::0 u 1:($2-$3*2) "%lf %*lf (%*lf) %lf (%lf) %*lf (%*lf)" t"" w l 0,"trbiaspar.txt" every :::0::0 u 1:($2+$3*2) "%lf %*lf (%*lf) %lf (%lf) %*lf (%*lf)" t"" w l 0,"trbiaspar.txt" every :::0::0 u 1:2 "%lf %*lf (%*lf) %*lf (%*lf) %lf (%lf)" t"LE in state (2)" w l ,"trbiaspar.txt" every :::0::0 u 1:($2-$3*2) "%lf %*lf (%*lf) %*lf (%*lf) %lf (%lf)" t"" w l 0,"trbiaspar.txt" every :::0::0 u 1:($2+$3*2) "%lf %*lf (%*lf) %*lf (%*lf) %lf (%lf)" t"" w l 0
                     19: set out "ebiaspar1.gif" 
                     20: replot
                     21: 
                     22: set ter gif small size 400,300
                     23: plot [70:95] "erbiaspar.txt" every :::0::0 u 1:2 t "e11" w l,"erbiaspar.txt" every :::0::0 u 1:3 t "e12" w l
                     24: set out "expbiaspar11.gif" 
                     25: replot
                     26: 
                     27: set ter gif small size 400,300
                     28: plot [70:95] "erbiaspar.txt" every :::0::0 u 1:4 t "e21" w l,"erbiaspar.txt" every :::0::0 u 1:5 t "e22" w l
                     29: set out "expbiaspar21.gif" 
                     30: replot
                     31: 
                     32: set xlabel "Age" 
                     33: set ylabel "Probability" 
                     34: set ter gif small size 400,300
                     35: plot [70:95] "pijrbiaspar.txt" u ($1==1 ? ($3):1/0):($5/($4+$5)) t"prev(1,2)" w l,"pijrbiaspar.txt" u ($1==1 ? ($3):1/0):($8/($7+$8)) t"prev(2,2)" w l
                     36: set out "pbiaspar11.gif" 
                     37: replot
                     38: 
                     39: p1=-12.290174 
                     40: p2=0.092161 
                     41: p3=-9.155590 
                     42: p4=0.046627 
                     43: p5=-2.629849 
                     44: p6=-0.022030 
                     45: p7=-7.958519 
                     46: p8=0.042614 
                     47: 
                     48: set ter gif small size 400,300
                     49: set log y
                     50: 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" 
                     51: set out "pebiaspar1.gif" 
                     52: replot
                     53: