From b0e8cfb7f4302af182c2b4f6e96a46815b0ed0e6 Mon Sep 17 00:00:00 2001 From: "N. Brouard" Date: Tue, 4 Aug 2015 07:15:27 +0000 Subject: [PATCH] 0.98q4 --- html/doc/biaspar-cov.htm | 30 - html/doc/biaspar.gp | 245 - html/doc/biaspar.htm | 47 +- html/doc/biaspar.imach | 4 +- html/doc/biaspar/ilkrbiaspar.txt | 41168 ++++++++++++++--------------- 5 files changed, 20592 insertions(+), 20902 deletions(-) diff --git a/html/doc/biaspar-cov.htm b/html/doc/biaspar-cov.htm index 4cc3b15..e69de29 100644 --- a/html/doc/biaspar-cov.htm +++ b/html/doc/biaspar-cov.htm @@ -1,30 +0,0 @@ - -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 fb51c34..e69de29 100644 --- a/html/doc/biaspar.gp +++ b/html/doc/biaspar.gp @@ -1,245 +0,0 @@ - -# 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 bc9a53e..2d47bd0 100644 --- a/html/doc/biaspar.htm +++ b/html/doc/biaspar.htm @@ -1,13 +1,14 @@ - -IMaCh 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=.
+ +IMaCh biaspar.htm + Imach version 0.98q2, April 2015,INED-EUROREVES-Institut de longevite-Japan Society for the Promotion of Science (Grant-in-Aid for Scientific Research 25293121), Intel Software 2015
$Revision$ $Date$
+Title=1st_example
Datafile=data1.txt Firstpass=1 Lastpass=4 Stepm=1 Weight=0 Model=

Total number of observations=8270
@@ -15,39 +16,3 @@ 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
-