version 1.13, 2002/02/20 17:02:08
|
version 1.16, 2002/02/20 17:12:32
|
Line 73 int npar=NPARMAX;
|
Line 73 int npar=NPARMAX;
|
int nlstate=2; /* Number of live states */
|
int nlstate=2; /* Number of live states */
|
int ndeath=1; /* Number of dead states */
|
int ndeath=1; /* Number of dead states */
|
int ncovmodel, ncov; /* Total number of covariables including constant a12*1 +b12*x ncovmodel=2 */
|
int ncovmodel, ncov; /* Total number of covariables including constant a12*1 +b12*x ncovmodel=2 */
|
|
int popbased=0;
|
|
|
int *wav; /* Number of waves for this individuual 0 is possible */
|
int *wav; /* Number of waves for this individuual 0 is possible */
|
int maxwav; /* Maxim number of waves */
|
int maxwav; /* Maxim number of waves */
|
Line 84 double jmean; /* Mean space between 2 wa
|
Line 85 double jmean; /* Mean space between 2 wa
|
double **oldm, **newm, **savm; /* Working pointers to matrices */
|
double **oldm, **newm, **savm; /* Working pointers to matrices */
|
double **oldms, **newms, **savms; /* Fixed working pointers to matrices */
|
double **oldms, **newms, **savms; /* Fixed working pointers to matrices */
|
FILE *fic,*ficpar, *ficparo,*ficres, *ficrespl, *ficrespij, *ficrest,*ficresf;
|
FILE *fic,*ficpar, *ficparo,*ficres, *ficrespl, *ficrespij, *ficrest,*ficresf;
|
FILE *ficgp, *fichtm,*ficresprob;
|
FILE *ficgp, *fichtm,*ficresprob,*ficpop;
|
FILE *ficreseij;
|
FILE *ficreseij;
|
char filerese[FILENAMELENGTH];
|
char filerese[FILENAMELENGTH];
|
FILE *ficresvij;
|
FILE *ficresvij;
|
Line 1149 void lubksb(double **a, int n, int *indx
|
Line 1150 void lubksb(double **a, int n, int *indx
|
}
|
}
|
|
|
/************ Frequencies ********************/
|
/************ Frequencies ********************/
|
void freqsummary(char fileres[], int agemin, int agemax, int **s, double **agev, int nlstate, int imx, int *Tvar, int **nbcode, int *ncodemax)
|
void freqsummary(char fileres[], int agemin, int agemax, int **s, double **agev, int nlstate, int imx, int *Tvar, int **nbcode, int *ncodemax, int fprev1,int lprev1)
|
{ /* Some frequencies */
|
{ /* Some frequencies */
|
|
|
int i, m, jk, k1, i1, j1, bool, z1,z2,j;
|
int i, m, jk, k1, i1, j1, bool, z1,z2,j;
|
Line 1191 void freqsummary(char fileres[], int ag
|
Line 1192 void freqsummary(char fileres[], int ag
|
bool=0;
|
bool=0;
|
}
|
}
|
if (bool==1) {
|
if (bool==1) {
|
for(m=firstpass; m<=lastpass-1; m++){
|
for(m=fprev1; m<=lprev1; m++){
|
if(agev[m][i]==0) agev[m][i]=agemax+1;
|
if(agev[m][i]==0) agev[m][i]=agemax+1;
|
if(agev[m][i]==1) agev[m][i]=agemax+2;
|
if(agev[m][i]==1) agev[m][i]=agemax+2;
|
freq[s[m][i]][s[m+1][i]][(int)agev[m][i]] += weight[i];
|
freq[s[m][i]][s[m+1][i]][(int)agev[m][i]] += weight[i];
|
Line 1215 void freqsummary(char fileres[], int ag
|
Line 1216 void freqsummary(char fileres[], int ag
|
printf("Age %d", i);
|
printf("Age %d", i);
|
for(jk=1; jk <=nlstate ; jk++){
|
for(jk=1; jk <=nlstate ; jk++){
|
for(m=-1, pp[jk]=0; m <=nlstate+ndeath ; m++)
|
for(m=-1, pp[jk]=0; m <=nlstate+ndeath ; m++)
|
pp[jk] += freq[jk][m][i];
|
pp[jk] += freq[jk][m][i];
|
}
|
}
|
for(jk=1; jk <=nlstate ; jk++){
|
for(jk=1; jk <=nlstate ; jk++){
|
for(m=-1, pos=0; m <=0 ; m++)
|
for(m=-1, pos=0; m <=0 ; m++)
|
Line 1225 void freqsummary(char fileres[], int ag
|
Line 1226 void freqsummary(char fileres[], int ag
|
else
|
else
|
printf(" %d.=%.0f loss[%d]=NaNQ%%",jk,pp[jk],jk);
|
printf(" %d.=%.0f loss[%d]=NaNQ%%",jk,pp[jk],jk);
|
}
|
}
|
for(jk=1; jk <=nlstate ; jk++){
|
|
for(m=1, pp[jk]=0; m <=nlstate+ndeath; m++)
|
for(jk=1; jk <=nlstate ; jk++){
|
|
for(m=0, pp[jk]=0; m <=nlstate+ndeath; m++)
|
pp[jk] += freq[jk][m][i];
|
pp[jk] += freq[jk][m][i];
|
}
|
}
|
|
|
for(jk=1,pos=0; jk <=nlstate ; jk++)
|
for(jk=1,pos=0; jk <=nlstate ; jk++)
|
pos += pp[jk];
|
pos += pp[jk];
|
for(jk=1; jk <=nlstate ; jk++){
|
for(jk=1; jk <=nlstate ; jk++){
|
Line 1262 void freqsummary(char fileres[], int ag
|
Line 1265 void freqsummary(char fileres[], int ag
|
|
|
} /* End of Freq */
|
} /* End of Freq */
|
|
|
|
/************ Prevalence ********************/
|
|
void prevalence(int agemin, int agemax, int **s, double **agev, int nlstate, int imx, int *Tvar, int **nbcode, int *ncodemax, int fprev1,int lprev1)
|
|
{ /* Some frequencies */
|
|
|
|
int i, m, jk, k1, i1, j1, bool, z1,z2,j;
|
|
double ***freq; /* Frequencies */
|
|
double *pp;
|
|
double pos;
|
|
|
|
pp=vector(1,nlstate);
|
|
probs= ma3x(1,130 ,1,8, 1,8);
|
|
|
|
freq=ma3x(-1,nlstate+ndeath,-1,nlstate+ndeath,agemin,agemax+3);
|
|
j1=0;
|
|
|
|
j=cptcoveff;
|
|
if (cptcovn<1) {j=1;ncodemax[1]=1;}
|
|
|
|
for(k1=1; k1<=j;k1++){
|
|
for(i1=1; i1<=ncodemax[k1];i1++){
|
|
j1++;
|
|
|
|
for (i=-1; i<=nlstate+ndeath; i++)
|
|
for (jk=-1; jk<=nlstate+ndeath; jk++)
|
|
for(m=agemin; m <= agemax+3; m++)
|
|
freq[i][jk][m]=0;
|
|
|
|
for (i=1; i<=imx; i++) {
|
|
bool=1;
|
|
if (cptcovn>0) {
|
|
for (z1=1; z1<=cptcoveff; z1++)
|
|
if (covar[Tvaraff[z1]][i]!= nbcode[Tvaraff[z1]][codtab[j1][z1]])
|
|
bool=0;
|
|
}
|
|
if (bool==1) {
|
|
for(m=fprev1; m<=lprev1; m++){
|
|
if(agev[m][i]==0) agev[m][i]=agemax+1;
|
|
if(agev[m][i]==1) agev[m][i]=agemax+2;
|
|
freq[s[m][i]][s[m+1][i]][(int)agev[m][i]] += weight[i];
|
|
freq[s[m][i]][s[m+1][i]][(int) agemax+3] += weight[i];
|
|
}
|
|
}
|
|
}
|
|
for(i=(int)agemin; i <= (int)agemax+3; i++){
|
|
for(jk=1; jk <=nlstate ; jk++){
|
|
for(m=-1, pp[jk]=0; m <=nlstate+ndeath ; m++)
|
|
pp[jk] += freq[jk][m][i];
|
|
}
|
|
for(jk=1; jk <=nlstate ; jk++){
|
|
for(m=-1, pos=0; m <=0 ; m++)
|
|
pos += freq[jk][m][i];
|
|
}
|
|
|
|
for(jk=1; jk <=nlstate ; jk++){
|
|
for(m=0, pp[jk]=0; m <=nlstate+ndeath; m++)
|
|
pp[jk] += freq[jk][m][i];
|
|
}
|
|
|
|
for(jk=1,pos=0; jk <=nlstate ; jk++) pos += pp[jk];
|
|
|
|
for(jk=1; jk <=nlstate ; jk++){
|
|
if( i <= (int) agemax){
|
|
if(pos>=1.e-5){
|
|
probs[i][jk][j1]= pp[jk]/pos;
|
|
}
|
|
}
|
|
}
|
|
|
|
}
|
|
}
|
|
}
|
|
|
|
|
|
free_ma3x(freq,-1,nlstate+ndeath,-1,nlstate+ndeath,(int) agemin,(int) agemax+3);
|
|
free_vector(pp,1,nlstate);
|
|
|
|
} /* End of Freq */
|
/************* Waves Concatenation ***************/
|
/************* Waves Concatenation ***************/
|
|
|
void concatwav(int wav[], int **dh, int **mw, int **s, double *agedc, double **agev, int firstpass, int lastpass, int imx, int nlstate, int stepm)
|
void concatwav(int wav[], int **dh, int **mw, int **s, double *agedc, double **agev, int firstpass, int lastpass, int imx, int nlstate, int stepm)
|
Line 1492 void varevsij(char fileres[], double ***
|
Line 1572 void varevsij(char fileres[], double ***
|
}
|
}
|
hpxij(p3mat,nhstepm,age,hstepm,xp,nlstate,stepm,oldm,savm, ij);
|
hpxij(p3mat,nhstepm,age,hstepm,xp,nlstate,stepm,oldm,savm, ij);
|
prevalim(prlim,nlstate,xp,age,oldm,savm,ftolpl,ij);
|
prevalim(prlim,nlstate,xp,age,oldm,savm,ftolpl,ij);
|
|
|
|
if (popbased==1) {
|
|
for(i=1; i<=nlstate;i++)
|
|
prlim[i][i]=probs[(int)age][i][ij];
|
|
}
|
|
|
for(j=1; j<= nlstate; j++){
|
for(j=1; j<= nlstate; j++){
|
for(h=0; h<=nhstepm; h++){
|
for(h=0; h<=nhstepm; h++){
|
for(i=1, gp[h][j]=0.;i<=nlstate;i++)
|
for(i=1, gp[h][j]=0.;i<=nlstate;i++)
|
Line 1503 void varevsij(char fileres[], double ***
|
Line 1589 void varevsij(char fileres[], double ***
|
xp[i] = x[i] - (i==theta ?delti[theta]:0);
|
xp[i] = x[i] - (i==theta ?delti[theta]:0);
|
hpxij(p3mat,nhstepm,age,hstepm,xp,nlstate,stepm,oldm,savm, ij);
|
hpxij(p3mat,nhstepm,age,hstepm,xp,nlstate,stepm,oldm,savm, ij);
|
prevalim(prlim,nlstate,xp,age,oldm,savm,ftolpl,ij);
|
prevalim(prlim,nlstate,xp,age,oldm,savm,ftolpl,ij);
|
|
|
|
if (popbased==1) {
|
|
for(i=1; i<=nlstate;i++)
|
|
prlim[i][i]=probs[(int)age][i][ij];
|
|
}
|
|
|
for(j=1; j<= nlstate; j++){
|
for(j=1; j<= nlstate; j++){
|
for(h=0; h<=nhstepm; h++){
|
for(h=0; h<=nhstepm; h++){
|
for(i=1, gm[h][j]=0.;i<=nlstate;i++)
|
for(i=1, gm[h][j]=0.;i<=nlstate;i++)
|
gm[h][j] += prlim[i][i]*p3mat[i][j][h];
|
gm[h][j] += prlim[i][i]*p3mat[i][j][h];
|
}
|
}
|
}
|
}
|
|
|
for(j=1; j<= nlstate; j++)
|
for(j=1; j<= nlstate; j++)
|
for(h=0; h<=nhstepm; h++){
|
for(h=0; h<=nhstepm; h++){
|
gradg[h][theta][j]= (gp[h][j]-gm[h][j])/2./delti[theta];
|
gradg[h][theta][j]= (gp[h][j]-gm[h][j])/2./delti[theta];
|
Line 1767 int main()
|
Line 1860 int main()
|
char fileres[FILENAMELENGTH], filerespij[FILENAMELENGTH], filereso[FILENAMELENGTH], fileresf[FILENAMELENGTH];
|
char fileres[FILENAMELENGTH], filerespij[FILENAMELENGTH], filereso[FILENAMELENGTH], fileresf[FILENAMELENGTH];
|
char filerest[FILENAMELENGTH];
|
char filerest[FILENAMELENGTH];
|
char fileregp[FILENAMELENGTH];
|
char fileregp[FILENAMELENGTH];
|
|
char popfile[FILENAMELENGTH];
|
char path[80],pathc[80],pathcd[80],pathtot[80],model[20];
|
char path[80],pathc[80],pathcd[80],pathtot[80],model[20];
|
int firstobs=1, lastobs=10;
|
int firstobs=1, lastobs=10;
|
int sdeb, sfin; /* Status at beginning and end */
|
int sdeb, sfin; /* Status at beginning and end */
|
int c, h , cpt,l;
|
int c, h , cpt,l;
|
int ju,jl, mi;
|
int ju,jl, mi;
|
int i1,j1, k1,k2,k3,jk,aa,bb, stepsize, ij;
|
int i1,j1, k1,k2,k3,jk,aa,bb, stepsize, ij;
|
int jnais,jdc,jint4,jint1,jint2,jint3,**outcome,**adl,*tab;
|
int jnais,jdc,jint4,jint1,jint2,jint3,**outcome,**adl,*tab;
|
|
int mobilav=0, fprev, lprev ,fprevfore=1, lprevfore=1,nforecast,popforecast=0;
|
int hstepm, nhstepm;
|
int hstepm, nhstepm;
|
|
int *popage;
|
|
|
double bage, fage, age, agelim, agebase;
|
double bage, fage, age, agelim, agebase;
|
double ftolpl=FTOL;
|
double ftolpl=FTOL;
|
double **prlim;
|
double **prlim;
|
Line 1788 int main()
|
Line 1884 int main()
|
double ***eij, ***vareij;
|
double ***eij, ***vareij;
|
double **varpl; /* Variances of prevalence limits by age */
|
double **varpl; /* Variances of prevalence limits by age */
|
double *epj, vepp;
|
double *epj, vepp;
|
double kk1;
|
double kk1, kk2;
|
|
double *popeffectif,*popcount;
|
|
|
char version[80]="Imach version 64b, May 2001, INED-EUROREVES ";
|
char version[80]="Imach version 64b, May 2001, INED-EUROREVES ";
|
char *alph[]={"a","a","b","c","d","e"}, str[4];
|
char *alph[]={"a","a","b","c","d","e"}, str[4];
|
Line 1804 int main()
|
Line 1901 int main()
|
gettimeofday(&start_time, (struct timezone*)0); */ /* at first time */
|
gettimeofday(&start_time, (struct timezone*)0); */ /* at first time */
|
|
|
|
|
printf("\nIMACH, Version 0.64b");
|
printf("\nIMACH, Version 0.7");
|
printf("\nEnter the parameter file name: ");
|
printf("\nEnter the parameter file name: ");
|
|
|
#ifdef windows
|
#ifdef windows
|
Line 1852 split(pathtot, path,optionfile);
|
Line 1949 split(pathtot, path,optionfile);
|
fscanf(ficpar,"title=%s datafile=%s lastobs=%d firstpass=%d lastpass=%d\nftol=%lf stepm=%d ncov=%d nlstate=%d ndeath=%d maxwav=%d mle=%d weight=%d\nmodel=%s\n",title, datafile, &lastobs, &firstpass,&lastpass,&ftol, &stepm, &ncov, &nlstate,&ndeath, &maxwav, &mle, &weightopt,model);
|
fscanf(ficpar,"title=%s datafile=%s lastobs=%d firstpass=%d lastpass=%d\nftol=%lf stepm=%d ncov=%d nlstate=%d ndeath=%d maxwav=%d mle=%d weight=%d\nmodel=%s\n",title, datafile, &lastobs, &firstpass,&lastpass,&ftol, &stepm, &ncov, &nlstate,&ndeath, &maxwav, &mle, &weightopt,model);
|
printf("title=%s datafile=%s lastobs=%d firstpass=%d lastpass=%d\nftol=%e stepm=%d ncov=%d nlstate=%d ndeath=%d maxwav=%d mle=%d weight=%d\nmodel=%s\n", title, datafile, lastobs, firstpass,lastpass,ftol, stepm, ncov, nlstate,ndeath, maxwav, mle, weightopt,model);
|
printf("title=%s datafile=%s lastobs=%d firstpass=%d lastpass=%d\nftol=%e stepm=%d ncov=%d nlstate=%d ndeath=%d maxwav=%d mle=%d weight=%d\nmodel=%s\n", title, datafile, lastobs, firstpass,lastpass,ftol, stepm, ncov, nlstate,ndeath, maxwav, mle, weightopt,model);
|
fprintf(ficparo,"title=%s datafile=%s lastobs=%d firstpass=%d lastpass=%d\nftol=%e stepm=%d ncov=%d nlstate=%d ndeath=%d maxwav=%d mle=%d weight=%d\nmodel=%s\n", title, datafile, lastobs, firstpass,lastpass,ftol,stepm,ncov,nlstate,ndeath,maxwav, mle, weightopt,model);
|
fprintf(ficparo,"title=%s datafile=%s lastobs=%d firstpass=%d lastpass=%d\nftol=%e stepm=%d ncov=%d nlstate=%d ndeath=%d maxwav=%d mle=%d weight=%d\nmodel=%s\n", title, datafile, lastobs, firstpass,lastpass,ftol,stepm,ncov,nlstate,ndeath,maxwav, mle, weightopt,model);
|
|
while((c=getc(ficpar))=='#' && c!= EOF){
|
|
ungetc(c,ficpar);
|
|
fgets(line, MAXLINE, ficpar);
|
|
puts(line);
|
|
fputs(line,ficparo);
|
|
}
|
|
ungetc(c,ficpar);
|
|
|
|
fscanf(ficpar,"fprevalence=%d lprevalence=%d pop_based=%d\n",&fprev,&lprev,&popbased);
|
|
fprintf(ficparo,"fprevalence=%d lprevalence=%d pop_based=%d\n",fprev,lprev,popbased);
|
|
|
|
while((c=getc(ficpar))=='#' && c!= EOF){
|
|
ungetc(c,ficpar);
|
|
fgets(line, MAXLINE, ficpar);
|
|
puts(line);
|
|
fputs(line,ficparo);
|
|
}
|
|
ungetc(c,ficpar);
|
|
|
|
fscanf(ficpar,"fprevalence=%d lprevalence=%d nforecast=%d mob_average=%d\n",&fprevfore,&lprevfore,&nforecast,&mobilav);
|
|
fprintf(ficparo,"fprevalence=%d lprevalence=%d nforecast=%d mob_average=%d\n",fprevfore,lprevfore,nforecast,mobilav);
|
|
|
|
|
|
while((c=getc(ficpar))=='#' && c!= EOF){
|
|
ungetc(c,ficpar);
|
|
fgets(line, MAXLINE, ficpar);
|
|
puts(line);
|
|
fputs(line,ficparo);
|
|
}
|
|
ungetc(c,ficpar);
|
|
|
|
fscanf(ficpar,"popforecast=%d popfile=%s\n",&popforecast,popfile);
|
|
|
covar=matrix(0,NCOVMAX,1,n);
|
covar=matrix(0,NCOVMAX,1,n);
|
cptcovn=0;
|
cptcovn=0;
|
if (strlen(model)>1) cptcovn=nbocc(model,'+')+1;
|
if (strlen(model)>1) cptcovn=nbocc(model,'+')+1;
|
Line 2014 split(pathtot, path,optionfile);
|
Line 2143 split(pathtot, path,optionfile);
|
if ((s[1][i]==3) && (s[2][i]==2)) s[2][i]=3;
|
if ((s[1][i]==3) && (s[2][i]==2)) s[2][i]=3;
|
if ((s[2][i]==3) && (s[3][i]==2)) s[3][i]=3;
|
if ((s[2][i]==3) && (s[3][i]==2)) s[3][i]=3;
|
if ((s[3][i]==3) && (s[4][i]==2)) s[4][i]=3;
|
if ((s[3][i]==3) && (s[4][i]==2)) s[4][i]=3;
|
}
|
}
|
for (i=1; i<=imx; i++) printf("%d %.lf %.lf %.lf %.lf/%.lf %.lf/%.lf %.lf/%.lf %d %.lf/%.lf %d %.lf/%.lf %d %.lf/%.lf %d\n",num[i],(covar[1][i]), (covar[2][i]), (weight[i]), (moisnais[i]), (annais[i]), (moisdc[i]), (andc[i]), (mint[1][i]), (anint[1][i]), (s[1][i]), (mint[2][i]), (anint[2][i]), (s[2][i]), (mint[3][i]), (anint[3][i]), (s[3][i]), (mint[4][i]), (anint[4][i]), (s[4][i]));*/
|
for (i=1; i<=imx; i++) printf("%d %.lf %.lf %.lf %.lf/%.lf %.lf/%.lf %.lf/%.lf %d %.lf/%.lf %d %.lf/%.lf %d %.lf/%.lf %d\n",num[i],(covar[1][i]), (covar[2][i]), (weight[i]), (moisnais[i]), (annais[i]), (moisdc[i]), (andc[i]), (mint[1][i]), (anint[1][i]), (s[1][i]), (mint[2][i]), (anint[2][i]), (s[2][i]), (mint[3][i]), (anint[3][i]), (s[3][i]), (mint[4][i]), (anint[4][i]), (s[4][i]));*/
|
|
|
/* Calculation of the number of parameter from char model*/
|
/* Calculation of the number of parameter from char model*/
|
Tvar=ivector(1,15);
|
Tvar=ivector(1,15);
|
Line 2206 printf("Total number of individuals= %d,
|
Line 2335 printf("Total number of individuals= %d,
|
|
|
/* Calculates basic frequencies. Computes observed prevalence at single age
|
/* Calculates basic frequencies. Computes observed prevalence at single age
|
and prints on file fileres'p'. */
|
and prints on file fileres'p'. */
|
freqsummary(fileres, agemin, agemax, s, agev, nlstate, imx,Tvar,nbcode, ncodemax);
|
freqsummary(fileres, agemin, agemax, s, agev, nlstate, imx,Tvar,nbcode, ncodemax, fprev, lprev);
|
|
|
pmmij= matrix(1,nlstate+ndeath,1,nlstate+ndeath); /* creation */
|
pmmij= matrix(1,nlstate+ndeath,1,nlstate+ndeath); /* creation */
|
oldms= matrix(1,nlstate+ndeath,1,nlstate+ndeath); /* creation */
|
oldms= matrix(1,nlstate+ndeath,1,nlstate+ndeath); /* creation */
|
Line 2223 printf("Total number of individuals= %d,
|
Line 2352 printf("Total number of individuals= %d,
|
}
|
}
|
|
|
/*--------- results files --------------*/
|
/*--------- results files --------------*/
|
fprintf(ficres,"\ntitle=%s datafile=%s lastobs=%d firstpass=%d lastpass=%d\nftol=%e stepm=%d ncov=%d nlstate=%d ndeath=%d maxwav=%d mle=%d weight=%d\nmodel=%s\n", title, datafile, lastobs, firstpass,lastpass,ftol, stepm, ncov, nlstate, ndeath, maxwav, mle,weightopt,model);
|
fprintf(ficres,"title=%s datafile=%s lastobs=%d firstpass=%d lastpass=%d\nftol=%e stepm=%d ncov=%d nlstate=%d ndeath=%d maxwav=%d mle=%d weight=%d\nmodel=%s\n", title, datafile, lastobs, firstpass,lastpass,ftol, stepm, ncov, nlstate, ndeath, maxwav, mle,weightopt,model);
|
|
fprintf(ficres,"fprevalence=%d lprevalence=%d pop_based=%d\n",fprev,lprev,popbased);
|
|
fprintf(ficres,"fprevalence=%d lprevalence=%d nforecast=%d mob_average=%d\n",fprevfore,lprevfore,nforecast,mobilav);
|
|
|
jk=1;
|
jk=1;
|
fprintf(ficres,"# Parameters\n");
|
fprintf(ficres,"# Parameters\n");
|
printf("# Parameters\n");
|
printf("# Parameters\n");
|
Line 2470 ij=1;
|
Line 2601 ij=1;
|
fclose(ficgp);
|
fclose(ficgp);
|
|
|
chdir(path);
|
chdir(path);
|
free_matrix(agev,1,maxwav,1,imx);
|
|
free_ivector(wav,1,imx);
|
free_ivector(wav,1,imx);
|
free_imatrix(dh,1,lastpass-firstpass+1,1,imx);
|
free_imatrix(dh,1,lastpass-firstpass+1,1,imx);
|
free_imatrix(mw,1,lastpass-firstpass+1,1,imx);
|
free_imatrix(mw,1,lastpass-firstpass+1,1,imx);
|
|
|
free_imatrix(s,1,maxwav+1,1,n);
|
|
|
|
|
|
free_ivector(num,1,n);
|
free_ivector(num,1,n);
|
free_vector(agedc,1,n);
|
free_vector(agedc,1,n);
|
free_vector(weight,1,n);
|
|
/*free_matrix(covar,1,NCOVMAX,1,n);*/
|
/*free_matrix(covar,1,NCOVMAX,1,n);*/
|
fclose(ficparo);
|
fclose(ficparo);
|
fclose(ficres);
|
fclose(ficres);
|
Line 2511 chdir(path);
|
Line 2637 chdir(path);
|
printf("Problem with %s \n",optionfilehtm);goto end;
|
printf("Problem with %s \n",optionfilehtm);goto end;
|
}
|
}
|
|
|
fprintf(fichtm,"<body><ul> <font size=\"6\">Imach, Version 0.64b </font> <hr size=\"2\" color=\"#EC5E5E\">
|
fprintf(fichtm,"<body><ul> <font size=\"6\">Imach, Version 0.7 </font> <hr size=\"2\" color=\"#EC5E5E\">
|
Titre=%s <br>Datafile=%s Firstpass=%d Lastpass=%d Stepm=%d Weight=%d Model=%s<br>
|
Titre=%s <br>Datafile=%s Firstpass=%d Lastpass=%d Stepm=%d Weight=%d Model=%s<br>
|
Total number of observations=%d <br>
|
Total number of observations=%d <br>
|
Interval (in months) between two waves: Min=%d Max=%d Mean=%.2lf<br>
|
Interval (in months) between two waves: Min=%d Max=%d Mean=%.2lf<br>
|
Line 2525 Interval (in months) between two waves:
|
Line 2651 Interval (in months) between two waves:
|
- Life expectancies by age and initial health status: <a href=\"e%s\">e%s</a> <br>
|
- Life expectancies by age and initial health status: <a href=\"e%s\">e%s</a> <br>
|
- Variances of life expectancies by age and initial health status: <a href=\"v%s\">v%s</a><br>
|
- Variances of life expectancies by age and initial health status: <a href=\"v%s\">v%s</a><br>
|
- Health expectancies with their variances: <a href=\"t%s\">t%s</a> <br>
|
- Health expectancies with their variances: <a href=\"t%s\">t%s</a> <br>
|
- Standard deviation of stationary prevalences: <a href=\"vpl%s\">vpl%s</a> <br><br>",title,datafile,firstpass,lastpass,stepm, weightopt,model,imx,jmin,jmax,jmean,fileres,fileres,fileres,fileres,fileres,fileres,fileres,fileres,fileres,fileres,fileres,fileres,fileres,fileres,fileres,fileres,fileres,fileres);
|
- Standard deviation of stationary prevalences: <a href=\"vpl%s\">vpl%s</a> <br>
|
|
- Prevalences forecasting: <a href=\"f%s\">f%s</a> <br>
|
|
<br>",title,datafile,firstpass,lastpass,stepm, weightopt,model,imx,jmin,jmax,jmean,fileres,fileres,fileres,fileres,fileres,fileres,fileres,fileres,fileres,fileres,fileres,fileres,fileres,fileres,fileres,fileres,fileres,fileres,fileres,fileres);
|
|
|
fprintf(fichtm," <li>Graphs</li><p>");
|
fprintf(fichtm," <li>Graphs</li><p>");
|
|
|
Line 2663 fclose(fichtm);
|
Line 2791 fclose(fichtm);
|
|
|
fclose(ficrespij);
|
fclose(ficrespij);
|
|
|
exit(0);
|
|
/*---------- Forecasting ------------------*/
|
/*---------- Forecasting ------------------*/
|
|
|
strcpy(fileresf,"f");
|
strcpy(fileresf,"f");
|
Line 2673 fclose(fichtm);
|
Line 2800 fclose(fichtm);
|
}
|
}
|
printf("Computing forecasting: result on file '%s' \n", fileresf);
|
printf("Computing forecasting: result on file '%s' \n", fileresf);
|
|
|
/* Mobile average */
|
prevalence(agemin, agemax, s, agev, nlstate, imx,Tvar,nbcode, ncodemax, fprevfore, lprevfore);
|
|
|
/* for (agedeb=bage; agedeb<=fage; agedeb++)
|
free_matrix(agev,1,maxwav,1,imx);
|
for (i=1; i<=nlstate;i++)
|
/* Mobile average */
|
for (cptcod=1;cptcod<=ncodemax[cptcoveff];cptcod++)
|
|
printf("%f %d i=%d j1=%d\n", probs[(int)agedeb][i][cptcod],(int) agedeb,i,cptcod);*/
|
|
|
|
if (cptcoveff==0) ncodemax[cptcoveff]=1;
|
if (cptcoveff==0) ncodemax[cptcoveff]=1;
|
|
|
mobaverage= ma3x(1,130 ,1,8, 1,8);
|
if (mobilav==1) {
|
for (agedeb=bage+3; agedeb<=fage-2; agedeb++)
|
mobaverage= ma3x(1, AGESUP,1,NCOVMAX, 1,NCOVMAX);
|
for (i=1; i<=nlstate;i++)
|
for (agedeb=bage+3; agedeb<=fage-2; agedeb++)
|
for (cptcod=1;cptcod<=ncodemax[cptcov];cptcod++)
|
for (i=1; i<=nlstate;i++)
|
mobaverage[(int)agedeb][i][cptcod]=0.;
|
for (cptcod=1;cptcod<=ncodemax[cptcov];cptcod++)
|
|
mobaverage[(int)agedeb][i][cptcod]=0.;
|
for (agedeb=bage+4; agedeb<=fage; agedeb++){
|
|
for (i=1; i<=nlstate;i++){
|
for (agedeb=bage+4; agedeb<=fage; agedeb++){
|
for (cptcod=1;cptcod<=ncodemax[cptcoveff];cptcod++){
|
for (i=1; i<=nlstate;i++){
|
for (cpt=0;cpt<=4;cpt++){
|
for (cptcod=1;cptcod<=ncodemax[cptcoveff];cptcod++){
|
mobaverage[(int)agedeb-2][i][cptcod]=mobaverage[(int)agedeb-2][i][cptcod]+probs[(int)agedeb-cpt][i][cptcod];
|
for (cpt=0;cpt<=4;cpt++){
|
|
mobaverage[(int)agedeb-2][i][cptcod]=mobaverage[(int)agedeb-2][i][cptcod]+probs[(int)agedeb-cpt][i][cptcod];
|
}
|
}
|
mobaverage[(int)agedeb-2][i][cptcod]=mobaverage[(int)agedeb-2][i][cptcod]/5;
|
mobaverage[(int)agedeb-2][i][cptcod]=mobaverage[(int)agedeb-2][i][cptcod]/5;
|
|
}
|
}
|
}
|
}
|
}
|
}
|
}
|
|
|
/* if (cptcod==2) printf("m=%f p=%f %d age=%d ",mobaverage[(int)agedeb-2][i][cptcod],probs[(int)agedeb-cpt][i][cptcod],cpt,(int)agedeb-2);*/
|
|
|
|
|
|
stepsize=(int) (stepm+YEARM-1)/YEARM;
|
stepsize=(int) (stepm+YEARM-1)/YEARM;
|
if (stepm<=24) stepsize=2;
|
if (stepm<=12) stepsize=1;
|
|
|
agelim=AGESUP;
|
agelim=AGESUP;
|
hstepm=stepsize*YEARM; /* Every year of age */
|
hstepm=stepsize*YEARM; /* Every year of age */
|
hstepm=hstepm/stepm; /* Typically 2 years, = 2/6 months = 4 */
|
hstepm=hstepm/stepm; /* Typically 2 years, = 2 years/6 months = 4 */
|
hstepm=12;
|
|
k=0;
|
if (popforecast==1) {
|
|
if((ficpop=fopen(popfile,"r"))==NULL) {
|
|
printf("Problem with population file : %s\n",popfile);goto end;
|
|
}
|
|
popage=ivector(0,AGESUP);
|
|
popeffectif=vector(0,AGESUP);
|
|
popcount=vector(0,AGESUP);
|
|
|
|
i=1;
|
|
while ((c=fscanf(ficpop,"%d %lf\n",&popage[i],&popcount[i])) != EOF)
|
|
{
|
|
i=i+1;
|
|
}
|
|
imx=i;
|
|
|
|
for (i=1; i<imx;i++) popeffectif[popage[i]]=popcount[i];
|
|
}
|
|
|
for(cptcov=1;cptcov<=i1;cptcov++){
|
for(cptcov=1;cptcov<=i1;cptcov++){
|
for(cptcod=1;cptcod<=ncodemax[cptcoveff];cptcod++){
|
for(cptcod=1;cptcod<=ncodemax[cptcoveff];cptcod++){
|
k=k+1;
|
k=k+1;
|
Line 2717 fclose(fichtm);
|
Line 2858 fclose(fichtm);
|
for(j=1;j<=cptcoveff;j++) {
|
for(j=1;j<=cptcoveff;j++) {
|
fprintf(ficresf,"V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtab[k][j]]);
|
fprintf(ficresf,"V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtab[k][j]]);
|
}
|
}
|
|
|
fprintf(ficresf,"******\n");
|
fprintf(ficresf,"******\n");
|
|
|
fprintf(ficresf,"# StartingAge FinalAge Horizon(in years)");
|
fprintf(ficresf,"# StartingAge FinalAge Horizon(in years)");
|
for(j=1; j<=nlstate+ndeath;j++) fprintf(ficresf," P.%d",j);
|
for(j=1; j<=nlstate+ndeath;j++) fprintf(ficresf," P.%d",j);
|
|
if (popforecast==1) fprintf(ficresf," [Population]");
|
|
|
for (agedeb=fage; agedeb>=bage; agedeb--){
|
for (agedeb=fage; agedeb>=bage; agedeb--){
|
fprintf(ficresf,"\n%d %.f %.f 0 ",k,agedeb, agedeb);
|
fprintf(ficresf,"\n%.f %.f 0",agedeb, agedeb);
|
|
if (mobilav==1) {
|
for(j=1; j<=nlstate;j++)
|
for(j=1; j<=nlstate;j++)
|
fprintf(ficresf,"%.3f ",mobaverage[(int)agedeb][j][cptcod]);
|
fprintf(ficresf," %.3f",mobaverage[(int)agedeb][j][cptcod]);
|
}
|
}
|
for(j=1; j<=ndeath;j++) fprintf(ficresf,"0.");
|
else {
|
|
for(j=1; j<=nlstate;j++)
|
|
fprintf(ficresf," %.3f",probs[(int)agedeb][j][cptcod]);
|
|
}
|
|
|
for (cpt=1; cpt<=8;cpt++)
|
for(j=1; j<=ndeath;j++) fprintf(ficresf," 0.00000");
|
|
if (popforecast==1) fprintf(ficresf," [%.f] ",popeffectif[(int)agedeb]);
|
|
}
|
|
|
|
for (cpt=1; cpt<=nforecast;cpt++) {
|
|
fprintf(ficresf,"\n");
|
for (agedeb=fage; agedeb>=bage; agedeb--){ /* If stepm=6 months */
|
for (agedeb=fage; agedeb>=bage; agedeb--){ /* If stepm=6 months */
|
|
nhstepm=(int) rint((agelim-agedeb)*YEARM/stepm);
|
nhstepm=(int) rint((agelim-agedeb)*YEARM/stepm); /* Typically 20 years = 20*12/6=40 */
|
nhstepm = nhstepm/hstepm;
|
nhstepm = nhstepm/hstepm; /* Typically 40/4=10 */
|
/*printf("agedeb=%.lf stepm=%d hstepm=%d nhstepm=%d \n",agedeb,stepm,hstepm,nhstepm);*/
|
/*printf("stepm=%d hstepm=%d nhstepm=%d \n",stepm,hstepm,nhstepm);*/
|
|
|
|
p3mat=ma3x(1,nlstate+ndeath,1, nlstate+ndeath, 0,nhstepm);
|
p3mat=ma3x(1,nlstate+ndeath,1, nlstate+ndeath, 0,nhstepm);
|
oldm=oldms;savm=savms;
|
oldm=oldms;savm=savms;
|
Line 2743 fclose(fichtm);
|
Line 2891 fclose(fichtm);
|
|
|
for (h=0; h<=nhstepm; h++){
|
for (h=0; h<=nhstepm; h++){
|
|
|
if (h*hstepm/YEARM*stepm==cpt)
|
if (h*hstepm/YEARM*stepm==cpt)
|
fprintf(ficresf,"\n%d %.f %.f %.f",k,agedeb, agedeb+ h*hstepm/YEARM*stepm, h*hstepm/YEARM*stepm);
|
fprintf(ficresf,"\n%.f %.f %.f",agedeb, agedeb+ h*hstepm/YEARM*stepm, h*hstepm/YEARM*stepm);
|
|
|
|
|
for(j=1; j<=nlstate+ndeath;j++) {
|
for(j=1; j<=nlstate+ndeath;j++) {
|
kk1=0.;
|
kk1=0.;kk2=0;
|
for(i=1; i<=nlstate;i++) {
|
for(i=1; i<=nlstate;i++) {
|
/* kk1=kk1+p3mat[i][j][h]*probs[(int)agedeb][i][cptcod];*/
|
if (mobilav==1)
|
kk1=kk1+p3mat[i][j][h]*mobaverage[(int)agedeb][i][cptcod];
|
kk1=kk1+p3mat[i][j][h]*mobaverage[(int)agedeb][i][cptcod];
|
|
else kk1=kk1+p3mat[i][j][h]*probs[(int)agedeb][i][cptcod];
|
|
if (popforecast==1) kk2=kk1*popeffectif[(int)agedeb];
|
}
|
}
|
|
if (h*hstepm/YEARM*stepm==cpt) {
|
if (h*hstepm/YEARM*stepm==cpt)
|
fprintf(ficresf," %.3f", kk1);
|
fprintf(ficresf," %.5f ", kk1);
|
if (popforecast==1) fprintf(ficresf," [%.f]", kk2);
|
}
|
}
|
}
|
}
|
free_ma3x(p3mat,1,nlstate+ndeath,1, nlstate+ndeath, 0,nhstepm);
|
|
}
|
}
|
|
free_ma3x(p3mat,1,nlstate+ndeath,1, nlstate+ndeath, 0,nhstepm);
|
|
|
|
}
|
}
|
}
|
}
|
}
|
|
}
|
|
if (mobilav==1) free_ma3x(mobaverage,1, AGESUP,1,NCOVMAX, 1,NCOVMAX);
|
|
if (popforecast==1) {
|
|
free_ivector(popage,0,AGESUP);
|
|
free_vector(popeffectif,0,AGESUP);
|
|
free_vector(popcount,0,AGESUP);
|
|
}
|
|
free_imatrix(s,1,maxwav+1,1,n);
|
|
free_vector(weight,1,n);
|
fclose(ficresf);
|
fclose(ficresf);
|
|
|
/*---------- Health expectancies and variances ------------*/
|
/*---------- Health expectancies and variances ------------*/
|
|
|
strcpy(filerest,"t");
|
strcpy(filerest,"t");
|
Line 2822 fclose(fichtm);
|
Line 2983 fclose(fichtm);
|
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);
|
|
if (popbased==1) {
|
|
for(i=1; i<=nlstate;i++)
|
|
prlim[i][i]=probs[(int)age][i][k];
|
|
}
|
|
|
fprintf(ficrest," %.0f",age);
|
fprintf(ficrest," %.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++) {
|