imach/src/imach.c - diff

Return to imach.c CVS log

Up to [Local Repository] / imach / src

Diff for /imach/src/imach.c between versions 1.197 and 1.201

version 1.197, 2015/09/01 18:24:39	version 1.201, 2015/09/15 17:34:58
Line 1	Line 1
/* $Id$	/* $Id$
$State$	$State$
$Log$	$Log$
	Revision 1.201 2015/09/15 17:34:58 brouard
	Summary: 0.98r0

	- Some new graphs like suvival functions
	- Some bugs fixed like model=1+age+V2.

	Revision 1.200 2015/09/09 16:53:55 brouard
	Summary: Big bug thanks to Flavia

	Even model=1+age+V2. did not work anymore

	Revision 1.199 2015/09/07 14:09:23 brouard
	Summary: 0.98q6 changing default small png format for graph to vectorized svg.

	Revision 1.198 2015/09/03 07:14:39 brouard
	Summary: 0.98q5 Flavia

Revision 1.197 2015/09/01 18:24:39 brouard	Revision 1.197 2015/09/01 18:24:39 brouard
* empty log message *	* empty log message *

Line 785 char command[FILENAMELENGTH];	Line 802 char command[FILENAMELENGTH];
int outcmd=0;	int outcmd=0;

char fileres[FILENAMELENGTH], filerespij[FILENAMELENGTH], filereso[FILENAMELENGTH], rfileres[FILENAMELENGTH];	char fileres[FILENAMELENGTH], filerespij[FILENAMELENGTH], filereso[FILENAMELENGTH], rfileres[FILENAMELENGTH];
	char fileresu[FILENAMELENGTH]; /* Without r in front */
char filelog[FILENAMELENGTH]; /* Log file */	char filelog[FILENAMELENGTH]; /* Log file */
char filerest[FILENAMELENGTH];	char filerest[FILENAMELENGTH];
char fileregp[FILENAMELENGTH];	char fileregp[FILENAMELENGTH];
Line 875 double idx;	Line 892 double idx;
int *nbcode, Tvar; /*< model=V2 => Tvar[1]= 2 /	int *nbcode, Tvar; /*< model=V2 => Tvar[1]= 2 /
int *Tage;	int *Tage;
int Ndum; /* Freq of modality (tricode */	int Ndum; /* Freq of modality (tricode */
int codtab; /< codtab=imatrix(1,100,1,10); */	/* int *codtab;/ /*< codtab=imatrix(1,100,1,10); /
int *Tvard, Tprod, cptcovprod, *Tvaraff;	int *Tvard, Tprod, cptcovprod, *Tvaraff;
double lsurv, lpop, *tpop;	double lsurv, lpop, *tpop;

Line 1911 double prevalim(double prlim, int nl	Line 1928 double prevalim(double prlim, int nl
if(nagesqr==1)	if(nagesqr==1)
cov[3]= agefin*agefin;;	cov[3]= agefin*agefin;;
for (k=1; k<=cptcovn;k++) {	for (k=1; k<=cptcovn;k++) {
cov[2+nagesqr+k]=nbcode[Tvar[k]][codtab[ij][Tvar[k]]];	/* cov[2+nagesqr+k]=nbcode[Tvar[k]][codtabm(ij,Tvar[k])]; */
/* printf("prevalim ij=%d k=%d Tvar[%d]=%d nbcode=%d cov=%lf codtab[%d][Tvar[%d]]=%d \n",ij,k, k, Tvar[k],nbcode[Tvar[k]][codtab[ij][Tvar[k]]],cov[2+k], ij, k, codtab[ij][Tvar[k]]);*/	cov[2+nagesqr+k]=nbcode[Tvar[k]][codtabm(ij,k)];
	/* printf("prevalim ij=%d k=%d Tvar[%d]=%d nbcode=%d cov=%lf codtabm(%d,Tvar[%d])=%d \n",ij,k, k, Tvar[k],nbcode[Tvar[k]][codtabm(ij,Tvar[k])],cov[2+k], ij, k, codtabm(ij,Tvar[k])]); */
}	}
/wrong? for (k=1; k<=cptcovage;k++) cov[2+Tage[k]]=cov[2+Tage[k]]cov[2]; */	/wrong? for (k=1; k<=cptcovage;k++) cov[2+Tage[k]]=cov[2+Tage[k]]cov[2]; */
for (k=1; k<=cptcovage;k++) cov[2+nagesqr+Tage[k]]=nbcode[Tvar[k]][codtab[ij][Tvar[k]]]*cov[2];	/* for (k=1; k<=cptcovage;k++) cov[2+nagesqr+Tage[k]]=nbcode[Tvar[k]][codtabm(ij,Tvar[k])]cov[2]; /
	for (k=1; k<=cptcovage;k++) cov[2+nagesqr+Tage[k]]=nbcode[Tvar[k]][codtabm(ij,k)]*cov[2];
for (k=1; k<=cptcovprod;k++) /* Useless */	for (k=1; k<=cptcovprod;k++) /* Useless */
cov[2+nagesqr+Tprod[k]]=nbcode[Tvard[k][1]][codtab[ij][Tvard[k][1]]] * nbcode[Tvard[k][2]][codtab[ij][Tvard[k][2]]];	/* cov[2+nagesqr+Tprod[k]]=nbcode[Tvard[k][1]][codtabm(ij,Tvard[k][1])] * nbcode[Tvard[k][2]][codtabm(ij,Tvard[k][2])]; */
	cov[2+nagesqr+Tprod[k]]=nbcode[Tvard[k][1]][codtabm(ij,k)] * nbcode[Tvard[k][2]][codtabm(ij,k)];

/printf("ij=%d cptcovprod=%d tvar=%d ", ij, cptcovprod, Tvar[1]);/	/printf("ij=%d cptcovprod=%d tvar=%d ", ij, cptcovprod, Tvar[1]);/
/printf("ij=%d cov[3]=%lf cov[4]=%lf \n",ij, cov[3],cov[4]);/	/printf("ij=%d cov[3]=%lf cov[4]=%lf \n",ij, cov[3],cov[4]);/
Line 2090 double *hpxij(double *po, int nhstep	Line 2110 double *hpxij(double *po, int nhstep
if(nagesqr==1)	if(nagesqr==1)
cov[3]= agexact*agexact;	cov[3]= agexact*agexact;
for (k=1; k<=cptcovn;k++)	for (k=1; k<=cptcovn;k++)
cov[2+nagesqr+k]=nbcode[Tvar[k]][codtab[ij][Tvar[k]]];	cov[2+nagesqr+k]=nbcode[Tvar[k]][codtabm(ij,k)];
	/* cov[2+nagesqr+k]=nbcode[Tvar[k]][codtabm(ij,Tvar[k])]; */
for (k=1; k<=cptcovage;k++) /* Should start at cptcovn+1 */	for (k=1; k<=cptcovage;k++) /* Should start at cptcovn+1 */
/* cov[2+Tage[k]]=cov[2+Tage[k]]cov[2]; /	/* cov[2+Tage[k]]=cov[2+Tage[k]]cov[2]; /
cov[2+nagesqr+Tage[k]]=nbcode[Tvar[Tage[k]]][codtab[ij][Tvar[Tage[k]]]]*cov[2];	cov[2+nagesqr+Tage[k]]=nbcode[Tvar[Tage[k]]][codtabm(ij,k)]*cov[2];
	/* cov[2+nagesqr+Tage[k]]=nbcode[Tvar[Tage[k]]][codtabm(ij,Tvar[Tage[k]])]cov[2]; /
for (k=1; k<=cptcovprod;k++) /* Useless because included in cptcovn */	for (k=1; k<=cptcovprod;k++) /* Useless because included in cptcovn */
cov[2+nagesqr+Tprod[k]]=nbcode[Tvard[k][1]][codtab[ij][Tvard[k][1]]]*nbcode[Tvard[k][2]][codtab[ij][Tvard[k][2]]];	cov[2+nagesqr+Tprod[k]]=nbcode[Tvard[k][1]][codtabm(ij,k)]*nbcode[Tvard[k][2]][codtabm(ij,k)];
	/* cov[2+nagesqr+Tprod[k]]=nbcode[Tvard[k][1]][codtabm(ij,Tvard[k][1])]nbcode[Tvard[k][2]][codtabm(ij,Tvard[k][2])]; /


/printf("hxi cptcov=%d cptcode=%d\n",cptcov,cptcode);/	/printf("hxi cptcov=%d cptcode=%d\n",cptcov,cptcode);/
Line 2553 void likelione(FILE *ficres,double p[],	Line 2576 void likelione(FILE *ficres,double p[],
int k;	int k;

if(globpri !=0){ / Just counts and sums, no printings */	if(globpri !=0){ / Just counts and sums, no printings */
strcpy(fileresilk,"ilk");	strcpy(fileresilk,"ILK_");
strcat(fileresilk,fileres);	strcat(fileresilk,fileres);
if((ficresilk=fopen(fileresilk,"w"))==NULL) {	if((ficresilk=fopen(fileresilk,"w"))==NULL) {
printf("Problem with resultfile: %s\n", fileresilk);	printf("Problem with resultfile: %s\n", fileresilk);
Line 2570 void likelione(FILE *ficres,double p[],	Line 2593 void likelione(FILE *ficres,double p[],
fretone=(funcone)(p);	fretone=(funcone)(p);
if(*globpri !=0){	if(*globpri !=0){
fclose(ficresilk);	fclose(ficresilk);
fprintf(fichtm,"\n<br>File of contributions to the likelihood: <a href=\"%s\">%s</a><br>\n",subdirf(fileresilk),subdirf(fileresilk));	fprintf(fichtm,"\n<br>File of contributions to the likelihood (if mle=1): <a href=\"%s\">%s</a><br>\n",subdirf(fileresilk),subdirf(fileresilk));
fflush(fichtm);	fflush(fichtm);
}	}
return;	return;
Line 2603 void mlikeli(FILE *ficres,double p[], in	Line 2626 void mlikeli(FILE *ficres,double p[], in
for (j=1;j<=npar;j++)	for (j=1;j<=npar;j++)
xi[i][j]=(i==j ? 1.0 : 0.0);	xi[i][j]=(i==j ? 1.0 : 0.0);
printf("Powell\n"); fprintf(ficlog,"Powell\n");	printf("Powell\n"); fprintf(ficlog,"Powell\n");
strcpy(filerespow,"pow");	strcpy(filerespow,"POW_");
strcat(filerespow,fileres);	strcat(filerespow,fileres);
if((ficrespow=fopen(filerespow,"w"))==NULL) {	if((ficrespow=fopen(filerespow,"w"))==NULL) {
printf("Problem with resultfile: %s\n", filerespow);	printf("Problem with resultfile: %s\n", filerespow);
Line 2930 void freqsummary(char fileres[], int ia	Line 2953 void freqsummary(char fileres[], int ia

pp=vector(1,nlstate);	pp=vector(1,nlstate);
prop=matrix(1,nlstate,iagemin,iagemax+3);	prop=matrix(1,nlstate,iagemin,iagemax+3);
strcpy(fileresp,"p");	strcpy(fileresp,"P_");
strcat(fileresp,fileres);	strcat(fileresp,fileresu);
if((ficresp=fopen(fileresp,"w"))==NULL) {	if((ficresp=fopen(fileresp,"w"))==NULL) {
printf("Problem with prevalence resultfile: %s\n", fileresp);	printf("Problem with prevalence resultfile: %s\n", fileresp);
fprintf(ficlog,"Problem with prevalence resultfile: %s\n", fileresp);	fprintf(ficlog,"Problem with prevalence resultfile: %s\n", fileresp);
Line 2966 void freqsummary(char fileres[], int ia	Line 2989 void freqsummary(char fileres[], int ia
bool=1;	bool=1;
if (cptcovn>0) { /* Filter is here: Must be looked at for model=V1+V2+V3+V4 */	if (cptcovn>0) { /* Filter is here: Must be looked at for model=V1+V2+V3+V4 */
for (z1=1; z1<=cptcoveff; z1++)	for (z1=1; z1<=cptcoveff; z1++)
if (covar[Tvaraff[z1]][i]!= nbcode[Tvaraff[z1]][codtab[j1][z1]]){	if (covar[Tvaraff[z1]][i]!= nbcode[Tvaraff[z1]][codtabm(j1,z1)]){
/* Tests if the value of each of the covariates of i is equal to filter j1 */	/* Tests if the value of each of the covariates of i is equal to filter j1 */
bool=0;	bool=0;
/* printf("bool=%d i=%d, z1=%d, Tvaraff[%d]=%d, covar[Tvarff][%d]=%2f, codtab[%d][%d]=%d, nbcode[Tvaraff][codtab[%d][%d]=%d, j1=%d\n",	/* printf("bool=%d i=%d, z1=%d, Tvaraff[%d]=%d, covar[Tvarff][%d]=%2f, codtabm(%d,%d)=%d, nbcode[Tvaraff][codtabm(%d,%d)=%d, j1=%d\n",
bool,i,z1, z1, Tvaraff[z1],i,covar[Tvaraff[z1]][i],j1,z1,codtab[j1][z1],	bool,i,z1, z1, Tvaraff[z1],i,covar[Tvaraff[z1]][i],j1,z1,codtabm(j1,z1),
j1,z1,nbcode[Tvaraff[z1]][codtab[j1][z1]],j1);*/	j1,z1,nbcode[Tvaraff[z1]][codtabm(j1,z1)],j1);*/
/* For j1=7 in V1+V2+V3+V4 = 0 1 1 0 and codtab[7][3]=1 and nbcde[3][?]=1*/	/* For j1=7 in V1+V2+V3+V4 = 0 1 1 0 and codtabm(7,3)=1 and nbcde[3][?]=1*/
}	}
}	}

Line 3001 void freqsummary(char fileres[], int ia	Line 3024 void freqsummary(char fileres[], int ia
pstamp(ficresp);	pstamp(ficresp);
if (cptcovn>0) {	if (cptcovn>0) {
fprintf(ficresp, "\n#********** Variable ");	fprintf(ficresp, "\n#********** Variable ");
for (z1=1; z1<=cptcoveff; z1++) fprintf(ficresp, "V%d=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtab[j1][z1]]);	for (z1=1; z1<=cptcoveff; z1++) fprintf(ficresp, "V%d=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtabm(j1,z1)]);
fprintf(ficresp, "**********\n#");	fprintf(ficresp, "**********\n#");
fprintf(ficlog, "\n#********** Variable ");	fprintf(ficlog, "\n#********** Variable ");
for (z1=1; z1<=cptcoveff; z1++) fprintf(ficlog, "V%d=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtab[j1][z1]]);	for (z1=1; z1<=cptcoveff; z1++) fprintf(ficlog, "V%d=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtabm(j1,z1)]);
fprintf(ficlog, "**********\n#");	fprintf(ficlog, "**********\n#");
}	}
for(i=1; i<=nlstate;i++)	for(i=1; i<=nlstate;i++)
Line 3132 void prevalence(double ***probs, double	Line 3155 void prevalence(double ***probs, double
bool=1;	bool=1;
if (cptcovn>0) {	if (cptcovn>0) {
for (z1=1; z1<=cptcoveff; z1++)	for (z1=1; z1<=cptcoveff; z1++)
if (covar[Tvaraff[z1]][i]!= nbcode[Tvaraff[z1]][codtab[j1][z1]])	if (covar[Tvaraff[z1]][i]!= nbcode[Tvaraff[z1]][codtabm(j1,z1)])
bool=0;	bool=0;
}	}
if (bool==1) {	if (bool==1) {
Line 3799 void varevsij(char optionfilefiname[], d	Line 3822 void varevsij(char optionfilefiname[], d

if(popbased==1){	if(popbased==1){
if(mobilav!=0)	if(mobilav!=0)
strcpy(digitp,"-populbased-mobilav-");	strcpy(digitp,"-POPULBASED-MOBILAV_");
else strcpy(digitp,"-populbased-nomobil-");	else strcpy(digitp,"-POPULBASED-NOMOBIL_");
}	}
else	else
strcpy(digitp,"-stablbased-");	strcpy(digitp,"-STABLBASED_");

if (mobilav!=0) {	if (mobilav!=0) {
mobaverage= ma3x(1, AGESUP,1,NCOVMAX, 1,NCOVMAX);	mobaverage= ma3x(1, AGESUP,1,NCOVMAX, 1,NCOVMAX);
Line 3813 void varevsij(char optionfilefiname[], d	Line 3836 void varevsij(char optionfilefiname[], d
}	}
}	}

strcpy(fileresprobmorprev,"prmorprev");	strcpy(fileresprobmorprev,"PRMORPREV-");
sprintf(digit,"%-d",ij);	sprintf(digit,"%-d",ij);
/printf("DIGIT=%s, ij=%d ijr=%-d\|\n",digit, ij,ij);/	/printf("DIGIT=%s, ij=%d ijr=%-d\|\n",digit, ij,ij);/
strcat(fileresprobmorprev,digit); /* Tvar to be done */	strcat(fileresprobmorprev,digit); /* Tvar to be done */
Line 3836 void varevsij(char optionfilefiname[], d	Line 3859 void varevsij(char optionfilefiname[], d
}	}
fprintf(ficresprobmorprev,"\n");	fprintf(ficresprobmorprev,"\n");
fprintf(ficgp,"\n# Routine varevsij");	fprintf(ficgp,"\n# Routine varevsij");
/* fprintf(fichtm, "#Local time at start: %s", strstart);*/	fprintf(ficgp,"\nunset title \n");
	/* fprintf(fichtm, "#Local time at start: %s", strstart);*/
fprintf(fichtm,"\n<li><h4> Computing probabilities of dying over estepm months as a weighted average (i.e global mortality independent of initial healh state)</h4></li>\n");	fprintf(fichtm,"\n<li><h4> Computing probabilities of dying over estepm months as a weighted average (i.e global mortality independent of initial healh state)</h4></li>\n");
fprintf(fichtm,"\n<br>%s <br>\n",digitp);	fprintf(fichtm,"\n<br>%s <br>\n",digitp);
/* } */	/* } */
Line 4048 void varevsij(char optionfilefiname[], d	Line 4072 void varevsij(char optionfilefiname[], d
free_vector(gmp,nlstate+1,nlstate+ndeath);	free_vector(gmp,nlstate+1,nlstate+ndeath);
free_matrix(gradgp,1,npar,nlstate+1,nlstate+ndeath);	free_matrix(gradgp,1,npar,nlstate+1,nlstate+ndeath);
free_matrix(trgradgp,nlstate+1,nlstate+ndeath,1,npar); /* mu or p point j*/	free_matrix(trgradgp,nlstate+1,nlstate+ndeath,1,npar); /* mu or p point j*/
fprintf(ficgp,"\nunset parametric;unset label; set ter png small size 320, 240");	/* fprintf(ficgp,"\nunset parametric;unset label; set ter png small size 320, 240"); */
	fprintf(ficgp,"\nunset parametric;unset label; set ter svg size 640, 480");
/* for(j=nlstate+1; j<= nlstate+ndeath; j++){ // Only the first actually */	/* for(j=nlstate+1; j<= nlstate+ndeath; j++){ // Only the first actually */
fprintf(ficgp,"\n set log y; unset log x;set xlabel \"Age\"; set ylabel \"Force of mortality (year-1)\";");	fprintf(ficgp,"\n set log y; unset log x;set xlabel \"Age\"; set ylabel \"Force of mortality (year-1)\";");
	fprintf(ficgp,"\nset out \"%s%s.svg\";",subdirf3(optionfilefiname,"VARMUPTJGR-",digitp),digit);
/* fprintf(ficgp,"\n plot \"%s\" u 1:($3%6.3f) not w l 1 ",fileresprobmorprev,YEARM/estepm); /	/* fprintf(ficgp,"\n plot \"%s\" u 1:($3%6.3f) not w l 1 ",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) 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); */
Line 4058 void varevsij(char optionfilefiname[], d	Line 4084 void varevsij(char optionfilefiname[], d
fprintf(ficgp,"\n replot \"%s\" u 1:(($3+1.96*$4)) t \"95%% interval\" w l lt 2 ",subdirf(fileresprobmorprev));	fprintf(ficgp,"\n replot \"%s\" u 1:(($3+1.96*$4)) t \"95%% interval\" w l lt 2 ",subdirf(fileresprobmorprev));
fprintf(ficgp,"\n replot \"%s\" u 1:(($3-1.96*$4)) not w l lt 2 ",subdirf(fileresprobmorprev));	fprintf(ficgp,"\n replot \"%s\" u 1:(($3-1.96*$4)) not w l lt 2 ",subdirf(fileresprobmorprev));
fprintf(fichtm,"\n<br> File (multiple files are possible if covariates are present): <A href=\"%s\">%s</a>\n",subdirf(fileresprobmorprev),subdirf(fileresprobmorprev));	fprintf(fichtm,"\n<br> File (multiple files are possible if covariates are present): <A href=\"%s\">%s</a>\n",subdirf(fileresprobmorprev),subdirf(fileresprobmorprev));
fprintf(fichtm,"\n<br> Probability is computed over estepm=%d months. <br> <img src=\"%s%s.png\"> <br>\n", estepm,subdirf3(optionfilefiname,"varmuptjgr",digitp),digit);	fprintf(fichtm,"\n<br> Probability is computed over estepm=%d months. <br> <img src=\"%s%s.svg\"> <br>\n", estepm,subdirf3(optionfilefiname,"VARMUPTJGR-",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(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.svg\"> <br>\n", stepm,YEARM,digitp,digit);
*/	*/
/* fprintf(ficgp,"\nset out \"varmuptjgr%s%s%s.png\";replot;",digitp,optionfilefiname,digit); */	/* fprintf(ficgp,"\nset out \"varmuptjgr%s%s%s.svg\";replot;",digitp,optionfilefiname,digit); */
fprintf(ficgp,"\nset out \"%s%s.png\";replot;\n",subdirf3(optionfilefiname,"varmuptjgr",digitp),digit);	fprintf(ficgp,"\nset out;\nset out \"%s%s.svg\";replot;set out;\n",subdirf3(optionfilefiname,"VARMUPTJGR-",digitp),digit);

free_vector(xp,1,npar);	free_vector(xp,1,npar);
free_matrix(doldm,1,nlstate,1,nlstate);	free_matrix(doldm,1,nlstate,1,nlstate);
Line 4180 void varprob(char optionfilefiname[], do	Line 4206 void varprob(char optionfilefiname[], do
char fileresprobcor[FILENAMELENGTH];	char fileresprobcor[FILENAMELENGTH];
double ***varpij;	double ***varpij;

strcpy(fileresprob,"prob");	strcpy(fileresprob,"PROB_");
strcat(fileresprob,fileres);	strcat(fileresprob,fileres);
if((ficresprob=fopen(fileresprob,"w"))==NULL) {	if((ficresprob=fopen(fileresprob,"w"))==NULL) {
printf("Problem with resultfile: %s\n", fileresprob);	printf("Problem with resultfile: %s\n", fileresprob);
fprintf(ficlog,"Problem with resultfile: %s\n", fileresprob);	fprintf(ficlog,"Problem with resultfile: %s\n", fileresprob);
}	}
strcpy(fileresprobcov,"probcov");	strcpy(fileresprobcov,"PROBCOV_");
strcat(fileresprobcov,fileres);	strcat(fileresprobcov,fileres);
if((ficresprobcov=fopen(fileresprobcov,"w"))==NULL) {	if((ficresprobcov=fopen(fileresprobcov,"w"))==NULL) {
printf("Problem with resultfile: %s\n", fileresprobcov);	printf("Problem with resultfile: %s\n", fileresprobcov);
fprintf(ficlog,"Problem with resultfile: %s\n", fileresprobcov);	fprintf(ficlog,"Problem with resultfile: %s\n", fileresprobcov);
}	}
strcpy(fileresprobcor,"probcor");	strcpy(fileresprobcor,"PROBCOR_");
strcat(fileresprobcor,fileres);	strcat(fileresprobcor,fileres);
if((ficresprobcor=fopen(fileresprobcor,"w"))==NULL) {	if((ficresprobcor=fopen(fileresprobcor,"w"))==NULL) {
printf("Problem with resultfile: %s\n", fileresprobcor);	printf("Problem with resultfile: %s\n", fileresprobcor);
Line 4235 void varprob(char optionfilefiname[], do	Line 4261 void varprob(char optionfilefiname[], do
fprintf(fichtm,"\n<li><h4> Computing and drawing one step probabilities with their confidence intervals</h4></li>\n");	fprintf(fichtm,"\n<li><h4> Computing and drawing one step probabilities with their confidence intervals</h4></li>\n");
fprintf(fichtm,"\n");	fprintf(fichtm,"\n");

fprintf(fichtm,"\n<li><h4> <a href=\"%s\">Matrix of variance-covariance of one-step probabilities (drawings)</a></h4></br>this page is important in order to visualize confidence intervals and especially correlation between disability and recovery</li>\n",optionfilehtmcov);	fprintf(fichtm,"\n<li><h4> <a href=\"%s\">Matrix of variance-covariance of one-step probabilities (drawings)</a></h4> this page is important in order to visualize confidence intervals and especially correlation between disability and recovery, or more generally, way in and way back.</li>\n",optionfilehtmcov);
fprintf(fichtmcov,"Current page is file <a href=\"%s\">%s</a><br>\n\n<h4>Matrix of variance-covariance of pairs of step probabilities</h4>\n",optionfilehtmcov, optionfilehtmcov);	fprintf(fichtmcov,"Current page is file <a href=\"%s\">%s</a><br>\n\n<h4>Matrix of variance-covariance of pairs of step probabilities</h4>\n",optionfilehtmcov, optionfilehtmcov);
fprintf(fichtmcov,"\nEllipsoids of confidence centered on point (p<inf>ij</inf>, p<inf>kl</inf>) are estimated \	fprintf(fichtmcov,"\nEllipsoids of confidence centered on point (p<inf>ij</inf>, p<inf>kl</inf>) are estimated \
and drawn. It helps understanding how is the covariance between two incidences.\	and drawn. It helps understanding how is the covariance between two incidences.\
Line 4258 To be simple, these graphs help to under	Line 4284 To be simple, these graphs help to under
/j1++;/	/j1++;/
if (cptcovn>0) {	if (cptcovn>0) {
fprintf(ficresprob, "\n#********** Variable ");	fprintf(ficresprob, "\n#********** Variable ");
for (z1=1; z1<=cptcoveff; z1++) fprintf(ficresprob, "V%d=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtab[j1][z1]]);	for (z1=1; z1<=cptcoveff; z1++) fprintf(ficresprob, "V%d=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtabm(j1,z1)]);
fprintf(ficresprob, "**********\n#\n");	fprintf(ficresprob, "**********\n#\n");
fprintf(ficresprobcov, "\n#********** Variable ");	fprintf(ficresprobcov, "\n#********** Variable ");
for (z1=1; z1<=cptcoveff; z1++) fprintf(ficresprobcov, "V%d=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtab[j1][z1]]);	for (z1=1; z1<=cptcoveff; z1++) fprintf(ficresprobcov, "V%d=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtabm(j1,z1)]);
fprintf(ficresprobcov, "**********\n#\n");	fprintf(ficresprobcov, "**********\n#\n");

fprintf(ficgp, "\n#********** Variable ");	fprintf(ficgp, "\n#********** Variable ");
for (z1=1; z1<=cptcoveff; z1++) fprintf(ficgp, " V%d=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtab[j1][z1]]);	for (z1=1; z1<=cptcoveff; z1++) fprintf(ficgp, " V%d=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtabm(j1,z1)]);
fprintf(ficgp, "**********\n#\n");	fprintf(ficgp, "**********\n#\n");


fprintf(fichtmcov, "\n<hr size=\"2\" color=\"#EC5E5E\">********** Variable ");	fprintf(fichtmcov, "\n<hr size=\"2\" color=\"#EC5E5E\">********** Variable ");
for (z1=1; z1<=cptcoveff; z1++) fprintf(fichtm, "V%d=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtab[j1][z1]]);	for (z1=1; z1<=cptcoveff; z1++) fprintf(fichtm, "V%d=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtabm(j1,z1)]);
fprintf(fichtmcov, "**********\n<hr size=\"2\" color=\"#EC5E5E\">");	fprintf(fichtmcov, "**********\n<hr size=\"2\" color=\"#EC5E5E\">");

fprintf(ficresprobcor, "\n#********** Variable ");	fprintf(ficresprobcor, "\n#********** Variable ");
for (z1=1; z1<=cptcoveff; z1++) fprintf(ficresprobcor, "V%d=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtab[j1][z1]]);	for (z1=1; z1<=cptcoveff; z1++) fprintf(ficresprobcor, "V%d=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtabm(j1,z1)]);
fprintf(ficresprobcor, "**********\n#");	fprintf(ficresprobcor, "**********\n#");
}	}

Line 4287 To be simple, these graphs help to under	Line 4313 To be simple, these graphs help to under
if(nagesqr==1)	if(nagesqr==1)
cov[3]= age*age;	cov[3]= age*age;
for (k=1; k<=cptcovn;k++) {	for (k=1; k<=cptcovn;k++) {
cov[2+nagesqr+k]=nbcode[Tvar[k]][codtab[j1][Tvar[k]]];/* j1 1 2 3 4	cov[2+nagesqr+k]=nbcode[Tvar[k]][codtabm(j1,k)];
	/cov[2+nagesqr+k]=nbcode[Tvar[k]][codtabm(j1,Tvar[k])];//* j1 1 2 3 4
* 1 1 1 1 1	* 1 1 1 1 1
* 2 2 1 1 1	* 2 2 1 1 1
* 3 1 2 1 1	* 3 1 2 1 1
Line 4295 To be simple, these graphs help to under	Line 4322 To be simple, these graphs help to under
/* nbcode[1][1]=0 nbcode[1][2]=1;*/	/* nbcode[1][1]=0 nbcode[1][2]=1;*/
}	}
/* for (k=1; k<=cptcovage;k++) cov[2+Tage[k]]=cov[2+Tage[k]]cov[2]; /	/* for (k=1; k<=cptcovage;k++) cov[2+Tage[k]]=cov[2+Tage[k]]cov[2]; /
for (k=1; k<=cptcovage;k++) cov[2+Tage[k]]=nbcode[Tvar[Tage[k]]][codtab[ij][Tvar[Tage[k]]]]*cov[2];	for (k=1; k<=cptcovage;k++) cov[2+Tage[k]]=nbcode[Tvar[Tage[k]]][codtabm(ij,k)]*cov[2];
for (k=1; k<=cptcovprod;k++)	for (k=1; k<=cptcovprod;k++)
cov[2+nagesqr+Tprod[k]]=nbcode[Tvard[k][1]][codtab[ij][Tvard[k][1]]]*nbcode[Tvard[k][2]][codtab[ij][Tvard[k][2]]];	cov[2+nagesqr+Tprod[k]]=nbcode[Tvard[k][1]][codtabm(ij,k)]*nbcode[Tvard[k][2]][codtabm(ij,k)];


for(theta=1; theta <=npar; theta++){	for(theta=1; theta <=npar; theta++){
Line 4445 To be simple, these graphs help to under	Line 4472 To be simple, these graphs help to under
/* mu2+ v21lc1cost + v22lc2sin(t) */	/* mu2+ v21lc1cost + v22lc2sin(t) */
if(first==1){	if(first==1){
first=0;	first=0;
	fprintf(ficgp,"\n# Ellipsoids of confidence\n#\n");
fprintf(ficgp,"\nset parametric;unset label");	fprintf(ficgp,"\nset parametric;unset label");
fprintf(ficgp,"\nset log y;set log x; set xlabel \"p%1d%1d (year-1)\";set ylabel \"p%1d%1d (year-1)\"",k1,l1,k2,l2);	fprintf(ficgp,"\nset log y;set log x; set xlabel \"p%1d%1d (year-1)\";set ylabel \"p%1d%1d (year-1)\"",k1,l1,k2,l2);
fprintf(ficgp,"\nset ter png small size 320, 240");	fprintf(ficgp,"\nset ter svg size 640, 480");
fprintf(fichtmcov,"\n<br>Ellipsoids of confidence cov(p%1d%1d,p%1d%1d) expressed in year<sup>-1</sup>\	fprintf(fichtmcov,"\n<br>Ellipsoids of confidence cov(p%1d%1d,p%1d%1d) expressed in year<sup>-1</sup>\
:<a href=\"%s%d%1d%1d-%1d%1d.png\">\	:<a href=\"%s_%d%1d%1d-%1d%1d.svg\">\
%s%d%1d%1d-%1d%1d.png</A>, ",k1,l1,k2,l2,\	%s_%d%1d%1d-%1d%1d.svg</A>, ",k1,l1,k2,l2,\
subdirf2(optionfilefiname,"varpijgr"), j1,k1,l1,k2,l2,\	subdirf2(optionfilefiname,"VARPIJGR_"), j1,k1,l1,k2,l2,\
subdirf2(optionfilefiname,"varpijgr"), j1,k1,l1,k2,l2);	subdirf2(optionfilefiname,"VARPIJGR_"), j1,k1,l1,k2,l2);
fprintf(fichtmcov,"\n<br><img src=\"%s%d%1d%1d-%1d%1d.png\"> ",subdirf2(optionfilefiname,"varpijgr"), j1,k1,l1,k2,l2);	fprintf(fichtmcov,"\n<br><img src=\"%s_%d%1d%1d-%1d%1d.svg\"> ",subdirf2(optionfilefiname,"VARPIJGR_"), j1,k1,l1,k2,l2);
fprintf(fichtmcov,"\n<br> Correlation at age %d (%.3f),",(int) age, c12);	fprintf(fichtmcov,"\n<br> Correlation at age %d (%.3f),",(int) age, c12);
fprintf(ficgp,"\nset out \"%s%d%1d%1d-%1d%1d.png\"",subdirf2(optionfilefiname,"varpijgr"), j1,k1,l1,k2,l2);	fprintf(ficgp,"\nset out \"%s_%d%1d%1d-%1d%1d.svg\"",subdirf2(optionfilefiname,"VARPIJGR_"), 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);
fprintf(ficgp,"\nplot [-pi:pi] %11.3e+ %.3f(%11.3e%11.3ecos(t)+%11.3e%11.3esin(t)), %11.3e +%.3f(%11.3e%11.3ecos(t)+%11.3e%11.3esin(t)) not",\	fprintf(ficgp,"\nplot [-pi:pi] %11.3e+ %.3f(%11.3e%11.3ecos(t)+%11.3e%11.3esin(t)), %11.3e +%.3f(%11.3e%11.3ecos(t)+%11.3e%11.3esin(t)) not",\
Line 4472 To be simple, these graphs help to under	Line 4500 To be simple, these graphs help to under
}/* if first */	}/* if first */
} /* age mod 5 */	} /* age mod 5 */
} /* end loop age */	} /* end loop age */
fprintf(ficgp,"\nset out \"%s%d%1d%1d-%1d%1d.png\";replot;",subdirf2(optionfilefiname,"varpijgr"), j1,k1,l1,k2,l2);	fprintf(ficgp,"\nset out;\nset out \"%s_%d%1d%1d-%1d%1d.svg\";replot;set out;",subdirf2(optionfilefiname,"VARPIJGR_"), j1,k1,l1,k2,l2);
first=1;	first=1;
} /l12 /	} /l12 /
} /* k12 */	} /* k12 */
Line 4494 To be simple, these graphs help to under	Line 4522 To be simple, these graphs help to under


/***************** Printing html file *********/	/***************** Printing html file *********/
void printinghtml(char fileres[], char title[], char datafile[], int firstpass, \	void printinghtml(char fileresu[], char title[], char datafile[], int firstpass, \
int lastpass, int stepm, int weightopt, char model[],\	int lastpass, int stepm, int weightopt, char model[],\
int imx,int jmin, int jmax, double jmeanint,char rfileres[],\	int imx,int jmin, int jmax, double jmeanint,char rfileres[],\
int popforecast, int estepm ,\	int popforecast, int estepm ,\
Line 4507 void printinghtml(char fileres[], char t	Line 4535 void printinghtml(char fileres[], char t
</ul>");	</ul>");
fprintf(fichtm,"<ul><li><h4><a name='firstorder'>Result files (first order: no variance)</a></h4>\n \	fprintf(fichtm,"<ul><li><h4><a name='firstorder'>Result files (first order: no variance)</a></h4>\n \
- Observed prevalence in each state (during the period defined between %.lf/%.lf/%.lf and %.lf/%.lf/%.lf): <a href=\"%s\">%s</a> <br>\n ",	- Observed prevalence in each state (during the period defined between %.lf/%.lf/%.lf and %.lf/%.lf/%.lf): <a href=\"%s\">%s</a> <br>\n ",
jprev1, mprev1,anprev1,jprev2, mprev2,anprev2,subdirf2(fileres,"p"),subdirf2(fileres,"p"));	jprev1, mprev1,anprev1,jprev2, mprev2,anprev2,subdirf2(fileresu,"P_"),subdirf2(fileresu,"P_"));
fprintf(fichtm,"\	fprintf(fichtm,"\
- Estimated transition probabilities over %d (stepm) months: <a href=\"%s\">%s</a><br>\n ",	- Estimated transition probabilities over %d (stepm) months: <a href=\"%s\">%s</a><br>\n ",
stepm,subdirf2(fileres,"pij"),subdirf2(fileres,"pij"));	stepm,subdirf2(fileresu,"PIJ_"),subdirf2(fileresu,"PIJ_"));
fprintf(fichtm,"\	fprintf(fichtm,"\
- Period (stable) prevalence in each health state: <a href=\"%s\">%s</a> <br>\n",	- Period (stable) prevalence in each health state: <a href=\"%s\">%s</a> <br>\n",
subdirf2(fileres,"pl"),subdirf2(fileres,"pl"));	subdirf2(fileresu,"PL_"),subdirf2(fileresu,"PL_"));
fprintf(fichtm,"\	fprintf(fichtm,"\
- (a) Life expectancies by health status at initial age, ei. (b) health expectancies by health status at initial age, eij . If one or more covariates are included, specific tables for each value of the covariate are output in sequences within the same file (estepm=%2d months): \	- (a) Life expectancies by health status at initial age, ei. (b) health expectancies by health status at initial age, eij . If one or more covariates are included, specific tables for each value of the covariate are output in sequences within the same file (estepm=%2d months): \
<a href=\"%s\">%s</a> <br>\n",	<a href=\"%s\">%s</a> <br>\n",
estepm,subdirf2(fileres,"e"),subdirf2(fileres,"e"));	estepm,subdirf2(fileresu,"E_"),subdirf2(fileresu,"E_"));
fprintf(fichtm,"\	fprintf(fichtm,"\
- Population projections by age and states: \	- Population projections by age and states: \
<a href=\"%s\">%s</a> <br>\n</li>", subdirf2(fileres,"f"),subdirf2(fileres,"f"));	<a href=\"%s\">%s</a> <br>\n</li>", subdirf2(fileresu,"F_"),subdirf2(fileresu,"F_"));

fprintf(fichtm," \n<ul><li><b>Graphs</b></li><p>");	fprintf(fichtm," \n<ul><li><b>Graphs</b></li><p>");

Line 4534 fprintf(fichtm," \n<ul><li><b>Graphs</b>	Line 4562 fprintf(fichtm," \n<ul><li><b>Graphs</b>
if (cptcovn > 0) {	if (cptcovn > 0) {
fprintf(fichtm,"<hr size=\"2\" color=\"#EC5E5E\">************ Results for covariates");	fprintf(fichtm,"<hr size=\"2\" color=\"#EC5E5E\">************ Results for covariates");
for (cpt=1; cpt<=cptcoveff;cpt++){	for (cpt=1; cpt<=cptcoveff;cpt++){
fprintf(fichtm," V%d=%d ",Tvaraff[cpt],nbcode[Tvaraff[cpt]][codtab[jj1][cpt]]);	fprintf(fichtm," V%d=%d ",Tvaraff[cpt],nbcode[Tvaraff[cpt]][codtabm(jj1,cpt)]);
printf(" V%d=%d ",Tvaraff[cpt],nbcode[Tvaraff[cpt]][codtab[jj1][cpt]]);fflush(stdout);	printf(" V%d=%d ",Tvaraff[cpt],nbcode[Tvaraff[cpt]][codtabm(jj1,cpt)]);fflush(stdout);
}	}
fprintf(fichtm," ************\n<hr size=\"2\" color=\"#EC5E5E\">");	fprintf(fichtm," ************\n<hr size=\"2\" color=\"#EC5E5E\">");
}	}
	/* aij, bij */
	fprintf(fichtm,"<br>- Logit model, for example: logit(pij)=log(pij/pii)= aij+ bij age + V1 age + etc. as a function of age: <a href=\"%s_%d-1.svg\">%s_%d-1.svg</a><br> \
	<img src=\"%s_%d-1.svg\">",subdirf2(optionfilefiname,"PE_"),jj1,subdirf2(optionfilefiname,"PE_"),jj1,subdirf2(optionfilefiname,"PE_"),jj1);
/* Pij */	/* Pij */
fprintf(fichtm,"<br>- Pij or Conditional probabilities to be observed in state j being in state i, %d (stepm) months before: <a href=\"%s%d_1.png\">%s%d_1.png</a><br> \	fprintf(fichtm,"<br>\n- Pij or Conditional probabilities to be observed in state j being in state i, %d (stepm) months before: <a href=\"%s_%d-2.svg\">%s_%d-2.svg</a><br> \
<img src=\"%s%d_1.png\">",stepm,subdirf2(optionfilefiname,"pe"),jj1,subdirf2(optionfilefiname,"pe"),jj1,subdirf2(optionfilefiname,"pe"),jj1);	<img src=\"%s_%d-2.svg\">",stepm,subdirf2(optionfilefiname,"PE_"),jj1,subdirf2(optionfilefiname,"PE_"),jj1,subdirf2(optionfilefiname,"PE_"),jj1);
/* Quasi-incidences */	/* Quasi-incidences */
fprintf(fichtm,"<br>- Pij or Conditional probabilities to be observed in state j being in state i %d (stepm) months\	fprintf(fichtm,"<br>\n- Iij 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: <a href=\"%s%d_2.png\">%s%d_2.png</a><br> \	before but expressed in per year i.e. quasi incidences if stepm is small and probabilities too,\
<img src=\"%s%d_2.png\">",stepm,subdirf2(optionfilefiname,"pe"),jj1,subdirf2(optionfilefiname,"pe"),jj1,subdirf2(optionfilefiname,"pe"),jj1);	incidence (rates) are the limit when h tends to zero of the ratio of the probability hPij \
/* Period (stable) prevalence in each health state */	divided by h: hPij/h : <a href=\"%s_%d-3.svg\">%s_%d-3.svg</a><br> \
for(cpt=1; cpt<=nlstate;cpt++){	<img src=\"%s_%d-3.svg\">",stepm,subdirf2(optionfilefiname,"PE_"),jj1,subdirf2(optionfilefiname,"PE_"),jj1,subdirf2(optionfilefiname,"PE_"),jj1);
fprintf(fichtm,"<br>- Convergence to period (stable) prevalence in state %d. Or probability to be in state %d being in state (1 to %d) at different ages. <a href=\"%s%d_%d.png\">%s%d_%d.png</a><br> \	/* Survival functions (period) in state j */
<img src=\"%s%d_%d.png\">", cpt, cpt, nlstate, subdirf2(optionfilefiname,"p"),cpt,jj1,subdirf2(optionfilefiname,"p"),cpt,jj1,subdirf2(optionfilefiname,"p"),cpt,jj1);	for(cpt=1; cpt<=nlstate;cpt++){
}	fprintf(fichtm,"<br>\n- Survival functions in state %d. Or probability to survive in state %d being in state (1 to %d) at different ages. <a href=\"%s%d_%d.svg\">%s%d_%d.svg</a><br> \
	<img src=\"%s_%d-%d.svg\">", cpt, cpt, nlstate, subdirf2(optionfilefiname,"LIJ_"),cpt,jj1,subdirf2(optionfilefiname,"LIJ_"),cpt,jj1,subdirf2(optionfilefiname,"LIJ_"),cpt,jj1);
	}
	/* State specific survival functions (period) */
	for(cpt=1; cpt<=nlstate;cpt++){
	fprintf(fichtm,"<br>\n- Survival functions from state %d in any different live states and total.\
	Or probability to survive in various states (1 to %d) being in state %d at different ages.\
	<a href=\"%s%d_%d.svg\">%s%d_%d.svg</a><br> <img src=\"%s_%d-%d.svg\">", cpt, nlstate, cpt, subdirf2(optionfilefiname,"LIJT_"),cpt,jj1,subdirf2(optionfilefiname,"LIJT_"),cpt,jj1,subdirf2(optionfilefiname,"LIJT_"),cpt,jj1);
	}
	/* Period (stable) prevalence in each health state */
	for(cpt=1; cpt<=nlstate;cpt++){
	fprintf(fichtm,"<br>\n- Convergence to period (stable) prevalence in state %d. Or probability to be in state %d being in state (1 to %d) at different ages. <a href=\"%s%d_%d.svg\">%s%d_%d.svg</a><br> \
	<img src=\"%s_%d-%d.svg\">", cpt, cpt, nlstate, subdirf2(optionfilefiname,"P_"),cpt,jj1,subdirf2(optionfilefiname,"P_"),cpt,jj1,subdirf2(optionfilefiname,"P_"),cpt,jj1);
	}
for(cpt=1; cpt<=nlstate;cpt++) {	for(cpt=1; cpt<=nlstate;cpt++) {
fprintf(fichtm,"\n<br>- Life expectancy by health state (%d) at initial age and its decomposition into health expectancies in each alive state (1 to %d) : <a href=\"%s%d%d.png\">%s%d%d.png</a> <br> \	fprintf(fichtm,"\n<br>- Life expectancy by health state (%d) at initial age and its decomposition into health expectancies in each alive state (1 to %d) (or area under each survival functions): <a href=\"%s%d%d.svg\">%s%d%d.svg</a> <br> \
<img src=\"%s%d%d.png\">",cpt,nlstate,subdirf2(optionfilefiname,"exp"),cpt,jj1,subdirf2(optionfilefiname,"exp"),cpt,jj1,subdirf2(optionfilefiname,"exp"),cpt,jj1);	<img src=\"%s_%d%d.svg\">",cpt,nlstate,subdirf2(optionfilefiname,"EXP_"),cpt,jj1,subdirf2(optionfilefiname,"EXP_"),cpt,jj1,subdirf2(optionfilefiname,"EXP_"),cpt,jj1);
}	}
/* } /\* end i1 \/ /	/* } /\* end i1 \/ /
}/* End k1 */	}/* End k1 */
Line 4572 covariance matrix of the one-step probab	Line 4616 covariance matrix of the one-step probab
See page 'Matrix of variance-covariance of one-step probabilities' below. \n", rfileres,rfileres);	See page 'Matrix of variance-covariance of one-step probabilities' below. \n", rfileres,rfileres);

fprintf(fichtm," - Standard deviation of one-step probabilities: <a href=\"%s\">%s</a> <br>\n",	fprintf(fichtm," - Standard deviation of one-step probabilities: <a href=\"%s\">%s</a> <br>\n",
subdirf2(fileres,"prob"),subdirf2(fileres,"prob"));	subdirf2(fileresu,"PROB_"),subdirf2(fileresu,"PROB_"));
fprintf(fichtm,"\	fprintf(fichtm,"\
- Variance-covariance of one-step probabilities: <a href=\"%s\">%s</a> <br>\n",	- Variance-covariance of one-step probabilities: <a href=\"%s\">%s</a> <br>\n",
subdirf2(fileres,"probcov"),subdirf2(fileres,"probcov"));	subdirf2(fileresu,"PROBCOV_"),subdirf2(fileresu,"PROBCOV_"));

fprintf(fichtm,"\	fprintf(fichtm,"\
- Correlation matrix of one-step probabilities: <a href=\"%s\">%s</a> <br>\n",	- Correlation matrix of one-step probabilities: <a href=\"%s\">%s</a> <br>\n",
subdirf2(fileres,"probcor"),subdirf2(fileres,"probcor"));	subdirf2(fileresu,"PROBCOR_"),subdirf2(fileresu,"PROBCOR_"));
fprintf(fichtm,"\	fprintf(fichtm,"\
- Variances and covariances of health expectancies by age and <b>initial health status</b> (cov(e<sup>ij</sup>,e<sup>kl</sup>)(estepm=%2d months): \	- Variances and covariances of health expectancies by age and <b>initial health status</b> (cov(e<sup>ij</sup>,e<sup>kl</sup>)(estepm=%2d months): \
<a href=\"%s\">%s</a> <br>\n</li>",	<a href=\"%s\">%s</a> <br>\n</li>",
estepm,subdirf2(fileres,"cve"),subdirf2(fileres,"cve"));	estepm,subdirf2(fileresu,"CVE_"),subdirf2(fileresu,"CVE_"));
fprintf(fichtm,"\	fprintf(fichtm,"\
- (a) Health expectancies by health status at initial age (e<sup>ij</sup>) and standard errors (in parentheses) (b) life expectancies and standard errors (e<sup>i.</sup>=e<sup>i1</sup>+e<sup>i2</sup>+...)(estepm=%2d months): \	- (a) Health expectancies by health status at initial age (e<sup>ij</sup>) and standard errors (in parentheses) (b) life expectancies and standard errors (e<sup>i.</sup>=e<sup>i1</sup>+e<sup>i2</sup>+...)(estepm=%2d months): \
<a href=\"%s\">%s</a> <br>\n</li>",	<a href=\"%s\">%s</a> <br>\n</li>",
estepm,subdirf2(fileres,"stde"),subdirf2(fileres,"stde"));	estepm,subdirf2(fileresu,"STDE_"),subdirf2(fileresu,"STDE_"));
fprintf(fichtm,"\	fprintf(fichtm,"\
- Variances and covariances of health expectancies by age. Status (i) based health expectancies (in state j), e<sup>ij</sup> are weighted by the period prevalences in each state i (if popbased=1, an additional computation is done using the cross-sectional prevalences, i.e population based) (estepm=%d months): <a href=\"%s\">%s</a><br>\n",	- Variances and covariances of health expectancies by age. Status (i) based health expectancies (in state j), e<sup>ij</sup> are weighted by the period prevalences in each state i (if popbased=1, an additional computation is done using the cross-sectional prevalences, i.e population based) (estepm=%d months): <a href=\"%s\">%s</a><br>\n",
estepm, subdirf2(fileres,"v"),subdirf2(fileres,"v"));	estepm, subdirf2(fileresu,"V_"),subdirf2(fileresu,"V_"));
fprintf(fichtm,"\	fprintf(fichtm,"\
- Total life expectancy and total health expectancies to be spent in each health state e<sup>.j</sup> with their standard errors (if popbased=1, an additional computation is done using the cross-sectional prevalences, i.e population based) (estepm=%d months): <a href=\"%s\">%s</a> <br>\n",	- Total life expectancy and total health expectancies to be spent in each health state e<sup>.j</sup> with their standard errors (if popbased=1, an additional computation is done using the cross-sectional prevalences, i.e population based) (estepm=%d months): <a href=\"%s\">%s</a> <br>\n",
estepm, subdirf2(fileres,"t"),subdirf2(fileres,"t"));	estepm, subdirf2(fileresu,"T_"),subdirf2(fileresu,"T_"));
fprintf(fichtm,"\	fprintf(fichtm,"\
- Standard deviation of period (stable) prevalences: <a href=\"%s\">%s</a> <br>\n",\	- Standard deviation of period (stable) prevalences: <a href=\"%s\">%s</a> <br>\n",\
subdirf2(fileres,"vpl"),subdirf2(fileres,"vpl"));	subdirf2(fileresu,"VPL_"),subdirf2(fileresu,"VPL_"));

/* if(popforecast==1) fprintf(fichtm,"\n */	/* if(popforecast==1) fprintf(fichtm,"\n */
/* - Prevalences forecasting: <a href=\"f%s\">f%s</a> <br>\n */	/* - Prevalences forecasting: <a href=\"f%s\">f%s</a> <br>\n */
Line 4617 See page 'Matrix of variance-covariance	Line 4661 See page 'Matrix of variance-covariance
if (cptcovn > 0) {	if (cptcovn > 0) {
fprintf(fichtm,"<hr size=\"2\" color=\"#EC5E5E\">************ Results for covariates");	fprintf(fichtm,"<hr size=\"2\" color=\"#EC5E5E\">************ Results for covariates");
for (cpt=1; cpt<=cptcoveff;cpt++)	for (cpt=1; cpt<=cptcoveff;cpt++)
fprintf(fichtm," V%d=%d ",Tvaraff[cpt],nbcode[Tvaraff[cpt]][codtab[jj1][cpt]]);	fprintf(fichtm," V%d=%d ",Tvaraff[cpt],nbcode[Tvaraff[cpt]][codtabm(jj1,cpt)]);
fprintf(fichtm," ************\n<hr size=\"2\" color=\"#EC5E5E\">");	fprintf(fichtm," ************\n<hr size=\"2\" color=\"#EC5E5E\">");
}	}
for(cpt=1; cpt<=nlstate;cpt++) {	for(cpt=1; cpt<=nlstate;cpt++) {
fprintf(fichtm,"<br>- Observed (cross-sectional) and period (incidence based) \	fprintf(fichtm,"<br>- Observed (cross-sectional) and period (incidence based) \
prevalence (with 95%% confidence interval) in state (%d): %s%d_%d.png <br>\	prevalence (with 95%% confidence interval) in state (%d): %s%d_%d.svg <br>\
<img src=\"%s%d_%d.png\">",cpt,subdirf2(optionfilefiname,"v"),cpt,jj1,subdirf2(optionfilefiname,"v"),cpt,jj1);	<img src=\"%s_%d-%d.svg\">",cpt,subdirf2(optionfilefiname,"V_"),cpt,jj1,subdirf2(optionfilefiname,"V_"),cpt,jj1);
}	}
fprintf(fichtm,"\n<br>- Total life expectancy by age and \	fprintf(fichtm,"\n<br>- Total life expectancy by age and \
health expectancies in states (1) and (2). If popbased=1 the smooth (due to the model) \	health expectancies in states (1) and (2). If popbased=1 the smooth (due to the model) \
true period expectancies (those weighted with period prevalences are also\	true period expectancies (those weighted with period prevalences are also\
drawn in addition to the population based expectancies computed using\	drawn in addition to the population based expectancies computed using\
observed and cahotic prevalences: %s%d.png<br>\	observed and cahotic prevalences: %s_%d.svg<br>\
<img src=\"%s%d.png\">",subdirf2(optionfilefiname,"e"),jj1,subdirf2(optionfilefiname,"e"),jj1);	<img src=\"%s_%d.svg\">",subdirf2(optionfilefiname,"E_"),jj1,subdirf2(optionfilefiname,"E_"),jj1);
/* } /\* end i1 \/ /	/* } /\* end i1 \/ /
}/* End k1 */	}/* End k1 */
fprintf(fichtm,"</ul>");	fprintf(fichtm,"</ul>");
Line 4638 true period expectancies (those weighted	Line 4682 true period expectancies (those weighted
}	}

/***************** Gnuplot file ************/	/***************** Gnuplot file ************/
void printinggnuplot(char fileres[], char optionfilefiname[], double ageminpar, double agemaxpar, double fage , char pathc[], double p[]){	void printinggnuplot(char fileresu[], char optionfilefiname[], double ageminpar, double agemaxpar, double fage , char pathc[], double p[]){

char dirfileres[132],optfileres[132];	char dirfileres[132],optfileres[132];
int cpt=0,k1=0,i=0,k=0,j=0,jk=0,k2=0,k3=0,ij=0,l=0;	int cpt=0,k1=0,i=0,k=0,j=0,jk=0,k2=0,k3=0,ij=0,l=0;
int ng=0;	int ng=0;
	int vpopbased;
/* if((ficgp=fopen(optionfilegnuplot,"a"))==NULL) { */	/* if((ficgp=fopen(optionfilegnuplot,"a"))==NULL) { */
/* printf("Problem with file %s",optionfilegnuplot); */	/* printf("Problem with file %s",optionfilegnuplot); */
/* fprintf(ficlog,"Problem with file %s",optionfilegnuplot); */	/* fprintf(ficlog,"Problem with file %s",optionfilegnuplot); */
Line 4656 void printinggnuplot(char fileres[], cha	Line 4701 void printinggnuplot(char fileres[], cha
strcpy(dirfileres,optionfilefiname);	strcpy(dirfileres,optionfilefiname);
strcpy(optfileres,"vpl");	strcpy(optfileres,"vpl");
/* 1eme*/	/* 1eme*/
fprintf(ficgp,"\n# 1st: Period (stable) prevalence with CI: 'vpl' files\n");	fprintf(ficgp,"\n# 1st: Period (stable) prevalence with CI: 'VPL_' files\n");
for (cpt=1; cpt<= nlstate ; cpt ++) {	for (cpt=1; cpt<= nlstate ; cpt ++) {
for (k1=1; k1<= m ; k1 ++) { /* plot [100000000000000000000:-100000000000000000000] "mysbiaspar/vplrmysbiaspar.txt to check */	for (k1=1; k1<= m ; k1 ++) { /* plot [100000000000000000000:-100000000000000000000] "mysbiaspar/vplrmysbiaspar.txt to check */
fprintf(ficgp,"\nset out \"%s%d_%d.png\" \n",subdirf2(optionfilefiname,"v"),cpt,k1);	fprintf(ficgp,"\nset out \"%s_%d-%d.svg\" \n",subdirf2(optionfilefiname,"V_"),cpt,k1);
fprintf(ficgp,"\n#set out \"v%s%d_%d.png\" \n",optionfilefiname,cpt,k1);	fprintf(ficgp,"\n#set out \"V_%s_%d-%d.svg\" \n",optionfilefiname,cpt,k1);
fprintf(ficgp,"set xlabel \"Age\" \n\	fprintf(ficgp,"set xlabel \"Age\" \n\
set ylabel \"Probability\" \n\	set ylabel \"Probability\" \n\
set ter png small size 320, 240\n\	set ter svg size 640, 480\n\
plot [%.f:%.f] \"%s\" every :::%d::%d u 1:2 \"%%lf",ageminpar,fage,subdirf2(fileres,"vpl"),k1-1,k1-1);	plot [%.f:%.f] \"%s\" every :::%d::%d u 1:2 \"%%lf",ageminpar,fage,subdirf2(fileresu,"VPL_"),k1-1,k1-1);

for (i=1; i<= nlstate ; i ++) {	for (i=1; i<= nlstate ; i ++) {
if (i==cpt) fprintf(ficgp," %%lf (%%lf)");	if (i==cpt) fprintf(ficgp," %%lf (%%lf)");
else fprintf(ficgp," %%lf (%%lf)");	else fprintf(ficgp," %%lf (%%lf)");
}	}
fprintf(ficgp,"\" t\"Period (stable) prevalence\" w l lt 0,\"%s\" every :::%d::%d u 1:($2+1.96*$3) \"%%lf",subdirf2(fileres,"vpl"),k1-1,k1-1);	fprintf(ficgp,"\" t\"Period (stable) prevalence\" w l lt 0,\"%s\" every :::%d::%d u 1:($2+1.96*$3) \"%%lf",subdirf2(fileresu,"VPL_"),k1-1,k1-1);
for (i=1; i<= nlstate ; i ++) {	for (i=1; i<= nlstate ; i ++) {
if (i==cpt) fprintf(ficgp," %%lf (%%lf)");	if (i==cpt) fprintf(ficgp," %%lf (%%lf)");
else fprintf(ficgp," %%lf (%%lf)");	else fprintf(ficgp," %%lf (%%lf)");
}	}
fprintf(ficgp,"\" t\"95%% CI\" w l lt 1,\"%s\" every :::%d::%d u 1:($2-1.96*$3) \"%%lf",subdirf2(fileres,"vpl"),k1-1,k1-1);	fprintf(ficgp,"\" t\"95%% CI\" w l lt 1,\"%s\" every :::%d::%d u 1:($2-1.96*$3) \"%%lf",subdirf2(fileresu,"VPL_"),k1-1,k1-1);
for (i=1; i<= nlstate ; i ++) {	for (i=1; i<= nlstate ; i ++) {
if (i==cpt) fprintf(ficgp," %%lf (%%lf)");	if (i==cpt) fprintf(ficgp," %%lf (%%lf)");
else fprintf(ficgp," %%lf (%%lf)");	else fprintf(ficgp," %%lf (%%lf)");
}	}
fprintf(ficgp,"\" t\"\" w l lt 1,\"%s\" every :::%d::%d u 1:($%d) t\"Observed prevalence \" w l lt 2",subdirf2(fileres,"p"),k1-1,k1-1,2+4*(cpt-1));	fprintf(ficgp,"\" t\"\" w l lt 1,\"%s\" every :::%d::%d u 1:($%d) t\"Observed prevalence \" w l lt 2",subdirf2(fileresu,"P_"),k1-1,k1-1,2+4*(cpt-1));
}	fprintf(ficgp,"\nset out \n");
}	} /* k1 */
	} /* cpt */
/2 eme/	/2 eme/
fprintf(ficgp,"\n# 2nd: Total life expectancy with CI: 't' files\n");	fprintf(ficgp,"\n# 2nd: Total life expectancy with CI: 't' files\n");
for (k1=1; k1<= m ; k1 ++) {	for (k1=1; k1<= m ; k1 ++) {
fprintf(ficgp,"\nset out \"%s%d.png\" \n",subdirf2(optionfilefiname,"e"),k1);	fprintf(ficgp,"\nset out \"%s_%d.svg\" \n",subdirf2(optionfilefiname,"E_"),k1);
fprintf(ficgp,"set ylabel \"Years\" \nset ter png small size 320, 240\nplot [%.f:%.f] ",ageminpar,fage);	for(vpopbased=0; vpopbased <= popbased; vpopbased++){ /* Done for vpopbased=0 and vpopbased=1 if popbased==1*/
	if(vpopbased==0)
for (i=1; i<= nlstate+1 ; i ++) {	fprintf(ficgp,"set ylabel \"Years\" \nset ter svg size 640, 480\nplot [%.f:%.f] ",ageminpar,fage);
k=2*i;	else
fprintf(ficgp,"\"%s\" every :::%d::%d u 1:2 \"%%lf",subdirf2(fileres,"t"),k1-1,k1-1);	fprintf(ficgp,"\nreplot ");
for (j=1; j<= nlstate+1 ; j ++) {	for (i=1; i<= nlstate+1 ; i ++) {
if (j==i) fprintf(ficgp," %%lf (%%lf)");	k=2*i;
else fprintf(ficgp," %%lf (%%lf)");	fprintf(ficgp,"\"%s\" every :::%d::%d u 1:($2==%d && $4!=0 ?$4 : 1/0) \"%%lf %%lf %%lf",subdirf2(fileresu,"T_"),k1-1,k1-1, vpopbased);
}	for (j=1; j<= nlstate+1 ; j ++) {
if (i== 1) fprintf(ficgp,"\" t\"TLE\" w l ,");	if (j==i) fprintf(ficgp," %%lf (%%lf)");
else fprintf(ficgp,"\" t\"LE in state (%d)\" w l ,",i-1);	else fprintf(ficgp," %%lf (%%lf)");
fprintf(ficgp,"\"%s\" every :::%d::%d u 1:($2-$3*2) \"%%lf",subdirf2(fileres,"t"),k1-1,k1-1);	}
for (j=1; j<= nlstate+1 ; j ++) {	if (i== 1) fprintf(ficgp,"\" t\"TLE\" w l lt %d, \\\n",i);
if (j==i) fprintf(ficgp," %%lf (%%lf)");	else fprintf(ficgp,"\" t\"LE in state (%d)\" w l lt %d, \\\n",i-1,i+1);
else fprintf(ficgp," %%lf (%%lf)");	fprintf(ficgp,"\"%s\" every :::%d::%d u 1:($2==%d && $4!=0 ? $4-$5*2 : 1/0) \"%%lf %%lf %%lf",subdirf2(fileresu,"T_"),k1-1,k1-1,vpopbased);
}	for (j=1; j<= nlstate+1 ; j ++) {
fprintf(ficgp,"\" t\"\" w l lt 0,");	if (j==i) fprintf(ficgp," %%lf (%%lf)");
fprintf(ficgp,"\"%s\" every :::%d::%d u 1:($2+$3*2) \"%%lf",subdirf2(fileres,"t"),k1-1,k1-1);	else fprintf(ficgp," %%lf (%%lf)");
for (j=1; j<= nlstate+1 ; j ++) {	}
if (j==i) fprintf(ficgp," %%lf (%%lf)");	fprintf(ficgp,"\" t\"\" w l lt 0,");
else fprintf(ficgp," %%lf (%%lf)");	fprintf(ficgp,"\"%s\" every :::%d::%d u 1:($2==%d && $4!=0 ? $4+$5*2 : 1/0) \"%%lf %%lf %%lf",subdirf2(fileresu,"T_"),k1-1,k1-1,vpopbased);
}	for (j=1; j<= nlstate+1 ; j ++) {
if (i== (nlstate+1)) fprintf(ficgp,"\" t\"\" w l lt 0");	if (j==i) fprintf(ficgp," %%lf (%%lf)");
else fprintf(ficgp,"\" t\"\" w l lt 0,");	else fprintf(ficgp," %%lf (%%lf)");
}	}
}	if (i== (nlstate+1)) fprintf(ficgp,"\" t\"\" w l lt 0");
	else fprintf(ficgp,"\" t\"\" w l lt 0,\\\n");
	} /* state */
	} /* vpopbased */
	fprintf(ficgp,"\nset out;set out \"%s_%d.svg\"; replot; set out; \n",subdirf2(optionfilefiname,"E_"),k1); /* Buggy gnuplot */
	} /* k1 */
/3eme/	/3eme/

for (k1=1; k1<= m ; k1 ++) {	for (k1=1; k1<= m ; k1 ++) {
for (cpt=1; cpt<= nlstate ; cpt ++) {	for (cpt=1; cpt<= nlstate ; cpt ++) {
/* k=2+nlstate(2cpt-2); */	/* k=2+nlstate(2cpt-2); */
k=2+(nlstate+1)*(cpt-1);	k=2+(nlstate+1)*(cpt-1);
fprintf(ficgp,"\nset out \"%s%d%d.png\" \n",subdirf2(optionfilefiname,"exp"),cpt,k1);	fprintf(ficgp,"\nset out \"%s_%d%d.svg\" \n",subdirf2(optionfilefiname,"EXP_"),cpt,k1);
fprintf(ficgp,"set ter png small size 320, 240\n\	fprintf(ficgp,"set ter svg size 640, 480\n\
plot [%.f:%.f] \"%s\" every :::%d::%d u 1:%d t \"e%d1\" w l",ageminpar,fage,subdirf2(fileres,"e"),k1-1,k1-1,k,cpt);	plot [%.f:%.f] \"%s\" every :::%d::%d u 1:%d t \"e%d1\" w l",ageminpar,fage,subdirf2(fileresu,"E_"),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);	/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) ");	for (i=1; i<= nlstate*2 ; i ++) fprintf(ficgp,"\%%lf (\%%lf) ");
fprintf(ficgp,"\" t \"e%d1\" w l",cpt);	fprintf(ficgp,"\" t \"e%d1\" w l",cpt);
Line 4732 plot [%.f:%.f] \"%s\" every :::%d::%d u	Line 4782 plot [%.f:%.f] \"%s\" every :::%d::%d u

*/	*/
for (i=1; i< nlstate ; i ++) {	for (i=1; i< nlstate ; i ++) {
fprintf(ficgp," ,\"%s\" every :::%d::%d u 1:%d t \"e%d%d\" w l",subdirf2(fileres,"e"),k1-1,k1-1,k+i,cpt,i+1);	fprintf(ficgp," ,\"%s\" every :::%d::%d u 1:%d t \"e%d%d\" w l",subdirf2(fileresu,"E_"),k1-1,k1-1,k+i,cpt,i+1);
/* fprintf(ficgp," ,\"%s\" every :::%d::%d u 1:%d t \"e%d%d\" w l",subdirf2(fileres,"e"),k1-1,k1-1,k+2i,cpt,i+1);/	/* fprintf(ficgp," ,\"%s\" every :::%d::%d u 1:%d t \"e%d%d\" w l",subdirf2(fileres,"e"),k1-1,k1-1,k+2i,cpt,i+1);/

}	}
fprintf(ficgp," ,\"%s\" every :::%d::%d u 1:%d t \"e%d.\" w l",subdirf2(fileres,"e"),k1-1,k1-1,k+nlstate,cpt);	fprintf(ficgp," ,\"%s\" every :::%d::%d u 1:%d t \"e%d.\" w l",subdirf2(fileresu,"E_"),k1-1,k1-1,k+nlstate,cpt);
}	}
}	}

	/* Survival functions (period) from state i in state j by initial state i */
	for (k1=1; k1<= m ; k1 ++) { /* For each multivariate if any */
	for (cpt=1; cpt<=nlstate ; cpt ++) { /* For each life state */
	k=3;
	fprintf(ficgp,"\n#\n#\n# Survival functions in state j : 'lij' files, cov=%d state=%d",k1, cpt);
	fprintf(ficgp,"\nset out \"%s_%d-%d.svg\" \n",subdirf2(optionfilefiname,"LIJ_"),cpt,k1);
	fprintf(ficgp,"set xlabel \"Age\" \nset ylabel \"Probability to be alive\" \n\
	set ter svg size 640, 480\n\
	unset log y\n\
	plot [%.f:%.f] ", ageminpar, agemaxpar);
	for (i=1; i<= nlstate ; i ++){
	if(i==1)
	fprintf(ficgp,"\"%s\"",subdirf2(fileresu,"PIJ_"));
	else
	fprintf(ficgp,", '' ");
	l=(nlstate+ndeath)*(i-1)+1;
	fprintf(ficgp," u ($1==%d ? ($3):1/0):($%d/($%d",k1,k+l+(cpt-1),k+l);
	for (j=2; j<= nlstate+ndeath ; j ++)
	fprintf(ficgp,"+$%d",k+l+j-1);
	fprintf(ficgp,")) t \"l(%d,%d)\" w l",i,cpt);
	} /* nlstate */
	fprintf(ficgp,"\nset out\n");
	} /* end cpt state*/
	} /* end covariate */

	/* Survival functions (period) from state i in state j by final state j */
	for (k1=1; k1<= m ; k1 ++) { /* For each multivariate if any */
	for (cpt=1; cpt<=nlstate ; cpt ++) { /* For each inital state */
	k=3;
	fprintf(ficgp,"\n#\n#\n# Survival functions in state j and all livestates from state i by final state j: 'lij' files, cov=%d state=%d",k1, cpt);
	fprintf(ficgp,"\nset out \"%s_%d-%d.svg\" \n",subdirf2(optionfilefiname,"LIJT_"),cpt,k1);
	fprintf(ficgp,"set xlabel \"Age\" \nset ylabel \"Probability to be alive\" \n\
	set ter svg size 640, 480\n\
	unset log y\n\
	plot [%.f:%.f] ", ageminpar, agemaxpar);
	for (j=1; j<= nlstate ; j ++){ /* Lived in state j */
	if(j==1)
	fprintf(ficgp,"\"%s\"",subdirf2(fileresu,"PIJ_"));
	else
	fprintf(ficgp,", '' ");
	l=(nlstate+ndeath)*(cpt-1) +j;
	fprintf(ficgp," u (($1==%d && (floor($2)%%5 == 0)) ? ($3):1/0):($%d",k1,k+l);
	/* for (i=2; i<= nlstate+ndeath ; i ++) */
	/* fprintf(ficgp,"+$%d",k+l+i-1); */
	fprintf(ficgp,") t \"l(%d,%d)\" w l",cpt,j);
	} /* nlstate */
	fprintf(ficgp,", '' ");
	fprintf(ficgp," u (($1==%d && (floor($2)%%5 == 0)) ? ($3):1/0):(",k1);
	for (j=1; j<= nlstate ; j ++){ /* Lived in state j */
	l=(nlstate+ndeath)*(cpt-1) +j;
	if(j < nlstate)
	fprintf(ficgp,"$%d +",k+l);
	else
	fprintf(ficgp,"$%d) t\"l(%d,.)\" w l",k+l,cpt);
	}
	fprintf(ficgp,"\nset out\n");
	} /* end cpt state*/
	} /* end covariate */

/* CV preval stable (period) */	/* CV preval stable (period) */
for (k1=1; k1<= m ; k1 ++) { /* For each multivariate if any */	for (k1=1; k1<= m ; k1 ++) { /* For each multivariate if any */
for (cpt=1; cpt<=nlstate ; cpt ++) { /* For each life state */	for (cpt=1; cpt<=nlstate ; cpt ++) { /* For each life state */
k=3;	k=3;
fprintf(ficgp,"\n#\n#\n#CV preval stable (period): 'pij' files, cov=%d state=%d",k1, cpt);	fprintf(ficgp,"\n#\n#\n#CV preval stable (period): 'pij' files, cov=%d state=%d",k1, cpt);
fprintf(ficgp,"\nset out \"%s%d_%d.png\" \n",subdirf2(optionfilefiname,"p"),cpt,k1);	fprintf(ficgp,"\nset out \"%s_%d-%d.svg\" \n",subdirf2(optionfilefiname,"P_"),cpt,k1);
fprintf(ficgp,"set xlabel \"Age\" \nset ylabel \"Probability\" \n\	fprintf(ficgp,"set xlabel \"Age\" \nset ylabel \"Probability\" \n\
set ter png small size 320, 240\n\	set ter svg size 640, 480\n\
unset log y\n\	unset log y\n\
plot [%.f:%.f] ", ageminpar, agemaxpar);	plot [%.f:%.f] ", ageminpar, agemaxpar);
for (i=1; i<= nlstate ; i ++){	for (i=1; i<= nlstate ; i ++){
if(i==1)	if(i==1)
fprintf(ficgp,"\"%s\"",subdirf2(fileres,"pij"));	fprintf(ficgp,"\"%s\"",subdirf2(fileresu,"PIJ_"));
else	else
fprintf(ficgp,", '' ");	fprintf(ficgp,", '' ");
l=(nlstate+ndeath)*(i-1)+1;	l=(nlstate+ndeath)*(i-1)+1;
fprintf(ficgp," u ($1==%d ? ($3):1/0):($%d/($%d",k1,k+l+(cpt-1),k+l);	fprintf(ficgp," u ($1==%d ? ($3):1/0):($%d/($%d",k1,k+l+(cpt-1),k+l);
for (j=1; j<= (nlstate-1) ; j ++)	for (j=2; j<= nlstate ; j ++)
fprintf(ficgp,"+$%d",k+l+j);	fprintf(ficgp,"+$%d",k+l+j-1);
fprintf(ficgp,")) t \"prev(%d,%d)\" w l",i,cpt);	fprintf(ficgp,")) t \"prev(%d,%d)\" w l",i,cpt);
} /* nlstate */	} /* nlstate */
fprintf(ficgp,"\n");	fprintf(ficgp,"\nset out\n");
} /* end cpt state*/	} /* end cpt state*/
} /* end covariate */	} /* end covariate */

/* proba elementaires */	/* proba elementaires */
fprintf(ficgp,"\n##############\n#MLE estimated parameters\n#############\n");	fprintf(ficgp,"\n##############\n#MLE estimated parameters\n#############\n");
for(i=1,jk=1; i <=nlstate; i++){	for(i=1,jk=1; i <=nlstate; i++){
Line 4783 plot [%.f:%.f] ", ageminpar, agemaxpar)	Line 4892 plot [%.f:%.f] ", ageminpar, agemaxpar)
fprintf(ficgp,"##############\n#\n");	fprintf(ficgp,"##############\n#\n");

/goto avoid;/	/goto avoid;/
fprintf(ficgp,"\n##############\n#Graphics of of probabilities or incidences\n#############\n");	fprintf(ficgp,"\n##############\n#Graphics of probabilities or incidences\n#############\n");
fprintf(ficgp,"# logi(p12/p11)=a12+b12age+c12ageage+d12V1+e12V1*age\n");	fprintf(ficgp,"# logi(p12/p11)=a12+b12age+c12ageage+d12V1+e12V1*age\n");
fprintf(ficgp,"# logi(p12/p11)=p1 +p2age +p3ageage+ p4V1+ p5V1age\n");	fprintf(ficgp,"# logi(p12/p11)=p1 +p2age +p3ageage+ p4V1+ p5V1age\n");
fprintf(ficgp,"# logi(p13/p11)=a13+b13age+c13ageage+d13V1+e13V1*age\n");	fprintf(ficgp,"# logi(p13/p11)=a13+b13age+c13ageage+d13V1+e13V1*age\n");
Line 4797 plot [%.f:%.f] ", ageminpar, agemaxpar)	Line 4906 plot [%.f:%.f] ", ageminpar, agemaxpar)
fprintf(ficgp,"# +exp(a13+b13age+c13ageage+d13V1+e13V1*age))\n");	fprintf(ficgp,"# +exp(a13+b13age+c13ageage+d13V1+e13V1*age))\n");
fprintf(ficgp,"# +exp(a14+b14age+c14ageage+d14V1+e14V1*age)+...)\n");	fprintf(ficgp,"# +exp(a14+b14age+c14ageage+d14V1+e14V1*age)+...)\n");
fprintf(ficgp,"#\n");	fprintf(ficgp,"#\n");
for(ng=1; ng<=2;ng++){ /* Number of graphics: first is probabilities second is incidence per year*/	for(ng=1; ng<=3;ng++){ /* Number of graphics: first is logit, 2nd is probabilities, third is incidences per year*/
fprintf(ficgp,"# ng=%d\n",ng);	fprintf(ficgp,"# ng=%d\n",ng);
fprintf(ficgp,"# jk=1 to 2^%d=%d\n",cptcoveff,m);	fprintf(ficgp,"# jk=1 to 2^%d=%d\n",cptcoveff,m);
for(jk=1; jk <=m; jk++) {	for(jk=1; jk <=m; jk++) {
fprintf(ficgp,"# jk=%d\n",jk);	fprintf(ficgp,"# jk=%d\n",jk);
fprintf(ficgp,"\nset out \"%s%d_%d.png\" \n",subdirf2(optionfilefiname,"pe"),jk,ng);	fprintf(ficgp,"\nset out \"%s_%d-%d.svg\" ",subdirf2(optionfilefiname,"PE_"),jk,ng);
if (ng==2)	fprintf(ficgp,"\nset ter svg size 640, 480 ");
	if (ng==1){
	fprintf(ficgp,"\nset ylabel \"Value of the logit of the model\"\n"); /* exp(a12+b12x) could be nice /
	fprintf(ficgp,"\nunset log y");
	}else if (ng==2){
	fprintf(ficgp,"\nset ylabel \"Probability\"\n");
	fprintf(ficgp,"\nset log y");
	}else if (ng==3){
fprintf(ficgp,"\nset ylabel \"Quasi-incidence per year\"\n");	fprintf(ficgp,"\nset ylabel \"Quasi-incidence per year\"\n");
else	fprintf(ficgp,"\nset log y");
fprintf(ficgp,"\nset title \"Probability\"\n");	}else
fprintf(ficgp,"\nset ter png small size 320, 240\nset log y\nplot [%.f:%.f] ",ageminpar,agemaxpar);	fprintf(ficgp,"\nunset title ");
	fprintf(ficgp,"\nplot [%.f:%.f] ",ageminpar,agemaxpar);
i=1;	i=1;
for(k2=1; k2<=nlstate; k2++) {	for(k2=1; k2<=nlstate; k2++) {
k3=i;	k3=i;
for(k=1; k<=(nlstate+ndeath); k++) {	for(k=1; k<=(nlstate+ndeath); k++) {
if (k != k2){	if (k != k2){
if(ng==2)	switch( ng) {
	case 1:
if(nagesqr==0)	if(nagesqr==0)
fprintf(ficgp," %fexp(p%d+p%dx",YEARM/stepm,i,i+1);	fprintf(ficgp," p%d+p%d*x",i,i+1);
else /* nagesqr =1 */	else /* nagesqr =1 */
fprintf(ficgp," %fexp(p%d+p%dx+p%dxx",YEARM/stepm,i,i+1,i+1+nagesqr);	fprintf(ficgp," p%d+p%dx+p%dx*x",i,i+1,i+1+nagesqr);
else	break;
	case 2: /* ng=2 */
if(nagesqr==0)	if(nagesqr==0)
fprintf(ficgp," exp(p%d+p%d*x",i,i+1);	fprintf(ficgp," exp(p%d+p%d*x",i,i+1);
else /* nagesqr =1 */	else /* nagesqr =1 */
fprintf(ficgp," exp(p%d+p%dx+p%dx*x",i,i+1,i+1+nagesqr);	fprintf(ficgp," exp(p%d+p%dx+p%dx*x",i,i+1,i+1+nagesqr);
	break;
	case 3:
	if(nagesqr==0)
	fprintf(ficgp," %fexp(p%d+p%dx",YEARM/stepm,i,i+1);
	else /* nagesqr =1 */
	fprintf(ficgp," %fexp(p%d+p%dx+p%dxx",YEARM/stepm,i,i+1,i+1+nagesqr);
	break;
	}
ij=1;/* To be checked else nbcode[0][0] wrong */	ij=1;/* To be checked else nbcode[0][0] wrong */
for(j=3; j <=ncovmodel-nagesqr; j++) {	for(j=3; j <=ncovmodel-nagesqr; j++) {
/* printf("Tage[%d]=%d, j=%d\n", ij, Tage[ij], j); */	/* printf("Tage[%d]=%d, j=%d\n", ij, Tage[ij], j); */
if(ij <=cptcovage) { /* Bug valgrind */	if(ij <=cptcovage) { /* Bug valgrind */
if((j-2)==Tage[ij]) { /* Bug valgrind */	if((j-2)==Tage[ij]) { /* Bug valgrind */
fprintf(ficgp,"+p%d%dx",i+j+nagesqr-1,nbcode[Tvar[j-2]][codtab[jk][Tvar[j-2]]]);	fprintf(ficgp,"+p%d%dx",i+j+nagesqr-1,nbcode[Tvar[j-2]][codtabm(jk,j-2)]);
	/* fprintf(ficgp,"+p%d%dx",i+j+nagesqr-1,nbcode[Tvar[j-2]][codtabm(jk,Tvar[j-2])]); */
ij++;	ij++;
}	}
}	}
else	else
fprintf(ficgp,"+p%d*%d",i+j+nagesqr-1,nbcode[Tvar[j-2]][codtab[jk][j-2]]);	fprintf(ficgp,"+p%d*%d",i+j+nagesqr-1,nbcode[Tvar[j-2]][codtabm(jk,j-2)]);
}	}
fprintf(ficgp,")/(1");	if(ng != 1){
	fprintf(ficgp,")/(1");

for(k1=1; k1 <=nlstate; k1++){	for(k1=1; k1 <=nlstate; k1++){
if(nagesqr==0)	if(nagesqr==0)
fprintf(ficgp,"+exp(p%d+p%dx",k3+(k1-1)ncovmodel,k3+(k1-1)*ncovmodel+1);	fprintf(ficgp,"+exp(p%d+p%dx",k3+(k1-1)ncovmodel,k3+(k1-1)*ncovmodel+1);
else /* nagesqr =1 */	else /* nagesqr =1 */
fprintf(ficgp,"+exp(p%d+p%dx+p%dxx",k3+(k1-1)ncovmodel,k3+(k1-1)ncovmodel+1,k3+(k1-1)ncovmodel+1+nagesqr);	fprintf(ficgp,"+exp(p%d+p%dx+p%dxx",k3+(k1-1)ncovmodel,k3+(k1-1)ncovmodel+1,k3+(k1-1)ncovmodel+1+nagesqr);

ij=1;	ij=1;
for(j=3; j <=ncovmodel-nagesqr; j++){	for(j=3; j <=ncovmodel-nagesqr; j++){
if(ij <=cptcovage) { /* Bug valgrind */	if(ij <=cptcovage) { /* Bug valgrind */
if((j-2)==Tage[ij]) { /* Bug valgrind */	if((j-2)==Tage[ij]) { /* Bug valgrind */
fprintf(ficgp,"+p%d%dx",k3+(k1-1)*ncovmodel+1+j-2+nagesqr,nbcode[Tvar[j-2]][codtab[jk][Tvar[j-2]]]);	fprintf(ficgp,"+p%d%dx",k3+(k1-1)*ncovmodel+1+j-2+nagesqr,nbcode[Tvar[j-2]][codtabm(jk,j-2)]);
ij++;	/* fprintf(ficgp,"+p%d%dx",k3+(k1-1)ncovmodel+1+j-2+nagesqr,nbcode[Tvar[j-2]][codtabm(jk,Tvar[j-2])]); /
	ij++;
	}
}	}
	else
	fprintf(ficgp,"+p%d%d",k3+(k1-1)ncovmodel+1+j-2+nagesqr,nbcode[Tvar[j-2]][codtabm(jk,j-2)]);
}	}
else	fprintf(ficgp,")");
fprintf(ficgp,"+p%d%d",k3+(k1-1)ncovmodel+1+j-2+nagesqr,nbcode[Tvar[j-2]][codtab[jk][j-2]]);
}	}
fprintf(ficgp,")");	fprintf(ficgp,")");
	if(ng ==2)
	fprintf(ficgp," t \"p%d%d\" ", k2,k);
	else /* ng= 3 */
	fprintf(ficgp," t \"i%d%d\" ", k2,k);
	}else{ /* end ng <> 1 */
	fprintf(ficgp," t \"logit(p%d%d)\" ", k2,k);
}	}
fprintf(ficgp,") t \"p%d%d\" ", k2,k);
if ((k+k2)!= (nlstate*2+ndeath)) fprintf(ficgp,",");	if ((k+k2)!= (nlstate*2+ndeath)) fprintf(ficgp,",");
i=i+ncovmodel;	i=i+ncovmodel;
}	}
} /* end k */	} /* end k */
} /* end k2 */	} /* end k2 */
	fprintf(ficgp,"\n set out\n");
} /* end jk */	} /* end jk */
} /* end ng */	} /* end ng */
/* avoid: */	/* avoid: */
Line 4929 void prevforecast(char fileres[], double	Line 5067 void prevforecast(char fileres[], double
agelim=AGESUP;	agelim=AGESUP;
prevalence(probs, ageminpar, agemax, s, agev, nlstate, imx, Tvar, nbcode, ncodemax, mint, anint, dateprev1, dateprev2, firstpass, lastpass);	prevalence(probs, ageminpar, agemax, s, agev, nlstate, imx, Tvar, nbcode, ncodemax, mint, anint, dateprev1, dateprev2, firstpass, lastpass);

strcpy(fileresf,"f");	strcpy(fileresf,"F_");
strcat(fileresf,fileres);	strcat(fileresf,fileresu);
if((ficresf=fopen(fileresf,"w"))==NULL) {	if((ficresf=fopen(fileresf,"w"))==NULL) {
printf("Problem with forecast resultfile: %s\n", fileresf);	printf("Problem with forecast resultfile: %s\n", fileresf);
fprintf(ficlog,"Problem with forecast resultfile: %s\n", fileresf);	fprintf(ficlog,"Problem with forecast resultfile: %s\n", fileresf);
Line 4979 void prevforecast(char fileres[], double	Line 5117 void prevforecast(char fileres[], double
k=k+1;	k=k+1;
fprintf(ficresf,"\n#******");	fprintf(ficresf,"\n#******");
for(j=1;j<=cptcoveff;j++) {	for(j=1;j<=cptcoveff;j++) {
fprintf(ficresf," V%d=%d, hpijx=probability over h years, hp.jx is weighted by observed prev ",Tvaraff[j],nbcode[Tvaraff[j]][codtab[k][j]]);	fprintf(ficresf," V%d=%d, hpijx=probability over h years, hp.jx is weighted by observed prev ",Tvaraff[j],nbcode[Tvaraff[j]][codtabm(k,j)]);
}	}
fprintf(ficresf,"******\n");	fprintf(ficresf,"******\n");
fprintf(ficresf,"# Covariate valuofcovar yearproj age");	fprintf(ficresf,"# Covariate valuofcovar yearproj age");
Line 5003 void prevforecast(char fileres[], double	Line 5141 void prevforecast(char fileres[], double
if (hhstepm/YEARMstepm ==yearp) {	if (hhstepm/YEARMstepm ==yearp) {
fprintf(ficresf,"\n");	fprintf(ficresf,"\n");
for(j=1;j<=cptcoveff;j++)	for(j=1;j<=cptcoveff;j++)
fprintf(ficresf,"%d %d ",Tvaraff[j],nbcode[Tvaraff[j]][codtab[k][j]]);	fprintf(ficresf,"%d %d ",Tvaraff[j],nbcode[Tvaraff[j]][codtabm(k,j)]);
fprintf(ficresf,"%.f %.f ",anproj1+yearp,agec+hhstepm/YEARMstepm);	fprintf(ficresf,"%.f %.f ",anproj1+yearp,agec+hhstepm/YEARMstepm);
}	}
for(j=1; j<=nlstate+ndeath;j++) {	for(j=1; j<=nlstate+ndeath;j++) {
Line 5053 void populforecast(char fileres[], doubl	Line 5191 void populforecast(char fileres[], doubl
prevalence(probs, ageminpar, agemax, s, agev, nlstate, imx, Tvar, nbcode, ncodemax, mint, anint, dateprev1, dateprev2, firstpass, lastpass);	prevalence(probs, ageminpar, agemax, s, agev, nlstate, imx, Tvar, nbcode, ncodemax, mint, anint, dateprev1, dateprev2, firstpass, lastpass);


strcpy(filerespop,"pop");	strcpy(filerespop,"POP_");
strcat(filerespop,fileres);	strcat(filerespop,fileresu);
if((ficrespop=fopen(filerespop,"w"))==NULL) {	if((ficrespop=fopen(filerespop,"w"))==NULL) {
printf("Problem with forecast resultfile: %s\n", filerespop);	printf("Problem with forecast resultfile: %s\n", filerespop);
fprintf(ficlog,"Problem with forecast resultfile: %s\n", filerespop);	fprintf(ficlog,"Problem with forecast resultfile: %s\n", filerespop);
Line 5101 void populforecast(char fileres[], doubl	Line 5239 void populforecast(char fileres[], doubl
k=k+1;	k=k+1;
fprintf(ficrespop,"\n#******");	fprintf(ficrespop,"\n#******");
for(j=1;j<=cptcoveff;j++) {	for(j=1;j<=cptcoveff;j++) {
fprintf(ficrespop," V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtab[k][j]]);	fprintf(ficrespop," V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtabm(k,j)]);
}	}
fprintf(ficrespop,"******\n");	fprintf(ficrespop,"******\n");
fprintf(ficrespop,"# Age");	fprintf(ficrespop,"# Age");
Line 5407 double gompertz_f(const gsl_vector *v, v	Line 5545 double gompertz_f(const gsl_vector *v, v
#endif	#endif

/***************** Printing html file *********/	/***************** Printing html file *********/
void printinghtmlmort(char fileres[], char title[], char datafile[], int firstpass, \	void printinghtmlmort(char fileresu[], char title[], char datafile[], int firstpass, \
int lastpass, int stepm, int weightopt, char model[],\	int lastpass, int stepm, int weightopt, char model[],\
int imx, double p[],double **matcov,double agemortsup){	int imx, double p[],double **matcov,double agemortsup){
int i,k;	int i,k;
Line 5416 void printinghtmlmort(char fileres[], ch	Line 5554 void printinghtmlmort(char fileres[], ch
fprintf(fichtm," mu(age) =%lfexp(%lf(age-%d)) per year<br><br>",p[1],p[2],agegomp);	fprintf(fichtm," mu(age) =%lfexp(%lf(age-%d)) per year<br><br>",p[1],p[2],agegomp);
for (i=1;i<=2;i++)	for (i=1;i<=2;i++)
fprintf(fichtm," p[%d] = %lf [%f ; %f]<br>\n",i,p[i],p[i]-2sqrt(matcov[i][i]),p[i]+2sqrt(matcov[i][i]));	fprintf(fichtm," p[%d] = %lf [%f ; %f]<br>\n",i,p[i],p[i]-2sqrt(matcov[i][i]),p[i]+2sqrt(matcov[i][i]));
fprintf(fichtm,"<br><br><img src=\"graphmort.png\">");	fprintf(fichtm,"<br><br><img src=\"graphmort.svg\">");
fprintf(fichtm,"</ul>");	fprintf(fichtm,"</ul>");

fprintf(fichtm,"<ul><li><h4>Life table</h4>\n <br>");	fprintf(fichtm,"<ul><li><h4>Life table</h4>\n <br>");
Line 5431 fprintf(fichtm,"<ul><li><h4>Life table</	Line 5569 fprintf(fichtm,"<ul><li><h4>Life table</
}	}

/***************** Gnuplot file ************/	/***************** Gnuplot file ************/
void printinggnuplotmort(char fileres[], char optionfilefiname[], double ageminpar, double agemaxpar, double fage , char pathc[], double p[]){	void printinggnuplotmort(char fileresu[], char optionfilefiname[], double ageminpar, double agemaxpar, double fage , char pathc[], double p[]){

char dirfileres[132],optfileres[132];	char dirfileres[132],optfileres[132];

Line 5445 void printinggnuplotmort(char fileres[],	Line 5583 void printinggnuplotmort(char fileres[],

strcpy(dirfileres,optionfilefiname);	strcpy(dirfileres,optionfilefiname);
strcpy(optfileres,"vpl");	strcpy(optfileres,"vpl");
fprintf(ficgp,"set out \"graphmort.png\"\n ");	fprintf(ficgp,"set out \"graphmort.svg\"\n ");
fprintf(ficgp,"set xlabel \"Age\"\n set ylabel \"Force of mortality (per year)\" \n ");	fprintf(ficgp,"set xlabel \"Age\"\n set ylabel \"Force of mortality (per year)\" \n ");
fprintf(ficgp, "set ter png small size 320, 240\n set log y\n");	fprintf(ficgp, "set ter svg size 640, 480\n set log y\n");
/* fprintf(ficgp, "set size 0.65,0.65\n"); */	/* fprintf(ficgp, "set size 0.65,0.65\n"); */
fprintf(ficgp,"plot [%d:100] %lfexp(%lf(x-%d))",agegomp,p[1],p[2],agegomp);	fprintf(ficgp,"plot [%d:100] %lfexp(%lf(x-%d))",agegomp,p[1],p[2],agegomp);

Line 5758 int decodemodel ( char model[], int last	Line 5896 int decodemodel ( char model[], int last
/* k=1 Tvar[1]=2 (from V2) */	/* k=1 Tvar[1]=2 (from V2) */
/* k=5 Tvar[5] */	/* k=5 Tvar[5] */
/* for (k=1; k<=cptcovn;k++) { */	/* for (k=1; k<=cptcovn;k++) { */
/* cov[2+k]=nbcode[Tvar[k]][codtab[ij][Tvar[k]]]; */	/* cov[2+k]=nbcode[Tvar[k]][codtabm(ij,Tvar[k])]; */
/* } */	/* } */
/* for (k=1; k<=cptcovage;k++) cov[2+Tage[k]]=nbcode[Tvar[Tage[k]]][codtab[ij][Tvar[Tage[k]]]]cov[2]; /	/* for (k=1; k<=cptcovage;k++) cov[2+Tage[k]]=nbcode[Tvar[Tage[k]]][codtabm(ij,Tvar[Tage[k])]]cov[2]; /
/*	/*
* Treating invertedly V2+V1+V3age+V2V4 is as if written V2V4 +V3age + V1 + V2 */	* Treating invertedly V2+V1+V3age+V2V4 is as if written V2V4 +V3age + V1 + V2 */
for(k=cptcovt; k>=1;k--) /*< Number of covariates /	for(k=cptcovt; k>=1;k--) /*< Number of covariates /
Line 6144 int prevalence_limit(double p, double	Line 6282 int prevalence_limit(double p, double
double ftolpl = 1.e-10;	double ftolpl = 1.e-10;
double age, agebase, agelim;	double age, agebase, agelim;

strcpy(filerespl,"pl");	strcpy(filerespl,"PL_");
strcat(filerespl,fileres);	strcat(filerespl,fileresu);
if((ficrespl=fopen(filerespl,"w"))==NULL) {	if((ficrespl=fopen(filerespl,"w"))==NULL) {
printf("Problem with period (stable) prevalence resultfile: %s\n", filerespl);return 1;	printf("Problem with period (stable) prevalence resultfile: %s\n", filerespl);return 1;
fprintf(ficlog,"Problem with period (stable) prevalence resultfile: %s\n", filerespl);return 1;	fprintf(ficlog,"Problem with period (stable) prevalence resultfile: %s\n", filerespl);return 1;
Line 6171 int prevalence_limit(double p, double	Line 6309 int prevalence_limit(double p, double
//for(cptcod=1;cptcod<=ncodemax[cptcov];cptcod++){	//for(cptcod=1;cptcod<=ncodemax[cptcov];cptcod++){
k=k+1;	k=k+1;
/* to clean */	/* to clean */
//printf("cptcov=%d cptcod=%d codtab=%d\n",cptcov, cptcod,codtab[cptcod][cptcov]);	//printf("cptcov=%d cptcod=%d codtab=%d\n",cptcov, cptcod,codtabm(cptcod,cptcov));
fprintf(ficrespl,"\n#******");	fprintf(ficrespl,"#******");
printf("\n#******");	printf("#******");
fprintf(ficlog,"\n#******");	fprintf(ficlog,"#******");
for(j=1;j<=cptcoveff;j++) {	for(j=1;j<=cptcoveff;j++) {
fprintf(ficrespl," V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtab[k][j]]);	fprintf(ficrespl," V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtabm(k,j)]);
printf(" V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtab[k][j]]);	printf(" V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtabm(k,j)]);
fprintf(ficlog," V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtab[k][j]]);	fprintf(ficlog," V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtabm(k,j)]);
}	}
fprintf(ficrespl,"******\n");	fprintf(ficrespl,"******\n");
printf("******\n");	printf("******\n");
Line 6186 int prevalence_limit(double p, double	Line 6324 int prevalence_limit(double p, double

fprintf(ficrespl,"#Age ");	fprintf(ficrespl,"#Age ");
for(j=1;j<=cptcoveff;j++) {	for(j=1;j<=cptcoveff;j++) {
fprintf(ficrespl,"V%d %d",Tvaraff[j],nbcode[Tvaraff[j]][codtab[k][j]]);	fprintf(ficrespl,"V%d %d ",Tvaraff[j],nbcode[Tvaraff[j]][codtabm(k,j)]);
}	}
for(i=1; i<=nlstate;i++) fprintf(ficrespl,"%d-%d ",i,i);	for(i=1; i<=nlstate;i++) fprintf(ficrespl,"%d-%d ",i,i);
fprintf(ficrespl,"\n");	fprintf(ficrespl,"\n");
Line 6196 int prevalence_limit(double p, double	Line 6334 int prevalence_limit(double p, double
prevalim(prlim, nlstate, p, age, oldm, savm,ftolpl,k);	prevalim(prlim, nlstate, p, age, oldm, savm,ftolpl,k);
fprintf(ficrespl,"%.0f ",age );	fprintf(ficrespl,"%.0f ",age );
for(j=1;j<=cptcoveff;j++)	for(j=1;j<=cptcoveff;j++)
fprintf(ficrespl,"%d %d ",Tvaraff[j],nbcode[Tvaraff[j]][codtab[k][j]]);	fprintf(ficrespl,"%d %d ",Tvaraff[j],nbcode[Tvaraff[j]][codtabm(k,j)]);
for(i=1; i<=nlstate;i++)	for(i=1; i<=nlstate;i++)
fprintf(ficrespl," %.5f", prlim[i][i]);	fprintf(ficrespl," %.5f", prlim[i][i]);
fprintf(ficrespl,"\n");	fprintf(ficrespl,"\n");
Line 6218 int hPijx(double *p, int bage, int fage)	Line 6356 int hPijx(double *p, int bage, int fage)
double agedeb;	double agedeb;
double ***p3mat;	double ***p3mat;

strcpy(filerespij,"pij"); strcat(filerespij,fileres);	strcpy(filerespij,"PIJ_"); strcat(filerespij,fileresu);
if((ficrespij=fopen(filerespij,"w"))==NULL) {	if((ficrespij=fopen(filerespij,"w"))==NULL) {
printf("Problem with Pij resultfile: %s\n", filerespij); return 1;	printf("Problem with Pij resultfile: %s\n", filerespij); return 1;
fprintf(ficlog,"Problem with Pij resultfile: %s\n", filerespij); return 1;	fprintf(ficlog,"Problem with Pij resultfile: %s\n", filerespij); return 1;
Line 6243 int hPijx(double *p, int bage, int fage)	Line 6381 int hPijx(double *p, int bage, int fage)
for (k=1; k <= (int) pow(2,cptcoveff); k++){	for (k=1; k <= (int) pow(2,cptcoveff); k++){
fprintf(ficrespij,"\n#****** ");	fprintf(ficrespij,"\n#****** ");
for(j=1;j<=cptcoveff;j++)	for(j=1;j<=cptcoveff;j++)
fprintf(ficrespij,"V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtab[k][j]]);	fprintf(ficrespij,"V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtabm(k,j)]);
fprintf(ficrespij,"******\n");	fprintf(ficrespij,"******\n");

for (agedeb=fage; agedeb>=bage; agedeb--){ /* If stepm=6 months */	for (agedeb=fage; agedeb>=bage; agedeb--){ /* If stepm=6 months */
Line 6479 int main(int argc, char *argv[])	Line 6617 int main(int argc, char *argv[])
/* */	/* */
strcpy(fileres,"r");	strcpy(fileres,"r");
strcat(fileres, optionfilefiname);	strcat(fileres, optionfilefiname);
	strcat(fileresu, optionfilefiname); /* Without r in front */
strcat(fileres,".txt"); /* Other files have txt extension */	strcat(fileres,".txt"); /* Other files have txt extension */
	strcat(fileresu,".txt"); /* Other files have txt extension */

/* Main ---------arguments file --------*/	/* Main ---------arguments file --------*/

Line 6494 int main(int argc, char *argv[])	Line 6634 int main(int argc, char *argv[])


strcpy(filereso,"o");	strcpy(filereso,"o");
strcat(filereso,fileres);	strcat(filereso,fileresu);
if((ficparo=fopen(filereso,"w"))==NULL) { /* opened on subdirectory */	if((ficparo=fopen(filereso,"w"))==NULL) { /* opened on subdirectory */
printf("Problem with Output resultfile: %s\n", filereso);	printf("Problem with Output resultfile: %s\n", filereso);
fprintf(ficlog,"Problem with Output resultfile: %s\n", filereso);	fprintf(ficlog,"Problem with Output resultfile: %s\n", filereso);
Line 6557 int main(int argc, char *argv[])	Line 6697 int main(int argc, char *argv[])
}else	}else
break;	break;
}	}
if((num_filled=sscanf(line,"model=1+age%[^.\n]\n", model)) !=EOF){	if((num_filled=sscanf(line,"model=1+age%[^.\n]", model)) !=EOF){
if (num_filled != 1) {	if (num_filled == 0)
	model[0]='\0';
	else if (num_filled != 1){
printf("ERROR %d: Model should be at minimum 'model=1+age.' %s\n",num_filled, line);	printf("ERROR %d: Model should be at minimum 'model=1+age.' %s\n",num_filled, line);
fprintf(ficlog,"ERROR %d: Model should be at minimum 'model=1+age.' %s\n",num_filled, line);	fprintf(ficlog,"ERROR %d: Model should be at minimum 'model=1+age.' %s\n",num_filled, line);
model[0]='\0';	model[0]='\0';
Line 6568 int main(int argc, char *argv[])	Line 6710 int main(int argc, char *argv[])
if (model[0]=='+'){	if (model[0]=='+'){
for(i=1; i<=strlen(model);i++)	for(i=1; i<=strlen(model);i++)
modeltemp[i-1]=model[i];	modeltemp[i-1]=model[i];
	strcpy(model,modeltemp);
}	}
strcpy(model,modeltemp);
}	}
printf(" model=1+age%s modeltemp= %s, model=%s\n",model, modeltemp, model);fflush(stdout);	/* printf(" model=1+age%s modeltemp= %s, model=%s\n",model, modeltemp, model);fflush(stdout); */
}	}
/* fscanf(ficpar,"title=%s datafile=%s lastobs=%d firstpass=%d lastpass=%d\nftol=%lf stepm=%d ncovcol=%d nlstate=%d ndeath=%d maxwav=%d mle=%d weight=%d model=1+age+%s\n",title, datafile, &lastobs, &firstpass,&lastpass,&ftol, &stepm, &ncovcol, &nlstate,&ndeath, &maxwav, &mle, &weightopt,model); */	/* fscanf(ficpar,"title=%s datafile=%s lastobs=%d firstpass=%d lastpass=%d\nftol=%lf stepm=%d ncovcol=%d nlstate=%d ndeath=%d maxwav=%d mle=%d weight=%d model=1+age+%s\n",title, datafile, &lastobs, &firstpass,&lastpass,&ftol, &stepm, &ncovcol, &nlstate,&ndeath, &maxwav, &mle, &weightopt,model); */
/* numlinepar=numlinepar+3; /\* In general \/ /	/* numlinepar=numlinepar+3; /\* In general \/ /
Line 6803 Please run with mle=-1 to get a correct	Line 6945 Please run with mle=-1 to get a correct
strcat(rfileres,"."); /* */	strcat(rfileres,"."); /* */
strcat(rfileres,optionfilext); /* Other files have txt extension */	strcat(rfileres,optionfilext); /* Other files have txt extension */
if((ficres =fopen(rfileres,"w"))==NULL) {	if((ficres =fopen(rfileres,"w"))==NULL) {
printf("Problem writing new parameter file: %s\n", fileres);goto end;	printf("Problem writing new parameter file: %s\n", rfileres);goto end;
fprintf(ficlog,"Problem writing new parameter file: %s\n", fileres);goto end;	fprintf(ficlog,"Problem writing new parameter file: %s\n", rfileres);goto end;
}	}
fprintf(ficres,"#%s\n",version);	fprintf(ficres,"#%s\n",version);
} /* End of mle != -3 */	} /* End of mle != -3 */
Line 6914 Please run with mle=-1 to get a correct	Line 7056 Please run with mle=-1 to get a correct
V2+V1age, there are 3 covariates Tvar[2]=1 (V1)./	V2+V1age, there are 3 covariates Tvar[2]=1 (V1)./
/* 1 to ncodemax[j] is the maximum value of this jth covariate */	/* 1 to ncodemax[j] is the maximum value of this jth covariate */

codtab=imatrix(1,100,1,10); /* codtab[h,k]=( (h-1) - mod(k-1,2(k-1) )/2(k-1) */	/* codtab=imatrix(1,100,1,10);/ / codtab[h,k]=( (h-1) - mod(k-1,2(k-1) )/2(k-1) */
/printf(" codtab[1,1],codtab[100,10]=%d,%d\n", codtab[1][1],codtab[100][10]);/	/printf(" codtab[1,1],codtab[100,10]=%d,%d\n", codtab[1][1],codtabm(100,10));/
/* codtab gives the value 1 or 2 of the hth combination of k covariates (1 or 2).*/	/* codtab gives the value 1 or 2 of the hth combination of k covariates (1 or 2).*/
h=0;	h=0;

Line 6949 Please run with mle=-1 to get a correct	Line 7091 Please run with mle=-1 to get a correct
*/	*/
for(h=1; h <=100 ;h++){	for(h=1; h <=100 ;h++){
/* printf("h=%2d ", h); */	/* printf("h=%2d ", h); */
for(k=1; k <=10; k++){	/* for(k=1; k <=10; k++){ */
/* printf("k=%d %d ",k,codtabm(h,k)); */	/* printf("k=%d %d ",k,codtabm(h,k)); */
codtab[h][k]=codtabm(h,k);	/* codtab[h][k]=codtabm(h,k); */
}	/* } */
/* printf("\n"); */	/* printf("\n"); */
}	}
/* for(k=1;k<=cptcoveff; k++){ /\* scans any effective covariate \/ /	/* for(k=1;k<=cptcoveff; k++){ /\* scans any effective covariate \/ /
Line 6987 Please run with mle=-1 to get a correct	Line 7129 Please run with mle=-1 to get a correct
/* Initialisation of ----------- gnuplot -------------*/	/* Initialisation of ----------- gnuplot -------------*/
strcpy(optionfilegnuplot,optionfilefiname);	strcpy(optionfilegnuplot,optionfilefiname);
if(mle==-3)	if(mle==-3)
strcat(optionfilegnuplot,"-mort");	strcat(optionfilegnuplot,"-MORT_");
strcat(optionfilegnuplot,".gp");	strcat(optionfilegnuplot,".gp");

if((ficgp=fopen(optionfilegnuplot,"w"))==NULL) {	if((ficgp=fopen(optionfilegnuplot,"w"))==NULL) {
Line 7006 Please run with mle=-1 to get a correct	Line 7148 Please run with mle=-1 to get a correct

strcpy(optionfilehtm,optionfilefiname); /* Main html file */	strcpy(optionfilehtm,optionfilefiname); /* Main html file */
if(mle==-3)	if(mle==-3)
strcat(optionfilehtm,"-mort");	strcat(optionfilehtm,"-MORT_");
strcat(optionfilehtm,".htm");	strcat(optionfilehtm,".htm");
if((fichtm=fopen(optionfilehtm,"w"))==NULL) {	if((fichtm=fopen(optionfilehtm,"w"))==NULL) {
printf("Problem with %s \n",optionfilehtm);	printf("Problem with %s \n",optionfilehtm);
Line 7122 Interval (in months) between two waves:	Line 7264 Interval (in months) between two waves:
#else	#else
printf("Powell\n"); fprintf(ficlog,"Powell\n");	printf("Powell\n"); fprintf(ficlog,"Powell\n");
#endif	#endif
strcpy(filerespow,"pow-mort");	strcpy(filerespow,"POW-MORT_");
strcat(filerespow,fileres);	strcat(filerespow,fileresu);
if((ficrespow=fopen(filerespow,"w"))==NULL) {	if((ficrespow=fopen(filerespow,"w"))==NULL) {
printf("Problem with resultfile: %s\n", filerespow);	printf("Problem with resultfile: %s\n", filerespow);
fprintf(ficlog,"Problem with resultfile: %s\n", filerespow);	fprintf(ficlog,"Problem with resultfile: %s\n", filerespow);
Line 7277 Please run with mle=-1 to get a correct	Line 7419 Please run with mle=-1 to get a correct
This is probably because your parameter file doesn't \n contain the exact number of lines (or columns) corresponding to your model line.\n\	This is probably because your parameter file doesn't \n contain the exact number of lines (or columns) corresponding to your model line.\n\
Please run with mle=-1 to get a correct covariance matrix.\n",ageminpar,agemaxpar);	Please run with mle=-1 to get a correct covariance matrix.\n",ageminpar,agemaxpar);
}else	}else
printinggnuplotmort(fileres, optionfilefiname,ageminpar,agemaxpar,fage, pathc,p);	printinggnuplotmort(fileresu, optionfilefiname,ageminpar,agemaxpar,fage, pathc,p);
printinghtmlmort(fileres,title,datafile, firstpass, lastpass, \	printinghtmlmort(fileresu,title,datafile, firstpass, lastpass, \
stepm, weightopt,\	stepm, weightopt,\
model,imx,p,matcov,agemortsup);	model,imx,p,matcov,agemortsup);

Line 7548 Please run with mle=-1 to get a correct	Line 7690 Please run with mle=-1 to get a correct
This is probably because your parameter file doesn't \n contain the exact number of lines (or columns) corresponding to your model line.\n\	This is probably because your parameter file doesn't \n contain the exact number of lines (or columns) corresponding to your model line.\n\
Please run with mle=-1 to get a correct covariance matrix.\n",ageminpar,agemaxpar);	Please run with mle=-1 to get a correct covariance matrix.\n",ageminpar,agemaxpar);
}else	}else
printinggnuplot(fileres, optionfilefiname,ageminpar,agemaxpar,fage, pathc,p);	printinggnuplot(fileresu, optionfilefiname,ageminpar,agemaxpar,fage, pathc,p);

printinghtml(fileres,title,datafile, firstpass, lastpass, stepm, weightopt,\	printinghtml(fileresu,title,datafile, firstpass, lastpass, stepm, weightopt,\
model,imx,jmin,jmax,jmean,rfileres,popforecast,estepm,\	model,imx,jmin,jmax,jmean,rfileres,popforecast,estepm,\
jprev1,mprev1,anprev1,jprev2,mprev2,anprev2);	jprev1,mprev1,anprev1,jprev2,mprev2,anprev2);

Line 7602 Please run with mle=-1 to get a correct	Line 7744 Please run with mle=-1 to get a correct
/if((stepm == 1) && (strcmp(model,".")==0)){/	/if((stepm == 1) && (strcmp(model,".")==0)){/
if(prevfcast==1){	if(prevfcast==1){
/* if(stepm ==1){*/	/* if(stepm ==1){*/
prevforecast(fileres, anproj1, mproj1, jproj1, agemin, agemax, dateprev1, dateprev2, mobilavproj, bage, fage, firstpass, lastpass, anproj2, p, cptcoveff);	prevforecast(fileresu, anproj1, mproj1, jproj1, agemin, agemax, dateprev1, dateprev2, mobilavproj, bage, fage, firstpass, lastpass, anproj2, p, cptcoveff);
/* (popforecast==1) populforecast(fileres, anpyram,mpyram,jpyram, agemin,agemax, dateprev1, dateprev2,mobilav, agedeb, fage, popforecast, popfile, anpyram1,p, i1);*/	/* (popforecast==1) populforecast(fileres, anpyram,mpyram,jpyram, agemin,agemax, dateprev1, dateprev2,mobilav, agedeb, fage, popforecast, popfile, anpyram1,p, i1);*/
/* } */	/* } */
/* else{ */	/* else{ */
Line 7632 Please run with mle=-1 to get a correct	Line 7774 Please run with mle=-1 to get a correct

/---------- Health expectancies, no variances ------------/	/---------- Health expectancies, no variances ------------/

strcpy(filerese,"e");	strcpy(filerese,"E_");
strcat(filerese,fileres);	strcat(filerese,fileresu);
if((ficreseij=fopen(filerese,"w"))==NULL) {	if((ficreseij=fopen(filerese,"w"))==NULL) {
printf("Problem with Health Exp. resultfile: %s\n", filerese); exit(0);	printf("Problem with Health Exp. resultfile: %s\n", filerese); exit(0);
fprintf(ficlog,"Problem with Health Exp. resultfile: %s\n", filerese); exit(0);	fprintf(ficlog,"Problem with Health Exp. resultfile: %s\n", filerese); exit(0);
Line 7646 Please run with mle=-1 to get a correct	Line 7788 Please run with mle=-1 to get a correct
for (k=1; k <= (int) pow(2,cptcoveff); k++){	for (k=1; k <= (int) pow(2,cptcoveff); k++){
fprintf(ficreseij,"\n#****** ");	fprintf(ficreseij,"\n#****** ");
for(j=1;j<=cptcoveff;j++) {	for(j=1;j<=cptcoveff;j++) {
fprintf(ficreseij,"V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtab[k][j]]);	fprintf(ficreseij,"V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtabm(k,j)]);
}	}
fprintf(ficreseij,"******\n");	fprintf(ficreseij,"******\n");

Line 7663 Please run with mle=-1 to get a correct	Line 7805 Please run with mle=-1 to get a correct
/---------- Health expectancies and variances ------------/	/---------- Health expectancies and variances ------------/


strcpy(filerest,"t");	strcpy(filerest,"T_");
strcat(filerest,fileres);	strcat(filerest,fileresu);
if((ficrest=fopen(filerest,"w"))==NULL) {	if((ficrest=fopen(filerest,"w"))==NULL) {
printf("Problem with total LE resultfile: %s\n", filerest);goto end;	printf("Problem with total LE resultfile: %s\n", filerest);goto end;
fprintf(ficlog,"Problem with total LE resultfile: %s\n", filerest);goto end;	fprintf(ficlog,"Problem with total LE resultfile: %s\n", filerest);goto end;
Line 7673 Please run with mle=-1 to get a correct	Line 7815 Please run with mle=-1 to get a correct
fprintf(ficlog,"Computing Total Life expectancies with their standard errors: file '%s' \n", filerest);	fprintf(ficlog,"Computing Total Life expectancies with their standard errors: file '%s' \n", filerest);


strcpy(fileresstde,"stde");	strcpy(fileresstde,"STDE_");
strcat(fileresstde,fileres);	strcat(fileresstde,fileresu);
if((ficresstdeij=fopen(fileresstde,"w"))==NULL) {	if((ficresstdeij=fopen(fileresstde,"w"))==NULL) {
printf("Problem with Health Exp. and std errors resultfile: %s\n", fileresstde); exit(0);	printf("Problem with Health Exp. and std errors resultfile: %s\n", fileresstde); exit(0);
fprintf(ficlog,"Problem with Health Exp. and std errors resultfile: %s\n", fileresstde); exit(0);	fprintf(ficlog,"Problem with Health Exp. and std errors resultfile: %s\n", fileresstde); exit(0);
Line 7682 Please run with mle=-1 to get a correct	Line 7824 Please run with mle=-1 to get a correct
printf("Computing Health Expectancies and standard errors: result on file '%s' \n", fileresstde);	printf("Computing Health Expectancies and standard errors: result on file '%s' \n", fileresstde);
fprintf(ficlog,"Computing Health Expectancies and standard errors: result on file '%s' \n", fileresstde);	fprintf(ficlog,"Computing Health Expectancies and standard errors: result on file '%s' \n", fileresstde);

strcpy(filerescve,"cve");	strcpy(filerescve,"CVE_");
strcat(filerescve,fileres);	strcat(filerescve,fileresu);
if((ficrescveij=fopen(filerescve,"w"))==NULL) {	if((ficrescveij=fopen(filerescve,"w"))==NULL) {
printf("Problem with Covar. Health Exp. resultfile: %s\n", filerescve); exit(0);	printf("Problem with Covar. Health Exp. resultfile: %s\n", filerescve); exit(0);
fprintf(ficlog,"Problem with Covar. Health Exp. resultfile: %s\n", filerescve); exit(0);	fprintf(ficlog,"Problem with Covar. Health Exp. resultfile: %s\n", filerescve); exit(0);
Line 7691 Please run with mle=-1 to get a correct	Line 7833 Please run with mle=-1 to get a correct
printf("Computing Covar. of Health Expectancies: result on file '%s' \n", filerescve);	printf("Computing Covar. of Health Expectancies: result on file '%s' \n", filerescve);
fprintf(ficlog,"Computing Covar. of Health Expectancies: result on file '%s' \n", filerescve);	fprintf(ficlog,"Computing Covar. of Health Expectancies: result on file '%s' \n", filerescve);

strcpy(fileresv,"v");	strcpy(fileresv,"V_");
strcat(fileresv,fileres);	strcat(fileresv,fileresu);
if((ficresvij=fopen(fileresv,"w"))==NULL) {	if((ficresvij=fopen(fileresv,"w"))==NULL) {
printf("Problem with variance resultfile: %s\n", fileresv);exit(0);	printf("Problem with variance resultfile: %s\n", fileresv);exit(0);
fprintf(ficlog,"Problem with variance resultfile: %s\n", fileresv);exit(0);	fprintf(ficlog,"Problem with variance resultfile: %s\n", fileresv);exit(0);
Line 7706 Please run with mle=-1 to get a correct	Line 7848 Please run with mle=-1 to get a correct
for (k=1; k <= (int) pow(2,cptcoveff); k++){	for (k=1; k <= (int) pow(2,cptcoveff); k++){
fprintf(ficrest,"\n#****** ");	fprintf(ficrest,"\n#****** ");
for(j=1;j<=cptcoveff;j++)	for(j=1;j<=cptcoveff;j++)
fprintf(ficrest,"V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtab[k][j]]);	fprintf(ficrest,"V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtabm(k,j)]);
fprintf(ficrest,"******\n");	fprintf(ficrest,"******\n");

fprintf(ficresstdeij,"\n#****** ");	fprintf(ficresstdeij,"\n#****** ");
fprintf(ficrescveij,"\n#****** ");	fprintf(ficrescveij,"\n#****** ");
for(j=1;j<=cptcoveff;j++) {	for(j=1;j<=cptcoveff;j++) {
fprintf(ficresstdeij,"V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtab[k][j]]);	fprintf(ficresstdeij,"V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtabm(k,j)]);
fprintf(ficrescveij,"V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtab[k][j]]);	fprintf(ficrescveij,"V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtabm(k,j)]);
}	}
fprintf(ficresstdeij,"******\n");	fprintf(ficresstdeij,"******\n");
fprintf(ficrescveij,"******\n");	fprintf(ficrescveij,"******\n");

fprintf(ficresvij,"\n#****** ");	fprintf(ficresvij,"\n#****** ");
for(j=1;j<=cptcoveff;j++)	for(j=1;j<=cptcoveff;j++)
fprintf(ficresvij,"V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtab[k][j]]);	fprintf(ficresvij,"V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtabm(k,j)]);
fprintf(ficresvij,"******\n");	fprintf(ficresvij,"******\n");

eij=ma3x(1,nlstate,1,nlstate,(int) bage, (int) fage);	eij=ma3x(1,nlstate,1,nlstate,(int) bage, (int) fage);
Line 7735 Please run with mle=-1 to get a correct	Line 7877 Please run with mle=-1 to get a correct


for(vpopbased=0; vpopbased <= popbased; vpopbased++){ /* Done for vpopbased=0 and vpopbased=1 if popbased==1*/	for(vpopbased=0; vpopbased <= popbased; vpopbased++){ /* Done for vpopbased=0 and vpopbased=1 if popbased==1*/
oldm=oldms;savm=savms; /* Segmentation fault */	oldm=oldms;savm=savms; /* ZZ Segmentation fault */
cptcod= 0; /* To be deleted */	cptcod= 0; /* To be deleted */
varevsij(optionfilefiname, vareij, matcov, p, delti, nlstate, stepm, (int) bage, (int) fage, oldm, savm, prlim, ftolpl,k, estepm, cptcov,cptcod,vpopbased,mobilav, strstart); /* cptcod not initialized Intel */	varevsij(optionfilefiname, vareij, matcov, p, delti, nlstate, stepm, (int) bage, (int) fage, oldm, savm, prlim, ftolpl,k, estepm, cptcov,cptcod,vpopbased,mobilav, strstart); /* cptcod not initialized Intel */
fprintf(ficrest,"# Total life expectancy with std error and decomposition into time to be expected in each health state\n# (weighted average of eij where weights are ");	fprintf(ficrest,"# Total life expectancy with std error and decomposition into time to be expected in each health state\n# (weighted average of eij where weights are ");
Line 7743 Please run with mle=-1 to get a correct	Line 7885 Please run with mle=-1 to get a correct
fprintf(ficrest,"the age specific prevalence observed (cross-sectionally) in the population i.e cross-sectionally\n in each health state (popbased=1) (mobilav=%d)\n",mobilav);	fprintf(ficrest,"the age specific prevalence observed (cross-sectionally) in the population i.e cross-sectionally\n in each health state (popbased=1) (mobilav=%d)\n",mobilav);
else	else
fprintf(ficrest,"the age specific period (stable) prevalences in each health state \n");	fprintf(ficrest,"the age specific period (stable) prevalences in each health state \n");
fprintf(ficrest,"# Age e.. (std) ");	fprintf(ficrest,"# Age popbased mobilav e.. (std) ");
for (i=1;i<=nlstate;i++) fprintf(ficrest,"e.%d (std) ",i);	for (i=1;i<=nlstate;i++) fprintf(ficrest,"e.%d (std) ",i);
fprintf(ficrest,"\n");	fprintf(ficrest,"\n");
	/* printf("Which p?\n"); for(i=1;i<=npar;i++)printf("p[i=%d]=%lf,",i,p[i]);printf("\n"); */
epj=vector(1,nlstate+1);	epj=vector(1,nlstate+1);
for(age=bage; age <=fage ;age++){	for(age=bage; age <=fage ;age++){
prevalim(prlim, nlstate, p, age, oldm, savm,ftolpl,k);	prevalim(prlim, nlstate, p, age, oldm, savm,ftolpl,k); /ZZ Is it the correct prevalim /
if (vpopbased==1) {	if (vpopbased==1) {
if(mobilav ==0){	if(mobilav ==0){
for(i=1; i<=nlstate;i++)	for(i=1; i<=nlstate;i++)
Line 7760 Please run with mle=-1 to get a correct	Line 7902 Please run with mle=-1 to get a correct
}	}
}	}

fprintf(ficrest," %4.0f",age);	fprintf(ficrest," %4.0f %d %d",age, vpopbased, mobilav);
	/* printf(" age %4.0f ",age); */
for(j=1, epj[nlstate+1]=0.;j <=nlstate;j++){	for(j=1, epj[nlstate+1]=0.;j <=nlstate;j++){
for(i=1, epj[j]=0.;i <=nlstate;i++) {	for(i=1, epj[j]=0.;i <=nlstate;i++) {
epj[j] += prlim[i][i]*eij[i][j][(int)age];	epj[j] += prlim[i][i]*eij[i][j][(int)age];
/* printf("%lf %lf ", prlim[i][i] ,eij[i][j][(int)age]);*/	/ZZZ printf("%lf %lf ", prlim[i][i] ,eij[i][j][(int)age]);/
	/* printf("%lf %lf ", prlim[i][i] ,eij[i][j][(int)age]); */
}	}
epj[nlstate+1] +=epj[j];	epj[nlstate+1] +=epj[j];
}	}
	/* printf(" age %4.0f \n",age); */

for(i=1, vepp=0.;i <=nlstate;i++)	for(i=1, vepp=0.;i <=nlstate;i++)
for(j=1;j <=nlstate;j++)	for(j=1;j <=nlstate;j++)
Line 7799 Please run with mle=-1 to get a correct	Line 7944 Please run with mle=-1 to get a correct

/------- Variance of period (stable) prevalence------/	/------- Variance of period (stable) prevalence------/

strcpy(fileresvpl,"vpl");	strcpy(fileresvpl,"VPL_");
strcat(fileresvpl,fileres);	strcat(fileresvpl,fileresu);
if((ficresvpl=fopen(fileresvpl,"w"))==NULL) {	if((ficresvpl=fopen(fileresvpl,"w"))==NULL) {
printf("Problem with variance of period (stable) prevalence resultfile: %s\n", fileresvpl);	printf("Problem with variance of period (stable) prevalence resultfile: %s\n", fileresvpl);
exit(0);	exit(0);
Line 7813 Please run with mle=-1 to get a correct	Line 7958 Please run with mle=-1 to get a correct
for (k=1; k <= (int) pow(2,cptcoveff); k++){	for (k=1; k <= (int) pow(2,cptcoveff); k++){
fprintf(ficresvpl,"\n#****** ");	fprintf(ficresvpl,"\n#****** ");
for(j=1;j<=cptcoveff;j++)	for(j=1;j<=cptcoveff;j++)
fprintf(ficresvpl,"V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtab[k][j]]);	fprintf(ficresvpl,"V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtabm(k,j)]);
fprintf(ficresvpl,"******\n");	fprintf(ficresvpl,"******\n");

varpl=matrix(1,nlstate,(int) bage, (int) fage);	varpl=matrix(1,nlstate,(int) bage, (int) fage);
Line 7850 Please run with mle=-1 to get a correct	Line 7995 Please run with mle=-1 to get a correct
free_ivector(Tage,1,NCOVMAX);	free_ivector(Tage,1,NCOVMAX);

free_imatrix(nbcode,0,NCOVMAX,0,NCOVMAX);	free_imatrix(nbcode,0,NCOVMAX,0,NCOVMAX);
free_imatrix(codtab,1,100,1,10);	/* free_imatrix(codtab,1,100,1,10); */
fflush(fichtm);	fflush(fichtm);
fflush(ficgp);	fflush(ficgp);

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>

Removed from v.1.197
changed lines
	Added in v.1.201