|
version 1.3, 2006/02/09 15:20:33
|
version 1.5, 2015/08/18 16:43:19
|
|
Line 1
|
Line 1
|
|
|
|
| # Imach version 0.98b, January 2006, INED-EUROREVES
|
# Imach version 0.98q5, August 2015,INED-EUROREVES-Institut de longevite-Japan Society for the Promotion of Science (Grant-in-Aid for Scientific Research 25293121), Intel Software 2015 |
| # biaspar.gp
|
# biaspar.gp |
| set missing 'NaNq'
|
set datafile missing 'NaNq' |
| cd "C:/Program Files/IMaCh/doc/"
|
cd "/Users/nbrouard/Documents/imach/cvs/imach/build/osx/html/doc" |
|
|
|
| set out "biaspar/vbiaspar11.png"
|
# 1st: Period (stable) prevalence with CI: 'vpl' files |
|
|
|
| #set out "vbiaspar11.png"
|
set out "biaspar/vbiaspar1_1.png" |
| set xlabel "Age"
|
|
| set ylabel "Probability"
|
#set out "vbiaspar1_1.png" |
| set ter png small
|
set xlabel "Age" |
| set size 0.65,0.65
|
set ylabel "Probability" |
| plot [70:95] "biaspar/vplrbiaspar.txt" every :::0::0 u 1:2 "%lf %lf (%lf) %*lf (%*lf)" t"Stable prevalence" w l 0,"biaspar/vplrbiaspar.txt" every :::0::0 u 1:($2+1.96*$3) "%lf %lf (%lf) %*lf (%*lf)" t"95% CI" w l 1,"biaspar/vplrbiaspar.txt" every :::0::0 u 1:($2-1.96*$3) "%lf %lf (%lf) %*lf (%*lf)" t"" w l 1,"biaspar/prbiaspar.txt" every :::0::0 u 1:($2) t"Observed prevalence " w l 2
|
set ter png small size 320, 240 |
| set out "biaspar/vbiaspar21.png"
|
plot [70:95] "biaspar/vplrbiaspar.txt" every :::0::0 u 1:2 "%lf %lf (%lf) %*lf (%*lf)" t"Period (stable) prevalence" w l lt 0,"biaspar/vplrbiaspar.txt" every :::0::0 u 1:($2+1.96*$3) "%lf %lf (%lf) %*lf (%*lf)" t"95% CI" w l lt 1,"biaspar/vplrbiaspar.txt" every :::0::0 u 1:($2-1.96*$3) "%lf %lf (%lf) %*lf (%*lf)" t"" w l lt 1,"biaspar/prbiaspar.txt" every :::0::0 u 1:($2) t"Observed prevalence " w l lt 2 |
|
|
set out "biaspar/vbiaspar2_1.png" |
| #set out "vbiaspar21.png"
|
|
| set xlabel "Age"
|
#set out "vbiaspar2_1.png" |
| set ylabel "Probability"
|
set xlabel "Age" |
| set ter png small
|
set ylabel "Probability" |
| set size 0.65,0.65
|
set ter png small size 320, 240 |
| plot [70:95] "biaspar/vplrbiaspar.txt" every :::0::0 u 1:2 "%lf %*lf (%*lf) %lf (%lf)" t"Stable prevalence" w l 0,"biaspar/vplrbiaspar.txt" every :::0::0 u 1:($2+1.96*$3) "%lf %*lf (%*lf) %lf (%lf)" t"95% CI" w l 1,"biaspar/vplrbiaspar.txt" every :::0::0 u 1:($2-1.96*$3) "%lf %*lf (%*lf) %lf (%lf)" t"" w l 1,"biaspar/prbiaspar.txt" every :::0::0 u 1:($6) t"Observed prevalence " w l 2
|
plot [70:95] "biaspar/vplrbiaspar.txt" every :::0::0 u 1:2 "%lf %*lf (%*lf) %lf (%lf)" t"Period (stable) prevalence" w l lt 0,"biaspar/vplrbiaspar.txt" every :::0::0 u 1:($2+1.96*$3) "%lf %*lf (%*lf) %lf (%lf)" t"95% CI" w l lt 1,"biaspar/vplrbiaspar.txt" every :::0::0 u 1:($2-1.96*$3) "%lf %*lf (%*lf) %lf (%lf)" t"" w l lt 1,"biaspar/prbiaspar.txt" every :::0::0 u 1:($6) t"Observed prevalence " w l lt 2 |
| set out "biaspar/ebiaspar1.png"
|
# 2nd: Total life expectancy with CI: 't' files |
| set ylabel "Years"
|
|
| set ter png small
|
set out "biaspar/ebiaspar1.png" |
| set size 0.65,0.65
|
set ylabel "Years" |
| plot [70:95] "biaspar/trbiaspar.txt" every :::0::0 u 1:2 "%lf %lf (%lf) %*lf (%*lf) %*lf (%*lf)" t"TLE" w l ,"biaspar/trbiaspar.txt" every :::0::0 u 1:($2-$3*2) "%lf %lf (%lf) %*lf (%*lf) %*lf (%*lf)" t"" w l 0,"biaspar/trbiaspar.txt" every :::0::0 u 1:($2+$3*2) "%lf %lf (%lf) %*lf (%*lf) %*lf (%*lf)" t"" w l 0,"biaspar/trbiaspar.txt" every :::0::0 u 1:2 "%lf %*lf (%*lf) %lf (%lf) %*lf (%*lf)" t"LE in state (1)" w l ,"biaspar/trbiaspar.txt" every :::0::0 u 1:($2-$3*2) "%lf %*lf (%*lf) %lf (%lf) %*lf (%*lf)" t"" w l 0,"biaspar/trbiaspar.txt" every :::0::0 u 1:($2+$3*2) "%lf %*lf (%*lf) %lf (%lf) %*lf (%*lf)" t"" w l 0,"biaspar/trbiaspar.txt" every :::0::0 u 1:2 "%lf %*lf (%*lf) %*lf (%*lf) %lf (%lf)" t"LE in state (2)" w l ,"biaspar/trbiaspar.txt" every :::0::0 u 1:($2-$3*2) "%lf %*lf (%*lf) %*lf (%*lf) %lf (%lf)" t"" w l 0,"biaspar/trbiaspar.txt" every :::0::0 u 1:($2+$3*2) "%lf %*lf (%*lf) %*lf (%*lf) %lf (%lf)" t"" w l 0
|
set ter png small size 320, 240 |
| set out "biaspar/expbiaspar11.png"
|
plot [70:95] "biaspar/trbiaspar.txt" every :::0::0 u 1:2 "%lf %lf (%lf) %*lf (%*lf) %*lf (%*lf)" t"TLE" w l ,"biaspar/trbiaspar.txt" every :::0::0 u 1:($2-$3*2) "%lf %lf (%lf) %*lf (%*lf) %*lf (%*lf)" t"" w l lt 0,"biaspar/trbiaspar.txt" every :::0::0 u 1:($2+$3*2) "%lf %lf (%lf) %*lf (%*lf) %*lf (%*lf)" t"" w l lt 0,"biaspar/trbiaspar.txt" every :::0::0 u 1:2 "%lf %*lf (%*lf) %lf (%lf) %*lf (%*lf)" t"LE in state (1)" w l ,"biaspar/trbiaspar.txt" every :::0::0 u 1:($2-$3*2) "%lf %*lf (%*lf) %lf (%lf) %*lf (%*lf)" t"" w l lt 0,"biaspar/trbiaspar.txt" every :::0::0 u 1:($2+$3*2) "%lf %*lf (%*lf) %lf (%lf) %*lf (%*lf)" t"" w l lt 0,"biaspar/trbiaspar.txt" every :::0::0 u 1:2 "%lf %*lf (%*lf) %*lf (%*lf) %lf (%lf)" t"LE in state (2)" w l ,"biaspar/trbiaspar.txt" every :::0::0 u 1:($2-$3*2) "%lf %*lf (%*lf) %*lf (%*lf) %lf (%lf)" t"" w l lt 0,"biaspar/trbiaspar.txt" every :::0::0 u 1:($2+$3*2) "%lf %*lf (%*lf) %*lf (%*lf) %lf (%lf)" t"" w l lt 0 |
| set ter png small
|
set out "biaspar/expbiaspar11.png" |
| set size 0.65,0.65
|
set ter png small size 320, 240 |
| plot [70:95] "biaspar/erbiaspar.txt" every :::0::0 u 1:2 t "e11" w l ,"biaspar/erbiaspar.txt" every :::0::0 u 1:4 t "e12" w l
|
plot [70:95] "biaspar/erbiaspar.txt" every :::0::0 u 1:2 t "e11" w l ,"biaspar/erbiaspar.txt" every :::0::0 u 1:3 t "e12" w l ,"biaspar/erbiaspar.txt" every :::0::0 u 1:4 t "e1." w l |
| set out "biaspar/expbiaspar21.png"
|
set out "biaspar/expbiaspar21.png" |
| set ter png small
|
set ter png small size 320, 240 |
| set size 0.65,0.65
|
plot [70:95] "biaspar/erbiaspar.txt" every :::0::0 u 1:5 t "e21" w l ,"biaspar/erbiaspar.txt" every :::0::0 u 1:6 t "e22" w l ,"biaspar/erbiaspar.txt" every :::0::0 u 1:7 t "e2." w l |
| plot [70:95] "biaspar/erbiaspar.txt" every :::0::0 u 1:6 t "e21" w l ,"biaspar/erbiaspar.txt" every :::0::0 u 1:8 t "e22" w l
|
# |
| set out "biaspar/pbiaspar11.png"
|
# |
| set xlabel "Age"
|
#CV preval stable (period): 'pij' files, cov=1 state=1 |
| set ylabel "Probability"
|
set out "biaspar/pbiaspar1_1.png" |
| set ter png small
|
set xlabel "Age" |
| set size 0.65,0.65
|
set ylabel "Probability" |
| unset log y
|
set ter png small size 320, 240 |
| plot [70:95] "biaspar/pijrbiaspar.txt" u ($1==1 ? ($3):1/0):($5/($4+$5)) t"prev(1,2)" w l,"biaspar/pijrbiaspar.txt" u ($1==1 ? ($3):1/0):($8/($7+$8)) t"prev(2,2)" w l
|
unset log y |
|
|
plot [70:95] "biaspar/pijrbiaspar.txt" u ($1==1 ? ($3):1/0):($4/($4+$5)) t "prev(1,1)" w l, '' u ($1==1 ? ($3):1/0):($7/($7+$8)) t "prev(2,1)" w l |
| set out "biaspar/pbiaspar21.png"
|
|
| set xlabel "Age"
|
# |
| set ylabel "Probability"
|
# |
| set ter png small
|
#CV preval stable (period): 'pij' files, cov=1 state=2 |
| set size 0.65,0.65
|
set out "biaspar/pbiaspar2_1.png" |
| unset log y
|
set xlabel "Age" |
| plot [70:95] "biaspar/pijrbiaspar.txt" u ($1==1 ? ($3):1/0):($6/($4+$5)) t"prev(2,3)" w l,"biaspar/pijrbiaspar.txt" u ($1==1 ? ($3):1/0):($12/($10+$11)) t"prev(3,3)" w l
|
set ylabel "Probability" |
| p1=-12.245240
|
set ter png small size 320, 240 |
| p2=0.092358
|
unset log y |
| p3=-10.671890
|
plot [70:95] "biaspar/pijrbiaspar.txt" u ($1==1 ? ($3):1/0):($5/($4+$5)) t "prev(1,2)" w l, '' u ($1==1 ? ($3):1/0):($8/($7+$8)) t "prev(2,2)" w l |
| p4=0.060969
|
|
| p5=-2.645345
|
############## |
| p6=-0.022325
|
#MLE estimated parameters |
| p7=-4.773317
|
############# |
| p8=0.007859
|
# initial state 1 |
|
|
# current state 2 |
| set out "biaspar/pebiaspar11.png"
|
p1=-12.202711; p2=0.091764; |
|
|
# current state 3 |
| set title "Probability"
|
p3=-10.299238; p4=0.056626; |
|
|
# initial state 2 |
| set ter png small
|
# current state 1 |
| set size 0.65,0.65
|
p5=-1.545619; p6=-0.035395; |
| set log y
|
# current state 3 |
| 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"
|
p7=-5.565963; p8=0.016951; |
| set out "biaspar/pebiaspar12.png"
|
############## |
|
|
# |
| set ylabel "Quasi-incidence per year"
|
|
|
|
############## |
| set ter png small
|
#Graphics of of probabilities or incidences |
| set size 0.65,0.65
|
############# |
| set log y
|
# logi(p12/p11)=a12+b12*age+c12age*age+d12*V1+e12*V1*age |
| plot [70:95] 12.000000*exp(p1+p2*x)/(1+exp(p1+p2*x)+exp(p3+p4*x)) t "p12" , 12.000000*exp(p3+p4*x)/(1+exp(p1+p2*x)+exp(p3+p4*x)) t "p13" , 12.000000*exp(p5+p6*x)/(1+exp(p5+p6*x)+exp(p7+p8*x)) t "p21" , 12.000000*exp(p7+p8*x)/(1+exp(p5+p6*x)+exp(p7+p8*x)) t "p23"
|
# logi(p12/p11)=p1 +p2*age +p3*age*age+ p4*V1+ p5*V1*age |
| # Routine varprob
|
# logi(p13/p11)=a13+b13*age+c13age*age+d13*V1+e13*V1*age |
| set parametric;unset label
|
# logi(p13/p11)=p6 +p7*age +p8*age*age+ p9*V1+ p10*V1*age |
| set log y;set log x; set xlabel "p13 (year-1)";set ylabel "p12 (year-1)"
|
# p12+p13+p14+p11=1=p11(1+exp(a12+b12*age+c12age*age+d12*V1+e12*V1*age) |
| set ter png small
|
# +exp(a13+b13*age+c13age*age+d13*V1+e13*V1*age)+...) |
| set size 0.65,0.65
|
# p11=1/(1+exp(a12+b12*age+c12age*age+d12*V1+e12*V1*age) |
| set out "biaspar/varpijgrbiaspar113-12.png"
|
# +exp(a13+b13*age+c13age*age+d13*V1+e13*V1*age)+...) |
| set label "65" at 1.459e-002, 2.328e-002 center
|
# p12=exp(a12+b12*age+c12age*age+d12*V1+e12*V1*age)/ |
| # Age 65, p13 - p12
|
# (1+exp(a12+b12*age+c12age*age+d12*V1+e12*V1*age) |
| plot [-pi:pi] 1.459e-002+ 2.000*( 7.081e-001* 2.074e-003*cos(t)+ 7.061e-001* 1.411e-003*sin(t)), 2.328e-002 +2.000*(-7.061e-001* 2.074e-003*cos(t)+ 7.081e-001* 1.411e-003*sin(t)) not
|
# +exp(a13+b13*age+c13age*age+d13*V1+e13*V1*age)) |
| # Age 70, p13 - p12
|
# +exp(a14+b14*age+c14age*age+d14*V1+e14*V1*age)+...) |
| set label "70" at 1.976e-002, 3.688e-002 center
|
# |
| replot 1.976e-002+ 2.000*( 5.375e-001* 2.333e-003*cos(t)+ 8.433e-001* 1.486e-003*sin(t)), 3.688e-002 +2.000*(-8.433e-001* 2.333e-003*cos(t)+ 5.375e-001* 1.486e-003*sin(t)) not
|
# ng=1 |
| # Age 75, p13 - p12
|
# jk=1 to 2^0=1 |
| set label "75" at 2.674e-002, 5.839e-002 center
|
# jk=1 |
| replot 2.674e-002+ 2.000*( 4.212e-001* 2.525e-003*cos(t)+ 9.070e-001* 1.399e-003*sin(t)), 5.839e-002 +2.000*(-9.070e-001* 2.525e-003*cos(t)+ 4.212e-001* 1.399e-003*sin(t)) not
|
|
| # Age 80, p13 - p12
|
set out "biaspar/pebiaspar1_1.png" |
| set label "80" at 3.614e-002, 9.233e-002 center
|
|
| replot 3.614e-002+ 2.000*( 3.981e-001* 2.739e-003*cos(t)+ 9.173e-001* 1.383e-003*sin(t)), 9.233e-002 +2.000*(-9.173e-001* 2.739e-003*cos(t)+ 3.981e-001* 1.383e-003*sin(t)) not
|
set title "Probability" |
| # Age 85, p13 - p12
|
|
| set label "85" at 4.875e-002, 1.457e-001 center
|
set ter png small size 320, 240 |
| replot 4.875e-002+ 2.000*( 3.975e-001* 4.130e-003*cos(t)+ 9.176e-001* 2.278e-003*sin(t)), 1.457e-001 +2.000*(-9.176e-001* 4.130e-003*cos(t)+ 3.975e-001* 2.278e-003*sin(t)) not
|
set log y |
| # Age 90, p13 - p12
|
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" # ng=2 |
| set label "90" at 6.557e-002, 2.292e-001 center
|
# jk=1 to 2^0=1 |
| replot 6.557e-002+ 2.000*( 2.928e-001* 8.973e-003*cos(t)+ 9.562e-001* 4.707e-003*sin(t)), 2.292e-001 +2.000*(-9.562e-001* 8.973e-003*cos(t)+ 2.928e-001* 4.707e-003*sin(t)) not
|
# jk=1 |
| # Age 95, p13 - p12
|
|
| set label "95" at 8.778e-002, 3.591e-001 center
|
set out "biaspar/pebiaspar1_2.png" |
| replot 8.778e-002+ 2.000*( 2.165e-001* 1.969e-002*cos(t)+ 9.763e-001* 8.913e-003*sin(t)), 3.591e-001 +2.000*(-9.763e-001* 1.969e-002*cos(t)+ 2.165e-001* 8.913e-003*sin(t)) not
|
|
| set out "biaspar/varpijgrbiaspar113-12.png";replot;
|
set ylabel "Quasi-incidence per year" |
| set parametric;unset label
|
|
| set log y;set log x; set xlabel "p21 (year-1)";set ylabel "p12 (year-1)"
|
set ter png small size 320, 240 |
| set ter png small
|
set log y |
| set size 0.65,0.65
|
plot [70:95] 12.000000*exp(p1+p2*x)/(1+exp(p1+p2*x)+exp(p3+p4*x)) t "p12" , 12.000000*exp(p3+p4*x)/(1+exp(p1+p2*x)+exp(p3+p4*x)) t "p13" , 12.000000*exp(p5+p6*x)/(1+exp(p5+p6*x)+exp(p7+p8*x)) t "p21" , 12.000000*exp(p7+p8*x)/(1+exp(p5+p6*x)+exp(p7+p8*x)) t "p23" |
| set out "biaspar/varpijgrbiaspar121-12.png"
|
# Routine varprob |
| set label "65" at 1.936e-001, 2.328e-002 center
|
set parametric;unset label |
| # Age 65, p21 - p12
|
set log y;set log x; set xlabel "p13 (year-1)";set ylabel "p12 (year-1)" |
| plot [-pi:pi] 1.936e-001+ 2.000*( 9.998e-001* 3.139e-002*cos(t)+-1.878e-002* 1.672e-003*sin(t)), 2.328e-002 +2.000*( 1.878e-002* 3.139e-002*cos(t)+ 9.998e-001* 1.672e-003*sin(t)) not
|
set ter png small size 320, 240 |
| # Age 70, p21 - p12
|
set out "biaspar/varpijgrbiaspar113-12.png" |
| set label "70" at 1.734e-001, 3.688e-002 center
|
set label "65" at 1.597e-02, 2.337e-02 center |
| replot 1.734e-001+ 2.000*( 9.995e-001* 2.153e-002*cos(t)+-3.254e-002* 2.005e-003*sin(t)), 3.688e-002 +2.000*( 3.254e-002* 2.153e-002*cos(t)+ 9.995e-001* 2.005e-003*sin(t)) not
|
# Age 65, p13 - p12 |
| # Age 75, p21 - p12
|
plot [-pi:pi] 1.597e-02+ 2.000*( 7.325e-01* 2.084e-03*cos(t)+ 6.807e-01* 1.465e-03*sin(t)), 2.337e-02 +2.000*( -6.807e-01* 2.084e-03*cos(t)+ 7.325e-01* 1.465e-03*sin(t)) not |
| set label "75" at 1.552e-001, 5.839e-002 center
|
# Age 70, p13 - p12 |
| replot 1.552e-001+ 2.000*( 9.984e-001* 1.374e-002*cos(t)+-5.600e-002* 2.239e-003*sin(t)), 5.839e-002 +2.000*( 5.600e-002* 1.374e-002*cos(t)+ 9.984e-001* 2.239e-003*sin(t)) not
|
set label "70" at 2.117e-02, 3.691e-02 center |
| # Age 80, p21 - p12
|
replot 2.117e-02+ 2.000*( 5.240e-01* 2.315e-03*cos(t)+ 8.517e-01* 1.521e-03*sin(t)), 3.691e-02 +2.000*( -8.517e-01* 2.315e-03*cos(t)+ 5.240e-01* 1.521e-03*sin(t)) not |
| set label "80" at 1.389e-001, 9.233e-002 center
|
# Age 75, p13 - p12 |
| replot 1.389e-001+ 2.000*( 9.950e-001* 8.279e-003*cos(t)+-9.943e-002* 2.449e-003*sin(t)), 9.233e-002 +2.000*( 9.943e-002* 8.279e-003*cos(t)+ 9.950e-001* 2.449e-003*sin(t)) not
|
set label "75" at 2.803e-02, 5.827e-02 center |
| # Age 85, p21 - p12
|
replot 2.803e-02+ 2.000*( 3.929e-01* 2.501e-03*cos(t)+ 9.196e-01* 1.403e-03*sin(t)), 5.827e-02 +2.000*( -9.196e-01* 2.501e-03*cos(t)+ 3.929e-01* 1.403e-03*sin(t)) not |
| set label "85" at 1.243e-001, 1.457e-001 center
|
# Age 80, p13 - p12 |
| replot 1.243e-001+ 2.000*( 9.670e-001* 6.440e-003*cos(t)+-2.549e-001* 3.655e-003*sin(t)), 1.457e-001 +2.000*( 2.549e-001* 6.440e-003*cos(t)+ 9.670e-001* 3.655e-003*sin(t)) not
|
set label "80" at 3.707e-02, 9.186e-02 center |
| # Age 90, p21 - p12
|
replot 3.707e-02+ 2.000*( 3.828e-01* 2.718e-03*cos(t)+ 9.238e-01* 1.387e-03*sin(t)), 9.186e-02 +2.000*( -9.238e-01* 2.718e-03*cos(t)+ 3.828e-01* 1.387e-03*sin(t)) not |
| set label "90" at 1.112e-001, 2.292e-001 center
|
# Age 85, p13 - p12 |
| replot 1.112e-001+ 2.000*( 5.855e-001* 9.526e-003*cos(t)+-8.106e-001* 6.807e-003*sin(t)), 2.292e-001 +2.000*( 8.106e-001* 9.526e-003*cos(t)+ 5.855e-001* 6.807e-003*sin(t)) not
|
set label "85" at 4.893e-02, 1.445e-01 center |
| # Age 95, p21 - p12
|
replot 4.893e-02+ 2.000*( 4.026e-01* 4.082e-03*cos(t)+ 9.154e-01* 2.265e-03*sin(t)), 1.445e-01 +2.000*( -9.154e-01* 4.082e-03*cos(t)+ 4.026e-01* 2.265e-03*sin(t)) not |
| set label "95" at 9.952e-002, 3.591e-001 center
|
# Age 90, p13 - p12 |
| replot 9.952e-002+ 2.000*( 2.226e-001* 1.969e-002*cos(t)+-9.749e-001* 9.591e-003*sin(t)), 3.591e-001 +2.000*( 9.749e-001* 1.969e-002*cos(t)+ 2.226e-001* 9.591e-003*sin(t)) not
|
set label "90" at 6.441e-02, 2.268e-01 center |
| set out "biaspar/varpijgrbiaspar121-12.png";replot;
|
replot 6.441e-02+ 2.000*( 2.874e-01* 8.787e-03*cos(t)+ 9.578e-01* 4.578e-03*sin(t)), 2.268e-01 +2.000*( -9.578e-01* 8.787e-03*cos(t)+ 2.874e-01* 4.578e-03*sin(t)) not |
| set parametric;unset label
|
# Age 95, p13 - p12 |
| set log y;set log x; set xlabel "p23 (year-1)";set ylabel "p12 (year-1)"
|
set label "95" at 8.441e-02, 3.543e-01 center |
| set ter png small
|
replot 8.441e-02+ 2.000*( 2.019e-01* 1.921e-02*cos(t)+ 9.794e-01* 8.453e-03*sin(t)), 3.543e-01 +2.000*( -9.794e-01* 1.921e-02*cos(t)+ 2.019e-01* 8.453e-03*sin(t)) not |
| set size 0.65,0.65
|
set out "biaspar/varpijgrbiaspar113-12.png";replot; |
| set out "biaspar/varpijgrbiaspar123-12.png"
|
set parametric;unset label |
| set label "65" at 1.640e-001, 2.328e-002 center
|
set log y;set log x; set xlabel "p21 (year-1)";set ylabel "p12 (year-1)" |
| # Age 65, p23 - p12
|
set ter png small size 320, 240 |
| plot [-pi:pi] 1.640e-001+ 2.000*( 9.995e-001* 2.122e-002*cos(t)+-3.110e-002* 1.646e-003*sin(t)), 2.328e-002 +2.000*( 3.110e-002* 2.122e-002*cos(t)+ 9.995e-001* 1.646e-003*sin(t)) not
|
set out "biaspar/varpijgrbiaspar121-12.png" |
| # Age 70, p23 - p12
|
set label "65" at 2.481e-01, 2.337e-02 center |
| set label "70" at 1.708e-001, 3.688e-002 center
|
# Age 65, p21 - p12 |
| replot 1.708e-001+ 2.000*( 9.988e-001* 1.755e-002*cos(t)+-4.814e-002* 1.950e-003*sin(t)), 3.688e-002 +2.000*( 4.814e-002* 1.755e-002*cos(t)+ 9.988e-001* 1.950e-003*sin(t)) not
|
plot [-pi:pi] 2.481e-01+ 2.000*( 9.999e-01* 3.884e-02*cos(t)+ -1.528e-02* 1.677e-03*sin(t)), 2.337e-02 +2.000*( 1.528e-02* 3.884e-02*cos(t)+ 9.999e-01* 1.677e-03*sin(t)) not |
| # Age 75, p23 - p12
|
# Age 70, p21 - p12 |
| set label "75" at 1.778e-001, 5.839e-002 center
|
set label "70" at 2.084e-01, 3.691e-02 center |
| replot 1.778e-001+ 2.000*( 9.972e-001* 1.375e-002*cos(t)+-7.534e-002* 2.131e-003*sin(t)), 5.839e-002 +2.000*( 7.534e-002* 1.375e-002*cos(t)+ 9.972e-001* 2.131e-003*sin(t)) not
|
replot 2.084e-01+ 2.000*( 9.996e-01* 2.493e-02*cos(t)+ -2.848e-02* 2.005e-03*sin(t)), 3.691e-02 +2.000*( 2.848e-02* 2.493e-02*cos(t)+ 9.996e-01* 2.005e-03*sin(t)) not |
| # Age 80, p23 - p12
|
# Age 75, p21 - p12 |
| set label "80" at 1.851e-001, 9.233e-002 center
|
set label "75" at 1.749e-01, 5.827e-02 center |
| replot 1.851e-001+ 2.000*( 9.926e-001* 1.011e-002*cos(t)+-1.215e-001* 2.277e-003*sin(t)), 9.233e-002 +2.000*( 1.215e-001* 1.011e-002*cos(t)+ 9.926e-001* 2.277e-003*sin(t)) not
|
replot 1.749e-01+ 2.000*( 9.986e-01* 1.487e-02*cos(t)+ -5.306e-02* 2.233e-03*sin(t)), 5.827e-02 +2.000*( 5.306e-02* 1.487e-02*cos(t)+ 9.986e-01* 2.233e-03*sin(t)) not |
| # Age 85, p23 - p12
|
# Age 80, p21 - p12 |
| set label "85" at 1.926e-001, 1.457e-001 center
|
set label "80" at 1.467e-01, 9.186e-02 center |
| replot 1.926e-001+ 2.000*( 9.723e-001* 7.665e-003*cos(t)+-2.338e-001* 3.559e-003*sin(t)), 1.457e-001 +2.000*( 2.338e-001* 7.665e-003*cos(t)+ 9.723e-001* 3.559e-003*sin(t)) not
|
replot 1.467e-01+ 2.000*( 9.949e-01* 8.452e-03*cos(t)+ -1.006e-01* 2.434e-03*sin(t)), 9.186e-02 +2.000*( 1.006e-01* 8.452e-03*cos(t)+ 9.949e-01* 2.434e-03*sin(t)) not |
| # Age 90, p23 - p12
|
# Age 85, p21 - p12 |
| set label "90" at 2.004e-001, 2.292e-001 center
|
set label "85" at 1.230e-01, 1.445e-01 center |
| replot 2.004e-001+ 2.000*( 6.216e-001* 9.497e-003*cos(t)+-7.833e-001* 7.226e-003*sin(t)), 2.292e-001 +2.000*( 7.833e-001* 9.497e-003*cos(t)+ 6.216e-001* 7.226e-003*sin(t)) not
|
replot 1.230e-01+ 2.000*( 9.691e-01* 6.397e-03*cos(t)+ -2.468e-01* 3.619e-03*sin(t)), 1.445e-01 +2.000*( 2.468e-01* 6.397e-03*cos(t)+ 9.691e-01* 3.619e-03*sin(t)) not |
| # Age 95, p23 - p12
|
# Age 90, p21 - p12 |
| set label "95" at 2.085e-001, 3.591e-001 center
|
set label "90" at 1.031e-01, 2.268e-01 center |
| replot 2.085e-001+ 2.000*( 2.331e-001* 1.967e-002*cos(t)+-9.725e-001* 1.154e-002*sin(t)), 3.591e-001 +2.000*( 9.725e-001* 1.967e-002*cos(t)+ 2.331e-001* 1.154e-002*sin(t)) not
|
replot 1.031e-01+ 2.000*( 5.313e-01* 9.166e-03*cos(t)+ -8.472e-01* 6.591e-03*sin(t)), 2.268e-01 +2.000*( 8.472e-01* 9.166e-03*cos(t)+ 5.313e-01* 6.591e-03*sin(t)) not |
| set out "biaspar/varpijgrbiaspar123-12.png";replot;
|
# Age 95, p21 - p12 |
| set parametric;unset label
|
set label "95" at 8.634e-02, 3.543e-01 center |
| set log y;set log x; set xlabel "p21 (year-1)";set ylabel "p13 (year-1)"
|
replot 8.634e-02+ 2.000*( 1.792e-01* 1.914e-02*cos(t)+ -9.838e-01* 8.537e-03*sin(t)), 3.543e-01 +2.000*( 9.838e-01* 1.914e-02*cos(t)+ 1.792e-01* 8.537e-03*sin(t)) not |
| set ter png small
|
set out "biaspar/varpijgrbiaspar121-12.png";replot; |
| set size 0.65,0.65
|
set parametric;unset label |
| set out "biaspar/varpijgrbiaspar121-13.png"
|
set log y;set log x; set xlabel "p23 (year-1)";set ylabel "p12 (year-1)" |
| set label "65" at 1.936e-001, 1.459e-002 center
|
set ter png small size 320, 240 |
| # Age 65, p21 - p13
|
set out "biaspar/varpijgrbiaspar123-12.png" |
| plot [-pi:pi] 1.936e-001+ 2.000*( 1.000e+000* 3.138e-002*cos(t)+-2.829e-003* 1.772e-003*sin(t)), 1.459e-002 +2.000*( 2.829e-003* 3.138e-002*cos(t)+ 1.000e+000* 1.772e-003*sin(t)) not
|
set label "65" at 1.338e-01, 2.337e-02 center |
| # Age 70, p21 - p13
|
# Age 65, p23 - p12 |
| set label "70" at 1.734e-001, 1.976e-002 center
|
plot [-pi:pi] 1.338e-01+ 2.000*( 9.993e-01* 1.768e-02*cos(t)+ -3.696e-02* 1.655e-03*sin(t)), 2.337e-02 +2.000*( 3.696e-02* 1.768e-02*cos(t)+ 9.993e-01* 1.655e-03*sin(t)) not |
| replot 1.734e-001+ 2.000*( 1.000e+000* 2.152e-002*cos(t)+-3.296e-003* 1.771e-003*sin(t)), 1.976e-002 +2.000*( 3.296e-003* 2.152e-002*cos(t)+ 1.000e+000* 1.771e-003*sin(t)) not
|
# Age 70, p23 - p12 |
| # Age 75, p21 - p13
|
set label "70" at 1.460e-01, 3.691e-02 center |
| set label "75" at 1.552e-001, 2.674e-002 center
|
replot 1.460e-01+ 2.000*( 9.985e-01* 1.541e-02*cos(t)+ -5.441e-02* 1.957e-03*sin(t)), 3.691e-02 +2.000*( 5.441e-02* 1.541e-02*cos(t)+ 9.985e-01* 1.957e-03*sin(t)) not |
| replot 1.552e-001+ 2.000*( 1.000e+000* 1.372e-002*cos(t)+-3.023e-003* 1.655e-003*sin(t)), 2.674e-002 +2.000*( 3.023e-003* 1.372e-002*cos(t)+ 1.000e+000* 1.655e-003*sin(t)) not
|
# Age 75, p23 - p12 |
| # Age 80, p21 - p13
|
set label "75" at 1.591e-01, 5.827e-02 center |
| set label "80" at 1.389e-001, 3.614e-002 center
|
replot 1.591e-01+ 2.000*( 9.967e-01* 1.275e-02*cos(t)+ -8.076e-02* 2.137e-03*sin(t)), 5.827e-02 +2.000*( 8.076e-02* 1.275e-02*cos(t)+ 9.967e-01* 2.137e-03*sin(t)) not |
| replot 1.389e-001+ 2.000*( 1.000e+000* 8.242e-003*cos(t)+-4.787e-003* 1.672e-003*sin(t)), 3.614e-002 +2.000*( 4.787e-003* 8.242e-003*cos(t)+ 1.000e+000* 1.672e-003*sin(t)) not
|
# Age 80, p23 - p12 |
| # Age 85, p21 - p13
|
set label "80" at 1.734e-01, 9.186e-02 center |
| set label "85" at 1.243e-001, 4.875e-002 center
|
replot 1.734e-01+ 2.000*( 9.924e-01* 9.884e-03*cos(t)+ -1.231e-01* 2.277e-03*sin(t)), 9.186e-02 +2.000*( 1.231e-01* 9.884e-03*cos(t)+ 9.924e-01* 2.277e-03*sin(t)) not |
| replot 1.243e-001+ 2.000*( 9.991e-001* 6.301e-003*cos(t)+-4.129e-002* 2.647e-003*sin(t)), 4.875e-002 +2.000*( 4.129e-002* 6.301e-003*cos(t)+ 9.991e-001* 2.647e-003*sin(t)) not
|
# Age 85, p23 - p12 |
| # Age 90, p21 - p13
|
set label "85" at 1.889e-01, 1.445e-01 center |
| set label "90" at 1.112e-001, 6.557e-002 center
|
replot 1.889e-01+ 2.000*( 9.741e-01* 7.696e-03*cos(t)+ -2.260e-01* 3.522e-03*sin(t)), 1.445e-01 +2.000*( 2.260e-01* 7.696e-03*cos(t)+ 9.741e-01* 3.522e-03*sin(t)) not |
| replot 1.112e-001+ 2.000*( 9.920e-001* 7.882e-003*cos(t)+-1.263e-001* 5.157e-003*sin(t)), 6.557e-002 +2.000*( 1.263e-001* 7.882e-003*cos(t)+ 9.920e-001* 5.157e-003*sin(t)) not
|
# Age 90, p23 - p12 |
| # Age 95, p21 - p13
|
set label "90" at 2.057e-01, 2.268e-01 center |
| set label "95" at 9.952e-002, 8.778e-002 center
|
replot 2.057e-01+ 2.000*( 6.276e-01* 9.309e-03*cos(t)+ -7.786e-01* 7.132e-03*sin(t)), 2.268e-01 +2.000*( 7.786e-01* 9.309e-03*cos(t)+ 6.276e-01* 7.132e-03*sin(t)) not |
| replot 9.952e-002+ 2.000*( 8.713e-001* 1.061e-002*cos(t)+-4.908e-001* 9.377e-003*sin(t)), 8.778e-002 +2.000*( 4.908e-001* 1.061e-002*cos(t)+ 8.713e-001* 9.377e-003*sin(t)) not
|
# Age 95, p23 - p12 |
| set out "biaspar/varpijgrbiaspar121-13.png";replot;
|
set label "95" at 2.238e-01, 3.543e-01 center |
| set parametric;unset label
|
replot 2.238e-01+ 2.000*( 2.443e-01* 1.925e-02*cos(t)+ -9.697e-01* 1.172e-02*sin(t)), 3.543e-01 +2.000*( 9.697e-01* 1.925e-02*cos(t)+ 2.443e-01* 1.172e-02*sin(t)) not |
| set log y;set log x; set xlabel "p23 (year-1)";set ylabel "p13 (year-1)"
|
set out "biaspar/varpijgrbiaspar123-12.png";replot; |
| set ter png small
|
set parametric;unset label |
| set size 0.65,0.65
|
set log y;set log x; set xlabel "p21 (year-1)";set ylabel "p13 (year-1)" |
| set out "biaspar/varpijgrbiaspar123-13.png"
|
set ter png small size 320, 240 |
| set label "65" at 1.640e-001, 1.459e-002 center
|
set out "biaspar/varpijgrbiaspar121-13.png" |
| # Age 65, p23 - p13
|
set label "65" at 2.481e-01, 1.597e-02 center |
| plot [-pi:pi] 1.640e-001+ 2.000*( 9.992e-001* 2.122e-002*cos(t)+ 3.961e-002* 1.564e-003*sin(t)), 1.459e-002 +2.000*(-3.961e-002* 2.122e-002*cos(t)+ 9.992e-001* 1.564e-003*sin(t)) not
|
# Age 65, p21 - p13 |
| # Age 70, p23 - p13
|
plot [-pi:pi] 2.481e-01+ 2.000*( 1.000e+00* 3.883e-02*cos(t)+ -2.210e-03* 1.821e-03*sin(t)), 1.597e-02 +2.000*( 2.210e-03* 3.883e-02*cos(t)+ 1.000e+00* 1.821e-03*sin(t)) not |
| set label "70" at 1.708e-001, 1.976e-002 center
|
# Age 70, p21 - p13 |
| replot 1.708e-001+ 2.000*( 9.987e-001* 1.755e-002*cos(t)+ 5.131e-002* 1.529e-003*sin(t)), 1.976e-002 +2.000*(-5.131e-002* 1.755e-002*cos(t)+ 9.987e-001* 1.529e-003*sin(t)) not
|
set label "70" at 2.084e-01, 2.117e-02 center |
| # Age 75, p23 - p13
|
replot 2.084e-01+ 2.000*( 1.000e+00* 2.492e-02*cos(t)+ -2.683e-03* 1.773e-03*sin(t)), 2.117e-02 +2.000*( 2.683e-03* 2.492e-02*cos(t)+ 1.000e+00* 1.773e-03*sin(t)) not |
| set label "75" at 1.778e-001, 2.674e-002 center
|
# Age 75, p21 - p13 |
| replot 1.778e-001+ 2.000*( 9.977e-001* 1.374e-002*cos(t)+ 6.811e-002* 1.369e-003*sin(t)), 2.674e-002 +2.000*(-6.811e-002* 1.374e-002*cos(t)+ 9.977e-001* 1.369e-003*sin(t)) not
|
set label "75" at 1.749e-01, 2.803e-02 center |
| # Age 80, p23 - p13
|
replot 1.749e-01+ 2.000*( 1.000e+00* 1.485e-02*cos(t)+ -2.605e-03* 1.622e-03*sin(t)), 2.803e-02 +2.000*( 2.605e-03* 1.485e-02*cos(t)+ 1.000e+00* 1.622e-03*sin(t)) not |
| set label "80" at 1.851e-001, 3.614e-002 center
|
# Age 80, p21 - p13 |
| replot 1.851e-001+ 2.000*( 9.950e-001* 1.009e-002*cos(t)+ 1.004e-001* 1.338e-003*sin(t)), 3.614e-002 +2.000*(-1.004e-001* 1.009e-002*cos(t)+ 9.950e-001* 1.338e-003*sin(t)) not
|
set label "80" at 1.467e-01, 3.707e-02 center |
| # Age 85, p23 - p13
|
replot 1.467e-01+ 2.000*( 1.000e+00* 8.412e-03*cos(t)+ -4.929e-03* 1.650e-03*sin(t)), 3.707e-02 +2.000*( 4.929e-03* 8.412e-03*cos(t)+ 1.000e+00* 1.650e-03*sin(t)) not |
| set label "85" at 1.926e-001, 4.875e-002 center
|
# Age 85, p21 - p13 |
| replot 1.926e-001+ 2.000*( 9.803e-001* 7.637e-003*cos(t)+ 1.977e-001* 2.231e-003*sin(t)), 4.875e-002 +2.000*(-1.977e-001* 7.637e-003*cos(t)+ 9.803e-001* 2.231e-003*sin(t)) not
|
set label "85" at 1.230e-01, 4.893e-02 center |
| # Age 90, p23 - p13
|
replot 1.230e-01+ 2.000*( 9.992e-01* 6.267e-03*cos(t)+ -3.981e-02* 2.636e-03*sin(t)), 4.893e-02 +2.000*( 3.981e-02* 6.267e-03*cos(t)+ 9.992e-01* 2.636e-03*sin(t)) not |
| set label "90" at 2.004e-001, 6.557e-002 center
|
# Age 90, p21 - p13 |
| replot 2.004e-001+ 2.000*( 9.247e-001* 8.660e-003*cos(t)+ 3.806e-001* 4.366e-003*sin(t)), 6.557e-002 +2.000*(-3.806e-001* 8.660e-003*cos(t)+ 9.247e-001* 4.366e-003*sin(t)) not
|
set label "90" at 1.031e-01, 6.441e-02 center |
| # Age 95, p23 - p13
|
replot 1.031e-01+ 2.000*( 9.921e-01* 7.441e-03*cos(t)+ -1.255e-01* 5.013e-03*sin(t)), 6.441e-02 +2.000*( 1.255e-01* 7.441e-03*cos(t)+ 9.921e-01* 5.013e-03*sin(t)) not |
| set label "95" at 2.085e-001, 8.778e-002 center
|
# Age 95, p21 - p13 |
| replot 2.085e-001+ 2.000*( 8.590e-001* 1.328e-002*cos(t)+ 5.121e-001* 8.035e-003*sin(t)), 8.778e-002 +2.000*(-5.121e-001* 1.328e-002*cos(t)+ 8.590e-001* 8.035e-003*sin(t)) not
|
set label "95" at 8.634e-02, 8.441e-02 center |
| set out "biaspar/varpijgrbiaspar123-13.png";replot;
|
replot 8.634e-02+ 2.000*( 6.785e-01* 9.541e-03*cos(t)+ -7.346e-01* 8.651e-03*sin(t)), 8.441e-02 +2.000*( 7.346e-01* 9.541e-03*cos(t)+ 6.785e-01* 8.651e-03*sin(t)) not |
| set parametric;unset label
|
set out "biaspar/varpijgrbiaspar121-13.png";replot; |
| set log y;set log x; set xlabel "p23 (year-1)";set ylabel "p21 (year-1)"
|
set parametric;unset label |
| set ter png small
|
set log y;set log x; set xlabel "p23 (year-1)";set ylabel "p13 (year-1)" |
| set size 0.65,0.65
|
set ter png small size 320, 240 |
| set out "biaspar/varpijgrbiaspar123-21.png"
|
set out "biaspar/varpijgrbiaspar123-13.png" |
| set label "65" at 1.640e-001, 1.936e-001 center
|
set label "65" at 1.338e-01, 1.597e-02 center |
| # Age 65, p23 - p21
|
# Age 65, p23 - p13 |
| plot [-pi:pi] 1.640e-001+ 2.000*( 5.483e-002* 3.141e-002*cos(t)+ 9.985e-001* 2.117e-002*sin(t)), 1.936e-001 +2.000*(-9.985e-001* 3.141e-002*cos(t)+ 5.483e-002* 2.117e-002*sin(t)) not
|
plot [-pi:pi] 1.338e-01+ 2.000*( 9.991e-01* 1.769e-02*cos(t)+ 4.256e-02* 1.662e-03*sin(t)), 1.597e-02 +2.000*( -4.256e-02* 1.769e-02*cos(t)+ 9.991e-01* 1.662e-03*sin(t)) not |
| # Age 70, p23 - p21
|
# Age 70, p23 - p13 |
| set label "70" at 1.708e-001, 1.734e-001 center
|
set label "70" at 1.460e-01, 2.117e-02 center |
| replot 1.708e-001+ 2.000*( 9.437e-002* 2.155e-002*cos(t)+ 9.955e-001* 1.749e-002*sin(t)), 1.734e-001 +2.000*(-9.955e-001* 2.155e-002*cos(t)+ 9.437e-002* 1.749e-002*sin(t)) not
|
replot 1.460e-01+ 2.000*( 9.986e-01* 1.541e-02*cos(t)+ 5.267e-02* 1.580e-03*sin(t)), 2.117e-02 +2.000*( -5.267e-02* 1.541e-02*cos(t)+ 9.986e-01* 1.580e-03*sin(t)) not |
| # Age 75, p23 - p21
|
# Age 75, p23 - p13 |
| set label "75" at 1.778e-001, 1.552e-001 center
|
set label "75" at 1.591e-01, 2.803e-02 center |
| replot 1.778e-001+ 2.000*( 7.002e-001* 1.395e-002*cos(t)+ 7.140e-001* 1.347e-002*sin(t)), 1.552e-001 +2.000*(-7.140e-001* 1.395e-002*cos(t)+ 7.002e-001* 1.347e-002*sin(t)) not
|
replot 1.591e-01+ 2.000*( 9.977e-01* 1.273e-02*cos(t)+ 6.800e-02* 1.375e-03*sin(t)), 2.803e-02 +2.000*( -6.800e-02* 1.273e-02*cos(t)+ 9.977e-01* 1.375e-03*sin(t)) not |
| # Age 80, p23 - p21
|
# Age 80, p23 - p13 |
| set label "80" at 1.851e-001, 1.389e-001 center
|
set label "80" at 1.734e-01, 3.707e-02 center |
| replot 1.851e-001+ 2.000*( 9.962e-001* 1.005e-002*cos(t)+ 8.694e-002* 8.227e-003*sin(t)), 1.389e-001 +2.000*(-8.694e-002* 1.005e-002*cos(t)+ 9.962e-001* 8.227e-003*sin(t)) not
|
replot 1.734e-01+ 2.000*( 9.950e-01* 9.861e-03*cos(t)+ 9.978e-02* 1.332e-03*sin(t)), 3.707e-02 +2.000*( -9.978e-02* 9.861e-03*cos(t)+ 9.950e-01* 1.332e-03*sin(t)) not |
| # Age 85, p23 - p21
|
# Age 85, p23 - p13 |
| set label "85" at 1.926e-001, 1.243e-001 center
|
set label "85" at 1.889e-01, 4.893e-02 center |
| replot 1.926e-001+ 2.000*( 9.869e-001* 7.529e-003*cos(t)+ 1.612e-001* 6.261e-003*sin(t)), 1.243e-001 +2.000*(-1.612e-001* 7.529e-003*cos(t)+ 9.869e-001* 6.261e-003*sin(t)) not
|
replot 1.889e-01+ 2.000*( 9.810e-01* 7.672e-03*cos(t)+ 1.938e-01* 2.231e-03*sin(t)), 4.893e-02 +2.000*( -1.938e-01* 7.672e-03*cos(t)+ 9.810e-01* 2.231e-03*sin(t)) not |
| # Age 90, p23 - p21
|
# Age 90, p23 - p13 |
| set label "90" at 2.004e-001, 1.112e-001 center
|
set label "90" at 2.057e-01, 6.441e-02 center |
| replot 2.004e-001+ 2.000*( 8.505e-001* 8.378e-003*cos(t)+ 5.260e-001* 7.632e-003*sin(t)), 1.112e-001 +2.000*(-5.260e-001* 8.378e-003*cos(t)+ 8.505e-001* 7.632e-003*sin(t)) not
|
replot 2.057e-01+ 2.000*( 9.336e-01* 8.471e-03*cos(t)+ 3.582e-01* 4.337e-03*sin(t)), 6.441e-02 +2.000*( -3.582e-01* 8.471e-03*cos(t)+ 9.336e-01* 4.337e-03*sin(t)) not |
| # Age 95, p23 - p21
|
# Age 95, p23 - p13 |
| set label "95" at 2.085e-001, 9.952e-002 center
|
set label "95" at 2.238e-01, 8.441e-02 center |
| replot 2.085e-001+ 2.000*( 9.743e-001* 1.222e-002*cos(t)+ 2.251e-001* 1.022e-002*sin(t)), 9.952e-002 +2.000*(-2.251e-001* 1.222e-002*cos(t)+ 9.743e-001* 1.022e-002*sin(t)) not
|
replot 2.238e-01+ 2.000*( 9.049e-01* 1.305e-02*cos(t)+ 4.257e-01* 8.023e-03*sin(t)), 8.441e-02 +2.000*( -4.257e-01* 1.305e-02*cos(t)+ 9.049e-01* 8.023e-03*sin(t)) not |
| set out "biaspar/varpijgrbiaspar123-21.png";replot;
|
set out "biaspar/varpijgrbiaspar123-13.png";replot; |
| # Routine varevsij
|
set parametric;unset label |
| set noparametric;set nolabel; set ter png small;set size 0.65, 0.65
|
set log y;set log x; set xlabel "p23 (year-1)";set ylabel "p21 (year-1)" |
| set log y; set nolog x;set xlabel "Age"; set ylabel "Force of mortality (year-1)";
|
set ter png small size 320, 240 |
| plot "biaspar/prmorprev1-stablbased-rbiaspar.txt" u 1:($3) not w l 1
|
set out "biaspar/varpijgrbiaspar123-21.png" |
| replot "biaspar/prmorprev1-stablbased-rbiaspar.txt" u 1:(($3+1.96*$4)) t "95% interval" w l 2
|
set label "65" at 1.338e-01, 2.481e-01 center |
| replot "biaspar/prmorprev1-stablbased-rbiaspar.txt" u 1:(($3-1.96*$4)) not w l 2
|
# Age 65, p23 - p21 |
| set out "biaspar/varmuptjgr-stablbased-biaspar1.png";replot;
|
plot [-pi:pi] 1.338e-01+ 2.000*( 8.647e-03* 3.884e-02*cos(t)+ 1.000e+00* 1.767e-02*sin(t)), 2.481e-01 +2.000*( -1.000e+00* 3.884e-02*cos(t)+ 8.647e-03* 1.767e-02*sin(t)) not |
| |
# Age 70, p23 - p21 |
| |
set label "70" at 1.460e-01, 2.084e-01 center |
| |
replot 1.460e-01+ 2.000*( 1.255e-02* 2.493e-02*cos(t)+ 9.999e-01* 1.539e-02*sin(t)), 2.084e-01 +2.000*( -9.999e-01* 2.493e-02*cos(t)+ 1.255e-02* 1.539e-02*sin(t)) not |
| |
# Age 75, p23 - p21 |
| |
set label "75" at 1.591e-01, 1.749e-01 center |
| |
replot 1.591e-01+ 2.000*( 4.580e-02* 1.486e-02*cos(t)+ 9.990e-01* 1.270e-02*sin(t)), 1.749e-01 +2.000*( -9.990e-01* 1.486e-02*cos(t)+ 4.580e-02* 1.270e-02*sin(t)) not |
| |
# Age 80, p23 - p21 |
| |
set label "80" at 1.734e-01, 1.467e-01 center |
| |
replot 1.734e-01+ 2.000*( 9.957e-01* 9.824e-03*cos(t)+ 9.255e-02* 8.399e-03*sin(t)), 1.467e-01 +2.000*( -9.255e-02* 9.824e-03*cos(t)+ 9.957e-01* 8.399e-03*sin(t)) not |
| |
# Age 85, p23 - p21 |
| |
set label "85" at 1.889e-01, 1.230e-01 center |
| |
replot 1.889e-01+ 2.000*( 9.875e-01* 7.569e-03*cos(t)+ 1.578e-01* 6.226e-03*sin(t)), 1.230e-01 +2.000*( -1.578e-01* 7.569e-03*cos(t)+ 9.875e-01* 6.226e-03*sin(t)) not |
| |
# Age 90, p23 - p21 |
| |
set label "90" at 2.057e-01, 1.031e-01 center |
| |
replot 2.057e-01+ 2.000*( 9.508e-01* 8.134e-03*cos(t)+ 3.099e-01* 7.328e-03*sin(t)), 1.031e-01 +2.000*( -3.099e-01* 8.134e-03*cos(t)+ 9.508e-01* 7.328e-03*sin(t)) not |
| |
# Age 95, p23 - p21 |
| |
set label "95" at 2.238e-01, 8.634e-02 center |
| |
replot 2.238e-01+ 2.000*( 9.979e-01* 1.231e-02*cos(t)+ 6.495e-02* 9.056e-03*sin(t)), 8.634e-02 +2.000*( -6.495e-02* 1.231e-02*cos(t)+ 9.979e-01* 9.056e-03*sin(t)) not |
| |
set out "biaspar/varpijgrbiaspar123-21.png";replot; |
| |
# Routine varevsij |
| |
unset parametric;unset label; set ter png small size 320, 240 |
| |
set log y; unset log x;set xlabel "Age"; set ylabel "Force of mortality (year-1)"; |
| |
plot "biaspar/prmorprev1-stablbased-rbiaspar.txt" u 1:($3) not w l lt 1 |
| |
replot "biaspar/prmorprev1-stablbased-rbiaspar.txt" u 1:(($3+1.96*$4)) t "95% interval" w l lt 2 |
| |
replot "biaspar/prmorprev1-stablbased-rbiaspar.txt" u 1:(($3-1.96*$4)) not w l lt 2 |
| |
set out "biaspar/varmuptjgr-stablbased-biaspar1.png";replot; |
![[BACK]](/icons/small/back.gif){kind=icon}
Return to biaspar.gp CVS log ![[TXT]](/icons/small/text.gif){kind=icon}
![[DIR]](/icons/small/dir.gif){kind=icon}
Up to [Local Repository] / imach / html / doc