Return to prevlim.h CVS log

Up to [Local Repository] / imach / src

* prevlim.h (Module):

    1:   /*--------------- Prevalence limit  (period or stable prevalence) --------------*/
    2:   
    3:     strcpy(filerespl,"pl");
    4:     strcat(filerespl,fileres);
    5:     if((ficrespl=fopen(filerespl,"w"))==NULL) {
    6:       printf("Problem with period (stable) prevalence resultfile: %s\n", filerespl);goto end;
    7:       fprintf(ficlog,"Problem with period (stable) prevalence resultfile: %s\n", filerespl);goto end;
    8:     }
    9:     printf("Computing period (stable) prevalence: result on file '%s' \n", filerespl);
   10:     fprintf(ficlog,"Computing period (stable) prevalence: result on file '%s' \n", filerespl);
   11:     pstamp(ficrespl);
   12:     fprintf(ficrespl,"# Period (stable) prevalence \n");
   13:     fprintf(ficrespl,"#Age ");
   14:     for(i=1; i<=nlstate;i++) fprintf(ficrespl,"%d-%d ",i,i);
   15:     fprintf(ficrespl,"\n");
   16:   
   17:     prlim=matrix(1,nlstate,1,nlstate);
   18: 
   19:     agebase=ageminpar;
   20:     agelim=agemaxpar;
   21:     ftolpl=1.e-10;
   22:     i1=pow(2,cptcoveff);
   23:     if (cptcovn < 1){i1=1;}
   24: 
   25:     for(cptcov=1,k=0;cptcov<=i1;cptcov++){
   26:     /* for(cptcov=1,k=0;cptcov<=1;cptcov++){ */
   27:       //for(cptcod=1;cptcod<=ncodemax[cptcov];cptcod++){
   28: 	k=k+1;
   29: 	/* to clean */
   30: 	//printf("cptcov=%d cptcod=%d codtab=%d\n",cptcov, cptcod,codtab[cptcod][cptcov]);
   31: 	fprintf(ficrespl,"\n#******");
   32: 	printf("\n#******");
   33: 	fprintf(ficlog,"\n#******");
   34: 	for(j=1;j<=cptcoveff;j++) {
   35: 	  fprintf(ficrespl," V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtab[k][j]]);
   36: 	  printf(" V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtab[k][j]]);
   37: 	  fprintf(ficlog," V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtab[k][j]]);
   38: 	}
   39: 	fprintf(ficrespl,"******\n");
   40: 	printf("******\n");
   41: 	fprintf(ficlog,"******\n");
   42: 
   43: 	fprintf(ficrespl,"#Age ");
   44: 	for(j=1;j<=cptcoveff;j++) {
   45: 	  fprintf(ficrespl,"V%d %d",Tvaraff[j],nbcode[Tvaraff[j]][codtab[k][j]]);
   46: 	}
   47: 	for(i=1; i<=nlstate;i++) fprintf(ficrespl,"%d-%d ",i,i);
   48: 	fprintf(ficrespl,"\n");
   49: 	
   50: 	for (age=agebase; age<=agelim; age++){
   51: 	/* for (age=agebase; age<=agebase; age++){ */
   52: 	  prevalim(prlim, nlstate, p, age, oldm, savm,ftolpl,k);
   53: 	  fprintf(ficrespl,"%.0f ",age );
   54: 	  for(j=1;j<=cptcoveff;j++)
   55: 	    fprintf(ficrespl,"%d %d ",Tvaraff[j],nbcode[Tvaraff[j]][codtab[k][j]]);
   56: 	  for(i=1; i<=nlstate;i++)
   57: 	    fprintf(ficrespl," %.5f", prlim[i][i]);
   58: 	  fprintf(ficrespl,"\n");
   59: 	} /* Age */
   60: 	/* was end of cptcod */
   61:     } /* cptcov */