From 5eb91f05d127e7f493ec827b8179e6b30999635b Mon Sep 17 00:00:00 2001 From: "N. Brouard" Date: Tue, 25 Oct 2005 11:53:26 +0000 Subject: [PATCH] Version 0.98 MacOS/X Panther --- html/doc/biaspar-cov.htm | 60 ++--- html/doc/biaspar.gp | 490 ++++++++++++++++++------------------- html/doc/biaspar.htm | 100 ++++---- html/doc/biaspar.log | 506 +++++++++++++++++++-------------------- html/doc/orbiaspar.txt | 46 ++-- html/doc/rbiaspar.imach | 82 +++---- 6 files changed, 642 insertions(+), 642 deletions(-) diff --git a/html/doc/biaspar-cov.htm b/html/doc/biaspar-cov.htm index 2ccf884..5a3b8cb 100644 --- a/html/doc/biaspar-cov.htm +++ b/html/doc/biaspar-cov.htm @@ -1,30 +1,30 @@ - -IMaCh Cov rbiaspar.txt - Imach version 0.07a, May 2004, 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.365), 70 (-0.388), 75 (-0.431), 80 (-0.475), 85 (-0.419), 90 (-0.360), 95 (-0.346), -
Ellipsoids of confidence cov(p21,p12) expressed in year-1 :biaspar/varpijgrbiaspar121-12.png, -
-
Correlation at age 65 (0.346), 70 (0.341), 75 (0.328), 80 (0.297), 85 (0.283), 90 (0.316), 95 (0.332), -
Ellipsoids of confidence cov(p23,p12) expressed in year-1 :biaspar/varpijgrbiaspar123-12.png, -
-
Correlation at age 65 (0.367), 70 (0.390), 75 (0.426), 80 (0.453), 85 (0.359), 90 (0.257), 95 (0.241), -
Ellipsoids of confidence cov(p21,p13) expressed in year-1 :biaspar/varpijgrbiaspar121-13.png, -
-
Correlation at age 65 (0.034), 70 (0.024), 75 (0.012), 80 (0.017), 85 (0.080), 90 (0.103), 95 (0.096), -
Ellipsoids of confidence cov(p23,p13) expressed in year-1 :biaspar/varpijgrbiaspar123-13.png, -
-
Correlation at age 65 (-0.456), 70 (-0.492), 75 (-0.552), 80 (-0.594), 85 (-0.523), 90 (-0.456), 95 (-0.404), -
Ellipsoids of confidence cov(p23,p21) expressed in year-1 :biaspar/varpijgrbiaspar123-21.png, -
-
Correlation at age 65 (-0.025), 70 (-0.022), 75 (-0.020), 80 (-0.027), 85 (-0.060), 90 (-0.077), 95 (-0.065), \ No newline at end of file + +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 diff --git a/html/doc/biaspar.gp b/html/doc/biaspar.gp index 1b37560..8f411e3 100644 --- a/html/doc/biaspar.gp +++ b/html/doc/biaspar.gp @@ -1,245 +1,245 @@ - -# Imach version 0.07a, May 2004, INED-EUROREVES -# biaspar.gp -set missing 'NaNq' -cd "D:/imachcvs/imach/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.245160 -p2=0.092357 -p3=-10.672078 -p4=0.060971 -p5=-2.645815 -p6=-0.022320 -p7=-4.773208 -p8=0.007857 - -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*( 6.866e-001* 2.057e-003*cos(t)+ 7.270e-001* 1.402e-003*sin(t)), 2.328e-002 +2.000*(-7.270e-001* 2.057e-003*cos(t)+ 6.866e-001* 1.402e-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.191e-001* 2.326e-003*cos(t)+ 8.547e-001* 1.471e-003*sin(t)), 3.688e-002 +2.000*(-8.547e-001* 2.326e-003*cos(t)+ 5.191e-001* 1.471e-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.119e-001* 2.525e-003*cos(t)+ 9.112e-001* 1.387e-003*sin(t)), 5.839e-002 +2.000*(-9.112e-001* 2.525e-003*cos(t)+ 4.119e-001* 1.387e-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.972e-001* 2.741e-003*cos(t)+ 9.177e-001* 1.382e-003*sin(t)), 9.233e-002 +2.000*(-9.177e-001* 2.741e-003*cos(t)+ 3.972e-001* 1.382e-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.964e-001* 4.133e-003*cos(t)+ 9.181e-001* 2.272e-003*sin(t)), 1.457e-001 +2.000*(-9.181e-001* 4.133e-003*cos(t)+ 3.964e-001* 2.272e-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.870e-001* 8.981e-003*cos(t)+ 9.579e-001* 4.671e-003*sin(t)), 2.292e-001 +2.000*(-9.579e-001* 8.981e-003*cos(t)+ 2.870e-001* 4.671e-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.106e-001* 1.972e-002*cos(t)+ 9.776e-001* 8.813e-003*sin(t)), 3.591e-001 +2.000*(-9.776e-001* 1.972e-002*cos(t)+ 2.106e-001* 8.813e-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.151e-002*cos(t)+-1.959e-002* 1.668e-003*sin(t)), 2.328e-002 +2.000*( 1.959e-002* 3.151e-002*cos(t)+ 9.998e-001* 1.668e-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.994e-001* 2.161e-002*cos(t)+-3.389e-002* 2.001e-003*sin(t)), 3.688e-002 +2.000*( 3.389e-002* 2.161e-002*cos(t)+ 9.994e-001* 2.001e-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.983e-001* 1.378e-002*cos(t)+-5.800e-002* 2.235e-003*sin(t)), 5.839e-002 +2.000*( 5.800e-002* 1.378e-002*cos(t)+ 9.983e-001* 2.235e-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.949e-001* 8.292e-003*cos(t)+-1.011e-001* 2.447e-003*sin(t)), 9.233e-002 +2.000*( 1.011e-001* 8.292e-003*cos(t)+ 9.949e-001* 2.447e-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.669e-001* 6.441e-003*cos(t)+-2.553e-001* 3.657e-003*sin(t)), 1.457e-001 +2.000*( 2.553e-001* 6.441e-003*cos(t)+ 9.669e-001* 3.657e-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.879e-001* 9.566e-003*cos(t)+-8.090e-001* 6.788e-003*sin(t)), 2.292e-001 +2.000*( 8.090e-001* 9.566e-003*cos(t)+ 5.879e-001* 6.788e-003*sin(t)) not -# Age 95, p21 - p12 -set label "95" at 9.953e-002, 3.591e-001 center -replot 9.953e-002+ 2.000*( 2.290e-001* 1.977e-002*cos(t)+-9.734e-001* 9.572e-003*sin(t)), 3.591e-001 +2.000*( 9.734e-001* 1.977e-002*cos(t)+ 2.290e-001* 9.572e-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.107e-002*cos(t)+-3.113e-002* 1.654e-003*sin(t)), 2.328e-002 +2.000*( 3.113e-002* 2.107e-002*cos(t)+ 9.995e-001* 1.654e-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.745e-002*cos(t)+-4.823e-002* 1.959e-003*sin(t)), 3.688e-002 +2.000*( 4.823e-002* 1.745e-002*cos(t)+ 9.988e-001* 1.959e-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.971e-001* 1.368e-002*cos(t)+-7.559e-002* 2.139e-003*sin(t)), 5.839e-002 +2.000*( 7.559e-002* 1.368e-002*cos(t)+ 9.971e-001* 2.139e-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.925e-001* 1.009e-002*cos(t)+-1.220e-001* 2.279e-003*sin(t)), 9.233e-002 +2.000*( 1.220e-001* 1.009e-002*cos(t)+ 9.925e-001* 2.279e-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.722e-001* 7.671e-003*cos(t)+-2.342e-001* 3.560e-003*sin(t)), 1.457e-001 +2.000*( 2.342e-001* 7.671e-003*cos(t)+ 9.722e-001* 3.560e-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.128e-001* 9.486e-003*cos(t)+-7.902e-001* 7.229e-003*sin(t)), 2.292e-001 +2.000*( 7.902e-001* 9.486e-003*cos(t)+ 6.128e-001* 7.229e-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.248e-001* 1.970e-002*cos(t)+-9.744e-001* 1.148e-002*sin(t)), 3.591e-001 +2.000*( 9.744e-001* 1.970e-002*cos(t)+ 2.248e-001* 1.148e-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.150e-002*cos(t)+-1.880e-003* 1.741e-003*sin(t)), 1.459e-002 +2.000*( 1.880e-003* 3.150e-002*cos(t)+ 1.000e+000* 1.741e-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.160e-002*cos(t)+-1.974e-003* 1.743e-003*sin(t)), 1.976e-002 +2.000*( 1.974e-003* 2.160e-002*cos(t)+ 1.000e+000* 1.743e-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.376e-002*cos(t)+-1.404e-003* 1.636e-003*sin(t)), 2.674e-002 +2.000*( 1.404e-003* 1.376e-002*cos(t)+ 1.000e+000* 1.636e-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.253e-003*cos(t)+-3.684e-003* 1.671e-003*sin(t)), 3.614e-002 +2.000*( 3.684e-003* 8.253e-003*cos(t)+ 1.000e+000* 1.671e-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.992e-001* 6.301e-003*cos(t)+-4.077e-002* 2.642e-003*sin(t)), 4.875e-002 +2.000*( 4.077e-002* 6.301e-003*cos(t)+ 9.992e-001* 2.642e-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.932e-001* 7.891e-003*cos(t)+-1.164e-001* 5.116e-003*sin(t)), 6.557e-002 +2.000*( 1.164e-001* 7.891e-003*cos(t)+ 9.932e-001* 5.116e-003*sin(t)) not -# Age 95, p21 - p13 -set label "95" at 9.953e-002, 8.778e-002 center -replot 9.953e-002+ 2.000*( 9.056e-001* 1.057e-002*cos(t)+-4.241e-001* 9.328e-003*sin(t)), 8.778e-002 +2.000*( 4.241e-001* 1.057e-002*cos(t)+ 9.056e-001* 9.328e-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.993e-001* 2.108e-002*cos(t)+ 3.790e-002* 1.549e-003*sin(t)), 1.459e-002 +2.000*(-3.790e-002* 2.108e-002*cos(t)+ 9.993e-001* 1.549e-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.988e-001* 1.745e-002*cos(t)+ 4.951e-002* 1.516e-003*sin(t)), 1.976e-002 +2.000*(-4.951e-002* 1.745e-002*cos(t)+ 9.988e-001* 1.516e-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.978e-001* 1.368e-002*cos(t)+ 6.667e-002* 1.362e-003*sin(t)), 2.674e-002 +2.000*(-6.667e-002* 1.368e-002*cos(t)+ 9.978e-001* 1.362e-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.006e-002*cos(t)+ 1.004e-001* 1.338e-003*sin(t)), 3.614e-002 +2.000*(-1.004e-001* 1.006e-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.801e-001* 7.644e-003*cos(t)+ 1.985e-001* 2.219e-003*sin(t)), 4.875e-002 +2.000*(-1.985e-001* 7.644e-003*cos(t)+ 9.801e-001* 2.219e-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.273e-001* 8.612e-003*cos(t)+ 3.743e-001* 4.351e-003*sin(t)), 6.557e-002 +2.000*(-3.743e-001* 8.612e-003*cos(t)+ 9.273e-001* 4.351e-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.652e-001* 1.310e-002*cos(t)+ 5.014e-001* 8.033e-003*sin(t)), 8.778e-002 +2.000*(-5.014e-001* 1.310e-002*cos(t)+ 8.652e-001* 8.033e-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*( 3.032e-002* 3.151e-002*cos(t)+ 9.995e-001* 2.105e-002*sin(t)), 1.936e-001 +2.000*(-9.995e-001* 3.151e-002*cos(t)+ 3.032e-002* 2.105e-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*( 4.964e-002* 2.161e-002*cos(t)+ 9.988e-001* 1.741e-002*sin(t)), 1.734e-001 +2.000*(-9.988e-001* 2.161e-002*cos(t)+ 4.964e-002* 1.741e-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*( 5.539e-001* 1.385e-002*cos(t)+ 8.326e-001* 1.356e-002*sin(t)), 1.552e-001 +2.000*(-8.326e-001* 1.385e-002*cos(t)+ 5.539e-001* 1.356e-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.977e-001* 1.002e-002*cos(t)+ 6.827e-002* 8.244e-003*sin(t)), 1.389e-001 +2.000*(-6.827e-002* 1.002e-002*cos(t)+ 9.977e-001* 8.244e-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.865e-001* 7.536e-003*cos(t)+ 1.638e-001* 6.259e-003*sin(t)), 1.243e-001 +2.000*(-1.638e-001* 7.536e-003*cos(t)+ 9.865e-001* 6.259e-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.440e-001* 8.341e-003*cos(t)+ 5.363e-001* 7.658e-003*sin(t)), 1.112e-001 +2.000*(-5.363e-001* 8.341e-003*cos(t)+ 8.440e-001* 7.658e-003*sin(t)) not -# Age 95, p23 - p21 -set label "95" at 2.085e-001, 9.953e-002 center -replot 2.085e-001+ 2.000*( 9.795e-001* 1.210e-002*cos(t)+ 2.015e-001* 1.028e-002*sin(t)), 9.953e-002 +2.000*(-2.015e-001* 1.210e-002*cos(t)+ 9.795e-001* 1.028e-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; + +# 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; diff --git a/html/doc/biaspar.htm b/html/doc/biaspar.htm index ccd3c10..f0f9984 100644 --- a/html/doc/biaspar.htm +++ b/html/doc/biaspar.htm @@ -1,51 +1,51 @@ - -IMaCh rbiaspar.txt - Imach version 0.07a, May 2004, 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
-