From c16aa449b41307bb4cc346ff1817bd742e3263a2 Mon Sep 17 00:00:00 2001 From: =?utf8?q?Agn=C3=A8s=20Li=C3=A8vre?= Date: Thu, 9 Feb 2006 15:20:33 +0000 Subject: [PATCH] Update to 0.98b --- html/doc/biaspar-cov.htm | 60 +- html/doc/biaspar.gp | 490 +- html/doc/biaspar.htm | 102 +- html/doc/biaspar.log | 506 +- html/doc/biaspar/ebiaspar1.png | Bin 3364 -> 3385 bytes html/doc/biaspar/erbiaspar.txt | 71 +- html/doc/biaspar/expbiaspar11.png | Bin 2417 -> 2350 bytes html/doc/biaspar/expbiaspar21.png | Bin 2521 -> 2460 bytes html/doc/biaspar/ilkrbiaspar.txt | 41172 ++++++++-------- html/doc/biaspar/pbiaspar11.png | Bin 6148 -> 6007 bytes html/doc/biaspar/pbiaspar21.png | Bin 4021 -> 4015 bytes html/doc/biaspar/pebiaspar11.png | Bin 2700 -> 2584 bytes html/doc/biaspar/pebiaspar12.png | Bin 2735 -> 2658 bytes html/doc/biaspar/pijrbiaspar.txt | 3291 +- html/doc/biaspar/plrbiaspar.txt | 61 +- html/doc/biaspar/powrbiaspar.txt | 72 +- html/doc/biaspar/prbiaspar.txt | 73 +- .../prmorprev1-stablbased-rbiaspar.txt | 67 +- html/doc/biaspar/probcorrbiaspar.txt | 439 +- html/doc/biaspar/probcovrbiaspar.txt | 439 +- html/doc/biaspar/probrbiaspar.txt | 67 +- html/doc/biaspar/trbiaspar.txt | 69 +- .../varmuptjgr-stablbased-biaspar1.png | Bin 2695 -> 2579 bytes html/doc/biaspar/varpijgrbiaspar113-12.png | Bin 2384 -> 2276 bytes html/doc/biaspar/varpijgrbiaspar121-12.png | Bin 2367 -> 2276 bytes html/doc/biaspar/varpijgrbiaspar121-13.png | Bin 2542 -> 2415 bytes html/doc/biaspar/varpijgrbiaspar123-12.png | Bin 2412 -> 2325 bytes html/doc/biaspar/varpijgrbiaspar123-13.png | Bin 2631 -> 2517 bytes html/doc/biaspar/varpijgrbiaspar123-21.png | Bin 10048 -> 9622 bytes html/doc/biaspar/vbiaspar11.png | Bin 3492 -> 3469 bytes html/doc/biaspar/vbiaspar21.png | Bin 3434 -> 3411 bytes html/doc/biaspar/vplrbiaspar.txt | 71 +- html/doc/biaspar/vrbiaspar.txt | 73 +- html/doc/orbiaspar.txt | 46 +- html/doc/rbiaspar.imach | 82 +- 35 files changed, 23632 insertions(+), 23619 deletions(-) diff --git a/html/doc/biaspar-cov.htm b/html/doc/biaspar-cov.htm index 5a3b8cb..4cc3b15 100644 --- a/html/doc/biaspar-cov.htm +++ b/html/doc/biaspar-cov.htm @@ -1,30 +1,30 @@ - -IMaCh Cov rbiaspar.txt - Imach version 0.98, September 2005, INED-EUROREVES
$Revision$ $Date$
-Title=1st_example
Datafile=data1.txt Firstpass=1 Lastpass=4 Stepm=1 Weight=0 Model=.
- -

Matrix of variance-covariance of pairs of step probabilities

- file biaspar-cov.htm
- -Ellipsoids of confidence centered on point (pij, pkl) are estimatedand drawn. It helps understanding how is the covariance between two incidences. They are expressed in year-1 in order to be less dependent of stepm.
- -
Contour plot corresponding to x'cov-1x = 4 (where x is the column vector (pij,pkl)) are drawn. It can be understood this way: if pij and pkl where uncorrelated the (2x2) matrix of covariance would have been (1/(var pij), 0 , 0, 1/(var pkl)), and the confidence interval would be 2 standard deviations wide on each axis.
Now, if both incidences are correlated (usual case) we diagonalised the inverse of the covariance matrix and made the appropriate rotation to look at the uncorrelated principal directions.
To be simple, these graphs help to understand the significativity of each parameter in relation to a second other one.
- -
Ellipsoids of confidence cov(p13,p12) expressed in year-1 :biaspar/varpijgrbiaspar113-12.png, -
-
Correlation at age 65 (-0.393), 70 (-0.413), 75 (-0.448), 80 (-0.477), 85 (-0.422), 90 (-0.376), 95 (-0.369), -
Ellipsoids of confidence cov(p21,p12) expressed in year-1 :biaspar/varpijgrbiaspar121-12.png, -
-
Correlation at age 65 (0.317), 70 (0.314), 75 (0.307), 80 (0.290), 85 (0.281), 90 (0.298), 95 (0.306), -
Ellipsoids of confidence cov(p23,p12) expressed in year-1 :biaspar/varpijgrbiaspar123-12.png, -
-
Correlation at age 65 (0.380), 70 (0.402), 75 (0.435), 80 (0.456), 85 (0.358), 90 (0.266), 95 (0.256), -
Ellipsoids of confidence cov(p21,p13) expressed in year-1 :biaspar/varpijgrbiaspar121-13.png, -
-
Correlation at age 65 (0.029), 70 (0.021), 75 (0.010), 80 (0.018), 85 (0.079), 90 (0.099), 95 (0.091), -
Ellipsoids of confidence cov(p23,p13) expressed in year-1 :biaspar/varpijgrbiaspar123-13.png, -
-
Correlation at age 65 (-0.473), 70 (-0.506), 75 (-0.561), 80 (-0.595), 85 (-0.519), 90 (-0.461), 95 (-0.418), -
Ellipsoids of confidence cov(p23,p21) expressed in year-1 :biaspar/varpijgrbiaspar123-21.png, -
-
Correlation at age 65 (-0.027), 70 (-0.023), 75 (-0.021), 80 (-0.028), 85 (-0.059), 90 (-0.076), 95 (-0.064), \ No newline at end of file + +IMaCh Cov rbiaspar.txt + Imach version 0.98b, January 2006, INED-EUROREVES
$Revision$ $Date$
+Title=1st_example
Datafile=data1.txt Firstpass=1 Lastpass=4 Stepm=1 Weight=0 Model=.
+ +

Matrix of variance-covariance of pairs of step probabilities

+ file biaspar-cov.htm
+ +Ellipsoids of confidence centered on point (pij, pkl) are estimatedand drawn. It helps understanding how is the covariance between two incidences. They are expressed in year-1 in order to be less dependent of stepm.
+ +
Contour plot corresponding to x'cov-1x = 4 (where x is the column vector (pij,pkl)) are drawn. It can be understood this way: if pij and pkl where uncorrelated the (2x2) matrix of covariance would have been (1/(var pij), 0 , 0, 1/(var pkl)), and the confidence interval would be 2 standard deviations wide on each axis.
Now, if both incidences are correlated (usual case) we diagonalised the inverse of the covariance matrix and made the appropriate rotation to look at the uncorrelated principal directions.
To be simple, these graphs help to understand the significativity of each parameter in relation to a second other one.
+ +
Ellipsoids of confidence cov(p13,p12) expressed in year-1 :biaspar/varpijgrbiaspar113-12.png, +
+
Correlation at age 65 (-0.367), 70 (-0.389), 75 (-0.431), 80 (-0.474), 85 (-0.418), 90 (-0.361), 95 (-0.348), +
Ellipsoids of confidence cov(p21,p12) expressed in year-1 :biaspar/varpijgrbiaspar121-12.png, +
+
Correlation at age 65 (0.332), 70 (0.327), 75 (0.317), 80 (0.292), 85 (0.282), 90 (0.309), 95 (0.322), +
Ellipsoids of confidence cov(p23,p12) expressed in year-1 :biaspar/varpijgrbiaspar123-12.png, +
+
Correlation at age 65 (0.370), 70 (0.393), 75 (0.427), 80 (0.453), 85 (0.359), 90 (0.260), 95 (0.246), +
Ellipsoids of confidence cov(p21,p13) expressed in year-1 :biaspar/varpijgrbiaspar121-13.png, +
+
Correlation at age 65 (0.050), 70 (0.040), 75 (0.025), 80 (0.023), 85 (0.081), 90 (0.109), 95 (0.105), +
Ellipsoids of confidence cov(p23,p13) expressed in year-1 :biaspar/varpijgrbiaspar123-13.png, +
+
Correlation at age 65 (-0.471), 70 (-0.504), 75 (-0.560), 80 (-0.594), 85 (-0.519), 90 (-0.462), 95 (-0.418), +
Ellipsoids of confidence cov(p23,p21) expressed in year-1 :biaspar/varpijgrbiaspar123-21.png, +
+
Correlation at age 65 (-0.044), 70 (-0.040), 75 (-0.035), 80 (-0.035), 85 (-0.059), 90 (-0.083), 95 (-0.079), \ No newline at end of file diff --git a/html/doc/biaspar.gp b/html/doc/biaspar.gp index 8f411e3..fb51c34 100644 --- a/html/doc/biaspar.gp +++ b/html/doc/biaspar.gp @@ -1,245 +1,245 @@ - -# Imach version 0.98, September 2005, INED-EUROREVES -# biaspar.gp -set missing 'NaNq' -cd "/Users/brouard/bin/html/doc" - -set out "biaspar/vbiaspar11.png" - -#set out "vbiaspar11.png" -set xlabel "Age" -set ylabel "Probability" -set ter png small -set size 0.65,0.65 -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 out "biaspar/vbiaspar21.png" - -#set out "vbiaspar21.png" -set xlabel "Age" -set ylabel "Probability" -set ter png small -set size 0.65,0.65 -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 -set out "biaspar/ebiaspar1.png" -set ylabel "Years" -set ter png small -set size 0.65,0.65 -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 out "biaspar/expbiaspar11.png" -set ter png small -set size 0.65,0.65 -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 -set out "biaspar/expbiaspar21.png" -set ter png small -set size 0.65,0.65 -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" -set ylabel "Probability" -set ter png small -set size 0.65,0.65 -unset log y -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 - -set out "biaspar/pbiaspar21.png" -set xlabel "Age" -set ylabel "Probability" -set ter png small -set size 0.65,0.65 -unset log y -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 -p1=-12.245504 -p2=0.092361 -p3=-10.670953 -p4=0.060958 -p5=-2.645034 -p6=-0.022329 -p7=-4.775501 -p8=0.007884 - -set out "biaspar/pebiaspar11.png" - -set title "Probability" - -set ter png small -set size 0.65,0.65 -set log y -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" -set out "biaspar/pebiaspar12.png" - -set ylabel "Quasi-incidence per year" - -set ter png small -set size 0.65,0.65 -set log y -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" -# Routine varprob -set parametric;unset label -set log y;set log x; set xlabel "p13 (year-1)";set ylabel "p12 (year-1)" -set ter png small -set size 0.65,0.65 -set out "biaspar/varpijgrbiaspar113-12.png" -set label "65" at 1.460e-02, 2.328e-02 center -# Age 65, p13 - p12 -plot [-pi:pi] 1.460e-02+ 2.000*( 7.019e-01* 2.099e-03*cos(t)+ 7.123e-01* 1.386e-03*sin(t)), 2.328e-02 +2.000*( -7.123e-01* 2.099e-03*cos(t)+ 7.019e-01* 1.386e-03*sin(t)) not -# Age 70, p13 - p12 -set label "70" at 1.977e-02, 3.688e-02 center -replot 1.977e-02+ 2.000*( 5.418e-01* 2.359e-03*cos(t)+ 8.405e-01* 1.461e-03*sin(t)), 3.688e-02 +2.000*( -8.405e-01* 2.359e-03*cos(t)+ 5.418e-01* 1.461e-03*sin(t)) not -# Age 75, p13 - p12 -set label "75" at 2.675e-02, 5.839e-02 center -replot 2.675e-02+ 2.000*( 4.274e-01* 2.544e-03*cos(t)+ 9.041e-01* 1.383e-03*sin(t)), 5.839e-02 +2.000*( -9.041e-01* 2.544e-03*cos(t)+ 4.274e-01* 1.383e-03*sin(t)) not -# Age 80, p13 - p12 -set label "80" at 3.615e-02, 9.232e-02 center -replot 3.615e-02+ 2.000*( 3.992e-01* 2.744e-03*cos(t)+ 9.169e-01* 1.381e-03*sin(t)), 9.232e-02 +2.000*( -9.169e-01* 2.744e-03*cos(t)+ 3.992e-01* 1.381e-03*sin(t)) not -# Age 85, p13 - p12 -set label "85" at 4.875e-02, 1.457e-01 center -replot 4.875e-02+ 2.000*( 3.981e-01* 4.136e-03*cos(t)+ 9.174e-01* 2.268e-03*sin(t)), 1.457e-01 +2.000*( -9.174e-01* 4.136e-03*cos(t)+ 3.981e-01* 2.268e-03*sin(t)) not -# Age 90, p13 - p12 -set label "90" at 6.556e-02, 2.292e-01 center -replot 6.556e-02+ 2.000*( 2.987e-01* 9.027e-03*cos(t)+ 9.543e-01* 4.655e-03*sin(t)), 2.292e-01 +2.000*( -9.543e-01* 9.027e-03*cos(t)+ 2.987e-01* 4.655e-03*sin(t)) not -# Age 95, p13 - p12 -set label "95" at 8.777e-02, 3.591e-01 center -replot 8.777e-02+ 2.000*( 2.247e-01* 1.984e-02*cos(t)+ 9.744e-01* 8.795e-03*sin(t)), 3.591e-01 +2.000*( -9.744e-01* 1.984e-02*cos(t)+ 2.247e-01* 8.795e-03*sin(t)) not -set out "biaspar/varpijgrbiaspar113-12.png";replot; -set parametric;unset label -set log y;set log x; set xlabel "p21 (year-1)";set ylabel "p12 (year-1)" -set ter png small -set size 0.65,0.65 -set out "biaspar/varpijgrbiaspar121-12.png" -set label "65" at 1.936e-01, 2.328e-02 center -# Age 65, p21 - p12 -plot [-pi:pi] 1.936e-01+ 2.000*( 9.998e-01* 3.109e-02*cos(t)+ -1.822e-02* 1.692e-03*sin(t)), 2.328e-02 +2.000*( 1.822e-02* 3.109e-02*cos(t)+ 9.998e-01* 1.692e-03*sin(t)) not -# Age 70, p21 - p12 -set label "70" at 1.734e-01, 3.688e-02 center -replot 1.734e-01+ 2.000*( 9.995e-01* 2.135e-02*cos(t)+ -3.169e-02* 2.026e-03*sin(t)), 3.688e-02 +2.000*( 3.169e-02* 2.135e-02*cos(t)+ 9.995e-01* 2.026e-03*sin(t)) not -# Age 75, p21 - p12 -set label "75" at 1.552e-01, 5.839e-02 center -replot 1.552e-01+ 2.000*( 9.985e-01* 1.365e-02*cos(t)+ -5.500e-02* 2.256e-03*sin(t)), 5.839e-02 +2.000*( 5.500e-02* 1.365e-02*cos(t)+ 9.985e-01* 2.256e-03*sin(t)) not -# Age 80, p21 - p12 -set label "80" at 1.389e-01, 9.232e-02 center -replot 1.389e-01+ 2.000*( 9.951e-01* 8.259e-03*cos(t)+ -9.930e-02* 2.453e-03*sin(t)), 9.232e-02 +2.000*( 9.930e-02* 8.259e-03*cos(t)+ 9.951e-01* 2.453e-03*sin(t)) not -# Age 85, p21 - p12 -set label "85" at 1.243e-01, 1.457e-01 center -replot 1.243e-01+ 2.000*( 9.672e-01* 6.435e-03*cos(t)+ -2.542e-01* 3.661e-03*sin(t)), 1.457e-01 +2.000*( 2.542e-01* 6.435e-03*cos(t)+ 9.672e-01* 3.661e-03*sin(t)) not -# Age 90, p21 - p12 -set label "90" at 1.112e-01, 2.292e-01 center -replot 1.112e-01+ 2.000*( 5.677e-01* 9.492e-03*cos(t)+ -8.232e-01* 6.841e-03*sin(t)), 2.292e-01 +2.000*( 8.232e-01* 9.492e-03*cos(t)+ 5.677e-01* 6.841e-03*sin(t)) not -# Age 95, p21 - p12 -set label "95" at 9.952e-02, 3.591e-01 center -replot 9.952e-02+ 2.000*( 2.083e-01* 1.976e-02*cos(t)+ -9.781e-01* 9.575e-03*sin(t)), 3.591e-01 +2.000*( 9.781e-01* 1.976e-02*cos(t)+ 2.083e-01* 9.575e-03*sin(t)) not -set out "biaspar/varpijgrbiaspar121-12.png";replot; -set parametric;unset label -set log y;set log x; set xlabel "p23 (year-1)";set ylabel "p12 (year-1)" -set ter png small -set size 0.65,0.65 -set out "biaspar/varpijgrbiaspar123-12.png" -set label "65" at 1.639e-01, 2.328e-02 center -# Age 65, p23 - p12 -plot [-pi:pi] 1.639e-01+ 2.000*( 9.995e-01* 2.122e-02*cos(t)+ -3.216e-02* 1.649e-03*sin(t)), 2.328e-02 +2.000*( 3.216e-02* 2.122e-02*cos(t)+ 9.995e-01* 1.649e-03*sin(t)) not -# Age 70, p23 - p12 -set label "70" at 1.707e-01, 3.688e-02 center -replot 1.707e-01+ 2.000*( 9.988e-01* 1.755e-02*cos(t)+ -4.952e-02* 1.952e-03*sin(t)), 3.688e-02 +2.000*( 4.952e-02* 1.755e-02*cos(t)+ 9.988e-01* 1.952e-03*sin(t)) not -# Age 75, p23 - p12 -set label "75" at 1.777e-01, 5.839e-02 center -replot 1.777e-01+ 2.000*( 9.970e-01* 1.375e-02*cos(t)+ -7.687e-02* 2.132e-03*sin(t)), 5.839e-02 +2.000*( 7.687e-02* 1.375e-02*cos(t)+ 9.970e-01* 2.132e-03*sin(t)) not -# Age 80, p23 - p12 -set label "80" at 1.850e-01, 9.232e-02 center -replot 1.850e-01+ 2.000*( 9.925e-01* 1.011e-02*cos(t)+ -1.224e-01* 2.275e-03*sin(t)), 9.232e-02 +2.000*( 1.224e-01* 1.011e-02*cos(t)+ 9.925e-01* 2.275e-03*sin(t)) not -# Age 85, p23 - p12 -set label "85" at 1.926e-01, 1.457e-01 center -replot 1.926e-01+ 2.000*( 9.723e-01* 7.666e-03*cos(t)+ -2.339e-01* 3.562e-03*sin(t)), 1.457e-01 +2.000*( 2.339e-01* 7.666e-03*cos(t)+ 9.723e-01* 3.562e-03*sin(t)) not -# Age 90, p23 - p12 -set label "90" at 2.004e-01, 2.292e-01 center -replot 2.004e-01+ 2.000*( 6.167e-01* 9.540e-03*cos(t)+ -7.872e-01* 7.206e-03*sin(t)), 2.292e-01 +2.000*( 7.872e-01* 9.540e-03*cos(t)+ 6.167e-01* 7.206e-03*sin(t)) not -# Age 95, p23 - p12 -set label "95" at 2.085e-01, 3.591e-01 center -replot 2.085e-01+ 2.000*( 2.384e-01* 1.981e-02*cos(t)+ -9.712e-01* 1.149e-02*sin(t)), 3.591e-01 +2.000*( 9.712e-01* 1.981e-02*cos(t)+ 2.384e-01* 1.149e-02*sin(t)) not -set out "biaspar/varpijgrbiaspar123-12.png";replot; -set parametric;unset label -set log y;set log x; set xlabel "p21 (year-1)";set ylabel "p13 (year-1)" -set ter png small -set size 0.65,0.65 -set out "biaspar/varpijgrbiaspar121-13.png" -set label "65" at 1.936e-01, 1.460e-02 center -# Age 65, p21 - p13 -plot [-pi:pi] 1.936e-01+ 2.000*( 1.000e+00* 3.109e-02*cos(t)+ -1.676e-03* 1.773e-03*sin(t)), 1.460e-02 +2.000*( 1.676e-03* 3.109e-02*cos(t)+ 1.000e+00* 1.773e-03*sin(t)) not -# Age 70, p21 - p13 -set label "70" at 1.734e-01, 1.977e-02 center -replot 1.734e-01+ 2.000*( 1.000e+00* 2.134e-02*cos(t)+ -1.717e-03* 1.772e-03*sin(t)), 1.977e-02 +2.000*( 1.717e-03* 2.134e-02*cos(t)+ 1.000e+00* 1.772e-03*sin(t)) not -# Age 75, p21 - p13 -set label "75" at 1.552e-01, 2.675e-02 center -replot 1.552e-01+ 2.000*( 1.000e+00* 1.363e-02*cos(t)+ -1.187e-03* 1.657e-03*sin(t)), 2.675e-02 +2.000*( 1.187e-03* 1.363e-02*cos(t)+ 1.000e+00* 1.657e-03*sin(t)) not -# Age 80, p21 - p13 -set label "80" at 1.389e-01, 3.615e-02 center -replot 1.389e-01+ 2.000*( 1.000e+00* 8.222e-03*cos(t)+ -3.844e-03* 1.674e-03*sin(t)), 3.615e-02 +2.000*( 3.844e-03* 8.222e-03*cos(t)+ 1.000e+00* 1.674e-03*sin(t)) not -# Age 85, p21 - p13 -set label "85" at 1.243e-01, 4.875e-02 center -replot 1.243e-01+ 2.000*( 9.992e-01* 6.297e-03*cos(t)+ -4.035e-02* 2.643e-03*sin(t)), 4.875e-02 +2.000*( 4.035e-02* 6.297e-03*cos(t)+ 9.992e-01* 2.643e-03*sin(t)) not -# Age 90, p21 - p13 -set label "90" at 1.112e-01, 6.556e-02 center -replot 1.112e-01+ 2.000*( 9.932e-01* 7.824e-03*cos(t)+ -1.163e-01* 5.151e-03*sin(t)), 6.556e-02 +2.000*( 1.163e-01* 7.824e-03*cos(t)+ 9.932e-01* 5.151e-03*sin(t)) not -# Age 95, p21 - p13 -set label "95" at 9.952e-02, 8.777e-02 center -replot 9.952e-02+ 2.000*( 8.754e-01* 1.047e-02*cos(t)+ -4.835e-01* 9.400e-03*sin(t)), 8.777e-02 +2.000*( 4.835e-01* 1.047e-02*cos(t)+ 8.754e-01* 9.400e-03*sin(t)) not -set out "biaspar/varpijgrbiaspar121-13.png";replot; -set parametric;unset label -set log y;set log x; set xlabel "p23 (year-1)";set ylabel "p13 (year-1)" -set ter png small -set size 0.65,0.65 -set out "biaspar/varpijgrbiaspar123-13.png" -set label "65" at 1.639e-01, 1.460e-02 center -# Age 65, p23 - p13 -plot [-pi:pi] 1.639e-01+ 2.000*( 9.992e-01* 2.122e-02*cos(t)+ 3.972e-02* 1.562e-03*sin(t)), 1.460e-02 +2.000*( -3.972e-02* 2.122e-02*cos(t)+ 9.992e-01* 1.562e-03*sin(t)) not -# Age 70, p23 - p13 -set label "70" at 1.707e-01, 1.977e-02 center -replot 1.707e-01+ 2.000*( 9.987e-01* 1.756e-02*cos(t)+ 5.147e-02* 1.527e-03*sin(t)), 1.977e-02 +2.000*( -5.147e-02* 1.756e-02*cos(t)+ 9.987e-01* 1.527e-03*sin(t)) not -# Age 75, p23 - p13 -set label "75" at 1.777e-01, 2.675e-02 center -replot 1.777e-01+ 2.000*( 9.977e-01* 1.375e-02*cos(t)+ 6.833e-02* 1.368e-03*sin(t)), 2.675e-02 +2.000*( -6.833e-02* 1.375e-02*cos(t)+ 9.977e-01* 1.368e-03*sin(t)) not -# Age 80, p23 - p13 -set label "80" at 1.850e-01, 3.615e-02 center -replot 1.850e-01+ 2.000*( 9.949e-01* 1.009e-02*cos(t)+ 1.005e-01* 1.339e-03*sin(t)), 3.615e-02 +2.000*( -1.005e-01* 1.009e-02*cos(t)+ 9.949e-01* 1.339e-03*sin(t)) not -# Age 85, p23 - p13 -set label "85" at 1.926e-01, 4.875e-02 center -replot 1.926e-01+ 2.000*( 9.804e-01* 7.637e-03*cos(t)+ 1.972e-01* 2.228e-03*sin(t)), 4.875e-02 +2.000*( -1.972e-01* 7.637e-03*cos(t)+ 9.804e-01* 2.228e-03*sin(t)) not -# Age 90, p23 - p13 -set label "90" at 2.004e-01, 6.556e-02 center -replot 2.004e-01+ 2.000*( 9.254e-01* 8.649e-03*cos(t)+ 3.789e-01* 4.358e-03*sin(t)), 6.556e-02 +2.000*( -3.789e-01* 8.649e-03*cos(t)+ 9.254e-01* 4.358e-03*sin(t)) not -# Age 95, p23 - p13 -set label "95" at 2.085e-01, 8.777e-02 center -replot 2.085e-01+ 2.000*( 8.603e-01* 1.326e-02*cos(t)+ 5.098e-01* 8.022e-03*sin(t)), 8.777e-02 +2.000*( -5.098e-01* 1.326e-02*cos(t)+ 8.603e-01* 8.022e-03*sin(t)) not -set out "biaspar/varpijgrbiaspar123-13.png";replot; -set parametric;unset label -set log y;set log x; set xlabel "p23 (year-1)";set ylabel "p21 (year-1)" -set ter png small -set size 0.65,0.65 -set out "biaspar/varpijgrbiaspar123-21.png" -set label "65" at 1.639e-01, 1.936e-01 center -# Age 65, p23 - p21 -plot [-pi:pi] 1.639e-01+ 2.000*( 3.376e-02* 3.110e-02*cos(t)+ 9.994e-01* 2.119e-02*sin(t)), 1.936e-01 +2.000*( -9.994e-01* 3.110e-02*cos(t)+ 3.376e-02* 2.119e-02*sin(t)) not -# Age 70, p23 - p21 -set label "70" at 1.707e-01, 1.734e-01 center -replot 1.707e-01+ 2.000*( 5.801e-02* 2.135e-02*cos(t)+ 9.983e-01* 1.752e-02*sin(t)), 1.734e-01 +2.000*( -9.983e-01* 2.135e-02*cos(t)+ 5.801e-02* 1.752e-02*sin(t)) not -# Age 75, p23 - p21 -set label "75" at 1.777e-01, 1.552e-01 center -replot 1.777e-01+ 2.000*( 8.046e-01* 1.382e-02*cos(t)+ 5.938e-01* 1.352e-02*sin(t)), 1.552e-01 +2.000*( -5.938e-01* 1.382e-02*cos(t)+ 8.046e-01* 1.352e-02*sin(t)) not -# Age 80, p23 - p21 -set label "80" at 1.850e-01, 1.389e-01 center -replot 1.850e-01+ 2.000*( 9.977e-01* 1.005e-02*cos(t)+ 6.821e-02* 8.213e-03*sin(t)), 1.389e-01 +2.000*( -6.821e-02* 1.005e-02*cos(t)+ 9.977e-01* 8.213e-03*sin(t)) not -# Age 85, p23 - p21 -set label "85" at 1.926e-01, 1.243e-01 center -replot 1.926e-01+ 2.000*( 9.869e-01* 7.531e-03*cos(t)+ 1.615e-01* 6.256e-03*sin(t)), 1.243e-01 +2.000*( -1.615e-01* 7.531e-03*cos(t)+ 9.869e-01* 6.256e-03*sin(t)) not -# Age 90, p23 - p21 -set label "90" at 2.004e-01, 1.112e-01 center -replot 2.004e-01+ 2.000*( 8.742e-01* 8.336e-03*cos(t)+ 4.856e-01* 7.619e-03*sin(t)), 1.112e-01 +2.000*( -4.856e-01* 8.336e-03*cos(t)+ 8.742e-01* 7.619e-03*sin(t)) not -# Age 95, p23 - p21 -set label "95" at 2.085e-01, 9.952e-02 center -replot 2.085e-01+ 2.000*( 9.838e-01* 1.218e-02*cos(t)+ 1.795e-01* 1.016e-02*sin(t)), 9.952e-02 +2.000*( -1.795e-01* 1.218e-02*cos(t)+ 9.838e-01* 1.016e-02*sin(t)) not -set out "biaspar/varpijgrbiaspar123-21.png";replot; -# Routine varevsij -set noparametric;set nolabel; set ter png small;set size 0.65, 0.65 - set log y; set nolog x;set xlabel "Age"; set ylabel "Force of mortality (year-1)"; - plot "biaspar/prmorprev1-stablbased-rbiaspar.txt" u 1:($3) not w l 1 - replot "biaspar/prmorprev1-stablbased-rbiaspar.txt" u 1:(($3+1.96*$4)) t "95% interval" w l 2 - replot "biaspar/prmorprev1-stablbased-rbiaspar.txt" u 1:(($3-1.96*$4)) not w l 2 -set out "biaspar/varmuptjgr-stablbased-biaspar1.png";replot; + +# Imach version 0.98b, January 2006, INED-EUROREVES +# biaspar.gp +set missing 'NaNq' +cd "C:/Program Files/IMaCh/doc/" + +set out "biaspar/vbiaspar11.png" + +#set out "vbiaspar11.png" +set xlabel "Age" +set ylabel "Probability" +set ter png small +set size 0.65,0.65 +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 out "biaspar/vbiaspar21.png" + +#set out "vbiaspar21.png" +set xlabel "Age" +set ylabel "Probability" +set ter png small +set size 0.65,0.65 +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 +set out "biaspar/ebiaspar1.png" +set ylabel "Years" +set ter png small +set size 0.65,0.65 +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 out "biaspar/expbiaspar11.png" +set ter png small +set size 0.65,0.65 +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 +set out "biaspar/expbiaspar21.png" +set ter png small +set size 0.65,0.65 +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" +set ylabel "Probability" +set ter png small +set size 0.65,0.65 +unset log y +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 + +set out "biaspar/pbiaspar21.png" +set xlabel "Age" +set ylabel "Probability" +set ter png small +set size 0.65,0.65 +unset log y +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 +p1=-12.245240 +p2=0.092358 +p3=-10.671890 +p4=0.060969 +p5=-2.645345 +p6=-0.022325 +p7=-4.773317 +p8=0.007859 + +set out "biaspar/pebiaspar11.png" + +set title "Probability" + +set ter png small +set size 0.65,0.65 +set log y +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" +set out "biaspar/pebiaspar12.png" + +set ylabel "Quasi-incidence per year" + +set ter png small +set size 0.65,0.65 +set log y +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" +# Routine varprob +set parametric;unset label +set log y;set log x; set xlabel "p13 (year-1)";set ylabel "p12 (year-1)" +set ter png small +set size 0.65,0.65 +set out "biaspar/varpijgrbiaspar113-12.png" +set label "65" at 1.459e-002, 2.328e-002 center +# Age 65, p13 - p12 +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 +# Age 70, p13 - p12 +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 +# Age 75, p13 - p12 +set label "75" at 2.674e-002, 5.839e-002 center +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 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 +# Age 85, p13 - p12 +set label "85" at 4.875e-002, 1.457e-001 center +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 +# Age 90, p13 - p12 +set label "90" at 6.557e-002, 2.292e-001 center +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 +# Age 95, p13 - p12 +set label "95" at 8.778e-002, 3.591e-001 center +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 parametric;unset label +set log y;set log x; set xlabel "p21 (year-1)";set ylabel "p12 (year-1)" +set ter png small +set size 0.65,0.65 +set out "biaspar/varpijgrbiaspar121-12.png" +set label "65" at 1.936e-001, 2.328e-002 center +# Age 65, p21 - p12 +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 +# Age 70, p21 - p12 +set label "70" at 1.734e-001, 3.688e-002 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 75, p21 - p12 +set label "75" at 1.552e-001, 5.839e-002 center +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 +# Age 80, p21 - p12 +set label "80" at 1.389e-001, 9.233e-002 center +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 +# Age 85, p21 - p12 +set label "85" at 1.243e-001, 1.457e-001 center +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 +# Age 90, p21 - p12 +set label "90" at 1.112e-001, 2.292e-001 center +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 +# Age 95, p21 - p12 +set label "95" at 9.952e-002, 3.591e-001 center +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 out "biaspar/varpijgrbiaspar121-12.png";replot; +set parametric;unset label +set log y;set log x; set xlabel "p23 (year-1)";set ylabel "p12 (year-1)" +set ter png small +set size 0.65,0.65 +set out "biaspar/varpijgrbiaspar123-12.png" +set label "65" at 1.640e-001, 2.328e-002 center +# Age 65, p23 - p12 +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 +# Age 70, p23 - p12 +set label "70" at 1.708e-001, 3.688e-002 center +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 +# Age 75, p23 - p12 +set label "75" at 1.778e-001, 5.839e-002 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 +# Age 80, p23 - p12 +set label "80" at 1.851e-001, 9.233e-002 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 +# Age 85, p23 - p12 +set label "85" at 1.926e-001, 1.457e-001 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 +# Age 90, p23 - p12 +set label "90" at 2.004e-001, 2.292e-001 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 +# Age 95, p23 - p12 +set label "95" at 2.085e-001, 3.591e-001 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 +set out "biaspar/varpijgrbiaspar123-12.png";replot; +set parametric;unset label +set log y;set log x; set xlabel "p21 (year-1)";set ylabel "p13 (year-1)" +set ter png small +set size 0.65,0.65 +set out "biaspar/varpijgrbiaspar121-13.png" +set label "65" at 1.936e-001, 1.459e-002 center +# Age 65, p21 - p13 +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 +# Age 70, p21 - p13 +set label "70" at 1.734e-001, 1.976e-002 center +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 75, p21 - p13 +set label "75" at 1.552e-001, 2.674e-002 center +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 80, p21 - p13 +set label "80" at 1.389e-001, 3.614e-002 center +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 85, p21 - p13 +set label "85" at 1.243e-001, 4.875e-002 center +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 90, p21 - p13 +set label "90" at 1.112e-001, 6.557e-002 center +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 95, p21 - p13 +set label "95" at 9.952e-002, 8.778e-002 center +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 +set out "biaspar/varpijgrbiaspar121-13.png";replot; +set parametric;unset label +set log y;set log x; set xlabel "p23 (year-1)";set ylabel "p13 (year-1)" +set ter png small +set size 0.65,0.65 +set out "biaspar/varpijgrbiaspar123-13.png" +set label "65" at 1.640e-001, 1.459e-002 center +# Age 65, p23 - p13 +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 70, p23 - p13 +set label "70" at 1.708e-001, 1.976e-002 center +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 +# Age 75, p23 - p13 +set label "75" at 1.778e-001, 2.674e-002 center +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 +# Age 80, p23 - p13 +set label "80" at 1.851e-001, 3.614e-002 center +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 +# Age 85, p23 - p13 +set label "85" at 1.926e-001, 4.875e-002 center +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 +# Age 90, p23 - p13 +set label "90" at 2.004e-001, 6.557e-002 center +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 +# Age 95, p23 - p13 +set label "95" at 2.085e-001, 8.778e-002 center +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 out "biaspar/varpijgrbiaspar123-13.png";replot; +set parametric;unset label +set log y;set log x; set xlabel "p23 (year-1)";set ylabel "p21 (year-1)" +set ter png small +set size 0.65,0.65 +set out "biaspar/varpijgrbiaspar123-21.png" +set label "65" at 1.640e-001, 1.936e-001 center +# Age 65, p23 - p21 +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 +# Age 70, p23 - p21 +set label "70" at 1.708e-001, 1.734e-001 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 +# Age 75, p23 - p21 +set label "75" at 1.778e-001, 1.552e-001 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 +# Age 80, p23 - p21 +set label "80" at 1.851e-001, 1.389e-001 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 +# Age 85, p23 - p21 +set label "85" at 1.926e-001, 1.243e-001 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 +# Age 90, p23 - p21 +set label "90" at 2.004e-001, 1.112e-001 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 +# Age 95, p23 - p21 +set label "95" at 2.085e-001, 9.952e-002 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 +set out "biaspar/varpijgrbiaspar123-21.png";replot; +# Routine varevsij +set noparametric;set nolabel; set ter png small;set size 0.65, 0.65 + set log y; set nolog x;set xlabel "Age"; set ylabel "Force of mortality (year-1)"; + plot "biaspar/prmorprev1-stablbased-rbiaspar.txt" u 1:($3) not w l 1 + replot "biaspar/prmorprev1-stablbased-rbiaspar.txt" u 1:(($3+1.96*$4)) t "95% interval" w l 2 + replot "biaspar/prmorprev1-stablbased-rbiaspar.txt" u 1:(($3-1.96*$4)) not w l 2 +set out "biaspar/varmuptjgr-stablbased-biaspar1.png";replot; diff --git a/html/doc/biaspar.htm b/html/doc/biaspar.htm index f0f9984..bc9a53e 100644 --- a/html/doc/biaspar.htm +++ b/html/doc/biaspar.htm @@ -1,51 +1,53 @@ - -IMaCh rbiaspar.txt - Imach version 0.98, September 2005, INED-EUROREVES
$Revision$ $Date$
-Title=1st_example
Datafile=data1.txt Firstpass=1 Lastpass=4 Stepm=1 Weight=0 Model=.
- -
- -
Total number of observations=8270
-Youngest age at first (selected) pass 70.00, oldest age 104.17
-Interval (in months) between two waves: Min=1 Max=74 Mean=24.04
- -
File of contributions to the likelihood: biaspar/ilkrbiaspar.txt
-