From 997875c1d6824f0441c03316db33cdf2f467aa91 Mon Sep 17 00:00:00 2001 From: "N. Brouard" Date: Fri, 24 May 2002 13:00:54 +0000 Subject: [PATCH] Change of gif to png. Clarify some of the outputs. Add a new graph: pij*YEARM/stepm in order to get pseudo-incidence per year. --- src/imach.c | 141 +++++++++++++++++++++++++++++----------------------- 1 file changed, 80 insertions(+), 61 deletions(-) diff --git a/src/imach.c b/src/imach.c index e00b7e6..f542bab 100644 --- a/src/imach.c +++ b/src/imach.c @@ -1665,7 +1665,7 @@ void evsij(char fileres[], double ***eij, double x[], int nlstate, int stepm, in for(j=1;j<=nlstate*2;j++) varhe[i][j][(int)age] =0.; - printf("%d||",(int)age);fflush(stdout); + printf("%d|",(int)age);fflush(stdout); for(h=0;h<=nhstepm-1;h++){ for(k=0;k<=nhstepm-1;k++){ matprod2(dnewm,trgradg[h],1,nlstate*2,1,npar,1,npar,matcov); @@ -1724,7 +1724,7 @@ void varevsij(char fileres[], double ***vareij, double **matcov, double x[], dou double age,agelim, hf; int theta; - fprintf(ficresvij,"# Covariances of life expectancies\n"); + fprintf(ficresvij,"# Variance and covariance of health expectancies e.j \n# (weighted average of eij where weights are the stable prevalence in health states i\n"); fprintf(ficresvij,"# Age"); for(i=1; i<=nlstate;i++) for(j=1; j<=nlstate;j++) @@ -1859,7 +1859,7 @@ void varprevlim(char fileres[], double **varpl, double **matcov, double x[], dou double age,agelim; int theta; - fprintf(ficresvpl,"# Standard deviation of prevalences limit\n"); + fprintf(ficresvpl,"# Standard deviation of prevalence's limit\n"); fprintf(ficresvpl,"# Age"); for(i=1; i<=nlstate;i++) fprintf(ficresvpl," %1d-%1d",i,i); @@ -2061,10 +2061,12 @@ fprintf(ficresprob,"#One-step probabilities and standard deviation in parenthese /******************* Printing html file ***********/ void printinghtml(char fileres[], char title[], char datafile[], int firstpass, \ - int lastpass, int stepm, int weightopt, char model[],\ - int imx,int jmin, int jmax, double jmeanint,char optionfile[], \ - char optionfilehtm[],char rfileres[], char optionfilegnuplot[],\ - char version[], int popforecast, int estepm ){ + int lastpass, int stepm, int weightopt, char model[],\ + int imx,int jmin, int jmax, double jmeanint,char optionfile[], \ + char optionfilehtm[],char rfileres[], char optionfilegnuplot[],\ + char version[], int popforecast, int estepm ,/* \ */ + double jprev1, double mprev1,double anprev1, \ + double jprev2, double mprev2,double anprev2){ int jj1, k1, i1, cpt; FILE *fichtm; /*char optionfilehtm[FILENAMELENGTH];*/ @@ -2075,26 +2077,30 @@ void printinghtml(char fileres[], char title[], char datafile[], int firstpass, printf("Problem with %s \n",optionfilehtm), exit(0); } - fprintf(fichtm," %s

\n + fprintf(fichtm," %s

\n Title=%s
Datafile=%s Firstpass=%d Lastpass=%d Stepm=%d Weight=%d Model=%s
\n \n Total number of observations=%d
\n Interval (in months) between two waves: Min=%d Max=%d Mean=%.2lf
\n -

Outputs files
\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
\n",version,title,datafile,firstpass,lastpass,stepm, weightopt,model,imx,jmin,jmax,jmean,fileres,fileres,optionfilegnuplot,optionfilegnuplot); + + fprintf(fichtm,"
- 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
- 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 @@ -2117,10 +2123,15 @@ fprintf(fichtm,"
- Graphs