Return to mypar.plt CVS log

Up to [Local Repository] / imach096d / doc

jkhfskj

    1: cd "D:/cygwin/home/agnes/imach/imachcvs/imach-0.71/mytry/" 
    2: set xlabel "Age" 
    3: set ylabel "Probability" 
    4: set ter gif small size 400,300
    5: plot [70:100] "vplrmypar.txt" every :::0::0 u 1:2 "%lf %lf (%lf) %*lf (%*lf)" t"Stationary prevalence" w l 0,"vplrmypar.txt" every :::0::0 u 1:($2+2*$3) "%lf %lf (%lf) %*lf (%*lf)" t"95% CI" w l 1,"vplrmypar.txt" every :::0::0 u 1:($2-2*$3) "%lf %lf (%lf) %*lf (%*lf)" t"" w l 1,"prmypar.txt" every :::0::0 u 1:($2) t"Observed prevalence " w l 2
    6: set out "vmypar11.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" every :::0::0 u 1:2 "%lf %*lf (%*lf) %lf (%lf)" t"Stationary prevalence" w l 0,"vplrmypar.txt" every :::0::0 u 1:($2+2*$3) "%lf %*lf (%*lf) %lf (%lf)" t"95% CI" w l 1,"vplrmypar.txt" every :::0::0 u 1:($2-2*$3) "%lf %*lf (%*lf) %lf (%lf)" t"" w l 1,"prmypar.txt" every :::0::0 u 1:($6) t"Observed prevalence " w l 2
   13: set out "vmypar21.gif" 
   14: replot
   15: 
   16: set ylabel "Years" 
   17: set ter gif small size 400,300
   18: plot [70:100] "trmypar.txt" every :::0::0 u 1:2 "%lf %lf (%lf) %*lf (%*lf) %*lf (%*lf)" t"TLE" w l ,"trmypar.txt" every :::0::0 u 1:($2-$3*2) "%lf %lf (%lf) %*lf (%*lf) %*lf (%*lf)" t"" w l 0,"trmypar.txt" every :::0::0 u 1:($2+$3*2) "%lf %lf (%lf) %*lf (%*lf) %*lf (%*lf)" t"" w l 0,"trmypar.txt" every :::0::0 u 1:2 "%lf %*lf (%*lf) %lf (%lf) %*lf (%*lf)" t"LE in state (1)" w l ,"trmypar.txt" every :::0::0 u 1:($2-$3*2) "%lf %*lf (%*lf) %lf (%lf) %*lf (%*lf)" t"" w l 0,"trmypar.txt" every :::0::0 u 1:($2+$3*2) "%lf %*lf (%*lf) %lf (%lf) %*lf (%*lf)" t"" w l 0,"trmypar.txt" every :::0::0 u 1:2 "%lf %*lf (%*lf) %*lf (%*lf) %lf (%lf)" t"LE in state (2)" w l ,"trmypar.txt" every :::0::0 u 1:($2-$3*2) "%lf %*lf (%*lf) %*lf (%*lf) %lf (%lf)" t"" w l 0,"trmypar.txt" every :::0::0 u 1:($2+$3*2) "%lf %*lf (%*lf) %*lf (%*lf) %lf (%lf)" t"" w l 0
   19: set out "emypar1.gif" 
   20: replot
   21: 
   22: set ter gif small size 400,300
   23: plot [70:100] "ermypar.txt" every :::0::0 u 1:2 t "e11" w l,"ermypar.txt" every :::0::0 u 1:3 t "e12" w l
   24: set out "expmypar11.gif" 
   25: replot
   26: 
   27: set ter gif small size 400,300
   28: plot [70:100] "ermypar.txt" every :::0::0 u 1:4 t "e21" w l,"ermypar.txt" every :::0::0 u 1:5 t "e22" w l
   29: set out "expmypar21.gif" 
   30: replot
   31: 
   32: set xlabel "Age" 
   33: set ylabel "Probability" 
   34: set ter gif small size 400,300
   35: plot [70:100] "pijrmypar.txt" u ($1==1 ? ($3):1/0):($5/($4+$5)) t"prev(1,2)" w l,"pijrmypar.txt" u ($1==1 ? ($3):1/0):($8/($7+$8)) t"prev(2,2)" w l
   36: set out "pmypar11.gif" 
   37: replot
   38: 
   39: p1=-12.965571 
   40: p2=0.135127 
   41: p3=-7.482528 
   42: p4=0.068843 
   43: p5=-0.381775 
   44: p6=-0.010177 
   45: p7=-4.796242 
   46: p8=0.048192 
   47: 
   48: set ter gif small size 400,300
   49: set log y
   50: plot  [70: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 "pemypar1.gif" 
   52: replot
   53: