Return to graph.plt CVS log

Up to [Local Repository] / imach096d / src

Annotation of imach096d/src/graph.plt, revision 1.1

1.1     ! lievre      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: