version 1.51, 2002/07/19 12:22:25
|
version 1.52, 2002/07/19 18:49:30
|
Line 1862 void varevsij(char optionfilefiname[], d
|
Line 1862 void varevsij(char optionfilefiname[], d
|
exit(0);
|
exit(0);
|
}
|
}
|
else{
|
else{
|
fprintf(fichtm,"\n<li><h4> Computing step probabilities of dying and weighted average (i.e global mortality independent of initial healh state)</h4></li>\n");
|
fprintf(fichtm,"\n<li><h4> Computing probabilities of dying as a weighted average (i.e global mortality independent of initial healh state)</h4></li>\n");
|
}
|
}
|
varppt = matrix(nlstate+1,nlstate+ndeath,nlstate+1,nlstate+ndeath);
|
varppt = matrix(nlstate+1,nlstate+ndeath,nlstate+1,nlstate+ndeath);
|
|
|
Line 2013 void varevsij(char optionfilefiname[], d
|
Line 2013 void varevsij(char optionfilefiname[], d
|
/* This for computing force of mortality (h=1)as a weighted average */
|
/* This for computing force of mortality (h=1)as a weighted average */
|
for(j=nlstate+1,gmp[j]=0.;j<=nlstate+ndeath;j++){
|
for(j=nlstate+1,gmp[j]=0.;j<=nlstate+ndeath;j++){
|
for(i=1; i<= nlstate; i++)
|
for(i=1; i<= nlstate; i++)
|
gmp[j] += prlim[i][i]*p3mat[i][j][1];
|
gmp[j] += prlim[i][i]*p3mat[i][j][1];
|
}
|
}
|
/* end force of mortality */
|
/* end force of mortality */
|
|
|
Line 2049 void varevsij(char optionfilefiname[], d
|
Line 2049 void varevsij(char optionfilefiname[], d
|
fprintf(ficgp,"\n replot \"%s\" u 1:(($3+1.96*$4)*%6.3f) t \"95\%% interval\" w l 2 ",fileresprobmorprev,YEARM/estepm);
|
fprintf(ficgp,"\n replot \"%s\" u 1:(($3+1.96*$4)*%6.3f) t \"95\%% interval\" w l 2 ",fileresprobmorprev,YEARM/estepm);
|
fprintf(ficgp,"\n replot \"%s\" u 1:(($3-1.96*$4)*%6.3f) not w l 2 ",fileresprobmorprev,YEARM/estepm);
|
fprintf(ficgp,"\n replot \"%s\" u 1:(($3-1.96*$4)*%6.3f) not w l 2 ",fileresprobmorprev,YEARM/estepm);
|
fprintf(fichtm,"\n<br> File (multiple files are possible if covariates are present): <A href=\"%s\">%s</a>\n",fileresprobmorprev,fileresprobmorprev);
|
fprintf(fichtm,"\n<br> File (multiple files are possible if covariates are present): <A href=\"%s\">%s</a>\n",fileresprobmorprev,fileresprobmorprev);
|
fprintf(fichtm,"\n<br> Probability is computed over estepm=%d months and then divided by estepm and multiplied by %.0f in order to have the probability to die over a year <br> <img src=\"varmuptjgr%s%s.png\"> <br>\n", stepm,YEARM,digitp,digit);
|
fprintf(fichtm,"\n<br> Probability is computed over estepm=%d months. <br> <img src=\"varmuptjgr%s%s.png\"> <br>\n", stepm,digitp,digit);
|
|
/* fprintf(fichtm,"\n<br> Probability is computed over estepm=%d months and then divided by estepm and multiplied by %.0f in order to have the probability to die over a year <br> <img src=\"varmuptjgr%s%s.png\"> <br>\n", stepm,YEARM,digitp,digit);
|
|
*/
|
fprintf(ficgp,"\nset out \"varmuptjgr%s%s.png\";replot;",digitp,digit);
|
fprintf(ficgp,"\nset out \"varmuptjgr%s%s.png\";replot;",digitp,digit);
|
|
|
free_vector(xp,1,npar);
|
free_vector(xp,1,npar);
|
Line 2420 void varprob(char optionfilefiname[], do
|
Line 2422 void varprob(char optionfilefiname[], do
|
fprintf(ficgp,"\nset ter png small\nset size 0.65,0.65");
|
fprintf(ficgp,"\nset ter png small\nset size 0.65,0.65");
|
fprintf(fichtm,"\n<br>Ellipsoids of confidence cov(p%1d%1d,p%1d%1d) expressed in year<sup>-1</sup> :<a href=\"varpijgr%s%d%1d%1d-%1d%1d.png\">varpijgr%s%d%1d%1d-%1d%1d.png</A>, ",k1,l1,k2,l2,optionfilefiname, j1,k1,l1,k2,l2,optionfilefiname, j1,k1,l1,k2,l2);
|
fprintf(fichtm,"\n<br>Ellipsoids of confidence cov(p%1d%1d,p%1d%1d) expressed in year<sup>-1</sup> :<a href=\"varpijgr%s%d%1d%1d-%1d%1d.png\">varpijgr%s%d%1d%1d-%1d%1d.png</A>, ",k1,l1,k2,l2,optionfilefiname, j1,k1,l1,k2,l2,optionfilefiname, j1,k1,l1,k2,l2);
|
fprintf(fichtm,"\n<br><img src=\"varpijgr%s%d%1d%1d-%1d%1d.png\"> ",optionfilefiname, j1,k1,l1,k2,l2);
|
fprintf(fichtm,"\n<br><img src=\"varpijgr%s%d%1d%1d-%1d%1d.png\"> ",optionfilefiname, j1,k1,l1,k2,l2);
|
|
fprintf(fichtm,"\n<br> Correlation at age %d (%.3f),",(int) age, c12);
|
fprintf(ficgp,"\nset out \"varpijgr%s%d%1d%1d-%1d%1d.png\"",optionfilefiname, j1,k1,l1,k2,l2);
|
fprintf(ficgp,"\nset out \"varpijgr%s%d%1d%1d-%1d%1d.png\"",optionfilefiname, j1,k1,l1,k2,l2);
|
fprintf(ficgp,"\nset label \"%d\" at %11.3e,%11.3e center",(int) age, mu1,mu2);
|
fprintf(ficgp,"\nset label \"%d\" at %11.3e,%11.3e center",(int) age, mu1,mu2);
|
fprintf(ficgp,"\n# Age %d, p%1d%1d - p%1d%1d",(int) age, k1,l1,k2,l2);
|
fprintf(ficgp,"\n# Age %d, p%1d%1d - p%1d%1d",(int) age, k1,l1,k2,l2);
|
Line 2428 void varprob(char optionfilefiname[], do
|
Line 2431 void varprob(char optionfilefiname[], do
|
mu2,std,v21,sqrt(lc1),v22,sqrt(lc2));
|
mu2,std,v21,sqrt(lc1),v22,sqrt(lc2));
|
}else{
|
}else{
|
first=0;
|
first=0;
|
|
fprintf(fichtm," %d (%.3f),",(int) age, c12);
|
fprintf(ficgp,"\n# Age %d, p%1d%1d - p%1d%1d",(int) age, k1,l1,k2,l2);
|
fprintf(ficgp,"\n# Age %d, p%1d%1d - p%1d%1d",(int) age, k1,l1,k2,l2);
|
fprintf(ficgp,"\nset label \"%d\" at %11.3e,%11.3e center",(int) age, mu1,mu2);
|
fprintf(ficgp,"\nset label \"%d\" at %11.3e,%11.3e center",(int) age, mu1,mu2);
|
fprintf(ficgp,"\nreplot %11.3e+ %.3f*(%11.3e*%11.3e*cos(t)+%11.3e*%11.3e*sin(t)), %11.3e +%.3f*(%11.3e*%11.3e*cos(t)+%11.3e*%11.3e*sin(t)) not",\
|
fprintf(ficgp,"\nreplot %11.3e+ %.3f*(%11.3e*%11.3e*cos(t)+%11.3e*%11.3e*sin(t)), %11.3e +%.3f*(%11.3e*%11.3e*cos(t)+%11.3e*%11.3e*sin(t)) not",\
|