|
version 1.247, 2016/09/02 11:11:21
|
version 1.256, 2017/03/27 05:50:23
|
|
Line 1
|
Line 1
|
| /* $Id$ |
/* $Id$ |
| $State$ |
$State$ |
| $Log$ |
$Log$ |
| |
Revision 1.256 2017/03/27 05:50:23 brouard |
| |
Summary: Temporary |
| |
|
| |
Revision 1.255 2017/03/08 16:02:28 brouard |
| |
Summary: IMaCh version 0.99r10 bugs in gnuplot fixed |
| |
|
| |
Revision 1.254 2017/03/08 07:13:00 brouard |
| |
Summary: Fixing data parameter line |
| |
|
| |
Revision 1.253 2016/12/15 11:59:41 brouard |
| |
Summary: 0.99 in progress |
| |
|
| |
Revision 1.252 2016/09/15 21:15:37 brouard |
| |
*** empty log message *** |
| |
|
| |
Revision 1.251 2016/09/15 15:01:13 brouard |
| |
Summary: not working |
| |
|
| |
Revision 1.250 2016/09/08 16:07:27 brouard |
| |
Summary: continue |
| |
|
| |
Revision 1.249 2016/09/07 17:14:18 brouard |
| |
Summary: Starting values from frequencies |
| |
|
| |
Revision 1.248 2016/09/07 14:10:18 brouard |
| |
*** empty log message *** |
| |
|
| Revision 1.247 2016/09/02 11:11:21 brouard |
Revision 1.247 2016/09/02 11:11:21 brouard |
| *** empty log message *** |
*** empty log message *** |
| |
|
|
Line 116
|
Line 143
|
| Author: Nicolas Brouard |
Author: Nicolas Brouard |
| |
|
| Revision 1.210 2015/11/18 17:41:20 brouard |
Revision 1.210 2015/11/18 17:41:20 brouard |
| Summary: Start working on projected prevalences |
Summary: Start working on projected prevalences Revision 1.209 2015/11/17 22:12:03 brouard |
| |
|
| Revision 1.209 2015/11/17 22:12:03 brouard |
|
| Summary: Adding ftolpl parameter |
Summary: Adding ftolpl parameter |
| Author: N Brouard |
Author: N Brouard |
| |
|
|
Line 906 typedef struct {
|
Line 931 typedef struct {
|
| /* #include <libintl.h> */ |
/* #include <libintl.h> */ |
| /* #define _(String) gettext (String) */ |
/* #define _(String) gettext (String) */ |
| |
|
| #define MAXLINE 1024 /* Was 256. Overflow with 312 with 2 states and 4 covariates. Should be ok */ |
#define MAXLINE 2048 /* Was 256 and 1024. Overflow with 312 with 2 states and 4 covariates. Should be ok */ |
| |
|
| #define GNUPLOTPROGRAM "gnuplot" |
#define GNUPLOTPROGRAM "gnuplot" |
| /*#define GNUPLOTPROGRAM "..\\gp37mgw\\wgnuplot"*/ |
/*#define GNUPLOTPROGRAM "..\\gp37mgw\\wgnuplot"*/ |
|
Line 2381 void powell(double p[], double **xi, int
|
Line 2406 void powell(double p[], double **xi, int
|
| flatd++; |
flatd++; |
| } |
} |
| if(flatd >0){ |
if(flatd >0){ |
| printf("%d flat directions\n",flatd); |
printf("%d flat directions: ",flatd); |
| fprintf(ficlog,"%d flat directions\n",flatd); |
fprintf(ficlog,"%d flat directions :",flatd); |
| for (j=1;j<=n;j++) { |
for (j=1;j<=n;j++) { |
| if(flatdir[j]>0){ |
if(flatdir[j]>0){ |
| printf("%d ",j); |
printf("%d ",j); |
|
Line 4110 void ludcmp(double **a, int n, int *indx
|
Line 4135 void ludcmp(double **a, int n, int *indx
|
| big=0.0; |
big=0.0; |
| for (j=1;j<=n;j++) |
for (j=1;j<=n;j++) |
| if ((temp=fabs(a[i][j])) > big) big=temp; |
if ((temp=fabs(a[i][j])) > big) big=temp; |
| if (big == 0.0) nrerror("Singular matrix in routine ludcmp"); |
if (big == 0.0){ |
| |
printf(" Singular Hessian matrix at row %d:\n",i); |
| |
for (j=1;j<=n;j++) { |
| |
printf(" a[%d][%d]=%f,",i,j,a[i][j]); |
| |
fprintf(ficlog," a[%d][%d]=%f,",i,j,a[i][j]); |
| |
} |
| |
fflush(ficlog); |
| |
fclose(ficlog); |
| |
nrerror("Singular matrix in routine ludcmp"); |
| |
} |
| vv[i]=1.0/big; |
vv[i]=1.0/big; |
| } |
} |
| for (j=1;j<=n;j++) { |
for (j=1;j<=n;j++) { |
|
Line 4176 void pstamp(FILE *fichier)
|
Line 4210 void pstamp(FILE *fichier)
|
| fprintf(fichier,"# %s.%s\n#IMaCh version %s, %s\n#%s\n# %s", optionfilefiname,optionfilext,version,copyright, fullversion, strstart); |
fprintf(fichier,"# %s.%s\n#IMaCh version %s, %s\n#%s\n# %s", optionfilefiname,optionfilext,version,copyright, fullversion, strstart); |
| } |
} |
| |
|
| |
int linreg(int ifi, int ila, int *no, const double x[], const double y[], double* a, double* b, double* r, double* sa, double * sb) { |
| |
|
| |
/* y=a+bx regression */ |
| |
double sumx = 0.0; /* sum of x */ |
| |
double sumx2 = 0.0; /* sum of x**2 */ |
| |
double sumxy = 0.0; /* sum of x * y */ |
| |
double sumy = 0.0; /* sum of y */ |
| |
double sumy2 = 0.0; /* sum of y**2 */ |
| |
double sume2; /* sum of square or residuals */ |
| |
double yhat; |
| |
|
| |
double denom=0; |
| |
int i; |
| |
int ne=*no; |
| |
|
| |
for ( i=ifi, ne=0;i<=ila;i++) { |
| |
if(!isfinite(x[i]) || !isfinite(y[i])){ |
| |
/* printf(" x[%d]=%f, y[%d]=%f\n",i,x[i],i,y[i]); */ |
| |
continue; |
| |
} |
| |
ne=ne+1; |
| |
sumx += x[i]; |
| |
sumx2 += x[i]*x[i]; |
| |
sumxy += x[i] * y[i]; |
| |
sumy += y[i]; |
| |
sumy2 += y[i]*y[i]; |
| |
denom = (ne * sumx2 - sumx*sumx); |
| |
/* printf("ne=%d, i=%d,x[%d]=%f, y[%d]=%f sumx=%f, sumx2=%f, sumxy=%f, sumy=%f, sumy2=%f, denom=%f\n",ne,i,i,x[i],i,y[i], sumx, sumx2,sumxy, sumy, sumy2,denom); */ |
| |
} |
| |
|
| |
denom = (ne * sumx2 - sumx*sumx); |
| |
if (denom == 0) { |
| |
// vertical, slope m is infinity |
| |
*b = INFINITY; |
| |
*a = 0; |
| |
if (r) *r = 0; |
| |
return 1; |
| |
} |
| |
|
| |
*b = (ne * sumxy - sumx * sumy) / denom; |
| |
*a = (sumy * sumx2 - sumx * sumxy) / denom; |
| |
if (r!=NULL) { |
| |
*r = (sumxy - sumx * sumy / ne) / /* compute correlation coeff */ |
| |
sqrt((sumx2 - sumx*sumx/ne) * |
| |
(sumy2 - sumy*sumy/ne)); |
| |
} |
| |
*no=ne; |
| |
for ( i=ifi, ne=0;i<=ila;i++) { |
| |
if(!isfinite(x[i]) || !isfinite(y[i])){ |
| |
/* printf(" x[%d]=%f, y[%d]=%f\n",i,x[i],i,y[i]); */ |
| |
continue; |
| |
} |
| |
ne=ne+1; |
| |
yhat = y[i] - *a -*b* x[i]; |
| |
sume2 += yhat * yhat ; |
| |
|
| |
denom = (ne * sumx2 - sumx*sumx); |
| |
/* printf("ne=%d, i=%d,x[%d]=%f, y[%d]=%f sumx=%f, sumx2=%f, sumxy=%f, sumy=%f, sumy2=%f, denom=%f\n",ne,i,i,x[i],i,y[i], sumx, sumx2,sumxy, sumy, sumy2,denom); */ |
| |
} |
| |
*sb = sqrt(sume2/(ne-2)/(sumx2 - sumx * sumx /ne)); |
| |
*sa= *sb * sqrt(sumx2/ne); |
| |
|
| |
return 0; |
| |
} |
| |
|
| /************ Frequencies ********************/ |
/************ Frequencies ********************/ |
| void freqsummary(char fileres[], int iagemin, int iagemax, int **s, double **agev, int nlstate, int imx, \ |
void freqsummary(char fileres[], double p[], double pstart[], int iagemin, int iagemax, int **s, double **agev, int nlstate, int imx, \ |
| int *Tvaraff, int *invalidvarcomb, int **nbcode, int *ncodemax,double **mint,double **anint, char strstart[], \ |
int *Tvaraff, int *invalidvarcomb, int **nbcode, int *ncodemax,double **mint,double **anint, char strstart[], \ |
| int firstpass, int lastpass, int stepm, int weightopt, char model[]) |
int firstpass, int lastpass, int stepm, int weightopt, char model[]) |
| { /* Some frequencies */ |
{ /* Some frequencies as well as proposing some starting values */ |
| |
|
| int i, m, jk, j1, bool, z1,j, k, iv; |
int i, m, jk, j1, bool, z1,j, nj, nl, k, iv, jj=0; |
| int iind=0, iage=0; |
int iind=0, iage=0; |
| int mi; /* Effective wave */ |
int mi; /* Effective wave */ |
| int first; |
int first; |
| double ***freq; /* Frequencies */ |
double ***freq; /* Frequencies */ |
| |
double *x, *y, a,b,r, sa, sb; /* for regression, y=b+m*x and r is the correlation coefficient */ |
| |
int no; |
| double *meanq; |
double *meanq; |
| double **meanqt; |
double **meanqt; |
| double *pp, **prop, *posprop, *pospropt; |
double *pp, **prop, *posprop, *pospropt; |
|
Line 4195 void freqsummary(char fileres[], int ia
|
Line 4296 void freqsummary(char fileres[], int ia
|
| double agebegin, ageend; |
double agebegin, ageend; |
| |
|
| pp=vector(1,nlstate); |
pp=vector(1,nlstate); |
| prop=matrix(1,nlstate,iagemin-AGEMARGE,iagemax+3+AGEMARGE); |
prop=matrix(1,nlstate,iagemin-AGEMARGE,iagemax+4+AGEMARGE); |
| posprop=vector(1,nlstate); /* Counting the number of transition starting from a live state per age */ |
posprop=vector(1,nlstate); /* Counting the number of transition starting from a live state per age */ |
| pospropt=vector(1,nlstate); /* Counting the number of transition starting from a live state */ |
pospropt=vector(1,nlstate); /* Counting the number of transition starting from a live state */ |
| /* prop=matrix(1,nlstate,iagemin,iagemax+3); */ |
/* prop=matrix(1,nlstate,iagemin,iagemax+3); */ |
|
Line 4239 Title=%s <br>Datafile=%s Firstpass=%d La
|
Line 4340 Title=%s <br>Datafile=%s Firstpass=%d La
|
| } |
} |
| fprintf(ficresphtmfr,"Current page is file <a href=\"%s\">%s</a><br>\n\n<h4>Frequencies of all effective transitions of the model, by age at begin of transition, and covariate value at the begin of transition (if the covariate is a varying covariate) </h4>Unknown status is -1<br/>\n",fileresphtmfr, fileresphtmfr); |
fprintf(ficresphtmfr,"Current page is file <a href=\"%s\">%s</a><br>\n\n<h4>Frequencies of all effective transitions of the model, by age at begin of transition, and covariate value at the begin of transition (if the covariate is a varying covariate) </h4>Unknown status is -1<br/>\n",fileresphtmfr, fileresphtmfr); |
| |
|
| freq= ma3x(-5,nlstate+ndeath,-5,nlstate+ndeath,iagemin-AGEMARGE,iagemax+3+AGEMARGE); |
y= vector(iagemin-AGEMARGE,iagemax+4+AGEMARGE); |
| |
x= vector(iagemin-AGEMARGE,iagemax+4+AGEMARGE); |
| |
freq= ma3x(-5,nlstate+ndeath,-5,nlstate+ndeath,iagemin-AGEMARGE,iagemax+4+AGEMARGE); |
| j1=0; |
j1=0; |
| |
|
| /* j=ncoveff; /\* Only fixed dummy covariates *\/ */ |
/* j=ncoveff; /\* Only fixed dummy covariates *\/ */ |
| j=cptcoveff; /* Only dummy covariates of the model */ |
j=cptcoveff; /* Only dummy covariates of the model */ |
| if (cptcovn<1) {j=1;ncodemax[1]=1;} |
if (cptcovn<1) {j=1;ncodemax[1]=1;} |
| |
|
| first=1; |
|
| |
|
| /* Detects if a combination j1 is empty: for a multinomial variable like 3 education levels: |
/* Detects if a combination j1 is empty: for a multinomial variable like 3 education levels: |
| reference=low_education V1=0,V2=0 |
reference=low_education V1=0,V2=0 |
|
Line 4254 Title=%s <br>Datafile=%s Firstpass=%d La
|
Line 4356 Title=%s <br>Datafile=%s Firstpass=%d La
|
| high_educ V1=0 V2=1 |
high_educ V1=0 V2=1 |
| Then V1=1 and V2=1 is a noisy combination that we want to exclude for the list 2**cptcoveff |
Then V1=1 and V2=1 is a noisy combination that we want to exclude for the list 2**cptcoveff |
| */ |
*/ |
| |
dateintsum=0; |
| for (j1 = 1; j1 <= (int) pow(2,j); j1++){ /* Loop on covariates combination in order of model, excluding quantitatives V4=0, V3=0 for example, fixed or varying covariates */ |
k2cpt=0; |
| posproptt=0.; |
|
| /*printf("cptcoveff=%d Tvaraff=%d", cptcoveff,Tvaraff[1]); |
if(cptcoveff == 0 ) |
| scanf("%d", i);*/ |
nl=1; /* Constant model only */ |
| for (i=-5; i<=nlstate+ndeath; i++) |
else |
| for (jk=-5; jk<=nlstate+ndeath; jk++) |
nl=2; |
| |
for (nj = 1; nj <= nl; nj++){ /* nj= 1 constant model, nl number of loops. */ |
| |
if(nj==1) |
| |
j=0; /* First pass for the constant */ |
| |
else |
| |
j=cptcoveff; /* Other passes for the covariate values */ |
| |
first=1; |
| |
for (j1 = 1; j1 <= (int) pow(2,j); j1++){ /* Loop on covariates combination in order of model, excluding quantitatives, V4=0, V3=0 for example, fixed or varying covariates */ |
| |
posproptt=0.; |
| |
/*printf("cptcoveff=%d Tvaraff=%d", cptcoveff,Tvaraff[1]); |
| |
scanf("%d", i);*/ |
| |
for (i=-5; i<=nlstate+ndeath; i++) |
| |
for (jk=-5; jk<=nlstate+ndeath; jk++) |
| |
for(m=iagemin; m <= iagemax+3; m++) |
| |
freq[i][jk][m]=0; |
| |
|
| |
for (i=1; i<=nlstate; i++) { |
| for(m=iagemin; m <= iagemax+3; m++) |
for(m=iagemin; m <= iagemax+3; m++) |
| freq[i][jk][m]=0; |
prop[i][m]=0; |
| |
posprop[i]=0; |
| for (i=1; i<=nlstate; i++) { |
pospropt[i]=0; |
| for(m=iagemin; m <= iagemax+3; m++) |
} |
| prop[i][m]=0; |
/* for (z1=1; z1<= nqfveff; z1++) { */ |
| posprop[i]=0; |
/* meanq[z1]+=0.; */ |
| pospropt[i]=0; |
/* for(m=1;m<=lastpass;m++){ */ |
| } |
/* meanqt[m][z1]=0.; */ |
| /* for (z1=1; z1<= nqfveff; z1++) { */ |
/* } */ |
| /* meanq[z1]+=0.; */ |
/* } */ |
| /* for(m=1;m<=lastpass;m++){ */ |
|
| /* meanqt[m][z1]=0.; */ |
/* dateintsum=0; */ |
| /* } */ |
/* k2cpt=0; */ |
| /* } */ |
|
| |
/* For that combination of covariate j1, we count and print the frequencies in one pass */ |
| dateintsum=0; |
for (iind=1; iind<=imx; iind++) { /* For each individual iind */ |
| k2cpt=0; |
bool=1; |
| /* For that combination of covariate j1, we count and print the frequencies in one pass */ |
if(j !=0){ |
| for (iind=1; iind<=imx; iind++) { /* For each individual iind */ |
if(anyvaryingduminmodel==0){ /* If All fixed covariates */ |
| bool=1; |
if (cptcoveff >0) { /* Filter is here: Must be looked at for model=V1+V2+V3+V4 */ |
| if(anyvaryingduminmodel==0){ /* If All fixed covariates */ |
/* for (z1=1; z1<= nqfveff; z1++) { */ |
| if (cptcoveff >0) { /* Filter is here: Must be looked at for model=V1+V2+V3+V4 */ |
/* meanq[z1]+=coqvar[Tvar[z1]][iind]; /\* Computes mean of quantitative with selected filter *\/ */ |
| /* for (z1=1; z1<= nqfveff; z1++) { */ |
/* } */ |
| /* meanq[z1]+=coqvar[Tvar[z1]][iind]; /\* Computes mean of quantitative with selected filter *\/ */ |
for (z1=1; z1<=cptcoveff; z1++) { /* loops on covariates in the model */ |
| /* } */ |
/* if(Tvaraff[z1] ==-20){ */ |
| for (z1=1; z1<=cptcoveff; z1++) { |
/* /\* sumnew+=cotvar[mw[mi][iind]][z1][iind]; *\/ */ |
| /* if(Tvaraff[z1] ==-20){ */ |
/* }else if(Tvaraff[z1] ==-10){ */ |
| /* /\* sumnew+=cotvar[mw[mi][iind]][z1][iind]; *\/ */ |
/* /\* sumnew+=coqvar[z1][iind]; *\/ */ |
| /* }else if(Tvaraff[z1] ==-10){ */ |
/* }else */ |
| /* /\* sumnew+=coqvar[z1][iind]; *\/ */ |
if (covar[Tvaraff[z1]][iind]!= nbcode[Tvaraff[z1]][codtabm(j1,z1)]){ /* for combination j1 of covariates */ |
| /* }else */ |
/* Tests if this individual iind responded to combination j1 (V4=1 V3=0) */ |
| if (covar[Tvaraff[z1]][iind]!= nbcode[Tvaraff[z1]][codtabm(j1,z1)]){ |
bool=0; /* bool should be equal to 1 to be selected, one covariate value failed */ |
| /* Tests if this individual iind responded to j1 (V4=1 V3=0) */ |
/* printf("bool=%d i=%d, z1=%d, Tvaraff[%d]=%d, covar[Tvarff][%d]=%2f, codtabm(%d,%d)=%d, nbcode[Tvaraff][codtabm(%d,%d)=%d, j1=%d\n", |
| bool=0; |
bool,i,z1, z1, Tvaraff[z1],i,covar[Tvaraff[z1]][i],j1,z1,codtabm(j1,z1), |
| /* printf("bool=%d i=%d, z1=%d, Tvaraff[%d]=%d, covar[Tvarff][%d]=%2f, codtabm(%d,%d)=%d, nbcode[Tvaraff][codtabm(%d,%d)=%d, j1=%d\n", |
j1,z1,nbcode[Tvaraff[z1]][codtabm(j1,z1)],j1);*/ |
| bool,i,z1, z1, Tvaraff[z1],i,covar[Tvaraff[z1]][i],j1,z1,codtabm(j1,z1), |
/* For j1=7 in V1+V2+V3+V4 = 0 1 1 0 and codtabm(7,3)=1 and nbcde[3][?]=1*/ |
| j1,z1,nbcode[Tvaraff[z1]][codtabm(j1,z1)],j1);*/ |
} /* Onlyf fixed */ |
| /* For j1=7 in V1+V2+V3+V4 = 0 1 1 0 and codtabm(7,3)=1 and nbcde[3][?]=1*/ |
} /* end z1 */ |
| } /* Onlyf fixed */ |
} /* cptcovn > 0 */ |
| } /* end z1 */ |
} /* end any */ |
| } /* cptcovn > 0 */ |
}/* end j==0 */ |
| } /* end any */ |
if (bool==1){ /* We selected an individual iind satisfying combination j1 or all fixed */ |
| if (bool==1){ /* We selected an individual iind satisfying combination j1 or all fixed */ |
/* for(m=firstpass; m<=lastpass; m++){ */ |
| /* for(m=firstpass; m<=lastpass; m++){ */ |
for(mi=1; mi<wav[iind];mi++){ /* For that wave */ |
| for(mi=1; mi<wav[iind];mi++){ /* For that wave */ |
m=mw[mi][iind]; |
| m=mw[mi][iind]; |
if(j!=0){ |
| if(anyvaryingduminmodel==1){ /* Some are varying covariates */ |
if(anyvaryingduminmodel==1){ /* Some are varying covariates */ |
| for (z1=1; z1<=cptcoveff; z1++) { |
for (z1=1; z1<=cptcoveff; z1++) { |
| if( Fixed[Tmodelind[z1]]==1){ |
if( Fixed[Tmodelind[z1]]==1){ |
| iv= Tvar[Tmodelind[z1]]-ncovcol-nqv; |
iv= Tvar[Tmodelind[z1]]-ncovcol-nqv; |
| if (cotvar[m][iv][iind]!= nbcode[Tvaraff[z1]][codtabm(j1,z1)]) /* iv=1 to ntv, right modality */ |
if (cotvar[m][iv][iind]!= nbcode[Tvaraff[z1]][codtabm(j1,z1)]) /* iv=1 to ntv, right modality. If covariate's |
| bool=0; |
value is -1, we don't select. It differs from the |
| }else if( Fixed[Tmodelind[z1]]== 0) { /* fixed */ |
constant and age model which counts them. */ |
| if (covar[Tvaraff[z1]][iind]!= nbcode[Tvaraff[z1]][codtabm(j1,z1)]) { |
bool=0; /* not selected */ |
| bool=0; |
}else if( Fixed[Tmodelind[z1]]== 0) { /* fixed */ |
| |
if (covar[Tvaraff[z1]][iind]!= nbcode[Tvaraff[z1]][codtabm(j1,z1)]) { |
| |
bool=0; |
| |
} |
| |
} |
| } |
} |
| |
}/* Some are varying covariates, we tried to speed up if all fixed covariates in the model, avoiding waves loop */ |
| |
} /* end j==0 */ |
| |
/* bool =0 we keep that guy which corresponds to the combination of dummy values */ |
| |
if(bool==1){ |
| |
/* dh[m][iind] or dh[mw[mi][iind]][iind] is the delay between two effective (mi) waves m=mw[mi][iind] |
| |
and mw[mi+1][iind]. dh depends on stepm. */ |
| |
agebegin=agev[m][iind]; /* Age at beginning of wave before transition*/ |
| |
ageend=agev[m][iind]+(dh[m][iind])*stepm/YEARM; /* Age at end of wave and transition */ |
| |
if(m >=firstpass && m <=lastpass){ |
| |
k2=anint[m][iind]+(mint[m][iind]/12.); |
| |
/*if ((k2>=dateprev1) && (k2<=dateprev2)) {*/ |
| |
if(agev[m][iind]==0) agev[m][iind]=iagemax+1; /* All ages equal to 0 are in iagemax+1 */ |
| |
if(agev[m][iind]==1) agev[m][iind]=iagemax+2; /* All ages equal to 1 are in iagemax+2 */ |
| |
if (s[m][iind]>0 && s[m][iind]<=nlstate) /* If status at wave m is known and a live state */ |
| |
prop[s[m][iind]][(int)agev[m][iind]] += weight[iind]; /* At age of beginning of transition, where status is known */ |
| |
if (m<lastpass) { |
| |
/* if(s[m][iind]==4 && s[m+1][iind]==4) */ |
| |
/* printf(" num=%ld m=%d, iind=%d s1=%d s2=%d agev at m=%d\n", num[iind], m, iind,s[m][iind],s[m+1][iind], (int)agev[m][iind]); */ |
| |
if(s[m][iind]==-1) |
| |
printf(" num=%ld m=%d, iind=%d s1=%d s2=%d agev at m=%d agebegin=%.2f ageend=%.2f, agemed=%d\n", num[iind], m, iind,s[m][iind],s[m+1][iind], (int)agev[m][iind],agebegin, ageend, (int)((agebegin+ageend)/2.)); |
| |
freq[s[m][iind]][s[m+1][iind]][(int)agev[m][iind]] += weight[iind]; /* At age of beginning of transition, where status is known */ |
| |
/* if((int)agev[m][iind] == 55) */ |
| |
/* printf("j=%d, j1=%d Age %d, iind=%d, num=%09ld m=%d\n",j,j1,(int)agev[m][iind],iind, num[iind],m); */ |
| |
/* freq[s[m][iind]][s[m+1][iind]][(int)((agebegin+ageend)/2.)] += weight[iind]; */ |
| |
freq[s[m][iind]][s[m+1][iind]][iagemax+3] += weight[iind]; /* Total is in iagemax+3 *//* At age of beginning of transition, where status is known */ |
| |
} |
| |
} /* end if between passes */ |
| |
if ((agev[m][iind]>1) && (agev[m][iind]< (iagemax+3)) && (anint[m][iind]!=9999) && (mint[m][iind]!=99) && (j==0)) { |
| |
dateintsum=dateintsum+k2; /* on all covariates ?*/ |
| |
k2cpt++; |
| |
/* printf("iind=%ld dateintmean = %lf dateintsum=%lf k2cpt=%lf k2=%lf\n",iind, dateintsum/k2cpt, dateintsum,k2cpt, k2); */ |
| } |
} |
| } |
}else{ |
| }/* Some are varying covariates, we tried to speed up if all fixed covariates in the model, avoiding waves loop */ |
bool=1; |
| /* bool =0 we keep that guy which corresponds to the combination of dummy values */ |
}/* end bool 2 */ |
| if(bool==1){ |
} /* end m */ |
| /* dh[m][iind] or dh[mw[mi][iind]][iind] is the delay between two effective (mi) waves m=mw[mi][iind] |
} /* end bool */ |
| and mw[mi+1][iind]. dh depends on stepm. */ |
} /* end iind = 1 to imx */ |
| agebegin=agev[m][iind]; /* Age at beginning of wave before transition*/ |
/* prop[s][age] is feeded for any initial and valid live state as well as |
| ageend=agev[m][iind]+(dh[m][iind])*stepm/YEARM; /* Age at end of wave and transition */ |
freq[s1][s2][age] at single age of beginning the transition, for a combination j1 */ |
| if(m >=firstpass && m <=lastpass){ |
|
| k2=anint[m][iind]+(mint[m][iind]/12.); |
|
| /*if ((k2>=dateprev1) && (k2<=dateprev2)) {*/ |
/* fprintf(ficresp, "#Count between %.lf/%.lf/%.lf and %.lf/%.lf/%.lf\n",jprev1, mprev1,anprev1,jprev2, mprev2,anprev2);*/ |
| if(agev[m][iind]==0) agev[m][iind]=iagemax+1; /* All ages equal to 0 are in iagemax+1 */ |
pstamp(ficresp); |
| if(agev[m][iind]==1) agev[m][iind]=iagemax+2; /* All ages equal to 1 are in iagemax+2 */ |
if (cptcoveff>0 && j!=0){ |
| if (s[m][iind]>0 && s[m][iind]<=nlstate) /* If status at wave m is known and a live state */ |
printf( "\n#********** Variable "); |
| prop[s[m][iind]][(int)agev[m][iind]] += weight[iind]; /* At age of beginning of transition, where status is known */ |
fprintf(ficresp, "\n#********** Variable "); |
| if (m<lastpass) { |
fprintf(ficresphtm, "\n<br/><br/><h3>********** Variable "); |
| /* if(s[m][iind]==4 && s[m+1][iind]==4) */ |
fprintf(ficresphtmfr, "\n<br/><br/><h3>********** Variable "); |
| /* printf(" num=%ld m=%d, iind=%d s1=%d s2=%d agev at m=%d\n", num[iind], m, iind,s[m][iind],s[m+1][iind], (int)agev[m][iind]); */ |
fprintf(ficlog, "\n#********** Variable "); |
| if(s[m][iind]==-1) |
for (z1=1; z1<=cptcoveff; z1++){ |
| printf(" num=%ld m=%d, iind=%d s1=%d s2=%d agev at m=%d agebegin=%.2f ageend=%.2f, agemed=%d\n", num[iind], m, iind,s[m][iind],s[m+1][iind], (int)agev[m][iind],agebegin, ageend, (int)((agebegin+ageend)/2.)); |
if(!FixedV[Tvaraff[z1]]){ |
| freq[s[m][iind]][s[m+1][iind]][(int)agev[m][iind]] += weight[iind]; /* At age of beginning of transition, where status is known */ |
printf( "V%d(fixed)=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtabm(j1,z1)]); |
| /* freq[s[m][iind]][s[m+1][iind]][(int)((agebegin+ageend)/2.)] += weight[iind]; */ |
fprintf(ficresp, "V%d(fixed)=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtabm(j1,z1)]); |
| freq[s[m][iind]][s[m+1][iind]][iagemax+3] += weight[iind]; /* Total is in iagemax+3 *//* At age of beginning of transition, where status is known */ |
fprintf(ficresphtm, "V%d(fixed)=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtabm(j1,z1)]); |
| } |
fprintf(ficresphtmfr, "V%d(fixed)=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtabm(j1,z1)]); |
| } /* end if between passes */ |
fprintf(ficlog, "V%d(fixed)=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtabm(j1,z1)]); |
| if ((agev[m][iind]>1) && (agev[m][iind]< (iagemax+3)) && (anint[m][iind]!=9999) && (mint[m][iind]!=99)) { |
}else{ |
| dateintsum=dateintsum+k2; |
printf( "V%d(varying)=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtabm(j1,z1)]); |
| k2cpt++; |
fprintf(ficresp, "V%d(varying)=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtabm(j1,z1)]); |
| /* printf("iind=%ld dateintmean = %lf dateintsum=%lf k2cpt=%lf k2=%lf\n",iind, dateintsum/k2cpt, dateintsum,k2cpt, k2); */ |
fprintf(ficresphtm, "V%d(varying)=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtabm(j1,z1)]); |
| } |
fprintf(ficresphtmfr, "V%d(varying)=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtabm(j1,z1)]); |
| } /* end bool 2 */ |
fprintf(ficlog, "V%d(varying)=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtabm(j1,z1)]); |
| } /* end m */ |
} |
| } /* end bool */ |
} |
| } /* end iind = 1 to imx */ |
printf( "**********\n#"); |
| /* prop[s][age] is feeded for any initial and valid live state as well as |
fprintf(ficresp, "**********\n#"); |
| freq[s1][s2][age] at single age of beginning the transition, for a combination j1 */ |
fprintf(ficresphtm, "**********</h3>\n"); |
| |
fprintf(ficresphtmfr, "**********</h3>\n"); |
| |
fprintf(ficlog, "**********\n"); |
| /* fprintf(ficresp, "#Count between %.lf/%.lf/%.lf and %.lf/%.lf/%.lf\n",jprev1, mprev1,anprev1,jprev2, mprev2,anprev2);*/ |
} |
| pstamp(ficresp); |
fprintf(ficresphtm,"<table style=\"text-align:center; border: 1px solid\">"); |
| if (cptcoveff>0){ |
for(i=1; i<=nlstate;i++) { |
| fprintf(ficresp, "\n#********** Variable "); |
fprintf(ficresp, " Age Prev(%d) N(%d) N ",i,i); |
| fprintf(ficresphtm, "\n<br/><br/><h3>********** Variable "); |
fprintf(ficresphtm, "<th>Age</th><th>Prev(%d)</th><th>N(%d)</th><th>N</th>",i,i); |
| fprintf(ficresphtmfr, "\n<br/><br/><h3>********** Variable "); |
} |
| fprintf(ficlog, "\n#********** Variable "); |
fprintf(ficresp, "\n"); |
| for (z1=1; z1<=cptcoveff; z1++){ |
fprintf(ficresphtm, "\n"); |
| if(DummyV[z1]){ |
|
| fprintf(ficresp, "V%d (fixed)=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtabm(j1,z1)]); |
/* Header of frequency table by age */ |
| fprintf(ficresphtm, "V%d (fixed)=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtabm(j1,z1)]); |
fprintf(ficresphtmfr,"<table style=\"text-align:center; border: 1px solid\">"); |
| fprintf(ficresphtmfr, "V%d (fixed)=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtabm(j1,z1)]); |
fprintf(ficresphtmfr,"<th>Age</th> "); |
| fprintf(ficlog, "V%d (fixed)=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtabm(j1,z1)]); |
for(jk=-1; jk <=nlstate+ndeath; jk++){ |
| }else{ |
for(m=-1; m <=nlstate+ndeath; m++){ |
| fprintf(ficresp, "V%d(varying)=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtabm(j1,z1)]); |
if(jk!=0 && m!=0) |
| fprintf(ficresphtm, "V%d(varying)=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtabm(j1,z1)]); |
fprintf(ficresphtmfr,"<th>%d%d</th> ",jk,m); |
| fprintf(ficresphtmfr, "V%d(varying)=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtabm(j1,z1)]); |
|
| fprintf(ficlog, "V%d(varying)=%d ",Tvaraff[z1],nbcode[Tvaraff[z1]][codtabm(j1,z1)]); |
|
| } |
|
| } |
|
| fprintf(ficresp, "**********\n#"); |
|
| fprintf(ficresphtm, "**********</h3>\n"); |
|
| fprintf(ficresphtmfr, "**********</h3>\n"); |
|
| fprintf(ficlog, "**********\n"); |
|
| } |
|
| fprintf(ficresphtm,"<table style=\"text-align:center; border: 1px solid\">"); |
|
| for(i=1; i<=nlstate;i++) { |
|
| fprintf(ficresp, " Age Prev(%d) N(%d) N ",i,i); |
|
| fprintf(ficresphtm, "<th>Age</th><th>Prev(%d)</th><th>N(%d)</th><th>N</th>",i,i); |
|
| } |
|
| fprintf(ficresp, "\n"); |
|
| fprintf(ficresphtm, "\n"); |
|
| |
|
| /* Header of frequency table by age */ |
|
| fprintf(ficresphtmfr,"<table style=\"text-align:center; border: 1px solid\">"); |
|
| fprintf(ficresphtmfr,"<th>Age</th> "); |
|
| for(jk=-1; jk <=nlstate+ndeath; jk++){ |
|
| for(m=-1; m <=nlstate+ndeath; m++){ |
|
| if(jk!=0 && m!=0) |
|
| fprintf(ficresphtmfr,"<th>%d%d</th> ",jk,m); |
|
| } |
|
| } |
|
| fprintf(ficresphtmfr, "\n"); |
|
| |
|
| /* For each age */ |
|
| for(iage=iagemin; iage <= iagemax+3; iage++){ |
|
| fprintf(ficresphtm,"<tr>"); |
|
| if(iage==iagemax+1){ |
|
| fprintf(ficlog,"1"); |
|
| fprintf(ficresphtmfr,"<tr><th>0</th> "); |
|
| }else if(iage==iagemax+2){ |
|
| fprintf(ficlog,"0"); |
|
| fprintf(ficresphtmfr,"<tr><th>Unknown</th> "); |
|
| }else if(iage==iagemax+3){ |
|
| fprintf(ficlog,"Total"); |
|
| fprintf(ficresphtmfr,"<tr><th>Total</th> "); |
|
| }else{ |
|
| if(first==1){ |
|
| first=0; |
|
| printf("See log file for details...\n"); |
|
| } |
} |
| fprintf(ficresphtmfr,"<tr><th>%d</th> ",iage); |
|
| fprintf(ficlog,"Age %d", iage); |
|
| } |
|
| for(jk=1; jk <=nlstate ; jk++){ |
|
| for(m=-1, pp[jk]=0; m <=nlstate+ndeath ; m++) |
|
| pp[jk] += freq[jk][m][iage]; |
|
| } |
} |
| for(jk=1; jk <=nlstate ; jk++){ |
fprintf(ficresphtmfr, "\n"); |
| for(m=-1, pos=0; m <=0 ; m++) |
|
| pos += freq[jk][m][iage]; |
/* For each age */ |
| if(pp[jk]>=1.e-10){ |
for(iage=iagemin; iage <= iagemax+3; iage++){ |
| |
fprintf(ficresphtm,"<tr>"); |
| |
if(iage==iagemax+1){ |
| |
fprintf(ficlog,"1"); |
| |
fprintf(ficresphtmfr,"<tr><th>0</th> "); |
| |
}else if(iage==iagemax+2){ |
| |
fprintf(ficlog,"0"); |
| |
fprintf(ficresphtmfr,"<tr><th>Unknown</th> "); |
| |
}else if(iage==iagemax+3){ |
| |
fprintf(ficlog,"Total"); |
| |
fprintf(ficresphtmfr,"<tr><th>Total</th> "); |
| |
}else{ |
| if(first==1){ |
if(first==1){ |
| printf(" %d.=%.0f loss[%d]=%.1f%%",jk,pp[jk],jk,100*pos/pp[jk]); |
first=0; |
| |
printf("See log file for details...\n"); |
| |
} |
| |
fprintf(ficresphtmfr,"<tr><th>%d</th> ",iage); |
| |
fprintf(ficlog,"Age %d", iage); |
| |
} |
| |
for(jk=1; jk <=nlstate ; jk++){ |
| |
for(m=-1, pp[jk]=0; m <=nlstate+ndeath ; m++) |
| |
pp[jk] += freq[jk][m][iage]; |
| |
} |
| |
for(jk=1; jk <=nlstate ; jk++){ |
| |
for(m=-1, pos=0; m <=0 ; m++) |
| |
pos += freq[jk][m][iage]; |
| |
if(pp[jk]>=1.e-10){ |
| |
if(first==1){ |
| |
printf(" %d.=%.0f loss[%d]=%.1f%%",jk,pp[jk],jk,100*pos/pp[jk]); |
| |
} |
| |
fprintf(ficlog," %d.=%.0f loss[%d]=%.1f%%",jk,pp[jk],jk,100*pos/pp[jk]); |
| |
}else{ |
| |
if(first==1) |
| |
printf(" %d.=%.0f loss[%d]=NaNQ%%",jk,pp[jk],jk); |
| |
fprintf(ficlog," %d.=%.0f loss[%d]=NaNQ%%",jk,pp[jk],jk); |
| } |
} |
| fprintf(ficlog," %d.=%.0f loss[%d]=%.1f%%",jk,pp[jk],jk,100*pos/pp[jk]); |
|
| }else{ |
|
| if(first==1) |
|
| printf(" %d.=%.0f loss[%d]=NaNQ%%",jk,pp[jk],jk); |
|
| fprintf(ficlog," %d.=%.0f loss[%d]=NaNQ%%",jk,pp[jk],jk); |
|
| } |
} |
| } |
|
| |
|
| for(jk=1; jk <=nlstate ; jk++){ |
|
| /* posprop[jk]=0; */ |
|
| for(m=0, pp[jk]=0; m <=nlstate+ndeath; m++)/* Summing on all ages */ |
|
| pp[jk] += freq[jk][m][iage]; |
|
| } /* pp[jk] is the total number of transitions starting from state jk and any ending status until this age */ |
|
| |
|
| for(jk=1,pos=0, pospropta=0.; jk <=nlstate ; jk++){ |
for(jk=1; jk <=nlstate ; jk++){ |
| pos += pp[jk]; /* pos is the total number of transitions until this age */ |
/* posprop[jk]=0; */ |
| posprop[jk] += prop[jk][iage]; /* prop is the number of transitions from a live state |
for(m=0, pp[jk]=0; m <=nlstate+ndeath; m++)/* Summing on all ages */ |
| |
pp[jk] += freq[jk][m][iage]; |
| |
} /* pp[jk] is the total number of transitions starting from state jk and any ending status until this age */ |
| |
|
| |
for(jk=1,pos=0, pospropta=0.; jk <=nlstate ; jk++){ |
| |
pos += pp[jk]; /* pos is the total number of transitions until this age */ |
| |
posprop[jk] += prop[jk][iage]; /* prop is the number of transitions from a live state |
| |
from jk at age iage prop[s[m][iind]][(int)agev[m][iind]] += weight[iind] */ |
| |
pospropta += prop[jk][iage]; /* prop is the number of transitions from a live state |
| from jk at age iage prop[s[m][iind]][(int)agev[m][iind]] += weight[iind] */ |
from jk at age iage prop[s[m][iind]][(int)agev[m][iind]] += weight[iind] */ |
| pospropta += prop[jk][iage]; /* prop is the number of transitions from a live state |
|
| from jk at age iage prop[s[m][iind]][(int)agev[m][iind]] += weight[iind] */ |
|
| } |
|
| for(jk=1; jk <=nlstate ; jk++){ |
|
| if(pos>=1.e-5){ |
|
| if(first==1) |
|
| printf(" %d.=%.0f prev[%d]=%.1f%%",jk,pp[jk],jk,100*pp[jk]/pos); |
|
| fprintf(ficlog," %d.=%.0f prev[%d]=%.1f%%",jk,pp[jk],jk,100*pp[jk]/pos); |
|
| }else{ |
|
| if(first==1) |
|
| printf(" %d.=%.0f prev[%d]=NaNQ%%",jk,pp[jk],jk); |
|
| fprintf(ficlog," %d.=%.0f prev[%d]=NaNQ%%",jk,pp[jk],jk); |
|
| } |
} |
| if( iage <= iagemax){ |
for(jk=1; jk <=nlstate ; jk++){ |
| if(pos>=1.e-5){ |
if(pos>=1.e-5){ |
| fprintf(ficresp," %d %.5f %.0f %.0f",iage,prop[jk][iage]/pospropta, prop[jk][iage],pospropta); |
if(first==1) |
| fprintf(ficresphtm,"<th>%d</th><td>%.5f</td><td>%.0f</td><td>%.0f</td>",iage,prop[jk][iage]/pospropta, prop[jk][iage],pospropta); |
printf(" %d.=%.0f prev[%d]=%.1f%%",jk,pp[jk],jk,100*pp[jk]/pos); |
| /*probs[iage][jk][j1]= pp[jk]/pos;*/ |
fprintf(ficlog," %d.=%.0f prev[%d]=%.1f%%",jk,pp[jk],jk,100*pp[jk]/pos); |
| /*printf("\niage=%d jk=%d j1=%d %.5f %.0f %.0f %f",iage,jk,j1,pp[jk]/pos, pp[jk],pos,probs[iage][jk][j1]);*/ |
}else{ |
| |
if(first==1) |
| |
printf(" %d.=%.0f prev[%d]=NaNQ%%",jk,pp[jk],jk); |
| |
fprintf(ficlog," %d.=%.0f prev[%d]=NaNQ%%",jk,pp[jk],jk); |
| |
} |
| |
if( iage <= iagemax){ |
| |
if(pos>=1.e-5){ |
| |
fprintf(ficresp," %d %.5f %.0f %.0f",iage,prop[jk][iage]/pospropta, prop[jk][iage],pospropta); |
| |
fprintf(ficresphtm,"<th>%d</th><td>%.5f</td><td>%.0f</td><td>%.0f</td>",iage,prop[jk][iage]/pospropta, prop[jk][iage],pospropta); |
| |
/*probs[iage][jk][j1]= pp[jk]/pos;*/ |
| |
/*printf("\niage=%d jk=%d j1=%d %.5f %.0f %.0f %f",iage,jk,j1,pp[jk]/pos, pp[jk],pos,probs[iage][jk][j1]);*/ |
| |
} |
| |
else{ |
| |
fprintf(ficresp," %d NaNq %.0f %.0f",iage,prop[jk][iage],pospropta); |
| |
fprintf(ficresphtm,"<th>%d</th><td>NaNq</td><td>%.0f</td><td>%.0f</td>",iage, prop[jk][iage],pospropta); |
| |
} |
| } |
} |
| else{ |
pospropt[jk] +=posprop[jk]; |
| fprintf(ficresp," %d NaNq %.0f %.0f",iage,prop[jk][iage],pospropta); |
} /* end loop jk */ |
| fprintf(ficresphtm,"<th>%d</th><td>NaNq</td><td>%.0f</td><td>%.0f</td>",iage, prop[jk][iage],pospropta); |
/* pospropt=0.; */ |
| |
for(jk=-1; jk <=nlstate+ndeath; jk++){ |
| |
for(m=-1; m <=nlstate+ndeath; m++){ |
| |
if(freq[jk][m][iage] !=0 ) { /* minimizing output */ |
| |
if(first==1){ |
| |
printf(" %d%d=%.0f",jk,m,freq[jk][m][iage]); |
| |
} |
| |
/* printf(" %d%d=%.0f",jk,m,freq[jk][m][iage]); */ |
| |
fprintf(ficlog," %d%d=%.0f",jk,m,freq[jk][m][iage]); |
| |
} |
| |
if(jk!=0 && m!=0) |
| |
fprintf(ficresphtmfr,"<td>%.0f</td> ",freq[jk][m][iage]); |
| } |
} |
| |
} /* end loop jk */ |
| |
posproptt=0.; |
| |
for(jk=1; jk <=nlstate; jk++){ |
| |
posproptt += pospropt[jk]; |
| |
} |
| |
fprintf(ficresphtmfr,"</tr>\n "); |
| |
if(iage <= iagemax){ |
| |
fprintf(ficresp,"\n"); |
| |
fprintf(ficresphtm,"</tr>\n"); |
| } |
} |
| pospropt[jk] +=posprop[jk]; |
if(first==1) |
| } /* end loop jk */ |
printf("Others in log...\n"); |
| /* pospropt=0.; */ |
fprintf(ficlog,"\n"); |
| for(jk=-1; jk <=nlstate+ndeath; jk++){ |
} /* end loop age iage */ |
| for(m=-1; m <=nlstate+ndeath; m++){ |
fprintf(ficresphtm,"<tr><th>Tot</th>"); |
| if(freq[jk][m][iage] !=0 ) { /* minimizing output */ |
for(jk=1; jk <=nlstate ; jk++){ |
| if(first==1){ |
if(posproptt < 1.e-5){ |
| printf(" %d%d=%.0f",jk,m,freq[jk][m][iage]); |
fprintf(ficresphtm,"<td>Nanq</td><td>%.0f</td><td>%.0f</td>",pospropt[jk],posproptt); |
| |
}else{ |
| |
fprintf(ficresphtm,"<td>%.5f</td><td>%.0f</td><td>%.0f</td>",pospropt[jk]/posproptt,pospropt[jk],posproptt); |
| |
} |
| |
} |
| |
fprintf(ficresphtm,"</tr>\n"); |
| |
fprintf(ficresphtm,"</table>\n"); |
| |
fprintf(ficresphtmfr,"</table>\n"); |
| |
if(posproptt < 1.e-5){ |
| |
fprintf(ficresphtm,"\n <p><b> This combination (%d) is not valid and no result will be produced</b></p>",j1); |
| |
fprintf(ficresphtmfr,"\n <p><b> This combination (%d) is not valid and no result will be produced</b></p>",j1); |
| |
fprintf(ficres,"\n This combination (%d) is not valid and no result will be produced\n\n",j1); |
| |
invalidvarcomb[j1]=1; |
| |
}else{ |
| |
fprintf(ficresphtm,"\n <p> This combination (%d) is valid and result will be produced.</p>",j1); |
| |
invalidvarcomb[j1]=0; |
| |
} |
| |
fprintf(ficresphtmfr,"</table>\n"); |
| |
fprintf(ficlog,"\n"); |
| |
if(j!=0){ |
| |
printf("#Freqsummary: Starting values for combination j1=%d:\n", j1); |
| |
for(i=1,jk=1; i <=nlstate; i++){ |
| |
for(k=1; k <=(nlstate+ndeath); k++){ |
| |
if (k != i) { |
| |
for(jj=1; jj <=ncovmodel; jj++){ /* For counting jk */ |
| |
if(jj==1){ /* Constant case (in fact cste + age) */ |
| |
if(j1==1){ /* All dummy covariates to zero */ |
| |
freq[i][k][iagemax+4]=freq[i][k][iagemax+3]; /* Stores case 0 0 0 */ |
| |
freq[i][i][iagemax+4]=freq[i][i][iagemax+3]; /* Stores case 0 0 0 */ |
| |
printf("%d%d ",i,k); |
| |
fprintf(ficlog,"%d%d ",i,k); |
| |
printf("%12.7f ln(%.0f/%.0f)= %f, OR=%f sd=%f \n",p[jk],freq[i][k][iagemax+3],freq[i][i][iagemax+3], log(freq[i][k][iagemax+3]/freq[i][i][iagemax+3]),freq[i][k][iagemax+3]/freq[i][i][iagemax+3], sqrt(1/freq[i][k][iagemax+3]+1/freq[i][i][iagemax+3])); |
| |
fprintf(ficlog,"%12.7f ln(%.0f/%.0f)= %12.7f \n",p[jk],freq[i][k][iagemax+3],freq[i][i][iagemax+3], log(freq[i][k][iagemax+3]/freq[i][i][iagemax+3])); |
| |
pstart[jk]= log(freq[i][k][iagemax+3]/freq[i][i][iagemax+3]); |
| |
} |
| |
}else if((j1==1) && (jj==2 || nagesqr==1)){ /* age or age*age parameter without covariate V4*age (to be done later) */ |
| |
for(iage=iagemin; iage <= iagemax+3; iage++){ |
| |
x[iage]= (double)iage; |
| |
y[iage]= log(freq[i][k][iage]/freq[i][i][iage]); |
| |
/* printf("i=%d, k=%d, jk=%d, j1=%d, jj=%d, y[%d]=%f\n",i,k,jk,j1,jj, iage, y[iage]); */ |
| |
} |
| |
linreg(iagemin,iagemax,&no,x,y,&a,&b,&r, &sa, &sb ); /* y= a+b*x with standard errors */ |
| |
pstart[jk]=b; |
| |
pstart[jk-1]=a; |
| |
}else if( j1!=1 && (j1==2 || (log(j1-1.)/log(2.)-(int)(log(j1-1.)/log(2.))) <0.010) && ( TvarsDind[(int)(log(j1-1.)/log(2.))+1]+2+nagesqr == jj) && Dummy[jj-2-nagesqr]==0){ /* We want only if the position, jj, in model corresponds to unique covariate equal to 1 in j1 combination */ |
| |
printf("j1=%d, jj=%d, (int)(log(j1-1.)/log(2.))+1=%d, TvarsDind[(int)(log(j1-1.)/log(2.))+1]=%d\n",j1, jj,(int)(log(j1-1.)/log(2.))+1,TvarsDind[(int)(log(j1-1.)/log(2.))+1]); |
| |
printf("j1=%d, jj=%d, (log(j1-1.)/log(2.))+1=%f, TvarsDind[(int)(log(j1-1.)/log(2.))+1]=%d\n",j1, jj,(log(j1-1.)/log(2.))+1,TvarsDind[(int)(log(j1-1.)/log(2.))+1]); |
| |
pstart[jk]= log((freq[i][k][iagemax+3]/freq[i][i][iagemax+3])/(freq[i][k][iagemax+4]/freq[i][i][iagemax+4])); |
| |
printf("%d%d ",i,k); |
| |
fprintf(ficlog,"%d%d ",i,k); |
| |
printf("jk=%d,i=%d,k=%d,p[%d]=%12.7f ln((%.0f/%.0f)/(%.0f/%.0f))= %f, OR=%f sd=%f \n",jk,i,k,jk,p[jk],freq[i][k][iagemax+3],freq[i][i][iagemax+3],freq[i][k][iagemax+4],freq[i][i][iagemax+4], log((freq[i][k][iagemax+3]/freq[i][i][iagemax+3])/(freq[i][k][iagemax+4]/freq[i][i][iagemax+4])),(freq[i][k][iagemax+3]/freq[i][i][iagemax+3])/(freq[i][k][iagemax+4]/freq[i][i][iagemax+4]), sqrt(1/freq[i][k][iagemax+3]+1/freq[i][i][iagemax+3]+1/freq[i][k][iagemax+4]+1/freq[i][i][iagemax+4])); |
| |
}else{ /* Other cases, like quantitative fixed or varying covariates */ |
| |
; |
| |
} |
| |
/* printf("%12.7f )", param[i][jj][k]); */ |
| |
/* fprintf(ficlog,"%12.7f )", param[i][jj][k]); */ |
| |
jk++; |
| |
} /* end jj */ |
| |
} /* end k!= i */ |
| |
} /* end k */ |
| |
} /* end i, jk */ |
| |
} /* end j !=0 */ |
| |
} /* end selected combination of covariate j1 */ |
| |
if(j==0){ /* We can estimate starting values from the occurences in each case */ |
| |
printf("#Freqsummary: Starting values for the constants:\n"); |
| |
fprintf(ficlog,"\n"); |
| |
for(i=1,jk=1; i <=nlstate; i++){ |
| |
for(k=1; k <=(nlstate+ndeath); k++){ |
| |
if (k != i) { |
| |
printf("%d%d ",i,k); |
| |
fprintf(ficlog,"%d%d ",i,k); |
| |
for(jj=1; jj <=ncovmodel; jj++){ |
| |
pstart[jk]=p[jk]; /* Setting pstart to p values by default */ |
| |
if(jj==1){ /* Age has to be done */ |
| |
pstart[jk]= log(freq[i][k][iagemax+3]/freq[i][i][iagemax+3]); |
| |
printf("%12.7f ln(%.0f/%.0f)= %12.7f ",p[jk],freq[i][k][iagemax+3],freq[i][i][iagemax+3], log(freq[i][k][iagemax+3]/freq[i][i][iagemax+3])); |
| |
fprintf(ficlog,"%12.7f ln(%.0f/%.0f)= %12.7f ",p[jk],freq[i][k][iagemax+3],freq[i][i][iagemax+3], log(freq[i][k][iagemax+3]/freq[i][i][iagemax+3])); |
| |
} |
| |
/* printf("%12.7f )", param[i][jj][k]); */ |
| |
/* fprintf(ficlog,"%12.7f )", param[i][jj][k]); */ |
| |
jk++; |
| } |
} |
| fprintf(ficlog," %d%d=%.0f",jk,m,freq[jk][m][iage]); |
printf("\n"); |
| |
fprintf(ficlog,"\n"); |
| } |
} |
| if(jk!=0 && m!=0) |
|
| fprintf(ficresphtmfr,"<td>%.0f</td> ",freq[jk][m][iage]); |
|
| } |
} |
| } /* end loop jk */ |
|
| posproptt=0.; |
|
| for(jk=1; jk <=nlstate; jk++){ |
|
| posproptt += pospropt[jk]; |
|
| } |
|
| fprintf(ficresphtmfr,"</tr>\n "); |
|
| if(iage <= iagemax){ |
|
| fprintf(ficresp,"\n"); |
|
| fprintf(ficresphtm,"</tr>\n"); |
|
| } |
} |
| if(first==1) |
printf("#Freqsummary\n"); |
| printf("Others in log...\n"); |
|
| fprintf(ficlog,"\n"); |
fprintf(ficlog,"\n"); |
| } /* end loop age iage */ |
for(jk=-1; jk <=nlstate+ndeath; jk++){ |
| fprintf(ficresphtm,"<tr><th>Tot</th>"); |
for(m=-1; m <=nlstate+ndeath; m++){ |
| for(jk=1; jk <=nlstate ; jk++){ |
/* param[i]|j][k]= freq[jk][m][iagemax+3] */ |
| if(posproptt < 1.e-5){ |
printf(" %d%d=%.0f",jk,m,freq[jk][m][iagemax+3]); |
| fprintf(ficresphtm,"<td>Nanq</td><td>%.0f</td><td>%.0f</td>",pospropt[jk],posproptt); |
fprintf(ficlog," %d%d=%.0f",jk,m,freq[jk][m][iagemax+3]); |
| }else{ |
/* if(freq[jk][m][iage] !=0 ) { /\* minimizing output *\/ */ |
| fprintf(ficresphtm,"<td>%.5f</td><td>%.0f</td><td>%.0f</td>",pospropt[jk]/posproptt,pospropt[jk],posproptt); |
/* printf(" %d%d=%.0f",jk,m,freq[jk][m][iagemax+3]); */ |
| |
/* fprintf(ficlog," %d%d=%.0f",jk,m,freq[jk][m][iagemax+3]); */ |
| |
/* } */ |
| |
} |
| |
} /* end loop jk */ |
| |
|
| |
printf("\n"); |
| |
fprintf(ficlog,"\n"); |
| |
} /* end j=0 */ |
| |
} /* end j */ |
| |
|
| |
if(mle == -2){ /* We want to use these values as starting values */ |
| |
for(i=1, jk=1; i <=nlstate; i++){ |
| |
for(j=1; j <=nlstate+ndeath; j++){ |
| |
if(j!=i){ |
| |
/*ca[0]= k+'a'-1;ca[1]='\0';*/ |
| |
printf("%1d%1d",i,j); |
| |
fprintf(ficparo,"%1d%1d",i,j); |
| |
for(k=1; k<=ncovmodel;k++){ |
| |
/* printf(" %lf",param[i][j][k]); */ |
| |
/* fprintf(ficparo," %lf",param[i][j][k]); */ |
| |
p[jk]=pstart[jk]; |
| |
printf(" %f ",pstart[jk]); |
| |
fprintf(ficparo," %f ",pstart[jk]); |
| |
jk++; |
| |
} |
| |
printf("\n"); |
| |
fprintf(ficparo,"\n"); |
| |
} |
| } |
} |
| } |
} |
| fprintf(ficresphtm,"</tr>\n"); |
} /* end mle=-2 */ |
| fprintf(ficresphtm,"</table>\n"); |
|
| fprintf(ficresphtmfr,"</table>\n"); |
|
| if(posproptt < 1.e-5){ |
|
| fprintf(ficresphtm,"\n <p><b> This combination (%d) is not valid and no result will be produced</b></p>",j1); |
|
| fprintf(ficresphtmfr,"\n <p><b> This combination (%d) is not valid and no result will be produced</b></p>",j1); |
|
| fprintf(ficres,"\n This combination (%d) is not valid and no result will be produced\n\n",j1); |
|
| invalidvarcomb[j1]=1; |
|
| }else{ |
|
| fprintf(ficresphtm,"\n <p> This combination (%d) is valid and result will be produced.</p>",j1); |
|
| invalidvarcomb[j1]=0; |
|
| } |
|
| fprintf(ficresphtmfr,"</table>\n"); |
|
| } /* end selected combination of covariate j1 */ |
|
| dateintmean=dateintsum/k2cpt; |
dateintmean=dateintsum/k2cpt; |
| |
|
| fclose(ficresp); |
fclose(ficresp); |
|
Line 4532 Title=%s <br>Datafile=%s Firstpass=%d La
|
Line 4770 Title=%s <br>Datafile=%s Firstpass=%d La
|
| fclose(ficresphtmfr); |
fclose(ficresphtmfr); |
| free_vector(meanq,1,nqfveff); |
free_vector(meanq,1,nqfveff); |
| free_matrix(meanqt,1,lastpass,1,nqtveff); |
free_matrix(meanqt,1,lastpass,1,nqtveff); |
| free_ma3x(freq,-5,nlstate+ndeath,-5,nlstate+ndeath, iagemin-AGEMARGE, iagemax+3+AGEMARGE); |
free_vector(x, iagemin-AGEMARGE, iagemax+4+AGEMARGE); |
| |
free_vector(y, iagemin-AGEMARGE, iagemax+4+AGEMARGE); |
| |
free_ma3x(freq,-5,nlstate+ndeath,-5,nlstate+ndeath, iagemin-AGEMARGE, iagemax+4+AGEMARGE); |
| free_vector(pospropt,1,nlstate); |
free_vector(pospropt,1,nlstate); |
| free_vector(posprop,1,nlstate); |
free_vector(posprop,1,nlstate); |
| free_matrix(prop,1,nlstate,iagemin-AGEMARGE, iagemax+3+AGEMARGE); |
free_matrix(prop,1,nlstate,iagemin-AGEMARGE, iagemax+4+AGEMARGE); |
| free_vector(pp,1,nlstate); |
free_vector(pp,1,nlstate); |
| /* End of freqsummary */ |
/* End of freqsummary */ |
| } |
} |
|
Line 4562 void prevalence(double ***probs, double
|
Line 4802 void prevalence(double ***probs, double
|
| iagemin= (int) agemin; |
iagemin= (int) agemin; |
| iagemax= (int) agemax; |
iagemax= (int) agemax; |
| /*pp=vector(1,nlstate);*/ |
/*pp=vector(1,nlstate);*/ |
| prop=matrix(1,nlstate,iagemin-AGEMARGE,iagemax+3+AGEMARGE); |
prop=matrix(1,nlstate,iagemin-AGEMARGE,iagemax+4+AGEMARGE); |
| /* freq=ma3x(-1,nlstate+ndeath,-1,nlstate+ndeath,iagemin,iagemax+3);*/ |
/* freq=ma3x(-1,nlstate+ndeath,-1,nlstate+ndeath,iagemin,iagemax+3);*/ |
| j1=0; |
j1=0; |
| |
|
|
Line 4572 void prevalence(double ***probs, double
|
Line 4812 void prevalence(double ***probs, double
|
| first=1; |
first=1; |
| for(j1=1; j1<= (int) pow(2,cptcoveff);j1++){ /* For each combination of covariate */ |
for(j1=1; j1<= (int) pow(2,cptcoveff);j1++){ /* For each combination of covariate */ |
| for (i=1; i<=nlstate; i++) |
for (i=1; i<=nlstate; i++) |
| for(iage=iagemin-AGEMARGE; iage <= iagemax+3+AGEMARGE; iage++) |
for(iage=iagemin-AGEMARGE; iage <= iagemax+4+AGEMARGE; iage++) |
| prop[i][iage]=0.0; |
prop[i][iage]=0.0; |
| printf("Prevalence combination of varying and fixed dummies %d\n",j1); |
printf("Prevalence combination of varying and fixed dummies %d\n",j1); |
| /* fprintf(ficlog," V%d=%d ",Tvaraff[j1],nbcode[Tvaraff[j1]][codtabm(k,j1)]); */ |
/* fprintf(ficlog," V%d=%d ",Tvaraff[j1],nbcode[Tvaraff[j1]][codtabm(k,j1)]); */ |
|
Line 4603 void prevalence(double ***probs, double
|
Line 4843 void prevalence(double ***probs, double
|
| if ((y2>=dateprev1) && (y2<=dateprev2)) { /* Here is the main selection (fractional years) */ |
if ((y2>=dateprev1) && (y2<=dateprev2)) { /* Here is the main selection (fractional years) */ |
| if(agev[m][i]==0) agev[m][i]=iagemax+1; |
if(agev[m][i]==0) agev[m][i]=iagemax+1; |
| if(agev[m][i]==1) agev[m][i]=iagemax+2; |
if(agev[m][i]==1) agev[m][i]=iagemax+2; |
| if((int)agev[m][i] <iagemin-AGEMARGE || (int)agev[m][i] >iagemax+3+AGEMARGE){ |
if((int)agev[m][i] <iagemin-AGEMARGE || (int)agev[m][i] >iagemax+4+AGEMARGE){ |
| printf("Error on individual # %d agev[m][i]=%f <%d-%d or > %d+3+%d m=%d; either change agemin or agemax or fix data\n",i, agev[m][i],iagemin,AGEMARGE, iagemax,AGEMARGE,m); |
printf("Error on individual # %d agev[m][i]=%f <%d-%d or > %d+3+%d m=%d; either change agemin or agemax or fix data\n",i, agev[m][i],iagemin,AGEMARGE, iagemax,AGEMARGE,m); |
| exit(1); |
exit(1); |
| } |
} |
|
Line 4640 void prevalence(double ***probs, double
|
Line 4880 void prevalence(double ***probs, double
|
| |
|
| /* free_ma3x(freq,-1,nlstate+ndeath,-1,nlstate+ndeath, iagemin, iagemax+3);*/ |
/* free_ma3x(freq,-1,nlstate+ndeath,-1,nlstate+ndeath, iagemin, iagemax+3);*/ |
| /*free_vector(pp,1,nlstate);*/ |
/*free_vector(pp,1,nlstate);*/ |
| free_matrix(prop,1,nlstate, iagemin-AGEMARGE,iagemax+3+AGEMARGE); |
free_matrix(prop,1,nlstate, iagemin-AGEMARGE,iagemax+4+AGEMARGE); |
| } /* End of prevalence */ |
} /* End of prevalence */ |
| |
|
| /************* Waves Concatenation ***************/ |
/************* Waves Concatenation ***************/ |
|
Line 6183 void printinghtml(char fileresu[], char
|
Line 6423 void printinghtml(char fileresu[], char
|
| |
|
| for(nres=1; nres <= nresult; nres++) /* For each resultline */ |
for(nres=1; nres <= nresult; nres++) /* For each resultline */ |
| for(k1=1; k1<=m;k1++){ /* For each combination of covariate */ |
for(k1=1; k1<=m;k1++){ /* For each combination of covariate */ |
| if(TKresult[nres]!= k1) |
if(m != 1 && TKresult[nres]!= k1) |
| continue; |
continue; |
| |
|
| /* for(i1=1; i1<=ncodemax[k1];i1++){ */ |
/* for(i1=1; i1<=ncodemax[k1];i1++){ */ |
|
Line 6234 divided by h: <sub>h</sub>P<sub>ij</sub>
|
Line 6474 divided by h: <sub>h</sub>P<sub>ij</sub>
|
| } |
} |
| /* Period (stable) prevalence in each health state */ |
/* Period (stable) prevalence in each health state */ |
| for(cpt=1; cpt<=nlstate;cpt++){ |
for(cpt=1; cpt<=nlstate;cpt++){ |
| fprintf(fichtm,"<br>\n- Convergence to period (stable) prevalence in state %d. Or probability to be in state %d being in state (1 to %d) at different ages. <a href=\"%s_%d-%d-%d.svg\">%s_%d-%d-%d.svg</a><br> \ |
fprintf(fichtm,"<br>\n- Convergence to period (stable) prevalence in state %d. Or probability to be in state %d some years earlier, knowing that we will be in state (1 to %d) at different ages. <a href=\"%s_%d-%d-%d.svg\">%s_%d-%d-%d.svg</a><br> \ |
| <img src=\"%s_%d-%d-%d.svg\">", cpt, cpt, nlstate, subdirf2(optionfilefiname,"P_"),cpt,k1,nres,subdirf2(optionfilefiname,"P_"),cpt,k1,nres,subdirf2(optionfilefiname,"P_"),cpt,k1,nres); |
<img src=\"%s_%d-%d-%d.svg\">", cpt, cpt, nlstate, subdirf2(optionfilefiname,"P_"),cpt,k1,nres,subdirf2(optionfilefiname,"P_"),cpt,k1,nres,subdirf2(optionfilefiname,"P_"),cpt,k1,nres); |
| } |
} |
| if(backcast==1){ |
if(backcast==1){ |
| /* Period (stable) back prevalence in each health state */ |
/* Period (stable) back prevalence in each health state */ |
| for(cpt=1; cpt<=nlstate;cpt++){ |
for(cpt=1; cpt<=nlstate;cpt++){ |
| fprintf(fichtm,"<br>\n- Convergence to period (stable) back prevalence in state %d. Or probability to be in state %d being in state (1 to %d) at different ages. <a href=\"%s_%d-%d-%d.svg\">%s_%d-%d-%d.svg</a><br> \ |
fprintf(fichtm,"<br>\n- Convergence to mixed (stable) back prevalence in state %d. Or probability to be in state %d at a younger age, knowing that we will be in state (1 to %d) at different older ages. <a href=\"%s_%d-%d-%d.svg\">%s_%d-%d-%d.svg</a><br> \ |
| <img src=\"%s_%d-%d-%d.svg\">", cpt, cpt, nlstate, subdirf2(optionfilefiname,"PB_"),cpt,k1,nres,subdirf2(optionfilefiname,"PB_"),cpt,k1,nres,subdirf2(optionfilefiname,"PB_"),cpt,k1,nres); |
<img src=\"%s_%d-%d-%d.svg\">", cpt, cpt, nlstate, subdirf2(optionfilefiname,"PB_"),cpt,k1,nres,subdirf2(optionfilefiname,"PB_"),cpt,k1,nres,subdirf2(optionfilefiname,"PB_"),cpt,k1,nres); |
| } |
} |
| } |
} |
|
Line 6315 See page 'Matrix of variance-covariance
|
Line 6555 See page 'Matrix of variance-covariance
|
| |
|
| for(nres=1; nres <= nresult; nres++){ /* For each resultline */ |
for(nres=1; nres <= nresult; nres++){ /* For each resultline */ |
| for(k1=1; k1<=m;k1++){ |
for(k1=1; k1<=m;k1++){ |
| if(TKresult[nres]!= k1) |
if(m != 1 && TKresult[nres]!= k1) |
| continue; |
continue; |
| /* for(i1=1; i1<=ncodemax[k1];i1++){ */ |
/* for(i1=1; i1<=ncodemax[k1];i1++){ */ |
| jj1++; |
jj1++; |
|
Line 6411 void printinggnuplot(char fileresu[], ch
|
Line 6651 void printinggnuplot(char fileresu[], ch
|
| for (k1=1; k1<= m ; k1 ++){ /* For each valid combination of covariate */ |
for (k1=1; k1<= m ; k1 ++){ /* For each valid combination of covariate */ |
| for(nres=1; nres <= nresult; nres++){ /* For each resultline */ |
for(nres=1; nres <= nresult; nres++){ /* For each resultline */ |
| /* plot [100000000000000000000:-100000000000000000000] "mysbiaspar/vplrmysbiaspar.txt to check */ |
/* plot [100000000000000000000:-100000000000000000000] "mysbiaspar/vplrmysbiaspar.txt to check */ |
| if(TKresult[nres]!= k1) |
if(m != 1 && TKresult[nres]!= k1) |
| continue; |
continue; |
| /* We are interested in selected combination by the resultline */ |
/* We are interested in selected combination by the resultline */ |
| /* printf("\n# 1st: Period (stable) prevalence with CI: 'VPL_' files and live state =%d ", cpt); */ |
/* printf("\n# 1st: Period (stable) prevalence with CI: 'VPL_' files and live state =%d ", cpt); */ |
|
Line 6493 void printinggnuplot(char fileresu[], ch
|
Line 6733 void printinggnuplot(char fileresu[], ch
|
| /*2 eme*/ |
/*2 eme*/ |
| for (k1=1; k1<= m ; k1 ++){ |
for (k1=1; k1<= m ; k1 ++){ |
| for(nres=1; nres <= nresult; nres++){ /* For each resultline */ |
for(nres=1; nres <= nresult; nres++){ /* For each resultline */ |
| if(TKresult[nres]!= k1) |
if(m != 1 && TKresult[nres]!= k1) |
| continue; |
continue; |
| fprintf(ficgp,"\n# 2nd: Total life expectancy with CI: 't' files "); |
fprintf(ficgp,"\n# 2nd: Total life expectancy with CI: 't' files "); |
| for (k=1; k<=cptcoveff; k++){ /* For each covariate and each value */ |
for (k=1; k<=cptcoveff; k++){ /* For each covariate and each value */ |
|
Line 6553 void printinggnuplot(char fileresu[], ch
|
Line 6793 void printinggnuplot(char fileresu[], ch
|
| /*3eme*/ |
/*3eme*/ |
| for (k1=1; k1<= m ; k1 ++){ |
for (k1=1; k1<= m ; k1 ++){ |
| for(nres=1; nres <= nresult; nres++){ /* For each resultline */ |
for(nres=1; nres <= nresult; nres++){ /* For each resultline */ |
| if(TKresult[nres]!= k1) |
if(m != 1 && TKresult[nres]!= k1) |
| continue; |
continue; |
| |
|
| for (cpt=1; cpt<= nlstate ; cpt ++) { |
for (cpt=1; cpt<= nlstate ; cpt ++) { |
|
Line 6602 plot [%.f:%.f] \"%s\" every :::%d::%d u
|
Line 6842 plot [%.f:%.f] \"%s\" every :::%d::%d u
|
| /* Survival functions (period) from state i in state j by initial state i */ |
/* Survival functions (period) from state i in state j by initial state i */ |
| for (k1=1; k1<=m; k1++){ /* For each covariate and each value */ |
for (k1=1; k1<=m; k1++){ /* For each covariate and each value */ |
| for(nres=1; nres <= nresult; nres++){ /* For each resultline */ |
for(nres=1; nres <= nresult; nres++){ /* For each resultline */ |
| if(TKresult[nres]!= k1) |
if(m != 1 && TKresult[nres]!= k1) |
| continue; |
continue; |
| for (cpt=1; cpt<=nlstate ; cpt ++) { /* For each life state cpt*/ |
for (cpt=1; cpt<=nlstate ; cpt ++) { /* For each life state cpt*/ |
| fprintf(ficgp,"\n#\n#\n# Survival functions in state j : 'LIJ_' files, cov=%d state=%d",k1, cpt); |
fprintf(ficgp,"\n#\n#\n# Survival functions in state j : 'LIJ_' files, cov=%d state=%d",k1, cpt); |
|
Line 6648 set ter svg size 640, 480\nunset log y\n
|
Line 6888 set ter svg size 640, 480\nunset log y\n
|
| /* Survival functions (period) from state i in state j by final state j */ |
/* Survival functions (period) from state i in state j by final state j */ |
| for (k1=1; k1<= m ; k1++){ /* For each covariate combination if any */ |
for (k1=1; k1<= m ; k1++){ /* For each covariate combination if any */ |
| for(nres=1; nres <= nresult; nres++){ /* For each resultline */ |
for(nres=1; nres <= nresult; nres++){ /* For each resultline */ |
| if(TKresult[nres]!= k1) |
if(m != 1 && TKresult[nres]!= k1) |
| continue; |
continue; |
| for (cpt=1; cpt<=nlstate ; cpt ++) { /* For each inital state */ |
for (cpt=1; cpt<=nlstate ; cpt ++) { /* For each inital state */ |
| fprintf(ficgp,"\n#\n#\n# Survival functions in state j and all livestates from state i by final state j: 'lij' files, cov=%d state=%d",k1, cpt); |
fprintf(ficgp,"\n#\n#\n# Survival functions in state j and all livestates from state i by final state j: 'lij' files, cov=%d state=%d",k1, cpt); |
|
Line 6702 set ter svg size 640, 480\nunset log y\n
|
Line 6942 set ter svg size 640, 480\nunset log y\n
|
| /* CV preval stable (period) for each covariate */ |
/* CV preval stable (period) for each covariate */ |
| for (k1=1; k1<= m ; k1 ++) /* For each covariate combination if any */ |
for (k1=1; k1<= m ; k1 ++) /* For each covariate combination if any */ |
| for(nres=1; nres <= nresult; nres++){ /* For each resultline */ |
for(nres=1; nres <= nresult; nres++){ /* For each resultline */ |
| if(TKresult[nres]!= k1) |
if(m != 1 && TKresult[nres]!= k1) |
| continue; |
continue; |
| for (cpt=1; cpt<=nlstate ; cpt ++) { /* For each life state */ |
for (cpt=1; cpt<=nlstate ; cpt ++) { /* For each life state of arrival */ |
| |
|
| fprintf(ficgp,"\n#\n#\n#CV preval stable (period): 'pij' files, covariatecombination#=%d state=%d",k1, cpt); |
fprintf(ficgp,"\n#\n#\n#CV preval stable (period): 'pij' files, covariatecombination#=%d state=%d",k1, cpt); |
| for (k=1; k<=cptcoveff; k++){ /* For each covariate and each value */ |
for (k=1; k<=cptcoveff; k++){ /* For each covariate and each value */ |
|
Line 6728 set ter svg size 640, 480\nunset log y\n
|
Line 6968 set ter svg size 640, 480\nunset log y\n
|
| fprintf(ficgp,"set xlabel \"Age\" \nset ylabel \"Probability\" \n\ |
fprintf(ficgp,"set xlabel \"Age\" \nset ylabel \"Probability\" \n\ |
| set ter svg size 640, 480\nunset log y\nplot [%.f:%.f] ", ageminpar, agemaxpar); |
set ter svg size 640, 480\nunset log y\nplot [%.f:%.f] ", ageminpar, agemaxpar); |
| k=3; /* Offset */ |
k=3; /* Offset */ |
| for (i=1; i<= nlstate ; i ++){ |
for (i=1; i<= nlstate ; i ++){ /* State of origin */ |
| if(i==1) |
if(i==1) |
| fprintf(ficgp,"\"%s\"",subdirf2(fileresu,"PIJ_")); |
fprintf(ficgp,"\"%s\"",subdirf2(fileresu,"PIJ_")); |
| else |
else |
| fprintf(ficgp,", '' "); |
fprintf(ficgp,", '' "); |
| l=(nlstate+ndeath)*(i-1)+1; |
l=(nlstate+ndeath)*(i-1)+1; /* 1, 1+ nlstate+ndeath, 1+2*(nlstate+ndeath) */ |
| fprintf(ficgp," u ($1==%d ? ($3):1/0):($%d/($%d",k1,k+l+(cpt-1),k+l); |
fprintf(ficgp," u ($1==%d ? ($3):1/0):($%d/($%d",k1,k+l+(cpt-1),k+l); |
| for (j=2; j<= nlstate ; j ++) |
for (j=2; j<= nlstate ; j ++) |
| fprintf(ficgp,"+$%d",k+l+j-1); |
fprintf(ficgp,"+$%d",k+l+j-1); |
|
Line 6749 set ter svg size 640, 480\nunset log y\n
|
Line 6989 set ter svg size 640, 480\nunset log y\n
|
| /* CV back preval stable (period) for each covariate */ |
/* CV back preval stable (period) for each covariate */ |
| for (k1=1; k1<= m ; k1 ++) /* For each covariate combination if any */ |
for (k1=1; k1<= m ; k1 ++) /* For each covariate combination if any */ |
| for(nres=1; nres <= nresult; nres++){ /* For each resultline */ |
for(nres=1; nres <= nresult; nres++){ /* For each resultline */ |
| if(TKresult[nres]!= k1) |
if(m != 1 && TKresult[nres]!= k1) |
| continue; |
continue; |
| for (cpt=1; cpt<=nlstate ; cpt ++) { /* For each life state */ |
for (cpt=1; cpt<=nlstate ; cpt ++) { /* For each life ending state */ |
| fprintf(ficgp,"\n#\n#\n#CV Back preval stable (period): 'pij' files, covariatecombination#=%d state=%d",k1, cpt); |
fprintf(ficgp,"\n#\n#\n#CV Back preval stable (period): 'pijb' files, covariatecombination#=%d state=%d",k1, cpt); |
| for (k=1; k<=cptcoveff; k++){ /* For each covariate and each value */ |
for (k=1; k<=cptcoveff; k++){ /* For each covariate and each value */ |
| lv= decodtabm(k1,k,cptcoveff); /* Should be the covariate number corresponding to k1 combination */ |
lv= decodtabm(k1,k,cptcoveff); /* Should be the covariate number corresponding to k1 combination */ |
| /* decodtabm(1,1,4) = 1 because h=1 k= (1) 1 1 1 */ |
/* decodtabm(1,1,4) = 1 because h=1 k= (1) 1 1 1 */ |
|
Line 6774 set ter svg size 640, 480\nunset log y\n
|
Line 7014 set ter svg size 640, 480\nunset log y\n
|
| fprintf(ficgp,"set xlabel \"Age\" \nset ylabel \"Probability\" \n\ |
fprintf(ficgp,"set xlabel \"Age\" \nset ylabel \"Probability\" \n\ |
| set ter svg size 640, 480\nunset log y\nplot [%.f:%.f] ", ageminpar, agemaxpar); |
set ter svg size 640, 480\nunset log y\nplot [%.f:%.f] ", ageminpar, agemaxpar); |
| k=3; /* Offset */ |
k=3; /* Offset */ |
| for (i=1; i<= nlstate ; i ++){ |
for (i=1; i<= nlstate ; i ++){ /* State of origin */ |
| if(i==1) |
if(i==1) |
| fprintf(ficgp,"\"%s\"",subdirf2(fileresu,"PIJB_")); |
fprintf(ficgp,"\"%s\"",subdirf2(fileresu,"PIJB_")); |
| else |
else |
| fprintf(ficgp,", '' "); |
fprintf(ficgp,", '' "); |
| /* l=(nlstate+ndeath)*(i-1)+1; */ |
/* l=(nlstate+ndeath)*(i-1)+1; */ |
| l=(nlstate+ndeath)*(cpt-1)+1; |
l=(nlstate+ndeath)*(cpt-1)+1; /* fixed for i; cpt=1 1, cpt=2 1+ nlstate+ndeath, 1+2*(nlstate+ndeath) */ |
| /* fprintf(ficgp," u ($1==%d ? ($3):1/0):($%d/($%d",k1,k+l+(cpt-1),k+l); /\* a vérifier *\/ */ |
/* fprintf(ficgp," u ($1==%d ? ($3):1/0):($%d/($%d",k1,k+l+(cpt-1),k+l); /\* a vérifier *\/ */ |
| /* fprintf(ficgp," u ($1==%d ? ($3):1/0):($%d/($%d",k1,k+l+(cpt-1),k+l+(cpt-1)+i-1); /\* a vérifier *\/ */ |
/* fprintf(ficgp," u ($1==%d ? ($3):1/0):($%d/($%d",k1,k+l+(cpt-1),k+l+(cpt-1)+i-1); /\* a vérifier *\/ */ |
| fprintf(ficgp," u ($1==%d ? ($3):1/0):($%d",k1,k+l+(cpt-1)+i-1); /* a vérifier */ |
fprintf(ficgp," u ($1==%d ? ($3):1/0):($%d",k1,k+l+i-1); /* To be verified */ |
| /* for (j=2; j<= nlstate ; j ++) */ |
/* for (j=2; j<= nlstate ; j ++) */ |
| /* fprintf(ficgp,"+$%d",k+l+j-1); */ |
/* fprintf(ficgp,"+$%d",k+l+j-1); */ |
| /* /\* fprintf(ficgp,"+$%d",k+l+j-1); *\/ */ |
/* /\* fprintf(ficgp,"+$%d",k+l+j-1); *\/ */ |
|
Line 6800 set ter svg size 640, 480\nunset log y\n
|
Line 7040 set ter svg size 640, 480\nunset log y\n
|
| |
|
| for (k1=1; k1<= m ; k1 ++) /* For each covariate combination if any */ |
for (k1=1; k1<= m ; k1 ++) /* For each covariate combination if any */ |
| for(nres=1; nres <= nresult; nres++){ /* For each resultline */ |
for(nres=1; nres <= nresult; nres++){ /* For each resultline */ |
| if(TKresult[nres]!= k1) |
if(m != 1 && TKresult[nres]!= k1) |
| continue; |
continue; |
| for (cpt=1; cpt<=nlstate ; cpt ++) { /* For each life state */ |
for (cpt=1; cpt<=nlstate ; cpt ++) { /* For each life state */ |
| fprintf(ficgp,"\n#\n#\n#Projection of prevalence to stable (period): 'PROJ_' files, covariatecombination#=%d state=%d",k1, cpt); |
fprintf(ficgp,"\n#\n#\n#Projection of prevalence to stable (period): 'PROJ_' files, covariatecombination#=%d state=%d",k1, cpt); |
|
Line 6931 set ter svg size 640, 480\nunset log y\n
|
Line 7171 set ter svg size 640, 480\nunset log y\n
|
| fprintf(ficgp,"# jk=1 to 2^%d=%d\n",cptcoveff,m);/* to be checked */ |
fprintf(ficgp,"# jk=1 to 2^%d=%d\n",cptcoveff,m);/* to be checked */ |
| for(jk=1; jk <=m; jk++) /* For each combination of covariate */ |
for(jk=1; jk <=m; jk++) /* For each combination of covariate */ |
| for(nres=1; nres <= nresult; nres++){ /* For each resultline */ |
for(nres=1; nres <= nresult; nres++){ /* For each resultline */ |
| if(TKresult[nres]!= jk) |
if(m != 1 && TKresult[nres]!= jk) |
| continue; |
continue; |
| fprintf(ficgp,"# Combination of dummy jk=%d and ",jk); |
fprintf(ficgp,"# Combination of dummy jk=%d and ",jk); |
| for (k4=1; k4<= nsq; k4++){ /* For each selected (single) quantitative value */ |
for (k4=1; k4<= nsq; k4++){ /* For each selected (single) quantitative value */ |
|
Line 7293 set ter svg size 640, 480\nunset log y\n
|
Line 7533 set ter svg size 640, 480\nunset log y\n
|
| /* if (h==(int)(YEARM*yearp)){ */ |
/* if (h==(int)(YEARM*yearp)){ */ |
| for(nres=1; nres <= nresult; nres++) /* For each resultline */ |
for(nres=1; nres <= nresult; nres++) /* For each resultline */ |
| for(k=1; k<=i1;k++){ |
for(k=1; k<=i1;k++){ |
| if(TKresult[nres]!= k) |
if(i1 != 1 && TKresult[nres]!= k) |
| continue; |
continue; |
| if(invalidvarcomb[k]){ |
if(invalidvarcomb[k]){ |
| printf("\nCombination (%d) projection ignored because no cases \n",k); |
printf("\nCombination (%d) projection ignored because no cases \n",k); |
|
Line 8244 int decoderesult ( char resultline[], in
|
Line 8484 int decoderesult ( char resultline[], in
|
| if (strlen(resultsav) >1){ |
if (strlen(resultsav) >1){ |
| j=nbocc(resultsav,'='); /**< j=Number of covariate values'=' */ |
j=nbocc(resultsav,'='); /**< j=Number of covariate values'=' */ |
| } |
} |
| |
if(j == 0){ /* Resultline but no = */ |
| |
TKresult[nres]=0; /* Combination for the nresult and the model */ |
| |
return (0); |
| |
} |
| |
|
| if( j != cptcovs ){ /* Be careful if a variable is in a product but not single */ |
if( j != cptcovs ){ /* Be careful if a variable is in a product but not single */ |
| printf("ERROR: the number of variable in the resultline, %d, differs from the number of variable used in the model line, %d.\n",j, cptcovs); |
printf("ERROR: the number of variable in the resultline, %d, differs from the number of variable used in the model line, %d.\n",j, cptcovs); |
| fprintf(ficlog,"ERROR: the number of variable in the resultline, %d, differs from the number of variable used in the model line, %d.\n",j, cptcovs); |
fprintf(ficlog,"ERROR: the number of variable in the resultline, %d, differs from the number of variable used in the model line, %d.\n",j, cptcovs); |
|
Line 8843 Dummy[k] 0=dummy (0 1), 1 quantitative (
|
Line 9088 Dummy[k] 0=dummy (0 1), 1 quantitative (
|
| } |
} |
| |
|
| int calandcheckages(int imx, int maxwav, double *agemin, double *agemax, int *nberr, int *nbwarn ) |
int calandcheckages(int imx, int maxwav, double *agemin, double *agemax, int *nberr, int *nbwarn ) |
| { |
{/* Check ages at death */ |
| int i, m; |
int i, m; |
| int firstone=0; |
int firstone=0; |
| |
|
|
Line 9180 int prevalence_limit(double *p, double *
|
Line 9425 int prevalence_limit(double *p, double *
|
| |
|
| for(k=1; k<=i1;k++){ /* For each combination k of dummy covariates in the model */ |
for(k=1; k<=i1;k++){ /* For each combination k of dummy covariates in the model */ |
| for(nres=1; nres <= nresult; nres++){ /* For each resultline */ |
for(nres=1; nres <= nresult; nres++){ /* For each resultline */ |
| if(TKresult[nres]!= k) |
if(i1 != 1 && TKresult[nres]!= k) |
| continue; |
continue; |
| |
|
| /* for(cptcov=1,k=0;cptcov<=i1;cptcov++){ */ |
/* for(cptcov=1,k=0;cptcov<=i1;cptcov++){ */ |
|
Line 9277 int back_prevalence_limit(double *p, dou
|
Line 9522 int back_prevalence_limit(double *p, dou
|
| |
|
| for(nres=1; nres <= nresult; nres++){ /* For each resultline */ |
for(nres=1; nres <= nresult; nres++){ /* For each resultline */ |
| for(k=1; k<=i1;k++){ /* For any combination of dummy covariates, fixed and varying */ |
for(k=1; k<=i1;k++){ /* For any combination of dummy covariates, fixed and varying */ |
| if(TKresult[nres]!= k) |
if(i1 != 1 && TKresult[nres]!= k) |
| continue; |
continue; |
| //printf("cptcov=%d cptcod=%d codtab=%d\n",cptcov, cptcod,codtabm(cptcod,cptcov)); |
//printf("cptcov=%d cptcod=%d codtab=%d\n",cptcov, cptcod,codtabm(cptcod,cptcov)); |
| fprintf(ficresplb,"#******"); |
fprintf(ficresplb,"#******"); |
|
Line 9336 int back_prevalence_limit(double *p, dou
|
Line 9581 int back_prevalence_limit(double *p, dou
|
| fprintf(ficresplb," %.3f %d\n", tot, *ncvyearp); |
fprintf(ficresplb," %.3f %d\n", tot, *ncvyearp); |
| } /* Age */ |
} /* Age */ |
| /* was end of cptcod */ |
/* was end of cptcod */ |
| |
/*fprintf(ficresplb,"\n");*/ /* Seems to be necessary for gnuplot only if two result lines and no covariate. */ |
| } /* end of any combination */ |
} /* end of any combination */ |
| } /* end of nres */ |
} /* end of nres */ |
| /* hBijx(p, bage, fage); */ |
/* hBijx(p, bage, fage); */ |
|
Line 9380 int hPijx(double *p, int bage, int fage)
|
Line 9626 int hPijx(double *p, int bage, int fage)
|
| /* k=k+1; */ |
/* k=k+1; */ |
| for(nres=1; nres <= nresult; nres++) /* For each resultline */ |
for(nres=1; nres <= nresult; nres++) /* For each resultline */ |
| for(k=1; k<=i1;k++){ |
for(k=1; k<=i1;k++){ |
| if(TKresult[nres]!= k) |
if(i1 != 1 && TKresult[nres]!= k) |
| continue; |
continue; |
| fprintf(ficrespij,"\n#****** "); |
fprintf(ficrespij,"\n#****** "); |
| for(j=1;j<=cptcoveff;j++) |
for(j=1;j<=cptcoveff;j++) |
|
Line 9452 int hPijx(double *p, int bage, int fage)
|
Line 9698 int hPijx(double *p, int bage, int fage)
|
| |
|
| /* hstepm=1; aff par mois*/ |
/* hstepm=1; aff par mois*/ |
| pstamp(ficrespijb); |
pstamp(ficrespijb); |
| fprintf(ficrespijb,"#****** h Pij x Back Probability to be in state i at age x-h being in j at x "); |
fprintf(ficrespijb,"#****** h Bij x Back probability to be in state i at age x-h being in j at x: B1j+B2j+...=1 "); |
| i1= pow(2,cptcoveff); |
i1= pow(2,cptcoveff); |
| /* for(cptcov=1,k=0;cptcov<=i1;cptcov++){ */ |
/* for(cptcov=1,k=0;cptcov<=i1;cptcov++){ */ |
| /* /\*for(cptcod=1;cptcod<=ncodemax[cptcov];cptcod++){*\/ */ |
/* /\*for(cptcod=1;cptcod<=ncodemax[cptcov];cptcod++){*\/ */ |
| /* k=k+1; */ |
/* k=k+1; */ |
| for(nres=1; nres <= nresult; nres++){ /* For each resultline */ |
for(nres=1; nres <= nresult; nres++){ /* For each resultline */ |
| for(k=1; k<=i1;k++){ /* For any combination of dummy covariates, fixed and varying */ |
for(k=1; k<=i1;k++){ /* For any combination of dummy covariates, fixed and varying */ |
| if(TKresult[nres]!= k) |
if(i1 != 1 && TKresult[nres]!= k) |
| continue; |
continue; |
| fprintf(ficrespijb,"\n#****** "); |
fprintf(ficrespijb,"\n#****** "); |
| for(j=1;j<=cptcoveff;j++) |
for(j=1;j<=cptcoveff;j++) |
|
Line 9486 int hPijx(double *p, int bage, int fage)
|
Line 9732 int hPijx(double *p, int bage, int fage)
|
| /* hbxij(p3mat,nhstepm,agedeb,hstepm,p,nlstate,stepm,oldm,savm, k); */ |
/* hbxij(p3mat,nhstepm,agedeb,hstepm,p,nlstate,stepm,oldm,savm, k); */ |
| hbxij(p3mat,nhstepm,agedeb,hstepm,p,prevacurrent,nlstate,stepm, k); |
hbxij(p3mat,nhstepm,agedeb,hstepm,p,prevacurrent,nlstate,stepm, k); |
| /* hbxij(p3mat,nhstepm,agedeb,hstepm,p,prevacurrent,nlstate,stepm,oldm,savm, dnewm, doldm, dsavm, k); */ |
/* hbxij(p3mat,nhstepm,agedeb,hstepm,p,prevacurrent,nlstate,stepm,oldm,savm, dnewm, doldm, dsavm, k); */ |
| fprintf(ficrespijb,"# Cov Agex agex-h hpijx with i,j="); |
fprintf(ficrespijb,"# Cov Agex agex-h hbijx with i,j="); |
| for(i=1; i<=nlstate;i++) |
for(i=1; i<=nlstate;i++) |
| for(j=1; j<=nlstate+ndeath;j++) |
for(j=1; j<=nlstate+ndeath;j++) |
| fprintf(ficrespijb," %1d-%1d",i,j); |
fprintf(ficrespijb," %1d-%1d",i,j); |
|
Line 9576 int main(int argc, char *argv[])
|
Line 9822 int main(int argc, char *argv[])
|
| double **prlim; |
double **prlim; |
| double **bprlim; |
double **bprlim; |
| double ***param; /* Matrix of parameters */ |
double ***param; /* Matrix of parameters */ |
| double *p; |
double ***paramstart; /* Matrix of starting parameter values */ |
| |
double *p, *pstart; /* p=param[1][1] pstart is for starting values guessed by freqsummary */ |
| double **matcov; /* Matrix of covariance */ |
double **matcov; /* Matrix of covariance */ |
| double **hess; /* Hessian matrix */ |
double **hess; /* Hessian matrix */ |
| double ***delti3; /* Scale */ |
double ***delti3; /* Scale */ |
|
Line 9922 int main(int argc, char *argv[])
|
Line 10169 int main(int argc, char *argv[])
|
| ungetc(c,ficpar); |
ungetc(c,ficpar); |
| |
|
| param= ma3x(1,nlstate,1,nlstate+ndeath-1,1,ncovmodel); |
param= ma3x(1,nlstate,1,nlstate+ndeath-1,1,ncovmodel); |
| |
paramstart= ma3x(1,nlstate,1,nlstate+ndeath-1,1,ncovmodel); |
| for(i=1; i <=nlstate; i++){ |
for(i=1; i <=nlstate; i++){ |
| j=0; |
j=0; |
| for(jj=1; jj <=nlstate+ndeath; jj++){ |
for(jj=1; jj <=nlstate+ndeath; jj++){ |
|
Line 9958 run imach with mle=-1 to get a correct t
|
Line 10206 run imach with mle=-1 to get a correct t
|
| } |
} |
| fflush(ficlog); |
fflush(ficlog); |
| |
|
| /* Reads scales values */ |
/* Reads parameters values */ |
| p=param[1][1]; |
p=param[1][1]; |
| |
pstart=paramstart[1][1]; |
| |
|
| /* Reads comments: lines beginning with '#' */ |
/* Reads comments: lines beginning with '#' */ |
| while((c=getc(ficpar))=='#' && c!= EOF){ |
while((c=getc(ficpar))=='#' && c!= EOF){ |
|
Line 10224 Please run with mle=-1 to get a correct
|
Line 10473 Please run with mle=-1 to get a correct
|
| */ |
*/ |
| |
|
| concatwav(wav, dh, bh, mw, s, agedc, agev, firstpass, lastpass, imx, nlstate, stepm); |
concatwav(wav, dh, bh, mw, s, agedc, agev, firstpass, lastpass, imx, nlstate, stepm); |
| /* */ |
/* Concatenates waves */ |
| |
|
| free_vector(moisdc,1,n); |
free_vector(moisdc,1,n); |
| free_vector(andc,1,n); |
free_vector(andc,1,n); |
|
Line 10405 Title=%s <br>Datafile=%s Firstpass=%d La
|
Line 10654 Title=%s <br>Datafile=%s Firstpass=%d La
|
| /* Calculates basic frequencies. Computes observed prevalence at single age |
/* Calculates basic frequencies. Computes observed prevalence at single age |
| and for any valid combination of covariates |
and for any valid combination of covariates |
| and prints on file fileres'p'. */ |
and prints on file fileres'p'. */ |
| freqsummary(fileres, agemin, agemax, s, agev, nlstate, imx, Tvaraff, invalidvarcomb, nbcode, ncodemax,mint,anint,strstart, \ |
freqsummary(fileres, p, pstart, agemin, agemax, s, agev, nlstate, imx, Tvaraff, invalidvarcomb, nbcode, ncodemax,mint,anint,strstart, \ |
| firstpass, lastpass, stepm, weightopt, model); |
firstpass, lastpass, stepm, weightopt, model); |
| |
|
| fprintf(fichtm,"\n"); |
fprintf(fichtm,"\n"); |
|
Line 10427 Interval (in months) between two waves:
|
Line 10676 Interval (in months) between two waves:
|
| /* For mortality only */ |
/* For mortality only */ |
| if (mle==-3){ |
if (mle==-3){ |
| ximort=matrix(1,NDIM,1,NDIM); |
ximort=matrix(1,NDIM,1,NDIM); |
| for(i=1;i<=NDIM;i++) |
for(i=1;i<=NDIM;i++) |
| for(j=1;j<=NDIM;j++) |
for(j=1;j<=NDIM;j++) |
| ximort[i][j]=0.; |
ximort[i][j]=0.; |
| /* ximort=gsl_matrix_alloc(1,NDIM,1,NDIM); */ |
/* ximort=gsl_matrix_alloc(1,NDIM,1,NDIM); */ |
| cens=ivector(1,n); |
cens=ivector(1,n); |
| ageexmed=vector(1,n); |
ageexmed=vector(1,n); |
|
Line 10884 Please run with mle=-1 to get a correct
|
Line 11133 Please run with mle=-1 to get a correct
|
| fprintf(ficparo,"agemin=%.0f agemax=%.0f bage=%.0f fage=%.0f estepm=%d, ftolpl=%e\n",ageminpar,agemaxpar,bage,fage, estepm, ftolpl); |
fprintf(ficparo,"agemin=%.0f agemax=%.0f bage=%.0f fage=%.0f estepm=%d, ftolpl=%e\n",ageminpar,agemaxpar,bage,fage, estepm, ftolpl); |
| |
|
| /* Other stuffs, more or less useful */ |
/* Other stuffs, more or less useful */ |
| while((c=getc(ficpar))=='#' && c!= EOF){ |
while(fgets(line, MAXLINE, ficpar)) { |
| ungetc(c,ficpar); |
/* If line starts with a # it is a comment */ |
| fgets(line, MAXLINE, ficpar); |
if (line[0] == '#') { |
| fputs(line,stdout); |
numlinepar++; |
| fputs(line,ficparo); |
fputs(line,stdout); |
| |
fputs(line,ficparo); |
| |
fputs(line,ficlog); |
| |
continue; |
| |
}else |
| |
break; |
| } |
} |
| ungetc(c,ficpar); |
|
| |
if((num_filled=sscanf(line,"begin-prev-date=%lf/%lf/%lf end-prev-date=%lf/%lf/%lf mov_average=%d\n",&jprev1, &mprev1,&anprev1,&jprev2, &mprev2,&anprev2,&mobilav)) !=EOF){ |
| fscanf(ficpar,"begin-prev-date=%lf/%lf/%lf end-prev-date=%lf/%lf/%lf mov_average=%d\n",&jprev1, &mprev1,&anprev1,&jprev2, &mprev2,&anprev2,&mobilav); |
|
| fprintf(ficparo,"begin-prev-date=%.lf/%.lf/%.lf end-prev-date=%.lf/%.lf/%.lf mov_average=%d\n",jprev1, mprev1,anprev1,jprev2, mprev2,anprev2,mobilav); |
if (num_filled != 7) { |
| fprintf(ficres,"begin-prev-date=%.lf/%.lf/%.lf end-prev-date=%.lf/%.lf/%.lf mov_average=%d\n",jprev1, mprev1,anprev1,jprev2, mprev2,anprev2,mobilav); |
printf("Error: Not 7 (data)parameters in line but %d, for example:begin-prev-date=1/1/1990 end-prev-date=1/6/2004 mov_average=0\n, your line=%s . Probably you are running an older format.\n",num_filled,line); |
| printf("begin-prev-date=%.lf/%.lf/%.lf end-prev-date=%.lf/%.lf/%.lf mov_average=%d\n",jprev1, mprev1,anprev1,jprev2, mprev2,anprev2,mobilav); |
fprintf(ficlog,"Error: Not 7 (data)parameters in line but %d, for example:begin-prev-date=1/1/1990 end-prev-date=1/6/2004 mov_average=0\n, your line=%s . Probably you are running an older format.\n",num_filled,line); |
| fprintf(ficlog,"begin-prev-date=%.lf/%.lf/%.lf end-prev-date=%.lf/%.lf/%.lf mov_average=%d\n",jprev1, mprev1,anprev1,jprev2, mprev2,anprev2,mobilav); |
goto end; |
| |
} |
| while((c=getc(ficpar))=='#' && c!= EOF){ |
/* fscanf(ficpar,"begin-prev-date=%lf/%lf/%lf end-prev-date=%lf/%lf/%lf mov_average=%d\n",&jprev1, &mprev1,&anprev1,&jprev2, &mprev2,&anprev2,&mobilav); */ |
| ungetc(c,ficpar); |
printf("begin-prev-date=%.lf/%.lf/%.lf end-prev-date=%.lf/%.lf/%.lf mov_average=%d\n",jprev1, mprev1,anprev1,jprev2, mprev2,anprev2,mobilav); |
| fgets(line, MAXLINE, ficpar); |
fprintf(ficparo,"begin-prev-date=%.lf/%.lf/%.lf end-prev-date=%.lf/%.lf/%.lf mov_average=%d\n",jprev1, mprev1,anprev1,jprev2, mprev2,anprev2,mobilav); |
| fputs(line,stdout); |
fprintf(ficres,"begin-prev-date=%.lf/%.lf/%.lf end-prev-date=%.lf/%.lf/%.lf mov_average=%d\n",jprev1, mprev1,anprev1,jprev2, mprev2,anprev2,mobilav); |
| fputs(line,ficparo); |
printf("begin-prev-date=%.lf/%.lf/%.lf end-prev-date=%.lf/%.lf/%.lf mov_average=%d\n",jprev1, mprev1,anprev1,jprev2, mprev2,anprev2,mobilav); |
| |
fprintf(ficlog,"begin-prev-date=%.lf/%.lf/%.lf end-prev-date=%.lf/%.lf/%.lf mov_average=%d\n",jprev1, mprev1,anprev1,jprev2, mprev2,anprev2,mobilav); |
| |
} |
| |
|
| |
while(fgets(line, MAXLINE, ficpar)) { |
| |
/* If line starts with a # it is a comment */ |
| |
if (line[0] == '#') { |
| |
numlinepar++; |
| |
fputs(line,stdout); |
| |
fputs(line,ficparo); |
| |
fputs(line,ficlog); |
| |
continue; |
| |
}else |
| |
break; |
| } |
} |
| ungetc(c,ficpar); |
|
| |
|
| |
|
| dateprev1=anprev1+(mprev1-1)/12.+(jprev1-1)/365.; |
dateprev1=anprev1+(mprev1-1)/12.+(jprev1-1)/365.; |
| dateprev2=anprev2+(mprev2-1)/12.+(jprev2-1)/365.; |
dateprev2=anprev2+(mprev2-1)/12.+(jprev2-1)/365.; |
| |
|
| fscanf(ficpar,"pop_based=%d\n",&popbased); |
if((num_filled=sscanf(line,"pop_based=%d\n",&popbased)) !=EOF){ |
| fprintf(ficlog,"pop_based=%d\n",popbased); |
if (num_filled != 1) { |
| fprintf(ficparo,"pop_based=%d\n",popbased); |
printf("Error: Not 1 (data)parameters in line but %d, for example:pop_based=0\n, your line=%s . Probably you are running an older format.\n",num_filled,line); |
| fprintf(ficres,"pop_based=%d\n",popbased); |
fprintf(ficlog,"Error: Not 1 (data)parameters in line but %d, for example: pop_based=1\n, your line=%s . Probably you are running an older format.\n",num_filled,line); |
| |
goto end; |
| while((c=getc(ficpar))=='#' && c!= EOF){ |
} |
| ungetc(c,ficpar); |
printf("pop_based=%d\n",popbased); |
| fgets(line, MAXLINE, ficpar); |
fprintf(ficlog,"pop_based=%d\n",popbased); |
| fputs(line,stdout); |
fprintf(ficparo,"pop_based=%d\n",popbased); |
| fputs(line,ficres); |
fprintf(ficres,"pop_based=%d\n",popbased); |
| fputs(line,ficparo); |
|
| } |
} |
| ungetc(c,ficpar); |
|
| |
while(fgets(line, MAXLINE, ficpar)) { |
| fscanf(ficpar,"prevforecast=%d starting-proj-date=%lf/%lf/%lf final-proj-date=%lf/%lf/%lf mobil_average=%d\n",&prevfcast,&jproj1,&mproj1,&anproj1,&jproj2,&mproj2,&anproj2,&mobilavproj); |
/* If line starts with a # it is a comment */ |
| fprintf(ficparo,"prevforecast=%d starting-proj-date=%.lf/%.lf/%.lf final-proj-date=%.lf/%.lf/%.lf mobil_average=%d\n",prevfcast,jproj1,mproj1,anproj1,jproj2,mproj2,anproj2,mobilavproj); |
if (line[0] == '#') { |
| printf("prevforecast=%d starting-proj-date=%.lf/%.lf/%.lf final-proj-date=%.lf/%.lf/%.lf mobil_average=%d\n",prevfcast,jproj1,mproj1,anproj1,jproj2,mproj2,anproj2,mobilavproj); |
numlinepar++; |
| fprintf(ficlog,"prevforecast=%d starting-proj-date=%.lf/%.lf/%.lf final-proj-date=%.lf/%.lf/%.lf mobil_average=%d\n",prevfcast,jproj1,mproj1,anproj1,jproj2,mproj2,anproj2,mobilavproj); |
fputs(line,stdout); |
| fprintf(ficres,"prevforecast=%d starting-proj-date=%.lf/%.lf/%.lf final-proj-date=%.lf/%.lf/%.lf mobil_average=%d\n",prevfcast,jproj1,mproj1,anproj1,jproj2,mproj2,anproj2,mobilavproj); |
fputs(line,ficparo); |
| |
fputs(line,ficlog); |
| |
continue; |
| |
}else |
| |
break; |
| |
} |
| |
/* while((c=getc(ficpar))=='#' && c!= EOF){ */ |
| |
/* ungetc(c,ficpar); */ |
| |
/* fgets(line, MAXLINE, ficpar); */ |
| |
/* fputs(line,stdout); */ |
| |
/* fputs(line,ficres); */ |
| |
/* fputs(line,ficparo); */ |
| |
/* } */ |
| |
/* ungetc(c,ficpar); */ |
| |
|
| |
/* fscanf(ficpar,"prevforecast=%d starting-proj-date=%lf/%lf/%lf final-proj-date=%lf/%lf/%lf mobil_average=%d\n",&prevfcast,&jproj1,&mproj1,&anproj1,&jproj2,&mproj2,&anproj2,&mobilavproj); */ |
| |
if((num_filled=sscanf(line,"prevforecast=%d starting-proj-date=%lf/%lf/%lf final-proj-date=%lf/%lf/%lf mobil_average=%d\n",&prevfcast,&jproj1,&mproj1,&anproj1,&jproj2,&mproj2,&anproj2,&mobilavproj)) !=EOF){ |
| |
if (num_filled != 8) { |
| |
printf("Error: Not 8 (data)parameters in line but %d, for example:prevforecast=1 starting-proj-date=1/1/1990 final-proj-date=1/1/2000 mobil_average=0\n, your line=%s . Probably you are running an older format.\n",num_filled,line); |
| |
fprintf(ficlog,"Error: Not 8 (data)parameters in line but %d, for example:prevforecast=1 starting-proj-date=1/1/1990 final-proj-date=1/1/2000 mov_average=0\n, your line=%s . Probably you are running an older format.\n",num_filled,line); |
| |
goto end; |
| |
} |
| |
fprintf(ficparo,"prevforecast=%d starting-proj-date=%.lf/%.lf/%.lf final-proj-date=%.lf/%.lf/%.lf mobil_average=%d\n",prevfcast,jproj1,mproj1,anproj1,jproj2,mproj2,anproj2,mobilavproj); |
| |
printf("prevforecast=%d starting-proj-date=%.lf/%.lf/%.lf final-proj-date=%.lf/%.lf/%.lf mobil_average=%d\n",prevfcast,jproj1,mproj1,anproj1,jproj2,mproj2,anproj2,mobilavproj); |
| |
fprintf(ficlog,"prevforecast=%d starting-proj-date=%.lf/%.lf/%.lf final-proj-date=%.lf/%.lf/%.lf mobil_average=%d\n",prevfcast,jproj1,mproj1,anproj1,jproj2,mproj2,anproj2,mobilavproj); |
| |
fprintf(ficres,"prevforecast=%d starting-proj-date=%.lf/%.lf/%.lf final-proj-date=%.lf/%.lf/%.lf mobil_average=%d\n",prevfcast,jproj1,mproj1,anproj1,jproj2,mproj2,anproj2,mobilavproj); |
| /* day and month of proj2 are not used but only year anproj2.*/ |
/* day and month of proj2 are not used but only year anproj2.*/ |
| |
|
| while((c=getc(ficpar))=='#' && c!= EOF){ |
|
| ungetc(c,ficpar); |
|
| fgets(line, MAXLINE, ficpar); |
|
| fputs(line,stdout); |
|
| fputs(line,ficparo); |
|
| fputs(line,ficres); |
|
| } |
} |
| ungetc(c,ficpar); |
while(fgets(line, MAXLINE, ficpar)) { |
| |
/* If line starts with a # it is a comment */ |
| |
if (line[0] == '#') { |
| |
numlinepar++; |
| |
fputs(line,stdout); |
| |
fputs(line,ficparo); |
| |
fputs(line,ficlog); |
| |
continue; |
| |
}else |
| |
break; |
| |
} |
| |
/* while((c=getc(ficpar))=='#' && c!= EOF){ */ |
| |
/* ungetc(c,ficpar); */ |
| |
/* fgets(line, MAXLINE, ficpar); */ |
| |
/* fputs(line,stdout); */ |
| |
/* fputs(line,ficparo); */ |
| |
/* fputs(line,ficres); */ |
| |
/* } */ |
| |
/* ungetc(c,ficpar); */ |
| |
|
| fscanf(ficpar,"backcast=%d starting-back-date=%lf/%lf/%lf final-back-date=%lf/%lf/%lf mobil_average=%d\n",&backcast,&jback1,&mback1,&anback1,&jback2,&mback2,&anback2,&mobilavproj); |
fscanf(ficpar,"backcast=%d starting-back-date=%lf/%lf/%lf final-back-date=%lf/%lf/%lf mobil_average=%d\n",&backcast,&jback1,&mback1,&anback1,&jback2,&mback2,&anback2,&mobilavproj); |
| fprintf(ficparo,"backcast=%d starting-back-date=%.lf/%.lf/%.lf final-back-date=%.lf/%.lf/%.lf mobil_average=%d\n",backcast,jback1,mback1,anback1,jback2,mback2,anback2,mobilavproj); |
if((num_filled=sscanf(line,"backcast=%d starting-back-date=%lf/%lf/%lf final-back-date=%lf/%lf/%lf mobil_average=%d\n",&backcast,&jback1,&mback1,&anback1,&jback2,&mback2,&anback2,&mobilavproj)) !=EOF){ |
| fprintf(ficlog,"backcast=%d starting-back-date=%.lf/%.lf/%.lf final-back-date=%.lf/%.lf/%.lf mobil_average=%d\n",backcast,jback1,mback1,anback1,jback2,mback2,anback2,mobilavproj); |
if (num_filled != 8) { |
| fprintf(ficres,"backcast=%d starting-back-date=%.lf/%.lf/%.lf final-back-date=%.lf/%.lf/%.lf mobil_average=%d\n",backcast,jback1,mback1,anback1,jback2,mback2,anback2,mobilavproj); |
printf("Error: Not 8 (data)parameters in line but %d, for example:backcast=1 starting-back-date=1/1/1990 finloal-back-date=1/1/1970 mobil_average=1\n, your line=%s . Probably you are running an older format.\n",num_filled,line); |
| |
fprintf(ficlog,"Error: Not 8 (data)parameters in line but %d, for example:backcast=1 starting-back-date=1/1/1990 finloal-back-date=1/1/1970 mobil_average=1\n, your line=%s . Probably you are running an older format.\n",num_filled,line); |
| |
goto end; |
| |
} |
| |
printf("backcast=%d starting-back-date=%.lf/%.lf/%.lf final-back-date=%.lf/%.lf/%.lf mobil_average=%d\n",backcast,jback1,mback1,anback1,jback2,mback2,anback2,mobilavproj); |
| |
fprintf(ficparo,"backcast=%d starting-back-date=%.lf/%.lf/%.lf final-back-date=%.lf/%.lf/%.lf mobil_average=%d\n",backcast,jback1,mback1,anback1,jback2,mback2,anback2,mobilavproj); |
| |
fprintf(ficlog,"backcast=%d starting-back-date=%.lf/%.lf/%.lf final-back-date=%.lf/%.lf/%.lf mobil_average=%d\n",backcast,jback1,mback1,anback1,jback2,mback2,anback2,mobilavproj); |
| |
fprintf(ficres,"backcast=%d starting-back-date=%.lf/%.lf/%.lf final-back-date=%.lf/%.lf/%.lf mobil_average=%d\n",backcast,jback1,mback1,anback1,jback2,mback2,anback2,mobilavproj); |
| /* day and month of proj2 are not used but only year anproj2.*/ |
/* day and month of proj2 are not used but only year anproj2.*/ |
| |
} |
| /* Results */ |
/* Results */ |
| nresult=0; |
nresult=0; |
| while(fgets(line, MAXLINE, ficpar)) { |
while(fgets(line, MAXLINE, ficpar)) { |
|
Line 10964 Please run with mle=-1 to get a correct
|
Line 11272 Please run with mle=-1 to get a correct
|
| while((num_filled=sscanf(line,"result:%[^\n]\n",resultline)) !=EOF){ |
while((num_filled=sscanf(line,"result:%[^\n]\n",resultline)) !=EOF){ |
| if (num_filled == 0){ |
if (num_filled == 0){ |
| resultline[0]='\0'; |
resultline[0]='\0'; |
| |
printf("Warning %d: no result line should be at minimum 'result: V2=0 V1=1 or result:.\n%s\n", num_filled, line); |
| break; |
break; |
| } else if (num_filled != 1){ |
} else if (num_filled != 1){ |
| printf("ERROR %d: result line should be at minimum 'result=' %s\n",num_filled, line); |
printf("ERROR %d: result line should be at minimum 'result: V2=0 V1=1 or result:.' %s\n",num_filled, line); |
| } |
} |
| nresult++; /* Sum of resultlines */ |
nresult++; /* Sum of resultlines */ |
| printf("Result %d: result=%s\n",nresult, resultline); |
printf("Result %d: result=%s\n",nresult, resultline); |
|
Line 11138 Please run with mle=-1 to get a correct
|
Line 11447 Please run with mle=-1 to get a correct
|
| |
|
| for(nres=1; nres <= nresult; nres++) /* For each resultline */ |
for(nres=1; nres <= nresult; nres++) /* For each resultline */ |
| for(k=1; k<=i1;k++){ /* For any combination of dummy covariates, fixed and varying */ |
for(k=1; k<=i1;k++){ /* For any combination of dummy covariates, fixed and varying */ |
| if(TKresult[nres]!= k) |
if(i1 != 1 && TKresult[nres]!= k) |
| continue; |
continue; |
| fprintf(ficreseij,"\n#****** "); |
fprintf(ficreseij,"\n#****** "); |
| printf("\n#****** "); |
printf("\n#****** "); |
|
Line 11211 Please run with mle=-1 to get a correct
|
Line 11520 Please run with mle=-1 to get a correct
|
| |
|
| for(nres=1; nres <= nresult; nres++) /* For each resultline */ |
for(nres=1; nres <= nresult; nres++) /* For each resultline */ |
| for(k=1; k<=i1;k++){ /* For any combination of dummy covariates, fixed and varying */ |
for(k=1; k<=i1;k++){ /* For any combination of dummy covariates, fixed and varying */ |
| if(TKresult[nres]!= k) |
if(i1 != 1 && TKresult[nres]!= k) |
| continue; |
continue; |
| printf("\n#****** Result for:"); |
printf("\n#****** Result for:"); |
| fprintf(ficrest,"\n#****** Result for:"); |
fprintf(ficrest,"\n#****** Result for:"); |
|
Line 11352 Please run with mle=-1 to get a correct
|
Line 11661 Please run with mle=-1 to get a correct
|
| |
|
| for(nres=1; nres <= nresult; nres++) /* For each resultline */ |
for(nres=1; nres <= nresult; nres++) /* For each resultline */ |
| for(k=1; k<=i1;k++){ |
for(k=1; k<=i1;k++){ |
| if(TKresult[nres]!= k) |
if(i1 != 1 && TKresult[nres]!= k) |
| continue; |
continue; |
| fprintf(ficresvpl,"\n#****** "); |
fprintf(ficresvpl,"\n#****** "); |
| printf("\n#****** "); |
printf("\n#****** "); |
![[BACK]](/icons/small/back.gif){kind=icon}
Return to imach.c CVS log ![[TXT]](/icons/small/text.gif){kind=icon}
![[DIR]](/icons/small/dir.gif){kind=icon}
Up to [Local Repository] / imach / src