set xlabel "Age" 
set ylabel "Probability" 
set ter gif small size 400,300
plot [70:100] "pijrbiaspar.txt" u 2:($4/($3+$4)) t"prev(1,2)" w l,"pijrbiaspar.txt" u 2:($7/($6+$7)) t"prev(2,2)" w l
set out "pbiaspar1.gif" 
replot

set ter gif small size 400,300
plot [70:100] "vplrbiaspar.txt" u 1:2 "%lf %lf (%lf) %lf (%lf)" t"Stationary prevalence" w l 0,"vplrbiaspar.txt" u 1:($2+2*$3) "%lf %lf (%lf) %lf (%lf)" t"95% CI" w l 1,"vplrbiaspar.txt" u 1:($2-2*$3) "%lf %lf (%lf) %lf (%lf)" t"" w l 1,"prbiaspar.txt" u 1:($2) t"Observed prevalence " w l 2
set out "vbiaspar1.gif" 
replot

set ter gif small size 400,300
plot [70:100] "vplrbiaspar.txt" u 1:4 "%lf %lf (%lf) %lf (%lf)" t"Stationary prevalence" w l 0,"vplrbiaspar.txt" u 1:($4+2*$5) "%lf %lf (%lf) %lf (%lf)" t"95% CI" w l 1,"vplrbiaspar.txt" u 1:($4-2*$5) "%lf %lf (%lf) %lf (%lf)" t"" w l 1,"prbiaspar.txt" u 1:($6) t"Observed prevalence " w l 2
set out "vbiaspar2.gif" 
replot

set ylabel "Years" 
set ter gif small size 400,300
plot [70:100] "trbiaspar.txt" u 1:2 "%lf %lf (%lf) %lf (%lf) %lf (%lf)" t"TLE" w l ,"trbiaspar.txt" u 1:($2-2*$3) "%lf %lf (%lf) %lf (%lf) %lf (%lf)" t"" w l 0,"trbiaspar.txt" u 1:($2+2*$3) "%lf %lf (%lf) %lf (%lf) %lf (%lf)" t"" w l 0,"trbiaspar.txt" u 1:4 "%lf %lf (%lf) %lf (%lf) %lf (%lf)" t"LE in state (1)" w l ,"trbiaspar.txt" u 1:($4-2*$5) "%lf %lf (%lf) %lf (%lf) %lf (%lf)" t"" w l 0,"trbiaspar.txt" u 1:($4+2*$5) "%lf %lf (%lf) %lf (%lf) %lf (%lf)" t"" w l 0,"trbiaspar.txt" u 1:6 "%lf %lf (%lf) %lf (%lf) %lf (%lf)" t"LE in state (2)" w l ,"trbiaspar.txt" u 1:($6-2*$7) "%lf %lf (%lf) %lf (%lf) %lf (%lf)" t"" w l 0,"trbiaspar.txt" u 1:($6+2*$7) "%lf %lf (%lf) %lf (%lf) %lf (%lf)" t"" w l 0
set out "ebiaspar.gif" 
replot