version 1.77, 2003/05/16 16:49:47
|
version 1.83, 2003/06/10 13:39:11
|
Line 1
|
Line 1
|
/* $Id$ |
/* $Id$ |
|
$State$ |
|
$Log$ |
|
Revision 1.83 2003/06/10 13:39:11 lievre |
|
*** empty log message *** |
|
|
|
Revision 1.82 2003/06/05 15:57:20 brouard |
|
Add log in imach.c and fullversion number is now printed. |
|
|
|
*/ |
|
/* |
Interpolated Markov Chain |
Interpolated Markov Chain |
|
|
Short summary of the programme: |
Short summary of the programme: |
Line 119
|
Line 129
|
#define ODIRSEPARATOR '\\' |
#define ODIRSEPARATOR '\\' |
#endif |
#endif |
|
|
char version[80]="Imach version 0.95a, May 2003, INED-EUROREVES "; |
/* $Id$ */ |
|
/* $State$ */ |
|
|
|
char version[]="Imach version 0.95a1, June 2003, INED-EUROREVES "; |
|
char fullversion[]="$Revision$ $Date$"; |
int erreur; /* Error number */ |
int erreur; /* Error number */ |
int nvar; |
int nvar; |
int cptcovn=0, cptcovage=0, cptcoveff=0,cptcov; |
int cptcovn=0, cptcovage=0, cptcoveff=0,cptcov; |
Line 1799 void concatwav(int wav[], int **dh, int
|
Line 1813 void concatwav(int wav[], int **dh, int
|
sum=sum+j; |
sum=sum+j; |
/*if (j<0) printf("j=%d num=%d \n",j,i);*/ |
/*if (j<0) printf("j=%d num=%d \n",j,i);*/ |
/* printf("%d %d %d %d\n", s[mw[mi][i]][i] ,s[mw[mi+1][i]][i],j,i);*/ |
/* printf("%d %d %d %d\n", s[mw[mi][i]][i] ,s[mw[mi+1][i]][i],j,i);*/ |
if(j<0)printf("Warning! Negative delay (%d) between waves %d and %d of individual at line %d who is aged %.1f with statuses %d %d\n ",j,mw[mi][i],mw[mi+1][i], i,agev[mw[mi][i]][i],s[mw[mi][i]][i] ,s[mw[mi+1][i]][i]); |
if(j<0)printf("Error! Negative delay (%d to death) between waves %d and %d of individual %d at line %d who is aged %.1f with statuses from %d to %d\n ",j,mw[mi][i],mw[mi+1][i],num[i], i,agev[mw[mi][i]][i],s[mw[mi][i]][i] ,s[mw[mi+1][i]][i]); |
} |
} |
} |
} |
else{ |
else{ |
Line 1810 void concatwav(int wav[], int **dh, int
|
Line 1824 void concatwav(int wav[], int **dh, int
|
else if (j <= jmin)jmin=j; |
else if (j <= jmin)jmin=j; |
/* if (j<10) printf("j=%d jmin=%d num=%d ",j,jmin,i); */ |
/* if (j<10) printf("j=%d jmin=%d num=%d ",j,jmin,i); */ |
/*printf("%d %lf %d %d %d\n", i,agev[mw[mi][i]][i],j,s[mw[mi][i]][i] ,s[mw[mi+1][i]][i]);*/ |
/*printf("%d %lf %d %d %d\n", i,agev[mw[mi][i]][i],j,s[mw[mi][i]][i] ,s[mw[mi+1][i]][i]);*/ |
if(j<0)printf("Warning! Negative delay (%d to death) between waves %d and %d of individual at line %d who is aged %.1f with statuses %d %d\n ",j,mw[mi][i],mw[mi+1][i], i,agev[mw[mi][i]][i],s[mw[mi][i]][i] ,s[mw[mi+1][i]][i]); |
if(j<0)printf("Error! Negative delay (%d) between waves %d and %d of individual %d at line %d who is aged %.1f with statuses from %d to %d\n ",j,mw[mi][i],mw[mi+1][i],num[i], i,agev[mw[mi][i]][i],s[mw[mi][i]][i] ,s[mw[mi+1][i]][i]); |
sum=sum+j; |
sum=sum+j; |
} |
} |
jk= j/stepm; |
jk= j/stepm; |
Line 2987 m=pow(2,cptcoveff);
|
Line 3001 m=pow(2,cptcoveff);
|
fprintf(ficgp,"\nset out \"p%s%d%d.png\" \n",strtok(optionfile, "."),cpt,k1); |
fprintf(ficgp,"\nset out \"p%s%d%d.png\" \n",strtok(optionfile, "."),cpt,k1); |
fprintf(ficgp,"set xlabel \"Age\" \nset ylabel \"Probability\" \nset ter png small\nset size 0.65,0.65\nplot [%.f:%.f] \"pij%s\" u ($1==%d ? ($3):1/0):($%d/($%d",ageminpar,agemaxpar,fileres,k1,k+cpt+1,k+1); |
fprintf(ficgp,"set xlabel \"Age\" \nset ylabel \"Probability\" \nset ter png small\nset size 0.65,0.65\nplot [%.f:%.f] \"pij%s\" u ($1==%d ? ($3):1/0):($%d/($%d",ageminpar,agemaxpar,fileres,k1,k+cpt+1,k+1); |
|
|
for (i=1; i<= nlstate ; i ++) |
for (i=1; i< nlstate ; i ++) |
fprintf(ficgp,"+$%d",k+i+1); |
fprintf(ficgp,"+$%d",k+i+1); |
fprintf(ficgp,")) t\"prev(%d,%d)\" w l",cpt,cpt+1); |
fprintf(ficgp,")) t\"prev(%d,%d)\" w l",cpt,cpt+1); |
|
|
Line 3454 int main(int argc, char *argv[])
|
Line 3468 int main(int argc, char *argv[])
|
gettimeofday(&start_time, (struct timezone*)0); */ /* at first time */ |
gettimeofday(&start_time, (struct timezone*)0); */ /* at first time */ |
getcwd(pathcd, size); |
getcwd(pathcd, size); |
|
|
printf("\n%s",version); |
printf("\n%s\n%s",version,fullversion); |
if(argc <=1){ |
if(argc <=1){ |
printf("\nEnter the parameter file name: "); |
printf("\nEnter the parameter file name: "); |
scanf("%s",pathtot); |
scanf("%s",pathtot); |