\n +
+
- Parameter files
\n - Copy of the parameter file: o%s
\n - - Gnuplot file name: %s
\n - - Observed prevalence in each state: p%s
\n - - Stationary prevalence in each state: pl%s
\n - - Transition probabilities: pij%s
\n - - Life expectancies by age and initial health status (estepm=%2d months): e%s
\n",version,title,datafile,firstpass,lastpass,stepm, weightopt,model,imx,jmin,jmax,jmean,fileres,fileres,optionfilegnuplot,optionfilegnuplot,fileres,fileres,fileres,fileres,fileres,fileres,estepm,fileres,fileres); - - fprintf(fichtm,"\n - - Parameter file with estimated parameters and the covariance matrix: %s
\n - - Variance of one-step probabilities: prob%s
\n - - Variances of life expectancies by age and initial health status (estepm=%d months): v%s
\n - - Health expectancies with their variances: t%s
\n - - Standard deviation of stationary prevalences: vpl%s
\n",rfileres,rfileres,fileres,fileres, estepm, fileres,fileres,fileres,fileres,fileres,fileres); + - Gnuplot file name: %s
- Result files (first order: no variance)
\n + - Observed prevalence in each state (during the period defined between %.lf/%.lf/%.lf and %.lf/%.lf/%.lf): p%s
\n + - Estimated transition probabilities over %d (stepm) months: pij%s
\n + - Stable prevalence in each health state: pl%s
\n + - Life expectancies by age and initial health status (estepm=%2d months): + e%s
\n ", \
+ jprev1, mprev1,anprev1,jprev2, mprev2,anprev2,fileres,fileres,stepm,fileres,fileres,fileres,fileres,estepm,fileres,fileres);
+
+ fprintf(fichtm,"\n - Result files (second order: variances)
\n + - Parameter file with estimated parameters and covariance matrix: %s
\n + - Variance of one-step probabilities: prob%s
\n + - Variances and covariances of life expectancies by age and initial health status (estepm=%d months): v%s
\n + - Health expectancies with their variances (no covariance): t%s
\n + - Standard deviation of stable prevalences: vpl%s
\n",rfileres,rfileres,fileres,fileres, estepm, fileres,fileres,fileres,fileres,fileres,fileres); if(popforecast==1) fprintf(fichtm,"\n - Prevalences forecasting: f%s
\n @@ -2110,17 +2117,22 @@ fprintf(fichtm," - Graphs
"); jj1=0; for(k1=1; k1<=m;k1++){ for(i1=1; i1<=ncodemax[k1];i1++){ - jj1++; - if (cptcovn > 0) { - fprintf(fichtm,"
************ Results for covariates"); - for (cpt=1; cpt<=cptcoveff;cpt++) - fprintf(fichtm," V%d=%d ",Tvaraff[cpt],nbcode[Tvaraff[cpt]][codtab[jj1][cpt]]); - fprintf(fichtm," ************\n
"); - } - fprintf(fichtm,"
- Probabilities: pe%s%d.png
-
![](\"pe%s%d.png\")
",strtok(optionfile, "."),jj1,strtok(optionfile, "."),jj1);
+ jj1++;
+ if (cptcovn > 0) {
+ fprintf(fichtm,"************ Results for covariates"); + for (cpt=1; cpt<=cptcoveff;cpt++) + fprintf(fichtm," V%d=%d ",Tvaraff[cpt],nbcode[Tvaraff[cpt]][codtab[jj1][cpt]]); + fprintf(fichtm," ************\n
"); + } + /* Pij */ + fprintf(fichtm,"
- Pij or Conditional probabilities to be observed in state j being in state i %d (stepm) months before: pe%s%d1.png
+
![](\"pe%s%d1.png\")
",stepm,strtok(optionfile, "."),jj1,strtok(optionfile, "."),jj1);
+ /* Quasi-incidences */
+ fprintf(fichtm,"- Pij or Conditional probabilities to be observed in state j being in state i %d (stepm) months before but expressed in per year i.e. quasi incidences if stepm is small and probabilities too: pe%s%d2.png
+
![](\"pe%s%d2.png\")
",stepm,strtok(optionfile, "."),jj1,strtok(optionfile, "."),jj1);
+ /* Stable prevalence in each health state */
for(cpt=1; cpt+ fprintf(fichtm,"
- Stable prevalence in each health state : p%s%d%d.png
![](\"p%s%d%d.png\")
",strtok(optionfile, "."),cpt,jj1,strtok(optionfile, "."),cpt,jj1);
}
for(cpt=1; cpt<=nlstate;cpt++) {
@@ -2137,7 +2149,7 @@ health expectancies in states (1) and (2
![](\"e%s%d.png\")
",strtok(optionfile, "."),jj1,strtok(optionfile, "."),jj1);
fprintf(fichtm,"\n");
}
- }
+ }
fclose(fichtm);
}
@@ -2145,7 +2157,7 @@ fclose(fichtm);
void printinggnuplot(char fileres[],char optionfilefiname[],char optionfile[],char optionfilegnuplot[], double ageminpar, double agemaxpar, double fage , char pathc[], double p[]){
int m,cpt,k1,i,k,j,jk,k2,k3,ij,l;
-
+ int ng;
strcpy(optionfilegnuplot,optionfilefiname);
strcat(optionfilegnuplot,".gp.txt");
if((ficgp=fopen(optionfilegnuplot,"w"))==NULL) {
@@ -2162,11 +2174,11 @@ m=pow(2,cptcoveff);
for (k1=1; k1<= m ; k1 ++) {
#ifdef windows
- fprintf(ficgp,"\nset out \"v%s%d%d.png\" \n\n",strtok(optionfile, "."),cpt,k1);
+ fprintf(ficgp,"\nset out \"v%s%d%d.png\" \n",strtok(optionfile, "."),cpt,k1);
fprintf(ficgp,"set xlabel \"Age\" \nset ylabel \"Probability\" \nset ter png small\nset size 0.65,0.65\nplot [%.f:%.f] \"vpl%s\" every :::%d::%d u 1:2 \"\%%lf",ageminpar,fage,fileres,k1-1,k1-1);
#endif
#ifdef unix
-fprintf(ficgp,"\nset out \"v%s%d%d.png\" \n\n",strtok(optionfile, "."),cpt,k1);
+fprintf(ficgp,"\nset out \"v%s%d%d.png\" \n",strtok(optionfile, "."),cpt,k1);
fprintf(ficgp,"set xlabel \"Age\" \nset ylabel \"Probability\" \nplot [%.f:%.f] \"vpl%s\" u 1:2 \"\%%lf",ageminpar,fage,fileres);
#endif
@@ -2193,7 +2205,7 @@ fprintf(ficgp,"\nset ter png small\nset
/*2 eme*/
for (k1=1; k1<= m ; k1 ++) {
- fprintf(ficgp,"\nset out \"e%s%d.png\" \n\n",strtok(optionfile, "."),k1);
+ fprintf(ficgp,"\nset out \"e%s%d.png\" \n",strtok(optionfile, "."),k1);
fprintf(ficgp,"set ylabel \"Years\" \nset ter png small\nset size 0.65,0.65\nplot [%.f:%.f] ",ageminpar,fage);
for (i=1; i<= nlstate+1 ; i ++) {
@@ -2226,7 +2238,7 @@ fprintf(ficgp,"\nset ter png small\nset
for (k1=1; k1<= m ; k1 ++) {
for (cpt=1; cpt<= nlstate ; cpt ++) {
k=2+nlstate*(2*cpt-2);
- fprintf(ficgp,"\nset out \"exp%s%d%d.png\" \n\n",strtok(optionfile, "."),cpt,k1);
+ fprintf(ficgp,"\nset out \"exp%s%d%d.png\" \n",strtok(optionfile, "."),cpt,k1);
fprintf(ficgp,"set ter png small\nset size 0.65,0.65\nplot [%.f:%.f] \"e%s\" every :::%d::%d u 1:%d t \"e%d1\" w l",ageminpar,fage,fileres,k1-1,k1-1,k,cpt);
/*fprintf(ficgp,",\"e%s\" every :::%d::%d u 1:($%d-2*$%d) \"\%%lf ",fileres,k1-1,k1-1,k,k+1);
for (i=1; i<= nlstate*2 ; i ++) fprintf(ficgp,"\%%lf (\%%lf) ");
@@ -2247,7 +2259,7 @@ fprintf(ficgp,"\" t \"e%d1\" w l",cpt);
for (k1=1; k1<= m ; k1 ++) {
for (cpt=1; cpt