Return to prevlim.h CVS log

Up to [Local Repository] / imach / src

Summary: prevlim.h

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