version 1.28, 2002/03/01 17:56:23
|
version 1.29, 2002/03/06 19:02:50
|
Line 1556 void evsij(char fileres[], double ***eij
|
Line 1556 void evsij(char fileres[], double ***eij
|
|
|
for(i=1; i<=nlstate;i++)
|
for(i=1; i<=nlstate;i++)
|
for(j=1; j<=nlstate;j++)
|
for(j=1; j<=nlstate;j++)
|
for (h=0, eij[i][j][(int)age]=0; h<=nhstepm; h++){
|
for (h=0, eij[i][j][(int)age]=0; h<=nhstepm-1; h++){
|
eij[i][j][(int)age] +=p3mat[i][j][h];
|
eij[i][j][(int)age] +=(p3mat[i][j][h]+p3mat[i][j][h+1])/2.0;
|
}
|
}
|
|
|
hf=1;
|
hf=1;
|
if (stepm >= YEARM) hf=stepm/YEARM;
|
if (stepm >= YEARM) hf=stepm/YEARM;
|
fprintf(ficreseij,"%.0f",age );
|
|
for(i=1; i<=nlstate;i++)
|
for(i=1; i<=nlstate;i++)
|
for(j=1; j<=nlstate;j++){
|
for(j=1; j<=nlstate;j++){
|
fprintf(ficreseij," %.4f", hf*eij[i][j][(int)age]);
|
if (j==i) eij[i][j][(int)age]= (eij[i][j][(int)age]-0.5*stepm/12./hf);
|
|
}
|
|
|
|
fprintf(ficreseij,"%3.0f",age );
|
|
for(i=1; i<=nlstate;i++)
|
|
for(j=1; j<=nlstate;j++){
|
|
fprintf(ficreseij," %9.4f", hf*eij[i][j][(int)age]);
|
}
|
}
|
fprintf(ficreseij,"\n");
|
fprintf(ficreseij,"\n");
|
free_ma3x(p3mat,1,nlstate+ndeath,1, nlstate+ndeath, 0,nhstepm);
|
free_ma3x(p3mat,1,nlstate+ndeath,1, nlstate+ndeath, 0,nhstepm);
|
Line 1882 fclose(ficresprob);
|
Line 1888 fclose(ficresprob);
|
}
|
}
|
|
|
/******************* Printing html file ***********/
|
/******************* Printing html file ***********/
|
void printinghtml(char fileres[], char title[], char datafile[], int firstpass, int lastpass, int stepm, int weightopt, char model[],int imx,int jmin, int jmax, double jmeanint,char optionfile[],char optionfilehtm[] ){
|
void printinghtml(char fileres[], char title[], char datafile[], int firstpass, int lastpass, int stepm, int weightopt, char model[],int imx,int jmin, int jmax, double jmeanint,char optionfile[],char optionfilehtm[],char rfileres[] ){
|
int jj1, k1, i1, cpt;
|
int jj1, k1, i1, cpt;
|
FILE *fichtm;
|
FILE *fichtm;
|
/*char optionfilehtm[FILENAMELENGTH];*/
|
/*char optionfilehtm[FILENAMELENGTH];*/
|
Line 1893 void printinghtml(char fileres[], char t
|
Line 1899 void printinghtml(char fileres[], char t
|
printf("Problem with %s \n",optionfilehtm), exit(0);
|
printf("Problem with %s \n",optionfilehtm), exit(0);
|
}
|
}
|
|
|
fprintf(fichtm,"<body><ul> <font size=\"6\">Imach, Version 0.7 </font> <hr size=\"2\" color=\"#EC5E5E\">
|
fprintf(fichtm,"<body><ul> <font size=\"6\">Imach, Version 0.71 </font> <hr size=\"2\" color=\"#EC5E5E\">
|
Title=%s <br>Datafile=%s Firstpass=%d Lastpass=%d Stepm=%d Weight=%d Model=%s<br>
|
Title=%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>
|
Line 1911 Interval (in months) between two waves:
|
Line 1917 Interval (in months) between two waves:
|
- Standard deviation of stationary prevalences: <a href=\"vpl%s\">vpl%s</a> <br>
|
- Standard deviation of stationary prevalences: <a href=\"vpl%s\">vpl%s</a> <br>
|
- Prevalences forecasting: <a href=\"f%s\">f%s</a> <br>
|
- Prevalences forecasting: <a href=\"f%s\">f%s</a> <br>
|
- Population forecasting (if popforecast=1): <a href=\"pop%s\">pop%s</a> <br>
|
- Population forecasting (if popforecast=1): <a href=\"pop%s\">pop%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,fileres,fileres);
|
<br>",title,datafile,firstpass,lastpass,stepm, weightopt,model,imx,jmin,jmax,jmean,fileres,fileres,rfileres,rfileres,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 2424 int main(int argc, char *argv[])
|
Line 2430 int main(int argc, char *argv[])
|
char optionfile[FILENAMELENGTH], datafile[FILENAMELENGTH], filerespl[FILENAMELENGTH];
|
char optionfile[FILENAMELENGTH], datafile[FILENAMELENGTH], filerespl[FILENAMELENGTH];
|
char optionfilext[10], optionfilefiname[FILENAMELENGTH], optionfilehtm[FILENAMELENGTH], optionfilegnuplot[FILENAMELENGTH], plotcmd[FILENAMELENGTH];
|
char optionfilext[10], optionfilefiname[FILENAMELENGTH], optionfilehtm[FILENAMELENGTH], optionfilegnuplot[FILENAMELENGTH], plotcmd[FILENAMELENGTH];
|
|
|
char fileres[FILENAMELENGTH], filerespij[FILENAMELENGTH], filereso[FILENAMELENGTH];
|
char fileres[FILENAMELENGTH], filerespij[FILENAMELENGTH], filereso[FILENAMELENGTH], rfileres[FILENAMELENGTH];
|
|
|
char filerest[FILENAMELENGTH];
|
char filerest[FILENAMELENGTH];
|
char fileregp[FILENAMELENGTH];
|
char fileregp[FILENAMELENGTH];
|
Line 2456 int main(int argc, char *argv[])
|
Line 2462 int main(int argc, char *argv[])
|
double dateprev1, dateprev2,jproj1,mproj1,anproj1,jproj2,mproj2,anproj2;
|
double dateprev1, dateprev2,jproj1,mproj1,anproj1,jproj2,mproj2,anproj2;
|
|
|
|
|
char version[80]="Imach version 0.7, February 2002, INED-EUROREVES ";
|
char version[80]="Imach version 0.71, February 2002, INED-EUROREVES ";
|
char *alph[]={"a","a","b","c","d","e"}, str[4];
|
char *alph[]={"a","a","b","c","d","e"}, str[4];
|
|
|
|
|
Line 2624 while((c=getc(ficpar))=='#' && c!= EOF){
|
Line 2630 while((c=getc(ficpar))=='#' && c!= EOF){
|
printf("\n");
|
printf("\n");
|
|
|
|
|
/*-------- data file ----------*/
|
/*-------- Rewriting paramater file ----------*/
|
if((ficres =fopen(fileres,"w"))==NULL) {
|
strcpy(rfileres,"r"); /* "Rparameterfile */
|
printf("Problem with resultfile: %s\n", fileres);goto end;
|
strcat(rfileres,optionfilefiname); /* Parameter file first name*/
|
|
strcat(rfileres,"."); /* */
|
|
strcat(rfileres,optionfilext); /* Other files have txt extension */
|
|
if((ficres =fopen(rfileres,"w"))==NULL) {
|
|
printf("Problem writing new parameter file: %s\n", fileres);goto end;
|
}
|
}
|
fprintf(ficres,"#%s\n",version);
|
fprintf(ficres,"#%s\n",version);
|
|
|
|
/*-------- data file ----------*/
|
if((fic=fopen(datafile,"r"))==NULL) {
|
if((fic=fopen(datafile,"r"))==NULL) {
|
printf("Problem with datafile: %s\n", datafile);goto end;
|
printf("Problem with datafile: %s\n", datafile);goto end;
|
}
|
}
|
Line 2902 printf("Total number of individuals= %d,
|
Line 2913 printf("Total number of individuals= %d,
|
}
|
}
|
|
|
/*--------- results files --------------*/
|
/*--------- results files --------------*/
|
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,"title=%s datafile=%s lastobs=%d firstpass=%d lastpass=%d\nftol=%e stepm=%d ncov=%d nlstate=%d ndeath=%d maxwav=%d mle= 0 weight=%d\nmodel=%s\n", title, datafile, lastobs, firstpass,lastpass,ftol, stepm, ncov, nlstate, ndeath, maxwav, weightopt,model);
|
|
|
|
|
jk=1;
|
jk=1;
|
Line 3057 while((c=getc(ficpar))=='#' && c!= EOF){
|
Line 3068 while((c=getc(ficpar))=='#' && c!= EOF){
|
|
|
/*--------- index.htm --------*/
|
/*--------- index.htm --------*/
|
|
|
printinghtml(fileres,title,datafile, firstpass, lastpass, stepm, weightopt,model,imx,jmin,jmax,jmean,optionfile,optionfilehtm);
|
printinghtml(fileres,title,datafile, firstpass, lastpass, stepm, weightopt,model,imx,jmin,jmax,jmean,optionfile,optionfilehtm,rfileres);
|
|
|
|
|
/*--------------- Prevalence limit --------------*/
|
/*--------------- Prevalence limit --------------*/
|