|
version 1.288, 2018/05/02 20:58:27
|
version 1.307, 2021/03/08 18:11:32
|
|
Line 1
|
Line 1
|
| /* $Id$ |
/* $Id$ |
| $State$ |
$State$ |
| $Log$ |
$Log$ |
| |
Revision 1.307 2021/03/08 18:11:32 brouard |
| |
Summary: 0.99r22 fixed bug on result: |
| |
|
| |
Revision 1.306 2021/02/20 15:44:02 brouard |
| |
Summary: Version 0.99r21 |
| |
|
| |
* imach.c (Module): Fix bug on quitting after result lines! |
| |
(Module): Version 0.99r21 |
| |
|
| |
Revision 1.305 2021/02/20 15:28:30 brouard |
| |
* imach.c (Module): Fix bug on quitting after result lines! |
| |
|
| |
Revision 1.304 2021/02/12 11:34:20 brouard |
| |
* imach.c (Module): The use of a Windows BOM (huge) file is now an error |
| |
|
| |
Revision 1.303 2021/02/11 19:50:15 brouard |
| |
* (Module): imach.c Someone entered 'results:' instead of 'result:'. Now it is an error which is printed. |
| |
|
| |
Revision 1.302 2020/02/22 21:00:05 brouard |
| |
* (Module): imach.c Update mle=-3 (for computing Life expectancy |
| |
and life table from the data without any state) |
| |
|
| |
Revision 1.301 2019/06/04 13:51:20 brouard |
| |
Summary: Error in 'r'parameter file backcast yearsbproj instead of yearsfproj |
| |
|
| |
Revision 1.300 2019/05/22 19:09:45 brouard |
| |
Summary: version 0.99r19 of May 2019 |
| |
|
| |
Revision 1.299 2019/05/22 18:37:08 brouard |
| |
Summary: Cleaned 0.99r19 |
| |
|
| |
Revision 1.298 2019/05/22 18:19:56 brouard |
| |
*** empty log message *** |
| |
|
| |
Revision 1.297 2019/05/22 17:56:10 brouard |
| |
Summary: Fix bug by moving date2dmy and nhstepm which gaefin=-1 |
| |
|
| |
Revision 1.296 2019/05/20 13:03:18 brouard |
| |
Summary: Projection syntax simplified |
| |
|
| |
|
| |
We can now start projections, forward or backward, from the mean date |
| |
of inteviews up to or down to a number of years of projection: |
| |
prevforecast=1 yearsfproj=15.3 mobil_average=0 |
| |
or |
| |
prevforecast=1 starting-proj-date=1/1/2007 final-proj-date=12/31/2017 mobil_average=0 |
| |
or |
| |
prevbackcast=1 yearsbproj=12.3 mobil_average=1 |
| |
or |
| |
prevbackcast=1 starting-back-date=1/10/1999 final-back-date=1/1/1985 mobil_average=1 |
| |
|
| |
Revision 1.295 2019/05/18 09:52:50 brouard |
| |
Summary: doxygen tex bug |
| |
|
| |
Revision 1.294 2019/05/16 14:54:33 brouard |
| |
Summary: There was some wrong lines added |
| |
|
| |
Revision 1.293 2019/05/09 15:17:34 brouard |
| |
*** empty log message *** |
| |
|
| |
Revision 1.292 2019/05/09 14:17:20 brouard |
| |
Summary: Some updates |
| |
|
| |
Revision 1.291 2019/05/09 13:44:18 brouard |
| |
Summary: Before ncovmax |
| |
|
| |
Revision 1.290 2019/05/09 13:39:37 brouard |
| |
Summary: 0.99r18 unlimited number of individuals |
| |
|
| |
The number n which was limited to 20,000 cases is now unlimited, from firstobs to lastobs. If the number is too for the virtual memory, probably an error will occur. |
| |
|
| |
Revision 1.289 2018/12/13 09:16:26 brouard |
| |
Summary: Bug for young ages (<-30) will be in r17 |
| |
|
| Revision 1.288 2018/05/02 20:58:27 brouard |
Revision 1.288 2018/05/02 20:58:27 brouard |
| Summary: Some bugs fixed |
Summary: Some bugs fixed |
| |
|
|
Line 1046 typedef struct {
|
Line 1120 typedef struct {
|
| #define NINTERVMAX 8 |
#define NINTERVMAX 8 |
| #define NLSTATEMAX 8 /**< Maximum number of live states (for func) */ |
#define NLSTATEMAX 8 /**< Maximum number of live states (for func) */ |
| #define NDEATHMAX 8 /**< Maximum number of dead states (for func) */ |
#define NDEATHMAX 8 /**< Maximum number of dead states (for func) */ |
| #define NCOVMAX 20 /**< Maximum number of covariates, including generated covariates V1*V2 */ |
#define NCOVMAX 20 /**< Maximum number of covariates, including generated covariates V1*V2 */ |
| #define codtabm(h,k) (1 & (h-1) >> (k-1))+1 |
#define codtabm(h,k) (1 & (h-1) >> (k-1))+1 |
| /*#define decodtabm(h,k,cptcoveff)= (h <= (1<<cptcoveff)?(((h-1) >> (k-1)) & 1) +1 : -1)*/ |
/*#define decodtabm(h,k,cptcoveff)= (h <= (1<<cptcoveff)?(((h-1) >> (k-1)) & 1) +1 : -1)*/ |
| #define decodtabm(h,k,cptcoveff) (((h-1) >> (k-1)) & 1) +1 |
#define decodtabm(h,k,cptcoveff) (((h-1) >> (k-1)) & 1) +1 |
| #define MAXN 20000 |
/*#define MAXN 20000 */ /* Should by replaced by nobs, real number of observations and unlimited */ |
| #define YEARM 12. /**< Number of months per year */ |
#define YEARM 12. /**< Number of months per year */ |
| /* #define AGESUP 130 */ |
/* #define AGESUP 130 */ |
| /* #define AGESUP 150 */ |
/* #define AGESUP 150 */ |
|
Line 1074 typedef struct {
|
Line 1148 typedef struct {
|
| /* $State$ */ |
/* $State$ */ |
| #include "version.h" |
#include "version.h" |
| char version[]=__IMACH_VERSION__; |
char version[]=__IMACH_VERSION__; |
| char copyright[]="April 2018,INED-EUROREVES-Institut de longevite-Japan Society for the Promotion of Science (Grant-in-Aid for Scientific Research 25293121), Intel Software 2015-2018"; |
char copyright[]="Feb 2021,INED-EUROREVES-Institut de longevite-Japan Society for the Promotion of Science (Grant-in-Aid for Scientific Research 25293121), Intel Software 2015-2020, INED 2000-2021"; |
| char fullversion[]="$Revision$ $Date$"; |
char fullversion[]="$Revision$ $Date$"; |
| char strstart[80]; |
char strstart[80]; |
| char optionfilext[10], optionfilefiname[FILENAMELENGTH]; |
char optionfilext[10], optionfilefiname[FILENAMELENGTH]; |
|
Line 1098 int nqfveff=0; /**< nqfveff Number of Qu
|
Line 1172 int nqfveff=0; /**< nqfveff Number of Qu
|
| int ntveff=0; /**< ntveff number of effective time varying variables */ |
int ntveff=0; /**< ntveff number of effective time varying variables */ |
| int nqtveff=0; /**< ntqveff number of effective time varying quantitative variables */ |
int nqtveff=0; /**< ntqveff number of effective time varying quantitative variables */ |
| int cptcov=0; /* Working variable */ |
int cptcov=0; /* Working variable */ |
| |
int nobs=10; /* Number of observations in the data lastobs-firstobs */ |
| int ncovcombmax=NCOVMAX; /* Maximum calculated number of covariate combination = pow(2, cptcoveff) */ |
int ncovcombmax=NCOVMAX; /* Maximum calculated number of covariate combination = pow(2, cptcoveff) */ |
| int npar=NPARMAX; |
int npar=NPARMAX; /* Number of parameters (nlstate+ndeath-1)*nlstate*ncovmodel; */ |
| int nlstate=2; /* Number of live states */ |
int nlstate=2; /* Number of live states */ |
| int ndeath=1; /* Number of dead states */ |
int ndeath=1; /* Number of dead states */ |
| int ncovmodel=0, ncovcol=0; /* Total number of covariables including constant a12*1 +b12*x ncovmodel=2 */ |
int ncovmodel=0, ncovcol=0; /* Total number of covariables including constant a12*1 +b12*x ncovmodel=2 */ |
|
Line 1238 double **pmmij, ***probs; /* Global poin
|
Line 1313 double **pmmij, ***probs; /* Global poin
|
| double ***mobaverage, ***mobaverages; /* New global variable */ |
double ***mobaverage, ***mobaverages; /* New global variable */ |
| double *ageexmed,*agecens; |
double *ageexmed,*agecens; |
| double dateintmean=0; |
double dateintmean=0; |
| |
double anprojd, mprojd, jprojd; /* For eventual projections */ |
| |
double anprojf, mprojf, jprojf; |
| |
|
| |
double anbackd, mbackd, jbackd; /* For eventual backprojections */ |
| |
double anbackf, mbackf, jbackf; |
| |
double jintmean,mintmean,aintmean; |
| double *weight; |
double *weight; |
| int **s; /* Status */ |
int **s; /* Status */ |
| double *agedc; |
double *agedc; |
|
Line 2357 void powell(double p[], double **xi, int
|
Line 2437 void powell(double p[], double **xi, int
|
| /* But p and xit have been updated at the end of linmin, *fret corresponds to new p, xit */ |
/* But p and xit have been updated at the end of linmin, *fret corresponds to new p, xit */ |
| /* New value of last point Pn is not computed, P(n-1) */ |
/* New value of last point Pn is not computed, P(n-1) */ |
| for(j=1;j<=n;j++) { |
for(j=1;j<=n;j++) { |
| if(flatdir[j] >0){ |
if(flatdir[j] >0){ |
| printf(" p(%d)=%lf flat=%d ",j,p[j],flatdir[j]); |
printf(" p(%d)=%lf flat=%d ",j,p[j],flatdir[j]); |
| fprintf(ficlog," p(%d)=%lf flat=%d ",j,p[j],flatdir[j]); |
fprintf(ficlog," p(%d)=%lf flat=%d ",j,p[j],flatdir[j]); |
| } |
} |
| /* printf("\n"); */ |
/* printf("\n"); */ |
| /* fprintf(ficlog,"\n"); */ |
/* fprintf(ficlog,"\n"); */ |
| } |
} |
| /* if (2.0*fabs(fp-(*fret)) <= ftol*(fabs(fp)+fabs(*fret))) { /\* Did we reach enough precision? *\/ */ |
/* if (2.0*fabs(fp-(*fret)) <= ftol*(fabs(fp)+fabs(*fret))) { /\* Did we reach enough precision? *\/ */ |
| if (2.0*fabs(fp-(*fret)) <= ftol) { /* Did we reach enough precision? */ |
if (2.0*fabs(fp-(*fret)) <= ftol) { /* Did we reach enough precision? */ |
| /* We could compare with a chi^2. chisquare(0.95,ddl=1)=3.84 */ |
/* We could compare with a chi^2. chisquare(0.95,ddl=1)=3.84 */ |
|
Line 2951 double **pmij(double **ps, double *cov,
|
Line 3031 double **pmij(double **ps, double *cov,
|
| ps[ii][ii]=1; |
ps[ii][ii]=1; |
| } |
} |
| } |
} |
| |
|
| |
|
| /* for(ii=1; ii<= nlstate+ndeath; ii++){ */ |
/* for(ii=1; ii<= nlstate+ndeath; ii++){ */ |
| /* for(jj=1; jj<= nlstate+ndeath; jj++){ */ |
/* for(jj=1; jj<= nlstate+ndeath; jj++){ */ |
| /* printf(" pmij ps[%d][%d]=%lf ",ii,jj,ps[ii][jj]); */ |
/* printf(" pmij ps[%d][%d]=%lf ",ii,jj,ps[ii][jj]); */ |
|
Line 2972 double **pmij(double **ps, double *cov,
|
Line 3052 double **pmij(double **ps, double *cov,
|
| /* double **bmij(double **ps, double *cov, int ncovmodel, double *x, int nlstate, double ***prevacurrent, double ***dnewm, double **doldm, double **dsavm, int ij ) */ |
/* double **bmij(double **ps, double *cov, int ncovmodel, double *x, int nlstate, double ***prevacurrent, double ***dnewm, double **doldm, double **dsavm, int ij ) */ |
| double **bmij(double **ps, double *cov, int ncovmodel, double *x, int nlstate, double ***prevacurrent, int ij ) |
double **bmij(double **ps, double *cov, int ncovmodel, double *x, int nlstate, double ***prevacurrent, int ij ) |
| { |
{ |
| /* Computes the backward probability at age agefin and covariate combination ij. In fact cov is already filled and x too. |
/* Computes the backward probability at age agefin, cov[2], and covariate combination 'ij'. In fact cov is already filled and x too. |
| * Call to pmij(cov and x), call to cross prevalence, sums and inverses, left multiply, and returns in **ps as well as **bmij. |
* Call to pmij(cov and x), call to cross prevalence, sums and inverses, left multiply, and returns in **ps as well as **bmij. |
| */ |
*/ |
| int i, ii, j,k; |
int i, ii, j,k; |
|
Line 2980 double **pmij(double **ps, double *cov,
|
Line 3060 double **pmij(double **ps, double *cov,
|
| double **out, **pmij(); |
double **out, **pmij(); |
| double sumnew=0.; |
double sumnew=0.; |
| double agefin; |
double agefin; |
| double k3=0.; /* constant of the w_x diagonal matrixe (in order for B to sum to 1 even for death state) */ |
double k3=0.; /* constant of the w_x diagonal matrix (in order for B to sum to 1 even for death state) */ |
| double **dnewm, **dsavm, **doldm; |
double **dnewm, **dsavm, **doldm; |
| double **bbmij; |
double **bbmij; |
| |
|
|
Line 2999 double **pmij(double **ps, double *cov,
|
Line 3079 double **pmij(double **ps, double *cov,
|
| /* outputs pmmij which is a stochastic matrix in row */ |
/* outputs pmmij which is a stochastic matrix in row */ |
| |
|
| /* Diag(w_x) */ |
/* Diag(w_x) */ |
| /* Problem with prevacurrent which can be zero */ |
/* Rescaling the cross-sectional prevalence: Problem with prevacurrent which can be zero */ |
| sumnew=0.; |
sumnew=0.; |
| /*for (ii=1;ii<=nlstate+ndeath;ii++){*/ |
/*for (ii=1;ii<=nlstate+ndeath;ii++){*/ |
| for (ii=1;ii<=nlstate;ii++){ /* Only on live states */ |
for (ii=1;ii<=nlstate;ii++){ /* Only on live states */ |
| /* printf(" agefin=%d, ii=%d, ij=%d, prev=%f\n",(int)agefin,ii, ij, prevacurrent[(int)agefin][ii][ij]); */ |
/* printf(" agefin=%d, ii=%d, ij=%d, prev=%f\n",(int)agefin,ii, ij, prevacurrent[(int)agefin][ii][ij]); */ |
| sumnew+=prevacurrent[(int)agefin][ii][ij]; |
sumnew+=prevacurrent[(int)agefin][ii][ij]; |
| } |
} |
| if(sumnew >0.01){ /* At least some value in the prevalence */ |
if(sumnew >0.01){ /* At least some value in the prevalence */ |
|
Line 3026 double **pmij(double **ps, double *cov,
|
Line 3106 double **pmij(double **ps, double *cov,
|
| } |
} |
| /* End doldm, At the end doldm is diag[(w_i)] */ |
/* End doldm, At the end doldm is diag[(w_i)] */ |
| |
|
| /* left Product of this diag matrix by pmmij=Px (dnewm=dsavm*doldm) */ |
/* Left product of this diag matrix by pmmij=Px (dnewm=dsavm*doldm): diag[(w_i)*Px */ |
| bbmij=matprod2(dnewm, doldm,1,nlstate+ndeath,1,nlstate+ndeath,1,nlstate+ndeath, pmmij); /* Bug Valgrind */ |
bbmij=matprod2(dnewm, doldm,1,nlstate+ndeath,1,nlstate+ndeath,1,nlstate+ndeath, pmmij); /* was a Bug Valgrind */ |
| |
|
| /* Diag(Sum_i w^i_x p^ij_x */ |
/* Diag(Sum_i w^i_x p^ij_x, should be the prevalence at age x+stepm */ |
| /* w1 p11 + w2 p21 only on live states N1./N..*N11/N1. + N2./N..*N21/N2.=(N11+N21)/N..=N.1/N.. */ |
/* w1 p11 + w2 p21 only on live states N1./N..*N11/N1. + N2./N..*N21/N2.=(N11+N21)/N..=N.1/N.. */ |
| for (j=1;j<=nlstate+ndeath;j++){ |
for (j=1;j<=nlstate+ndeath;j++){ |
| sumnew=0.; |
sumnew=0.; |
|
Line 3047 double **pmij(double **ps, double *cov,
|
Line 3127 double **pmij(double **ps, double *cov,
|
| } /*End ii */ |
} /*End ii */ |
| } /* End j, At the end dsavm is diag[1/(w_1p1i+w_2 p2i)] for ALL states even if the sum is only for live states */ |
} /* End j, At the end dsavm is diag[1/(w_1p1i+w_2 p2i)] for ALL states even if the sum is only for live states */ |
| |
|
| ps=matprod2(ps, dnewm,1,nlstate+ndeath,1,nlstate+ndeath,1,nlstate+ndeath, dsavm); /* Bug Valgrind */ |
ps=matprod2(ps, dnewm,1,nlstate+ndeath,1,nlstate+ndeath,1,nlstate+ndeath, dsavm); /* was a Bug Valgrind */ |
| /* ps is now diag[w_i] * Px * diag [1/(w_1p1i+w_2 p2i)] */ |
/* ps is now diag[w_i] * Px * diag [1/(w_1p1i+w_2 p2i)] */ |
| /* end bmij */ |
/* end bmij */ |
| return ps; /*pointer is unchanged */ |
return ps; /*pointer is unchanged */ |
|
Line 3119 double **bpmij(double **ps, double *cov,
|
Line 3199 double **bpmij(double **ps, double *cov,
|
| ps[ii][ii]=1; |
ps[ii][ii]=1; |
| } |
} |
| } |
} |
| /* Added for backcast */ /* Transposed matrix too */ |
/* Added for prevbcast */ /* Transposed matrix too */ |
| for(jj=1; jj<= nlstate+ndeath; jj++){ |
for(jj=1; jj<= nlstate+ndeath; jj++){ |
| s1=0.; |
s1=0.; |
| for(ii=1; ii<= nlstate+ndeath; ii++){ |
for(ii=1; ii<= nlstate+ndeath; ii++){ |
|
Line 3879 return -l;
|
Line 3959 return -l;
|
| |
|
| |
|
| /*************** function likelione ***********/ |
/*************** function likelione ***********/ |
| void likelione(FILE *ficres,double p[], int npar, int nlstate, int *globpri, long *ipmx, double *sw, double *fretone, double (*funcone)(double [])) |
void likelione(FILE *ficres,double p[], int npar, int nlstate, int *globpri, long *ipmx, double *sw, double *fretone, double (*func)(double [])) |
| { |
{ |
| /* This routine should help understanding what is done with |
/* This routine should help understanding what is done with |
| the selection of individuals/waves and |
the selection of individuals/waves and |
|
Line 3903 void likelione(FILE *ficres,double p[],
|
Line 3983 void likelione(FILE *ficres,double p[],
|
| fprintf(ficresilk," -2*gipw/gsw*weight*ll(total)\n"); |
fprintf(ficresilk," -2*gipw/gsw*weight*ll(total)\n"); |
| } |
} |
| |
|
| *fretone=(*funcone)(p); |
*fretone=(*func)(p); |
| if(*globpri !=0){ |
if(*globpri !=0){ |
| fclose(ficresilk); |
fclose(ficresilk); |
| if (mle ==0) |
if (mle ==0) |
|
Line 4378 void pstamp(FILE *fichier)
|
Line 4458 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); |
| } |
} |
| |
|
| |
void date2dmy(double date,double *day, double *month, double *year){ |
| |
double yp=0., yp1=0., yp2=0.; |
| |
|
| |
yp1=modf(date,&yp);/* extracts integral of date in yp and |
| |
fractional in yp1 */ |
| |
*year=yp; |
| |
yp2=modf((yp1*12),&yp); |
| |
*month=yp; |
| |
yp1=modf((yp2*30.5),&yp); |
| |
*day=yp; |
| |
if(*day==0) *day=1; |
| |
if(*month==0) *month=1; |
| |
} |
| |
|
| |
|
| |
|
| /************ Frequencies ********************/ |
/************ Frequencies ********************/ |
|
Line 4952 Title=%s <br>Datafile=%s Firstpass=%d La
|
Line 5046 Title=%s <br>Datafile=%s Firstpass=%d La
|
| } |
} |
| } /* end mle=-2 */ |
} /* end mle=-2 */ |
| dateintmean=dateintsum/k2cpt; |
dateintmean=dateintsum/k2cpt; |
| |
date2dmy(dateintmean,&jintmean,&mintmean,&aintmean); |
| |
|
| fclose(ficresp); |
fclose(ficresp); |
| fclose(ficresphtm); |
fclose(ficresphtm); |
|
Line 5145 void prevalence(double ***probs, double
|
Line 5240 void prevalence(double ***probs, double
|
| |
|
| void concatwav(int wav[], int **dh, int **bh, int **mw, int **s, double *agedc, double **agev, int firstpass, int lastpass, int imx, int nlstate, int stepm) |
void concatwav(int wav[], int **dh, int **bh, int **mw, int **s, double *agedc, double **agev, int firstpass, int lastpass, int imx, int nlstate, int stepm) |
| { |
{ |
| /* Concatenates waves: wav[i] is the number of effective (useful waves) of individual i. |
/* Concatenates waves: wav[i] is the number of effective (useful waves in the sense that a non interview is useless) of individual i. |
| Death is a valid wave (if date is known). |
Death is a valid wave (if date is known). |
| mw[mi][i] is the mi (mi=1 to wav[i]) effective wave of individual i |
mw[mi][i] is the mi (mi=1 to wav[i]) effective wave of individual i |
| dh[m][i] or dh[mw[mi][i]][i] is the delay between two effective waves m=mw[mi][i] |
dh[m][i] or dh[mw[mi][i]][i] is the delay between two effective waves m=mw[mi][i] |
| and mw[mi+1][i]. dh depends on stepm. |
and mw[mi+1][i]. dh depends on stepm. s[m][i] exists for any wave from firstpass to lastpass |
| */ |
*/ |
| |
|
| int i=0, mi=0, m=0, mli=0; |
int i=0, mi=0, m=0, mli=0; |
|
Line 5187 void concatwav(int wav[], int **dh, int
|
Line 5282 void concatwav(int wav[], int **dh, int
|
| #else |
#else |
| if(s[m][i]==-1 && (int) andc[i] == 9999 && (int)anint[m][i] != 9999){ |
if(s[m][i]==-1 && (int) andc[i] == 9999 && (int)anint[m][i] != 9999){ |
| if(firsthree == 0){ |
if(firsthree == 0){ |
| printf("Information! Unknown status for individual %ld line=%d occurred at last wave %d at known date %d/%d. Please, check if your unknown date of death %d/%d means a live state %d at wave %d. This case(%d)/wave(%d) contributes to the likelihood as 1-p%d%d .\nOthers in log file only\n",num[i],i,lastpass,(int)mint[m][i],(int)anint[m][i], (int) moisdc[i], (int) andc[i], s[m][i], m, i, m, s[m][i], nlstate+ndeath); |
printf("Information! Unknown status for individual %ld line=%d occurred at last wave %d at known date %d/%d. Please, check if your unknown date of death %d/%d means a live state %d at wave %d. This case(%d)/wave(%d) contributes to the likelihood as 1-p_{%d%d} .\nOthers in log file only\n",num[i],i,lastpass,(int)mint[m][i],(int)anint[m][i], (int) moisdc[i], (int) andc[i], s[m][i], m, i, m, s[m][i], nlstate+ndeath); |
| firsthree=1; |
firsthree=1; |
| } |
} |
| fprintf(ficlog,"Information! Unknown status for individual %ld line=%d occurred at last wave %d at known date %d/%d. Please, check if your unknown date of death %d/%d means a live state %d at wave %d. This case(%d)/wave(%d) contributes to the likelihood as 1-p%d%d .\n",num[i],i,lastpass,(int)mint[m][i],(int)anint[m][i], (int) moisdc[i], (int) andc[i], s[m][i], m, i, m, s[m][i], nlstate+ndeath); |
fprintf(ficlog,"Information! Unknown status for individual %ld line=%d occurred at last wave %d at known date %d/%d. Please, check if your unknown date of death %d/%d means a live state %d at wave %d. This case(%d)/wave(%d) contributes to the likelihood as 1-p_{%d%d} .\n",num[i],i,lastpass,(int)mint[m][i],(int)anint[m][i], (int) moisdc[i], (int) andc[i], s[m][i], m, i, m, s[m][i], nlstate+ndeath); |
| mw[++mi][i]=m; |
mw[++mi][i]=m; |
| mli=m; |
mli=m; |
| } |
} |
|
Line 5994 void concatwav(int wav[], int **dh, int
|
Line 6089 void concatwav(int wav[], int **dh, int
|
| prlim[i][i]=mobaverage[(int)age][i][ij]; |
prlim[i][i]=mobaverage[(int)age][i][ij]; |
| } |
} |
| } |
} |
| /**< Computes the shifted transition matrix \f$ {}{h}_p^{ij}_x\f$ at horizon h. |
/**< Computes the shifted transition matrix \f$ {}{h}_p^{ij}x\f$ at horizon h. |
| */ |
*/ |
| hpxij(p3mat,nhstepm,age,hstepm,xp,nlstate,stepm,oldm,savm, ij,nres); /* Returns p3mat[i][j][h] for h=0 to nhstepm */ |
hpxij(p3mat,nhstepm,age,hstepm,xp,nlstate,stepm,oldm,savm, ij,nres); /* Returns p3mat[i][j][h] for h=0 to nhstepm */ |
| /**< And for each alive state j, sums over i \f$ w^i_x {}{h}_p^{ij}_x\f$, which are the probability |
/**< And for each alive state j, sums over i \f$ w^i_x {}{h}_p^{ij}x\f$, which are the probability |
| * at horizon h in state j including mortality. |
* at horizon h in state j including mortality. |
| */ |
*/ |
| for(j=1; j<= nlstate; j++){ |
for(j=1; j<= nlstate; j++){ |
|
Line 6787 To be simple, these graphs help to under
|
Line 6882 To be simple, these graphs help to under
|
| void printinghtml(char fileresu[], char title[], char datafile[], int firstpass, \ |
void printinghtml(char fileresu[], char title[], char datafile[], int firstpass, \ |
| int lastpass, int stepm, int weightopt, char model[],\ |
int lastpass, int stepm, int weightopt, char model[],\ |
| int imx,int jmin, int jmax, double jmeanint,char rfileres[],\ |
int imx,int jmin, int jmax, double jmeanint,char rfileres[],\ |
| int popforecast, int mobilav, int prevfcast, int mobilavproj, int backcast, int estepm , \ |
int popforecast, int mobilav, int prevfcast, int mobilavproj, int prevbcast, int estepm , \ |
| double jprev1, double mprev1,double anprev1, double dateprev1, double dateproj1, double dateback1, \ |
double jprev1, double mprev1,double anprev1, double dateprev1, double dateprojd, double dateback1, \ |
| double jprev2, double mprev2,double anprev2, double dateprev2, double dateproj2, double dateback2){ |
double jprev2, double mprev2,double anprev2, double dateprev2, double dateprojf, double dateback2){ |
| int jj1, k1, i1, cpt, k4, nres; |
int jj1, k1, i1, cpt, k4, nres; |
| |
|
| fprintf(fichtm,"<ul><li><a href='#firstorder'>Result files (first order: no variance)</a>\n \ |
fprintf(fichtm,"<ul><li><a href='#firstorder'>Result files (first order: no variance)</a>\n \ |
|
Line 6919 divided by h: <sub>h</sub>P<sub>ij</sub>
|
Line 7014 divided by h: <sub>h</sub>P<sub>ij</sub>
|
| <img src=\"%s_%d-3-%d.svg\">",stepm,subdirf2(optionfilefiname,"PE_"),k1,nres,subdirf2(optionfilefiname,"PE_"),k1,nres,subdirf2(optionfilefiname,"PE_"),k1,nres); |
<img src=\"%s_%d-3-%d.svg\">",stepm,subdirf2(optionfilefiname,"PE_"),k1,nres,subdirf2(optionfilefiname,"PE_"),k1,nres,subdirf2(optionfilefiname,"PE_"),k1,nres); |
| /* Survival functions (period) in state j */ |
/* Survival functions (period) in state j */ |
| for(cpt=1; cpt<=nlstate;cpt++){ |
for(cpt=1; cpt<=nlstate;cpt++){ |
| fprintf(fichtm,"<br>\n- Survival functions in state %d. Or probability to survive 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- Survival functions in state %d. And probability to be observed 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> \ |
| <img src=\"%s_%d-%d-%d.svg\">", cpt, cpt, nlstate, subdirf2(optionfilefiname,"LIJ_"),cpt,k1,nres,subdirf2(optionfilefiname,"LIJ_"),cpt,k1,nres,subdirf2(optionfilefiname,"LIJ_"),cpt,k1,nres); |
<img src=\"%s_%d-%d-%d.svg\">", cpt, cpt, nlstate, subdirf2(optionfilefiname,"LIJ_"),cpt,k1,nres,subdirf2(optionfilefiname,"LIJ_"),cpt,k1,nres,subdirf2(optionfilefiname,"LIJ_"),cpt,k1,nres); |
| } |
} |
| /* State specific survival functions (period) */ |
/* State specific survival functions (period) */ |
| for(cpt=1; cpt<=nlstate;cpt++){ |
for(cpt=1; cpt<=nlstate;cpt++){ |
| fprintf(fichtm,"<br>\n- Survival functions from state %d in each live state and total.\ |
fprintf(fichtm,"<br>\n- Survival functions in state %d and in any other live state (total).\ |
| Or probability to survive in various states (1 to %d) being in state %d at different ages. \ |
And probability to be observed in various states (up to %d) being in state %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, nlstate, cpt, subdirf2(optionfilefiname,"LIJT_"),cpt,k1,nres,subdirf2(optionfilefiname,"LIJT_"),cpt,k1,nres,subdirf2(optionfilefiname,"LIJT_"),cpt,k1,nres); |
<a href=\"%s_%d-%d-%d.svg\">%s_%d-%d-%d.svg</a><br> <img src=\"%s_%d-%d-%d.svg\">", cpt, nlstate, cpt, subdirf2(optionfilefiname,"LIJT_"),cpt,k1,nres,subdirf2(optionfilefiname,"LIJT_"),cpt,k1,nres,subdirf2(optionfilefiname,"LIJT_"),cpt,k1,nres); |
| } |
} |
| /* Period (forward stable) prevalence in each health state */ |
/* Period (forward stable) prevalence in each health state */ |
|
Line 6933 divided by h: <sub>h</sub>P<sub>ij</sub>
|
Line 7028 divided by h: <sub>h</sub>P<sub>ij</sub>
|
| fprintf(fichtm,"<br>\n- Convergence to period (stable) prevalence in state %d. Or probability for a person being in state (1 to %d) at different ages, to be in state %d some years after. <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 for a person being in state (1 to %d) at different ages, to be in state %d some years after. <a href=\"%s_%d-%d-%d.svg\">%s_%d-%d-%d.svg</a><br> \ |
| <img src=\"%s_%d-%d-%d.svg\">", cpt, nlstate, cpt, 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, nlstate, cpt, subdirf2(optionfilefiname,"P_"),cpt,k1,nres,subdirf2(optionfilefiname,"P_"),cpt,k1,nres,subdirf2(optionfilefiname,"P_"),cpt,k1,nres); |
| } |
} |
| if(backcast==1){ |
if(prevbcast==1){ |
| /* Backward prevalence in each health state */ |
/* Backward prevalence in each health state */ |
| for(cpt=1; cpt<=nlstate;cpt++){ |
for(cpt=1; cpt<=nlstate;cpt++){ |
| fprintf(fichtm,"<br>\n- Convergence to mixed (stable) back prevalence in state %d. Or probability for a person to be in state %d at a younger age, knowing that she/he was in state (1 to %d) at different older 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 for a person to be in state %d at a younger age, knowing that she/he was in state (1 to %d) at different older ages. <a href=\"%s_%d-%d-%d.svg\">%s_%d-%d-%d.svg</a><br> \ |
|
Line 6944 divided by h: <sub>h</sub>P<sub>ij</sub>
|
Line 7039 divided by h: <sub>h</sub>P<sub>ij</sub>
|
| /* Projection of prevalence up to period (forward stable) prevalence in each health state */ |
/* Projection of prevalence up to period (forward stable) prevalence in each health state */ |
| for(cpt=1; cpt<=nlstate;cpt++){ |
for(cpt=1; cpt<=nlstate;cpt++){ |
| fprintf(fichtm,"<br>\n- Projection of cross-sectional prevalence (estimated with cases observed from %.1f to %.1f and mobil_average=%d), from year %.1f up to year %.1f tending to period (stable) forward prevalence in state %d. Or probability to be in state %d being in an observed weighted state (from 1 to %d). <a href=\"%s_%d-%d-%d.svg\">%s_%d-%d-%d.svg</a><br> \ |
fprintf(fichtm,"<br>\n- Projection of cross-sectional prevalence (estimated with cases observed from %.1f to %.1f and mobil_average=%d), from year %.1f up to year %.1f tending to period (stable) forward prevalence in state %d. Or probability to be in state %d being in an observed weighted state (from 1 to %d). <a href=\"%s_%d-%d-%d.svg\">%s_%d-%d-%d.svg</a><br> \ |
| <img src=\"%s_%d-%d-%d.svg\">", dateprev1, dateprev2, mobilavproj, dateproj1, dateproj2, cpt, cpt, nlstate, subdirf2(optionfilefiname,"PROJ_"),cpt,k1,nres,subdirf2(optionfilefiname,"PROJ_"),cpt,k1,nres,subdirf2(optionfilefiname,"PROJ_"),cpt,k1,nres); |
<img src=\"%s_%d-%d-%d.svg\">", dateprev1, dateprev2, mobilavproj, dateprojd, dateprojf, cpt, cpt, nlstate, subdirf2(optionfilefiname,"PROJ_"),cpt,k1,nres,subdirf2(optionfilefiname,"PROJ_"),cpt,k1,nres,subdirf2(optionfilefiname,"PROJ_"),cpt,k1,nres); |
| } |
} |
| } |
} |
| if(backcast==1){ |
if(prevbcast==1){ |
| /* Back projection of prevalence up to stable (mixed) back-prevalence in each health state */ |
/* Back projection of prevalence up to stable (mixed) back-prevalence in each health state */ |
| for(cpt=1; cpt<=nlstate;cpt++){ |
for(cpt=1; cpt<=nlstate;cpt++){ |
| fprintf(fichtm,"<br>\n- Back projection of cross-sectional prevalence (estimated with cases observed from %.1f to %.1f and mobil_average=%d), \ |
fprintf(fichtm,"<br>\n- Back projection of cross-sectional prevalence (estimated with cases observed from %.1f to %.1f and mobil_average=%d), \ |
|
Line 7060 true period expectancies (those weighted
|
Line 7155 true period expectancies (those weighted
|
| } |
} |
| |
|
| /******************* Gnuplot file **************/ |
/******************* Gnuplot file **************/ |
| void printinggnuplot(char fileresu[], char optionfilefiname[], double ageminpar, double agemaxpar, double bage, double fage , int prevfcast, int backcast, char pathc[], double p[], int offyear, int offbyear){ |
void printinggnuplot(char fileresu[], char optionfilefiname[], double ageminpar, double agemaxpar, double bage, double fage , int prevfcast, int prevbcast, char pathc[], double p[], int offyear, int offbyear){ |
| |
|
| char dirfileres[132],optfileres[132]; |
char dirfileres[132],optfileres[132]; |
| char gplotcondition[132], gplotlabel[132]; |
char gplotcondition[132], gplotlabel[132]; |
|
Line 7214 void printinggnuplot(char fileresu[], ch
|
Line 7309 void printinggnuplot(char fileresu[], ch
|
| } /* end covariate */ |
} /* end covariate */ |
| } /* end if no covariate */ |
} /* end if no covariate */ |
| |
|
| if(backcast==1){ /* We need to get the corresponding values of the covariates involved in this combination k1 */ |
if(prevbcast==1){ /* We need to get the corresponding values of the covariates involved in this combination k1 */ |
| /* fprintf(ficgp,",\"%s\" every :::%d::%d u 1:($%d) t\"Backward stable prevalence\" w l lt 3",subdirf2(fileresu,"PLB_"),k1-1,k1-1,1+cpt); */ |
/* fprintf(ficgp,",\"%s\" every :::%d::%d u 1:($%d) t\"Backward stable prevalence\" w l lt 3",subdirf2(fileresu,"PLB_"),k1-1,k1-1,1+cpt); */ |
| fprintf(ficgp,",\"%s\" u 1:((",subdirf2(fileresu,"PLB_")); /* Age is in 1, nres in 2 to be fixed */ |
fprintf(ficgp,",\"%s\" u 1:((",subdirf2(fileresu,"PLB_")); /* Age is in 1, nres in 2 to be fixed */ |
| if(cptcoveff ==0){ |
if(cptcoveff ==0){ |
|
Line 7241 void printinggnuplot(char fileresu[], ch
|
Line 7336 void printinggnuplot(char fileresu[], ch
|
| } |
} |
| } /* end covariate */ |
} /* end covariate */ |
| } /* end if no covariate */ |
} /* end if no covariate */ |
| if(backcast == 1){ |
if(prevbcast == 1){ |
| fprintf(ficgp,", \"%s\" every :::%d::%d u 1:($2==%d ? $3:1/0) \"%%lf %%lf",subdirf2(fileresu,"VBL_"),nres-1,nres-1,nres); |
fprintf(ficgp,", \"%s\" every :::%d::%d u 1:($2==%d ? $3:1/0) \"%%lf %%lf",subdirf2(fileresu,"VBL_"),nres-1,nres-1,nres); |
| /* k1-1 error should be nres-1*/ |
/* k1-1 error should be nres-1*/ |
| for (i=1; i<= nlstate ; i ++) { |
for (i=1; i<= nlstate ; i ++) { |
|
Line 7260 void printinggnuplot(char fileresu[], ch
|
Line 7355 void printinggnuplot(char fileresu[], ch
|
| } |
} |
| fprintf(ficgp,"\" t\"\" w l lt 4"); |
fprintf(ficgp,"\" t\"\" w l lt 4"); |
| } /* end if backprojcast */ |
} /* end if backprojcast */ |
| } /* end if backcast */ |
} /* end if prevbcast */ |
| /* fprintf(ficgp,"\nset out ;unset label;\n"); */ |
/* fprintf(ficgp,"\nset out ;unset label;\n"); */ |
| fprintf(ficgp,"\nset out ;unset title;\n"); |
fprintf(ficgp,"\nset out ;unset title;\n"); |
| } /* nres */ |
} /* nres */ |
|
Line 7548 set ter svg size 640, 480\nunset log y\n
|
Line 7643 set ter svg size 640, 480\nunset log y\n
|
| |
|
| |
|
| /* 7eme */ |
/* 7eme */ |
| if(backcast == 1){ |
if(prevbcast == 1){ |
| /* CV backward prevalence for each covariate */ |
/* CV backward prevalence 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 */ |
|
Line 7600 set ter svg size 640, 480\nunset log y\n
|
Line 7695 set ter svg size 640, 480\nunset log y\n
|
| fprintf(ficgp,"\nset out; unset label;\n"); |
fprintf(ficgp,"\nset out; unset label;\n"); |
| } /* end cpt state*/ |
} /* end cpt state*/ |
| } /* end covariate */ |
} /* end covariate */ |
| } /* End if backcast */ |
} /* End if prevbcast */ |
| |
|
| /* 8eme */ |
/* 8eme */ |
| if(prevfcast==1){ |
if(prevfcast==1){ |
|
Line 7716 set ter svg size 640, 480\nunset log y\n
|
Line 7811 set ter svg size 640, 480\nunset log y\n
|
| } /* end covariate */ |
} /* end covariate */ |
| } /* End if prevfcast */ |
} /* End if prevfcast */ |
| |
|
| if(backcast==1){ |
if(prevbcast==1){ |
| /* Back projection from cross-sectional to stable (mixed) for each covariate */ |
/* Back projection from cross-sectional to stable (mixed) 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 */ |
|
Line 7829 set ter svg size 640, 480\nunset log y\n
|
Line 7924 set ter svg size 640, 480\nunset log y\n
|
| fprintf(ficgp,"\nset out; unset label;\n"); |
fprintf(ficgp,"\nset out; unset label;\n"); |
| } /* end cpt state*/ |
} /* end cpt state*/ |
| } /* end covariate */ |
} /* end covariate */ |
| } /* End if backcast */ |
} /* End if prevbcast */ |
| |
|
| |
|
| /* 9eme writing MLE parameters */ |
/* 9eme writing MLE parameters */ |
|
Line 8277 set ter svg size 640, 480\nunset log y\n
|
Line 8372 set ter svg size 640, 480\nunset log y\n
|
| }/* End movingaverage */ |
}/* End movingaverage */ |
| |
|
| |
|
| |
|
| /************** Forecasting ******************/ |
/************** Forecasting ******************/ |
| void prevforecast(char fileres[], double anproj1, double mproj1, double jproj1, double ageminpar, double agemax, double dateprev1, double dateprev2, int mobilav, double ***prev, double bage, double fage, int firstpass, int lastpass, double anproj2, double p[], int cptcoveff){ |
/* void prevforecast(char fileres[], double dateintmean, double anprojd, double mprojd, double jprojd, double ageminpar, double agemax, double dateprev1, double dateprev2, int mobilav, double ***prev, double bage, double fage, int firstpass, int lastpass, double anprojf, double p[], int cptcoveff)*/ |
| /* proj1, year, month, day of starting projection |
void prevforecast(char fileres[], double dateintmean, double dateprojd, double dateprojf, double ageminpar, double agemax, double dateprev1, double dateprev2, int mobilav, double ***prev, double bage, double fage, int firstpass, int lastpass, double p[], int cptcoveff){ |
| |
/* dateintemean, mean date of interviews |
| |
dateprojd, year, month, day of starting projection |
| |
dateprojf date of end of projection;year of end of projection (same day and month as proj1). |
| agemin, agemax range of age |
agemin, agemax range of age |
| dateprev1 dateprev2 range of dates during which prevalence is computed |
dateprev1 dateprev2 range of dates during which prevalence is computed |
| anproj2 year of en of projection (same day and month as proj1). |
|
| */ |
*/ |
| |
/* double anprojd, mprojd, jprojd; */ |
| |
/* double anprojf, mprojf, jprojf; */ |
| int yearp, stepsize, hstepm, nhstepm, j, k, cptcod, i, h, i1, k4, nres=0; |
int yearp, stepsize, hstepm, nhstepm, j, k, cptcod, i, h, i1, k4, nres=0; |
| double agec; /* generic age */ |
double agec; /* generic age */ |
| double agelim, ppij, yp,yp1,yp2,jprojmean,mprojmean,anprojmean; |
double agelim, ppij, yp,yp1,yp2; |
| double *popeffectif,*popcount; |
double *popeffectif,*popcount; |
| double ***p3mat; |
double ***p3mat; |
| /* double ***mobaverage; */ |
/* double ***mobaverage; */ |
|
Line 8323 set ter svg size 640, 480\nunset log y\n
|
Line 8423 set ter svg size 640, 480\nunset log y\n
|
| if(estepm > stepm){ /* Yes every two year */ |
if(estepm > stepm){ /* Yes every two year */ |
| stepsize=2; |
stepsize=2; |
| } |
} |
| |
hstepm=hstepm/stepm; |
| |
|
| hstepm=hstepm/stepm; |
|
| yp1=modf(dateintmean,&yp);/* extracts integral of datemean in yp and |
/* yp1=modf(dateintmean,&yp);/\* extracts integral of datemean in yp and */ |
| fractional in yp1 */ |
/* fractional in yp1 *\/ */ |
| anprojmean=yp; |
/* aintmean=yp; */ |
| yp2=modf((yp1*12),&yp); |
/* yp2=modf((yp1*12),&yp); */ |
| mprojmean=yp; |
/* mintmean=yp; */ |
| yp1=modf((yp2*30.5),&yp); |
/* yp1=modf((yp2*30.5),&yp); */ |
| jprojmean=yp; |
/* jintmean=yp; */ |
| if(jprojmean==0) jprojmean=1; |
/* if(jintmean==0) jintmean=1; */ |
| if(mprojmean==0) jprojmean=1; |
/* if(mintmean==0) mintmean=1; */ |
| |
|
| |
|
| |
/* date2dmy(dateintmean,&jintmean,&mintmean,&aintmean); */ |
| |
/* date2dmy(dateprojd,&jprojd, &mprojd, &anprojd); */ |
| |
/* date2dmy(dateprojf,&jprojf, &mprojf, &anprojf); */ |
| i1=pow(2,cptcoveff); |
i1=pow(2,cptcoveff); |
| if (cptcovn < 1){i1=1;} |
if (cptcovn < 1){i1=1;} |
| |
|
| fprintf(ficresf,"# Mean day of interviews %.lf/%.lf/%.lf (%.2f) between %.2f and %.2f \n",jprojmean,mprojmean,anprojmean,dateintmean,dateprev1,dateprev2); |
fprintf(ficresf,"# Mean day of interviews %.lf/%.lf/%.lf (%.2f) between %.2f and %.2f \n",jintmean,mintmean,aintmean,dateintmean,dateprev1,dateprev2); |
| |
|
| fprintf(ficresf,"#****** Routine prevforecast **\n"); |
fprintf(ficresf,"#****** Routine prevforecast **\n"); |
| |
|
|
Line 8364 set ter svg size 640, 480\nunset log y\n
|
Line 8469 set ter svg size 640, 480\nunset log y\n
|
| fprintf(ficresf," p%d%d",i,j); |
fprintf(ficresf," p%d%d",i,j); |
| fprintf(ficresf," wp.%d",j); |
fprintf(ficresf," wp.%d",j); |
| } |
} |
| for (yearp=0; yearp<=(anproj2-anproj1);yearp +=stepsize) { |
for (yearp=0; yearp<=(anprojf-anprojd);yearp +=stepsize) { |
| fprintf(ficresf,"\n"); |
fprintf(ficresf,"\n"); |
| fprintf(ficresf,"\n# Forecasting at date %.lf/%.lf/%.lf ",jproj1,mproj1,anproj1+yearp); |
fprintf(ficresf,"\n# Forecasting at date %.lf/%.lf/%.lf ",jprojd,mprojd,anprojd+yearp); |
| /* for (agec=fage; agec>=(ageminpar-1); agec--){ */ |
/* for (agec=fage; agec>=(ageminpar-1); agec--){ */ |
| for (agec=fage; agec>=(bage); agec--){ |
for (agec=fage; agec>=(bage); agec--){ |
| nhstepm=(int) rint((agelim-agec)*YEARM/stepm); |
nhstepm=(int) rint((agelim-agec)*YEARM/stepm); |
|
Line 8384 set ter svg size 640, 480\nunset log y\n
|
Line 8489 set ter svg size 640, 480\nunset log y\n
|
| fprintf(ficresf,"\n"); |
fprintf(ficresf,"\n"); |
| for(j=1;j<=cptcoveff;j++) |
for(j=1;j<=cptcoveff;j++) |
| fprintf(ficresf,"%d %d ",Tvaraff[j],nbcode[Tvaraff[j]][codtabm(k,j)]); |
fprintf(ficresf,"%d %d ",Tvaraff[j],nbcode[Tvaraff[j]][codtabm(k,j)]); |
| fprintf(ficresf,"%.f %.f ",anproj1+yearp,agec+h*hstepm/YEARM*stepm); |
fprintf(ficresf,"%.f %.f ",anprojd+yearp,agec+h*hstepm/YEARM*stepm); |
| |
|
| for(j=1; j<=nlstate+ndeath;j++) { |
for(j=1; j<=nlstate+ndeath;j++) { |
| ppij=0.; |
ppij=0.; |
|
Line 8412 set ter svg size 640, 480\nunset log y\n
|
Line 8517 set ter svg size 640, 480\nunset log y\n
|
| } |
} |
| |
|
| /************** Back Forecasting ******************/ |
/************** Back Forecasting ******************/ |
| void prevbackforecast(char fileres[], double ***prevacurrent, double anback1, double mback1, double jback1, double ageminpar, double agemax, double dateprev1, double dateprev2, int mobilav, double bage, double fage, int firstpass, int lastpass, double anback2, double p[], int cptcoveff){ |
/* void prevbackforecast(char fileres[], double ***prevacurrent, double anback1, double mback1, double jback1, double ageminpar, double agemax, double dateprev1, double dateprev2, int mobilav, double bage, double fage, int firstpass, int lastpass, double anback2, double p[], int cptcoveff){ */ |
| /* back1, year, month, day of starting backection |
void prevbackforecast(char fileres[], double ***prevacurrent, double dateintmean, double dateprojd, double dateprojf, double ageminpar, double agemax, double dateprev1, double dateprev2, int mobilav, double bage, double fage, int firstpass, int lastpass, double p[], int cptcoveff){ |
| |
/* back1, year, month, day of starting backprojection |
| agemin, agemax range of age |
agemin, agemax range of age |
| dateprev1 dateprev2 range of dates during which prevalence is computed |
dateprev1 dateprev2 range of dates during which prevalence is computed |
| anback2 year of end of backprojection (same day and month as back1). |
anback2 year of end of backprojection (same day and month as back1). |
|
Line 8421 set ter svg size 640, 480\nunset log y\n
|
Line 8527 set ter svg size 640, 480\nunset log y\n
|
| */ |
*/ |
| int yearp, stepsize, hstepm, nhstepm, j, k, cptcod, i, h, i1, k4, nres=0; |
int yearp, stepsize, hstepm, nhstepm, j, k, cptcod, i, h, i1, k4, nres=0; |
| double agec; /* generic age */ |
double agec; /* generic age */ |
| double agelim, ppij, ppi, yp,yp1,yp2,jprojmean,mprojmean,anprojmean; |
double agelim, ppij, ppi, yp,yp1,yp2; /* ,jintmean,mintmean,aintmean;*/ |
| double *popeffectif,*popcount; |
double *popeffectif,*popcount; |
| double ***p3mat; |
double ***p3mat; |
| /* double ***mobaverage; */ |
/* double ***mobaverage; */ |
|
Line 8464 set ter svg size 640, 480\nunset log y\n
|
Line 8570 set ter svg size 640, 480\nunset log y\n
|
| } |
} |
| |
|
| hstepm=hstepm/stepm; |
hstepm=hstepm/stepm; |
| yp1=modf(dateintmean,&yp);/* extracts integral of datemean in yp and |
/* yp1=modf(dateintmean,&yp);/\* extracts integral of datemean in yp and */ |
| fractional in yp1 */ |
/* fractional in yp1 *\/ */ |
| anprojmean=yp; |
/* aintmean=yp; */ |
| yp2=modf((yp1*12),&yp); |
/* yp2=modf((yp1*12),&yp); */ |
| mprojmean=yp; |
/* mintmean=yp; */ |
| yp1=modf((yp2*30.5),&yp); |
/* yp1=modf((yp2*30.5),&yp); */ |
| jprojmean=yp; |
/* jintmean=yp; */ |
| if(jprojmean==0) jprojmean=1; |
/* if(jintmean==0) jintmean=1; */ |
| if(mprojmean==0) jprojmean=1; |
/* if(mintmean==0) jintmean=1; */ |
| |
|
| i1=pow(2,cptcoveff); |
i1=pow(2,cptcoveff); |
| if (cptcovn < 1){i1=1;} |
if (cptcovn < 1){i1=1;} |
| |
|
| fprintf(ficresfb,"# Mean day of interviews %.lf/%.lf/%.lf (%.2f) between %.2f and %.2f \n",jprojmean,mprojmean,anprojmean,dateintmean,dateprev1,dateprev2); |
fprintf(ficresfb,"# Mean day of interviews %.lf/%.lf/%.lf (%.2f) between %.2f and %.2f \n",jintmean,mintmean,aintmean,dateintmean,dateprev1,dateprev2); |
| printf("# Mean day of interviews %.lf/%.lf/%.lf (%.2f) between %.2f and %.2f \n",jprojmean,mprojmean,anprojmean,dateintmean,dateprev1,dateprev2); |
printf("# Mean day of interviews %.lf/%.lf/%.lf (%.2f) between %.2f and %.2f \n",jintmean,mintmean,aintmean,dateintmean,dateprev1,dateprev2); |
| |
|
| fprintf(ficresfb,"#****** Routine prevbackforecast **\n"); |
fprintf(ficresfb,"#****** Routine prevbackforecast **\n"); |
| |
|
|
Line 8503 set ter svg size 640, 480\nunset log y\n
|
Line 8609 set ter svg size 640, 480\nunset log y\n
|
| fprintf(ficresfb," b%d%d",i,j); |
fprintf(ficresfb," b%d%d",i,j); |
| fprintf(ficresfb," b.%d",j); |
fprintf(ficresfb," b.%d",j); |
| } |
} |
| for (yearp=0; yearp>=(anback2-anback1);yearp -=stepsize) { |
for (yearp=0; yearp>=(anbackf-anbackd);yearp -=stepsize) { |
| /* for (yearp=0; yearp<=(anproj2-anproj1);yearp +=stepsize) { */ |
/* for (yearp=0; yearp<=(anproj2-anproj1);yearp +=stepsize) { */ |
| fprintf(ficresfb,"\n"); |
fprintf(ficresfb,"\n"); |
| fprintf(ficresfb,"\n# Back Forecasting at date %.lf/%.lf/%.lf ",jback1,mback1,anback1+yearp); |
fprintf(ficresfb,"\n# Back Forecasting at date %.lf/%.lf/%.lf ",jbackd,mbackd,anbackd+yearp); |
| /* printf("\n# Back Forecasting at date %.lf/%.lf/%.lf ",jback1,mback1,anback1+yearp); */ |
/* printf("\n# Back Forecasting at date %.lf/%.lf/%.lf ",jback1,mback1,anback1+yearp); */ |
| /* for (agec=bage; agec<=agemax-1; agec++){ /\* testing *\/ */ |
/* for (agec=bage; agec<=agemax-1; agec++){ /\* testing *\/ */ |
| for (agec=bage; agec<=fage; agec++){ /* testing */ |
for (agec=bage; agec<=fage; agec++){ /* testing */ |
|
Line 8529 set ter svg size 640, 480\nunset log y\n
|
Line 8635 set ter svg size 640, 480\nunset log y\n
|
| fprintf(ficresfb,"\n"); |
fprintf(ficresfb,"\n"); |
| for(j=1;j<=cptcoveff;j++) |
for(j=1;j<=cptcoveff;j++) |
| fprintf(ficresfb,"%d %d ",Tvaraff[j],nbcode[Tvaraff[j]][codtabm(k,j)]); |
fprintf(ficresfb,"%d %d ",Tvaraff[j],nbcode[Tvaraff[j]][codtabm(k,j)]); |
| fprintf(ficresfb,"%.f %.f ",anback1+yearp,agec-h*hstepm/YEARM*stepm); |
fprintf(ficresfb,"%.f %.f ",anbackd+yearp,agec-h*hstepm/YEARM*stepm); |
| for(i=1; i<=nlstate+ndeath;i++) { |
for(i=1; i<=nlstate+ndeath;i++) { |
| ppij=0.;ppi=0.; |
ppij=0.;ppi=0.; |
| for(j=1; j<=nlstate;j++) { |
for(j=1; j<=nlstate;j++) { |
|
Line 8977 void prwizard(int ncovmodel, int nlstate
|
Line 9083 void prwizard(int ncovmodel, int nlstate
|
| /******************* Gompertz Likelihood ******************************/ |
/******************* Gompertz Likelihood ******************************/ |
| double gompertz(double x[]) |
double gompertz(double x[]) |
| { |
{ |
| double A,B,L=0.0,sump=0.,num=0.; |
double A=0.0,B=0.,L=0.0,sump=0.,num=0.; |
| int i,n=0; /* n is the size of the sample */ |
int i,n=0; /* n is the size of the sample */ |
| |
|
| for (i=1;i<=imx ; i++) { |
for (i=1;i<=imx ; i++) { |
|
Line 8985 double gompertz(double x[])
|
Line 9091 double gompertz(double x[])
|
| /* sump=sump+1;*/ |
/* sump=sump+1;*/ |
| num=num+1; |
num=num+1; |
| } |
} |
| |
L=0.0; |
| |
/* agegomp=AGEGOMP; */ |
| /* for (i=0; i<=imx; i++) |
/* for (i=0; i<=imx; i++) |
| if (wav[i]>0) printf("i=%d ageex=%lf agecens=%lf agedc=%lf cens=%d %d\n" ,i,ageexmed[i],agecens[i],agedc[i],cens[i],wav[i]);*/ |
if (wav[i]>0) printf("i=%d ageex=%lf agecens=%lf agedc=%lf cens=%d %d\n" ,i,ageexmed[i],agecens[i],agedc[i],cens[i],wav[i]);*/ |
| |
|
| for (i=1;i<=imx ; i++) |
for (i=1;i<=imx ; i++) { |
| { |
/* mu(a)=mu(agecomp)*exp(teta*(age-agegomp)) |
| if (cens[i] == 1 && wav[i]>1) |
mu(a)=x[1]*exp(x[2]*(age-agegomp)); x[1] and x[2] are per year. |
| A=-x[1]/(x[2])*(exp(x[2]*(agecens[i]-agegomp))-exp(x[2]*(ageexmed[i]-agegomp))); |
* L= Product mu(agedeces)exp(-\int_ageexam^agedc mu(u) du ) for a death between agedc (in month) |
| |
* and agedc +1 month, cens[i]=0: log(x[1]/YEARM) |
| if (cens[i] == 0 && wav[i]>1) |
* + |
| |
* exp(-\int_ageexam^agecens mu(u) du ) when censored, cens[i]=1 |
| |
*/ |
| |
if (wav[i] > 1 || agedc[i] < AGESUP) { |
| |
if (cens[i] == 1){ |
| |
A=-x[1]/(x[2])*(exp(x[2]*(agecens[i]-agegomp))-exp(x[2]*(ageexmed[i]-agegomp))); |
| |
} else if (cens[i] == 0){ |
| A=-x[1]/(x[2])*(exp(x[2]*(agedc[i]-agegomp))-exp(x[2]*(ageexmed[i]-agegomp))) |
A=-x[1]/(x[2])*(exp(x[2]*(agedc[i]-agegomp))-exp(x[2]*(ageexmed[i]-agegomp))) |
| +log(x[1]/YEARM)+x[2]*(agedc[i]-agegomp)+log(YEARM); |
+log(x[1]/YEARM) +x[2]*(agedc[i]-agegomp)+log(YEARM); |
| |
} else |
| |
printf("Gompertz cens[%d] neither 1 nor 0\n",i); |
| /*if (wav[i] > 1 && agecens[i] > 15) {*/ /* ??? */ |
/*if (wav[i] > 1 && agecens[i] > 15) {*/ /* ??? */ |
| if (wav[i] > 1 ) { /* ??? */ |
L=L+A*weight[i]; |
| L=L+A*weight[i]; |
|
| /* printf("\ni=%d A=%f L=%lf x[1]=%lf x[2]=%lf ageex=%lf agecens=%lf cens=%d agedc=%lf weight=%lf\n",i,A,L,x[1],x[2],ageexmed[i]*12,agecens[i]*12,cens[i],agedc[i]*12,weight[i]);*/ |
/* printf("\ni=%d A=%f L=%lf x[1]=%lf x[2]=%lf ageex=%lf agecens=%lf cens=%d agedc=%lf weight=%lf\n",i,A,L,x[1],x[2],ageexmed[i]*12,agecens[i]*12,cens[i],agedc[i]*12,weight[i]);*/ |
| } |
} |
| } |
} |
| |
|
| /*printf("x1=%2.9f x2=%2.9f x3=%2.9f L=%f\n",x[1],x[2],x[3],L);*/ |
/*printf("x1=%2.9f x2=%2.9f x3=%2.9f L=%f\n",x[1],x[2],x[3],L);*/ |
| |
|
| return -2*L*num/sump; |
return -2*L*num/sump; |
| } |
} |
|
Line 9015 double gompertz(double x[])
|
Line 9127 double gompertz(double x[])
|
| /******************* Gompertz_f Likelihood ******************************/ |
/******************* Gompertz_f Likelihood ******************************/ |
| double gompertz_f(const gsl_vector *v, void *params) |
double gompertz_f(const gsl_vector *v, void *params) |
| { |
{ |
| double A,B,LL=0.0,sump=0.,num=0.; |
double A=0.,B=0.,LL=0.0,sump=0.,num=0.; |
| double *x= (double *) v->data; |
double *x= (double *) v->data; |
| int i,n=0; /* n is the size of the sample */ |
int i,n=0; /* n is the size of the sample */ |
| |
|
|
Line 9108 int readdata(char datafile[], int firsto
|
Line 9220 int readdata(char datafile[], int firsto
|
| int i=0, j=0, n=0, iv=0, v; |
int i=0, j=0, n=0, iv=0, v; |
| int lstra; |
int lstra; |
| int linei, month, year,iout; |
int linei, month, year,iout; |
| |
int noffset=0; /* This is the offset if BOM data file */ |
| char line[MAXLINE], linetmp[MAXLINE]; |
char line[MAXLINE], linetmp[MAXLINE]; |
| char stra[MAXLINE], strb[MAXLINE]; |
char stra[MAXLINE], strb[MAXLINE]; |
| char *stratrunc; |
char *stratrunc; |
|
Line 9141 int readdata(char datafile[], int firsto
|
Line 9254 int readdata(char datafile[], int firsto
|
| fprintf(ficlog,"Problem while opening datafile: %s with errno='%s'\n", datafile,strerror(errno));fflush(ficlog);return 1; |
fprintf(ficlog,"Problem while opening datafile: %s with errno='%s'\n", datafile,strerror(errno));fflush(ficlog);return 1; |
| } |
} |
| |
|
| i=1; |
/* Is it a BOM UTF-8 Windows file? */ |
| |
/* First data line */ |
| linei=0; |
linei=0; |
| |
while(fgets(line, MAXLINE, fic)) { |
| |
noffset=0; |
| |
if( line[0] == (char)0xEF && line[1] == (char)0xBB) /* EF BB BF */ |
| |
{ |
| |
noffset=noffset+3; |
| |
printf("# Data file '%s' is an UTF8 BOM file, please convert to UTF8 or ascii file and rerun.\n",datafile);fflush(stdout); |
| |
fprintf(ficlog,"# Data file '%s' is an UTF8 BOM file, please convert to UTF8 or ascii file and rerun.\n",datafile); |
| |
fflush(ficlog); return 1; |
| |
} |
| |
/* else if( line[0] == (char)0xFE && line[1] == (char)0xFF)*/ |
| |
else if( line[0] == (char)0xFF && line[1] == (char)0xFE) |
| |
{ |
| |
noffset=noffset+2; |
| |
printf("# Error Data file '%s' is a huge UTF16BE BOM file, please convert to UTF8 or ascii file (for example with dos2unix) and rerun.\n",datafile);fflush(stdout); |
| |
fprintf(ficlog,"# Error Data file '%s' is a huge UTF16BE BOM file, please convert to UTF8 or ascii file (for example with dos2unix) and rerun.\n",datafile); |
| |
fflush(ficlog); return 1; |
| |
} |
| |
else if( line[0] == 0 && line[1] == 0) |
| |
{ |
| |
if( line[2] == (char)0xFE && line[3] == (char)0xFF){ |
| |
noffset=noffset+4; |
| |
printf("# Error Data file '%s' is a huge UTF16BE BOM file, please convert to UTF8 or ascii file (for example with dos2unix) and rerun.\n",datafile);fflush(stdout); |
| |
fprintf(ficlog,"# Error Data file '%s' is a huge UTF16BE BOM file, please convert to UTF8 or ascii file (for example with dos2unix) and rerun.\n",datafile); |
| |
fflush(ficlog); return 1; |
| |
} |
| |
} else{ |
| |
;/*printf(" Not a BOM file\n");*/ |
| |
} |
| |
/* If line starts with a # it is a comment */ |
| |
if (line[noffset] == '#') { |
| |
linei=linei+1; |
| |
break; |
| |
}else{ |
| |
break; |
| |
} |
| |
} |
| |
fclose(fic); |
| |
if((fic=fopen(datafile,"r"))==NULL) { |
| |
printf("Problem while opening datafile: %s with errno='%s'\n", datafile,strerror(errno));fflush(stdout); |
| |
fprintf(ficlog,"Problem while opening datafile: %s with errno='%s'\n", datafile,strerror(errno));fflush(ficlog);return 1; |
| |
} |
| |
/* Not a Bom file */ |
| |
|
| |
i=1; |
| while ((fgets(line, MAXLINE, fic) != NULL) &&((i >= firstobs) && (i <=lastobs))) { |
while ((fgets(line, MAXLINE, fic) != NULL) &&((i >= firstobs) && (i <=lastobs))) { |
| linei=linei+1; |
linei=linei+1; |
| for(j=strlen(line); j>=0;j--){ /* Untabifies line */ |
for(j=strlen(line); j>=0;j--){ /* Untabifies line */ |
|
Line 9263 int readdata(char datafile[], int firsto
|
Line 9421 int readdata(char datafile[], int firsto
|
| return 1; |
return 1; |
| } |
} |
| anint[j][i]= (double) year; |
anint[j][i]= (double) year; |
| mint[j][i]= (double)month; |
mint[j][i]= (double)month; |
| |
/* if( (int)anint[j][i]+ (int)(mint[j][i])/12. < (int) (moisnais[i]/12.+annais[i])){ */ |
| |
/* printf("Warning reading data around '%s' at line number %d for individual %d, '%s'\nThe date of interview (%2d/%4d) at wave %d occurred before the date of birth (%2d/%4d).\n",strb, linei,i, line, mint[j][i],anint[j][i], moisnais[i],annais[i]); */ |
| |
/* fprintf(ficlog,"Warning reading data around '%s' at line number %d for individual %d, '%s'\nThe date of interview (%2d/%4d) at wave %d occurred before the date of birth (%2d/%4d).\n",strb, linei,i, line, mint[j][i],anint[j][i], moisnais[i],annais[i]); */ |
| |
/* } */ |
| strcpy(line,stra); |
strcpy(line,stra); |
| } /* End loop on waves */ |
} /* End loop on waves */ |
| |
|
|
Line 9302 int readdata(char datafile[], int firsto
|
Line 9464 int readdata(char datafile[], int firsto
|
| |
|
| } |
} |
| annais[i]=(double)(year); |
annais[i]=(double)(year); |
| moisnais[i]=(double)(month); |
moisnais[i]=(double)(month); |
| |
for (j=1;j<=maxwav;j++){ |
| |
if( (int)anint[j][i]+ (int)(mint[j][i])/12. < (int) (moisnais[i]/12.+annais[i])){ |
| |
printf("Warning reading data around '%s' at line number %d for individual %d, '%s'\nThe date of interview (%2d/%4d) at wave %d occurred before the date of birth (%2d/%4d).\n",strb, linei,i, line, (int)mint[j][i],(int)anint[j][i], j,(int)moisnais[i],(int)annais[i]); |
| |
fprintf(ficlog,"Warning reading data around '%s' at line number %d for individual %d, '%s'\nThe date of interview (%2d/%4d) at wave %d occurred before the date of birth (%2d/%4d).\n",strb, linei,i, line, (int)mint[j][i],(int)anint[j][i], j, (int)moisnais[i],(int)annais[i]); |
| |
} |
| |
} |
| |
|
| strcpy(line,stra); |
strcpy(line,stra); |
| |
|
| /* Sample weight */ |
/* Sample weight */ |
|
Line 9425 int decoderesult ( char resultline[], in
|
Line 9594 int decoderesult ( char resultline[], in
|
| char stra[80], strb[80], strc[80], strd[80],stre[80]; |
char stra[80], strb[80], strc[80], strd[80],stre[80]; |
| |
|
| removefirstspace(&resultline); |
removefirstspace(&resultline); |
| printf("decoderesult:%s\n",resultline); |
|
| |
|
| if (strstr(resultline,"v") !=0){ |
if (strstr(resultline,"v") !=0){ |
| printf("Error. 'v' must be in upper case 'V' result: %s ",resultline); |
printf("Error. 'v' must be in upper case 'V' result: %s ",resultline); |
|
Line 9440 int decoderesult ( char resultline[], in
|
Line 9608 int decoderesult ( char resultline[], in
|
| TKresult[nres]=0; /* Combination for the nresult and the model */ |
TKresult[nres]=0; /* Combination for the nresult and the model */ |
| return (0); |
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 10187 BOOL IsWow64()
|
Line 10354 BOOL IsWow64()
|
| #endif |
#endif |
| |
|
| void syscompilerinfo(int logged) |
void syscompilerinfo(int logged) |
| { |
{ |
| /* #include "syscompilerinfo.h"*/ |
#include <stdint.h> |
| |
|
| |
/* #include "syscompilerinfo.h"*/ |
| /* command line Intel compiler 32bit windows, XP compatible:*/ |
/* command line Intel compiler 32bit windows, XP compatible:*/ |
| /* /GS /W3 /Gy |
/* /GS /W3 /Gy |
| /Zc:wchar_t /Zi /O2 /Fd"Release\vc120.pdb" /D "WIN32" /D "NDEBUG" /D |
/Zc:wchar_t /Zi /O2 /Fd"Release\vc120.pdb" /D "WIN32" /D "NDEBUG" /D |
|
Line 10223 void syscompilerinfo(int logged)
|
Line 10392 void syscompilerinfo(int logged)
|
| /ManifestFile:"Release\IMaCh.exe.intermediate.manifest" /OPT:ICF |
/ManifestFile:"Release\IMaCh.exe.intermediate.manifest" /OPT:ICF |
| /NOLOGO /TLBID:1 |
/NOLOGO /TLBID:1 |
| */ |
*/ |
| |
|
| |
|
| #if defined __INTEL_COMPILER |
#if defined __INTEL_COMPILER |
| #if defined(__GNUC__) |
#if defined(__GNUC__) |
| struct utsname sysInfo; /* For Intel on Linux and OS/X */ |
struct utsname sysInfo; /* For Intel on Linux and OS/X */ |
|
Line 10239 void syscompilerinfo(int logged)
|
Line 10410 void syscompilerinfo(int logged)
|
| } |
} |
| #endif |
#endif |
| |
|
| #include <stdint.h> |
|
| |
|
| printf("Compiled with:");if(logged)fprintf(ficlog,"Compiled with:"); |
printf("Compiled with:");if(logged)fprintf(ficlog,"Compiled with:"); |
| #if defined(__clang__) |
#if defined(__clang__) |
| printf(" Clang/LLVM");if(logged)fprintf(ficlog," Clang/LLVM"); /* Clang/LLVM. ---------------------------------------------- */ |
printf(" Clang/LLVM");if(logged)fprintf(ficlog," Clang/LLVM"); /* Clang/LLVM. ---------------------------------------------- */ |
|
Line 10647 int hPijx(double *p, int bage, int fage)
|
Line 10816 int hPijx(double *p, int bage, int fage)
|
| /*if (stepm<=24) stepsize=2;*/ |
/*if (stepm<=24) stepsize=2;*/ |
| |
|
| /* agelim=AGESUP; */ |
/* agelim=AGESUP; */ |
| ageminl=30; |
ageminl=AGEINF; /* was 30 */ |
| hstepm=stepsize*YEARM; /* Every year of age */ |
hstepm=stepsize*YEARM; /* Every year of age */ |
| hstepm=hstepm/stepm; /* Typically 2 years, = 2/6 months = 4 */ |
hstepm=hstepm/stepm; /* Typically 2 years, = 2/6 months = 4 */ |
| |
|
|
Line 10677 int hPijx(double *p, int bage, int fage)
|
Line 10846 int hPijx(double *p, int bage, int fage)
|
| /* for (agedeb=fage; agedeb>=bage; agedeb--){ /\* If stepm=6 months *\/ */ |
/* for (agedeb=fage; agedeb>=bage; agedeb--){ /\* If stepm=6 months *\/ */ |
| for (agedeb=bage; agedeb<=fage; agedeb++){ /* If stepm=6 months and estepm=24 (2 years) */ |
for (agedeb=bage; agedeb<=fage; agedeb++){ /* If stepm=6 months and estepm=24 (2 years) */ |
| /* nhstepm=(int) rint((agelim-agedeb)*YEARM/stepm); /\* Typically 20 years = 20*12/6=40 *\/ */ |
/* nhstepm=(int) rint((agelim-agedeb)*YEARM/stepm); /\* Typically 20 years = 20*12/6=40 *\/ */ |
| nhstepm=(int) rint((agedeb-ageminl)*YEARM/stepm); /* Typically 20 years = 20*12/6=40 */ |
nhstepm=(int) rint((agedeb-ageminl)*YEARM/stepm+0.1)-1; /* Typically 20 years = 20*12/6=40 or 55*12/24=27.5-1.1=>27 */ |
| nhstepm = nhstepm/hstepm; /* Typically 40/4=10, because estepm=24 stepm=6 => hstepm=24/6=4 */ |
nhstepm = nhstepm/hstepm; /* Typically 40/4=10, because estepm=24 stepm=6 => hstepm=24/6=4 or 28*/ |
| |
|
| /* nhstepm=nhstepm*YEARM; aff par mois*/ |
/* nhstepm=nhstepm*YEARM; aff par mois*/ |
| |
|
|
Line 10726 int main(int argc, char *argv[])
|
Line 10895 int main(int argc, char *argv[])
|
| double ssval; |
double ssval; |
| #endif |
#endif |
| int movingaverage(double ***probs, double bage,double fage, double ***mobaverage, int mobilav); |
int movingaverage(double ***probs, double bage,double fage, double ***mobaverage, int mobilav); |
| int i,j, k, n=MAXN,iter=0,m,size=100, cptcod; |
int i,j, k, iter=0,m,size=100, cptcod; /* Suppressing because nobs */ |
| |
/* int i,j, k, n=MAXN,iter=0,m,size=100, cptcod; */ |
| int ncvyear=0; /* Number of years needed for the period prevalence to converge */ |
int ncvyear=0; /* Number of years needed for the period prevalence to converge */ |
| int jj, ll, li, lj, lk; |
int jj, ll, li, lj, lk; |
| int numlinepar=0; /* Current linenumber of parameter file */ |
int numlinepar=0; /* Current linenumber of parameter file */ |
|
Line 10761 int main(int argc, char *argv[])
|
Line 10931 int main(int argc, char *argv[])
|
| |
|
| char pathr[MAXLINE], pathimach[MAXLINE]; |
char pathr[MAXLINE], pathimach[MAXLINE]; |
| char *tok, *val; /* pathtot */ |
char *tok, *val; /* pathtot */ |
| int firstobs=1, lastobs=10; |
int firstobs=1, lastobs=10; /* nobs = lastobs-firstobs declared globally ;*/ |
| int c, h , cpt, c2; |
int c, h , cpt, c2; |
| int jl=0; |
int jl=0; |
| int i1, j1, jk, stepsize=0; |
int i1, j1, jk, stepsize=0; |
|
Line 10769 int main(int argc, char *argv[])
|
Line 10939 int main(int argc, char *argv[])
|
| |
|
| int *tab; |
int *tab; |
| int mobilavproj=0 , prevfcast=0 ; /* moving average of prev, If prevfcast=1 prevalence projection */ |
int mobilavproj=0 , prevfcast=0 ; /* moving average of prev, If prevfcast=1 prevalence projection */ |
| int backcast=0; |
/* double anprojd, mprojd, jprojd; /\* For eventual projections *\/ */ |
| |
/* double anprojf, mprojf, jprojf; */ |
| |
/* double jintmean,mintmean,aintmean; */ |
| |
int prvforecast = 0; /* Might be 1 (date of beginning of projection is a choice or 2 is the dateintmean */ |
| |
int prvbackcast = 0; /* Might be 1 (date of beginning of projection is a choice or 2 is the dateintmean */ |
| |
double yrfproj= 10.0; /* Number of years of forward projections */ |
| |
double yrbproj= 10.0; /* Number of years of backward projections */ |
| |
int prevbcast=0; /* defined as global for mlikeli and mle, replacing backcast */ |
| int mobilav=0,popforecast=0; |
int mobilav=0,popforecast=0; |
| int hstepm=0, nhstepm=0; |
int hstepm=0, nhstepm=0; |
| int agemortsup; |
int agemortsup; |
|
Line 10794 int main(int argc, char *argv[])
|
Line 10971 int main(int argc, char *argv[])
|
| double *epj, vepp; |
double *epj, vepp; |
| |
|
| double dateprev1, dateprev2; |
double dateprev1, dateprev2; |
| double jproj1=1,mproj1=1,anproj1=2000,jproj2=1,mproj2=1,anproj2=2000, dateproj1=0, dateproj2=0; |
double jproj1=1,mproj1=1,anproj1=2000,jproj2=1,mproj2=1,anproj2=2000, dateproj1=0, dateproj2=0, dateprojd=0, dateprojf=0; |
| double jback1=1,mback1=1,anback1=2000,jback2=1,mback2=1,anback2=2000, dateback1=0, dateback2=0; |
double jback1=1,mback1=1,anback1=2000,jback2=1,mback2=1,anback2=2000, dateback1=0, dateback2=0, datebackd=0, datebackf=0; |
| |
|
| |
|
| double **ximort; |
double **ximort; |
| char *alph[]={"a","a","b","c","d","e"}, str[4]="1234"; |
char *alph[]={"a","a","b","c","d","e"}, str[4]="1234"; |
|
Line 10990 int main(int argc, char *argv[])
|
Line 11168 int main(int argc, char *argv[])
|
| noffset=noffset+3; |
noffset=noffset+3; |
| printf("# File is an UTF8 Bom.\n"); // 0xBF |
printf("# File is an UTF8 Bom.\n"); // 0xBF |
| } |
} |
| else if( line[0] == (char)0xFE && line[1] == (char)0xFF) |
/* else if( line[0] == (char)0xFE && line[1] == (char)0xFF)*/ |
| |
else if( line[0] == (char)0xFF && line[1] == (char)0xFE) |
| { |
{ |
| noffset=noffset+2; |
noffset=noffset+2; |
| printf("# File is an UTF16BE BOM file\n"); |
printf("# File is an UTF16BE BOM file\n"); |
|
Line 11078 int main(int argc, char *argv[])
|
Line 11257 int main(int argc, char *argv[])
|
| } |
} |
| if((num_filled=sscanf(line,"model=1+age%[^.\n]", model)) !=EOF){ |
if((num_filled=sscanf(line,"model=1+age%[^.\n]", model)) !=EOF){ |
| if (num_filled != 1){ |
if (num_filled != 1){ |
| printf("ERROR %d: Model should be at minimum 'model=1+age' %s\n",num_filled, line); |
printf("ERROR %d: Model should be at minimum 'model=1+age+' instead of '%s'\n",num_filled, line); |
| fprintf(ficlog,"ERROR %d: Model should be at minimum 'model=1+age' %s\n",num_filled, line); |
fprintf(ficlog,"ERROR %d: Model should be at minimum 'model=1+age+' instead of '%s'\n",num_filled, line); |
| model[0]='\0'; |
model[0]='\0'; |
| goto end; |
goto end; |
| } |
} |
|
Line 11128 int main(int argc, char *argv[])
|
Line 11307 int main(int argc, char *argv[])
|
| ungetc(c,ficpar); |
ungetc(c,ficpar); |
| |
|
| |
|
| covar=matrix(0,NCOVMAX,1,n); /**< used in readdata */ |
covar=matrix(0,NCOVMAX,firstobs,lastobs); /**< used in readdata */ |
| if(nqv>=1)coqvar=matrix(1,nqv,1,n); /**< Fixed quantitative covariate */ |
if(nqv>=1)coqvar=matrix(1,nqv,firstobs,lastobs); /**< Fixed quantitative covariate */ |
| if(nqtv>=1)cotqvar=ma3x(1,maxwav,1,nqtv,1,n); /**< Time varying quantitative covariate */ |
if(nqtv>=1)cotqvar=ma3x(1,maxwav,1,nqtv,firstobs,lastobs); /**< Time varying quantitative covariate */ |
| if(ntv+nqtv>=1)cotvar=ma3x(1,maxwav,1,ntv+nqtv,1,n); /**< Time varying covariate (dummy and quantitative)*/ |
if(ntv+nqtv>=1)cotvar=ma3x(1,maxwav,1,ntv+nqtv,firstobs,lastobs); /**< Time varying covariate (dummy and quantitative)*/ |
| cptcovn=0; /*Number of covariates, i.e. number of '+' in model statement plus one, indepently of n in Vn*/ |
cptcovn=0; /*Number of covariates, i.e. number of '+' in model statement plus one, indepently of n in Vn*/ |
| /* v1+v2+v3+v2*v4+v5*age makes cptcovn = 5 |
/* v1+v2+v3+v2*v4+v5*age makes cptcovn = 5 |
| v1+v2*age+v2*v3 makes cptcovn = 3 |
v1+v2*age+v2*v3 makes cptcovn = 3 |
|
Line 11194 int main(int argc, char *argv[])
|
Line 11373 int main(int argc, char *argv[])
|
| for(jj=1; jj <=nlstate+ndeath; jj++){ |
for(jj=1; jj <=nlstate+ndeath; jj++){ |
| if(jj==i) continue; |
if(jj==i) continue; |
| j++; |
j++; |
| |
while((c=getc(ficpar))=='#' && c!= EOF){ |
| |
ungetc(c,ficpar); |
| |
fgets(line, MAXLINE, ficpar); |
| |
numlinepar++; |
| |
fputs(line,stdout); |
| |
fputs(line,ficparo); |
| |
fputs(line,ficlog); |
| |
} |
| |
ungetc(c,ficpar); |
| fscanf(ficpar,"%1d%1d",&i1,&j1); |
fscanf(ficpar,"%1d%1d",&i1,&j1); |
| if ((i1 != i) || (j1 != jj)){ |
if ((i1 != i) || (j1 != jj)){ |
| printf("Error in line parameters number %d, %1d%1d instead of %1d%1d \n \ |
printf("Error in line parameters number %d, %1d%1d instead of %1d%1d \n \ |
|
Line 11334 Please run with mle=-1 to get a correct
|
Line 11522 Please run with mle=-1 to get a correct
|
| |
|
| /* Main data |
/* Main data |
| */ |
*/ |
| n= lastobs; |
nobs=lastobs-firstobs+1; /* was = lastobs;*/ |
| num=lvector(1,n); |
/* num=lvector(1,n); */ |
| moisnais=vector(1,n); |
/* moisnais=vector(1,n); */ |
| annais=vector(1,n); |
/* annais=vector(1,n); */ |
| moisdc=vector(1,n); |
/* moisdc=vector(1,n); */ |
| andc=vector(1,n); |
/* andc=vector(1,n); */ |
| weight=vector(1,n); |
/* weight=vector(1,n); */ |
| agedc=vector(1,n); |
/* agedc=vector(1,n); */ |
| cod=ivector(1,n); |
/* cod=ivector(1,n); */ |
| for(i=1;i<=n;i++){ |
/* for(i=1;i<=n;i++){ */ |
| |
num=lvector(firstobs,lastobs); |
| |
moisnais=vector(firstobs,lastobs); |
| |
annais=vector(firstobs,lastobs); |
| |
moisdc=vector(firstobs,lastobs); |
| |
andc=vector(firstobs,lastobs); |
| |
weight=vector(firstobs,lastobs); |
| |
agedc=vector(firstobs,lastobs); |
| |
cod=ivector(firstobs,lastobs); |
| |
for(i=firstobs;i<=lastobs;i++){ |
| num[i]=0; |
num[i]=0; |
| moisnais[i]=0; |
moisnais[i]=0; |
| annais[i]=0; |
annais[i]=0; |
|
Line 11353 Please run with mle=-1 to get a correct
|
Line 11550 Please run with mle=-1 to get a correct
|
| cod[i]=0; |
cod[i]=0; |
| weight[i]=1.0; /* Equal weights, 1 by default */ |
weight[i]=1.0; /* Equal weights, 1 by default */ |
| } |
} |
| mint=matrix(1,maxwav,1,n); |
mint=matrix(1,maxwav,firstobs,lastobs); |
| anint=matrix(1,maxwav,1,n); |
anint=matrix(1,maxwav,firstobs,lastobs); |
| s=imatrix(1,maxwav+1,1,n); /* s[i][j] health state for wave i and individual j */ |
s=imatrix(1,maxwav+1,firstobs,lastobs); /* s[i][j] health state for wave i and individual j */ |
| tab=ivector(1,NCOVMAX); |
tab=ivector(1,NCOVMAX); |
| ncodemax=ivector(1,NCOVMAX); /* Number of code per covariate; if O and 1 only, 2**ncov; V1+V2+V3+V4=>16 */ |
ncodemax=ivector(1,NCOVMAX); /* Number of code per covariate; if O and 1 only, 2**ncov; V1+V2+V3+V4=>16 */ |
| ncodemaxwundef=ivector(1,NCOVMAX); /* Number of code per covariate; if - 1 O and 1 only, 2**ncov; V1+V2+V3+V4=>16 */ |
ncodemaxwundef=ivector(1,NCOVMAX); /* Number of code per covariate; if - 1 O and 1 only, 2**ncov; V1+V2+V3+V4=>16 */ |
|
Line 11457 Please run with mle=-1 to get a correct
|
Line 11654 Please run with mle=-1 to get a correct
|
| |
|
| |
|
| agegomp=(int)agemin; |
agegomp=(int)agemin; |
| free_vector(moisnais,1,n); |
free_vector(moisnais,firstobs,lastobs); |
| free_vector(annais,1,n); |
free_vector(annais,firstobs,lastobs); |
| /* free_matrix(mint,1,maxwav,1,n); |
/* free_matrix(mint,1,maxwav,1,n); |
| free_matrix(anint,1,maxwav,1,n);*/ |
free_matrix(anint,1,maxwav,1,n);*/ |
| /* free_vector(moisdc,1,n); */ |
/* free_vector(moisdc,1,n); */ |
|
Line 11484 Please run with mle=-1 to get a correct
|
Line 11681 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 */ |
/* Concatenates waves */ |
| |
|
| free_vector(moisdc,1,n); |
free_vector(moisdc,firstobs,lastobs); |
| free_vector(andc,1,n); |
free_vector(andc,firstobs,lastobs); |
| |
|
| /* Routine tricode is to calculate cptcoveff (real number of unique covariates) and to associate covariable number and modality */ |
/* Routine tricode is to calculate cptcoveff (real number of unique covariates) and to associate covariable number and modality */ |
| nbcode=imatrix(0,NCOVMAX,0,NCOVMAX); |
nbcode=imatrix(0,NCOVMAX,0,NCOVMAX); |
|
Line 11675 Title=%s <br>Datafile=%s Firstpass=%d La
|
Line 11872 Title=%s <br>Datafile=%s Firstpass=%d La
|
| fprintf(fichtm,"\n<li> Number of fixed quantitative variables: nqv=%d ", nqv); |
fprintf(fichtm,"\n<li> Number of fixed quantitative variables: nqv=%d ", nqv); |
| ncurrv=i; |
ncurrv=i; |
| for(i=ncurrv; i <=ncurrv-1+nqv; i++) fprintf(fichtm,"V%d ", i); |
for(i=ncurrv; i <=ncurrv-1+nqv; i++) fprintf(fichtm,"V%d ", i); |
| fprintf(fichtm,"\n<li> Number of time varying (wave varying) covariates: ntv=%d ", ntv); |
fprintf(fichtm,"\n<li> Number of time varying (wave varying) dummy covariates: ntv=%d ", ntv); |
| ncurrv=i; |
ncurrv=i; |
| for(i=ncurrv; i <=ncurrv-1+ntv; i++) fprintf(fichtm,"V%d ", i); |
for(i=ncurrv; i <=ncurrv-1+ntv; i++) fprintf(fichtm,"V%d ", i); |
| fprintf(fichtm,"\n<li>Number of quantitative time varying covariates: nqtv=%d ", nqtv); |
fprintf(fichtm,"\n<li>Number of time varying quantitative covariates: nqtv=%d ", nqtv); |
| ncurrv=i; |
ncurrv=i; |
| for(i=ncurrv; i <=ncurrv-1+nqtv; i++) fprintf(fichtm,"V%d ", i); |
for(i=ncurrv; i <=ncurrv-1+nqtv; i++) fprintf(fichtm,"V%d ", i); |
| fprintf(fichtm,"\n<li>Weights column \n<br>Number of alive states: nlstate=%d <br>Number of death states (not really implemented): ndeath=%d \n<li>Number of waves: maxwav=%d \n<li>Parameter for maximization (1), using parameter values (0), for design of parameters and variance-covariance matrix: mle=%d \n<li>Does the weight column be taken into account (1), or not (0): weight=%d</ul>\n", \ |
fprintf(fichtm,"\n<li>Weights column \n<br>Number of alive states: nlstate=%d <br>Number of death states (not really implemented): ndeath=%d \n<li>Number of waves: maxwav=%d \n<li>Parameter for maximization (1), using parameter values (0), for design of parameters and variance-covariance matrix: mle=%d \n<li>Does the weight column be taken into account (1), or not (0): weight=%d</ul>\n", \ |
|
Line 11710 Interval (in months) between two waves:
|
Line 11907 Interval (in months) between two waves:
|
| 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(firstobs,lastobs); |
| ageexmed=vector(1,n); |
ageexmed=vector(firstobs,lastobs); |
| agecens=vector(1,n); |
agecens=vector(firstobs,lastobs); |
| dcwave=ivector(1,n); |
dcwave=ivector(firstobs,lastobs); |
| |
|
| for (i=1; i<=imx; i++){ |
for (i=1; i<=imx; i++){ |
| dcwave[i]=-1; |
dcwave[i]=-1; |
|
Line 11747 Interval (in months) between two waves:
|
Line 11944 Interval (in months) between two waves:
|
| ximort[i][j]=(i == j ? 1.0 : 0.0); |
ximort[i][j]=(i == j ? 1.0 : 0.0); |
| } |
} |
| |
|
| /*p[1]=0.0268; p[NDIM]=0.083;*/ |
p[1]=0.0268; p[NDIM]=0.083; |
| /*printf("%lf %lf", p[1], p[2]);*/ |
/* printf("%lf %lf", p[1], p[2]); */ |
| |
|
| |
|
| #ifdef GSL |
#ifdef GSL |
|
Line 11874 Interval (in months) between two waves:
|
Line 12071 Interval (in months) between two waves:
|
| printf("%f [%f ; %f]\n",p[i],p[i]-2*sqrt(matcov[i][i]),p[i]+2*sqrt(matcov[i][i])); |
printf("%f [%f ; %f]\n",p[i],p[i]-2*sqrt(matcov[i][i]),p[i]+2*sqrt(matcov[i][i])); |
| fprintf(ficlog,"%f [%f ; %f]\n",p[i],p[i]-2*sqrt(matcov[i][i]),p[i]+2*sqrt(matcov[i][i])); |
fprintf(ficlog,"%f [%f ; %f]\n",p[i],p[i]-2*sqrt(matcov[i][i]),p[i]+2*sqrt(matcov[i][i])); |
| } |
} |
| lsurv=vector(1,AGESUP); |
lsurv=vector(agegomp,AGESUP); |
| lpop=vector(1,AGESUP); |
lpop=vector(agegomp,AGESUP); |
| tpop=vector(1,AGESUP); |
tpop=vector(agegomp,AGESUP); |
| lsurv[agegomp]=100000; |
lsurv[agegomp]=100000; |
| |
|
| for (k=agegomp;k<=AGESUP;k++) { |
for (k=agegomp;k<=AGESUP;k++) { |
|
Line 11923 Please run with mle=-1 to get a correct
|
Line 12120 Please run with mle=-1 to get a correct
|
| stepm, weightopt,\ |
stepm, weightopt,\ |
| model,imx,p,matcov,agemortsup); |
model,imx,p,matcov,agemortsup); |
| |
|
| free_vector(lsurv,1,AGESUP); |
free_vector(lsurv,agegomp,AGESUP); |
| free_vector(lpop,1,AGESUP); |
free_vector(lpop,agegomp,AGESUP); |
| free_vector(tpop,1,AGESUP); |
free_vector(tpop,agegomp,AGESUP); |
| free_matrix(ximort,1,NDIM,1,NDIM); |
free_matrix(ximort,1,NDIM,1,NDIM); |
| free_ivector(cens,1,n); |
free_ivector(dcwave,firstobs,lastobs); |
| free_vector(agecens,1,n); |
free_vector(agecens,firstobs,lastobs); |
| free_ivector(dcwave,1,n); |
free_vector(ageexmed,firstobs,lastobs); |
| |
free_ivector(cens,firstobs,lastobs); |
| #ifdef GSL |
#ifdef GSL |
| #endif |
#endif |
| } /* Endof if mle==-3 mortality only */ |
} /* Endof if mle==-3 mortality only */ |
|
Line 12125 Please run with mle=-1 to get a correct
|
Line 12323 Please run with mle=-1 to get a correct
|
| fputs(line,stdout); |
fputs(line,stdout); |
| fputs(line,ficparo); |
fputs(line,ficparo); |
| fputs(line,ficlog); |
fputs(line,ficlog); |
| |
fputs(line,ficres); |
| continue; |
continue; |
| }else |
}else |
| break; |
break; |
|
Line 12170 Please run with mle=-1 to get a correct
|
Line 12369 Please run with mle=-1 to get a correct
|
| fputs(line,stdout); |
fputs(line,stdout); |
| fputs(line,ficparo); |
fputs(line,ficparo); |
| fputs(line,ficlog); |
fputs(line,ficlog); |
| |
fputs(line,ficres); |
| continue; |
continue; |
| }else |
}else |
| break; |
break; |
|
Line 12195 Please run with mle=-1 to get a correct
|
Line 12395 Please run with mle=-1 to get a correct
|
| fputs(line,stdout); |
fputs(line,stdout); |
| fputs(line,ficparo); |
fputs(line,ficparo); |
| fputs(line,ficlog); |
fputs(line,ficlog); |
| |
fputs(line,ficres); |
| continue; |
continue; |
| }else |
}else |
| break; |
break; |
|
Line 12217 Please run with mle=-1 to get a correct
|
Line 12418 Please run with mle=-1 to get a correct
|
| } |
} |
| |
|
| /* Results */ |
/* Results */ |
| |
endishere=0; |
| nresult=0; |
nresult=0; |
| do{ |
do{ |
| if(!fgets(line, MAXLINE, ficpar)){ |
if(!fgets(line, MAXLINE, ficpar)){ |
| endishere=1; |
endishere=1; |
| parameterline=14; |
|
| }else if (line[0] == '#') { |
}else if (line[0] == '#') { |
| /* If line starts with a # it is a comment */ |
/* If line starts with a # it is a comment */ |
| numlinepar++; |
numlinepar++; |
| fputs(line,stdout); |
fputs(line,stdout); |
| fputs(line,ficparo); |
fputs(line,ficparo); |
| fputs(line,ficlog); |
fputs(line,ficlog); |
| |
fputs(line,ficres); |
| continue; |
continue; |
| }else if(sscanf(line,"prevforecast=%[^\n]\n",modeltemp)) |
}else if(sscanf(line,"prevforecast=%[^\n]\n",modeltemp)) |
| parameterline=11; |
parameterline=11; |
| else if(sscanf(line,"backcast=%[^\n]\n",modeltemp)) |
else if(sscanf(line,"prevbackcast=%[^\n]\n",modeltemp)) |
| parameterline=12; |
parameterline=12; |
| else if(sscanf(line,"result:%[^\n]\n",modeltemp)) |
else if(sscanf(line,"result:%[^\n]\n",modeltemp)){ |
| parameterline=13; |
parameterline=13; |
| |
} |
| else{ |
else{ |
| parameterline=14; |
parameterline=14; |
| } |
} |
| switch (parameterline){ |
switch (parameterline){ |
| case 11: |
case 11: |
| 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=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 && (num_filled == 8)){ |
| if (num_filled != 8) { |
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("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); |
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(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); |
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.*/ |
| dateproj1=anproj1+(mproj1-1)/12.+(jproj1-1)/365.; |
dateproj1=anproj1+(mproj1-1)/12.+(jproj1-1)/365.; |
| dateproj2=anproj2+(mproj2-1)/12.+(jproj2-1)/365.; |
dateproj2=anproj2+(mproj2-1)/12.+(jproj2-1)/365.; |
| |
prvforecast = 1; |
| |
} |
| |
else if((num_filled=sscanf(line,"prevforecast=%d yearsfproj=%lf mobil_average=%d\n",&prevfcast,&yrfproj,&mobilavproj)) !=EOF){/* && (num_filled == 3))*/ |
| |
printf("prevforecast=%d yearsfproj=%lf.2 mobil_average=%d\n",prevfcast,yrfproj,mobilavproj); |
| |
fprintf(ficlog,"prevforecast=%d yearsfproj=%lf.2 mobil_average=%d\n",prevfcast,yrfproj,mobilavproj); |
| |
fprintf(ficres,"prevforecast=%d yearsfproj=%lf.2 mobil_average=%d\n",prevfcast,yrfproj,mobilavproj); |
| |
prvforecast = 2; |
| |
} |
| |
else { |
| |
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\nnor 3 (data)parameters, for example:prevforecast=1 yearsfproj=10 mobil_average=0. Your line=%s . You are running probably 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 mobil_average=0\nnor 3 (data)parameters, for example:prevforecast=1 yearproj=10 mobil_average=0. Your line=%s . You are running probably an older format.\n, ",num_filled,line); |
| |
goto end; |
| } |
} |
| break; |
break; |
| case 12: |
case 12: |
| /*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);*/ |
if((num_filled=sscanf(line,"prevbackcast=%d starting-back-date=%lf/%lf/%lf final-back-date=%lf/%lf/%lf mobil_average=%d\n",&prevbcast,&jback1,&mback1,&anback1,&jback2,&mback2,&anback2,&mobilavproj)) !=EOF && (num_filled == 8)){ |
| 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(ficparo,"prevbackcast=%d starting-back-date=%.lf/%.lf/%.lf final-back-date=%.lf/%.lf/%.lf mobil_average=%d\n",prevbcast,jback1,mback1,anback1,jback2,mback2,anback2,mobilavproj); |
| if (num_filled != 8) { |
printf("prevbackcast=%d starting-back-date=%.lf/%.lf/%.lf final-back-date=%.lf/%.lf/%.lf mobil_average=%d\n",prevbcast,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 final-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,"prevbackcast=%d starting-back-date=%.lf/%.lf/%.lf final-back-date=%.lf/%.lf/%.lf mobil_average=%d\n",prevbcast,jback1,mback1,anback1,jback2,mback2,anback2,mobilavproj); |
| fprintf(ficlog,"Error: Not 8 (data)parameters in line but %d, for example:backcast=1 starting-back-date=1/1/1990 final-back-date=1/1/1970 mobil_average=1\n, your line=%s . Probably you are running an older format.\n",num_filled,line); |
fprintf(ficres,"prevbackcast=%d starting-back-date=%.lf/%.lf/%.lf final-back-date=%.lf/%.lf/%.lf mobil_average=%d\n",prevbcast,jback1,mback1,anback1,jback2,mback2,anback2,mobilavproj); |
| goto end; |
/* day and month of back2 are not used but only year anback2.*/ |
| } |
|
| 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.*/ |
|
| dateback1=anback1+(mback1-1)/12.+(jback1-1)/365.; |
dateback1=anback1+(mback1-1)/12.+(jback1-1)/365.; |
| dateback2=anback2+(mback2-1)/12.+(jback2-1)/365.; |
dateback2=anback2+(mback2-1)/12.+(jback2-1)/365.; |
| |
prvbackcast = 1; |
| |
} |
| |
else if((num_filled=sscanf(line,"prevbackcast=%d yearsbproj=%lf mobil_average=%d\n",&prevbcast,&yrbproj,&mobilavproj)) ==3){/* && (num_filled == 3))*/ |
| |
printf("prevbackcast=%d yearsbproj=%lf.2 mobil_average=%d\n",prevbcast,yrbproj,mobilavproj); |
| |
fprintf(ficlog,"prevbackcast=%d yearsbproj=%lf.2 mobil_average=%d\n",prevbcast,yrbproj,mobilavproj); |
| |
fprintf(ficres,"prevbackcast=%d yearsbproj=%lf.2 mobil_average=%d\n",prevbcast,yrbproj,mobilavproj); |
| |
prvbackcast = 2; |
| |
} |
| |
else { |
| |
printf("Error: Not 8 (data)parameters in line but %d, for example:prevbackcast=1 starting-back-date=1/1/1990 final-back-date=1/1/2000 mobil_average=0\nnor 3 (data)parameters, for example:prevbackcast=1 yearsbproj=10 mobil_average=0. Your line=%s . You are running probably an older format.\n, ",num_filled,line); |
| |
fprintf(ficlog,"Error: Not 8 (data)parameters in line but %d, for example:prevbackcast=1 starting-back-date=1/1/1990 final-back-date=1/1/2000 mobil_average=0\nnor 3 (data)parameters, for example:prevbackcast=1 yearbproj=10 mobil_average=0. Your line=%s . You are running probably an older format.\n, ",num_filled,line); |
| |
goto end; |
| } |
} |
| break; |
break; |
| case 13: |
case 13: |
| if((num_filled=sscanf(line,"result:%[^\n]\n",resultline)) !=EOF){ |
num_filled=sscanf(line,"result:%[^\n]\n",resultline); |
| if (num_filled == 0){ |
nresult++; /* Sum of resultlines */ |
| resultline[0]='\0'; |
printf("Result %d: result:%s\n",nresult, resultline); |
| printf("Warning %d: no result line! It should be at minimum 'result: V2=0 V1=1 or result:.\n%s\n", num_filled, line); |
if(nresult > MAXRESULTLINES){ |
| fprintf(ficlog,"Warning %d: no result line! It should be at minimum 'result: V2=0 V1=1 or result:.\n%s\n", num_filled, line); |
printf("ERROR: Current version of IMaCh limits the number of resultlines to %d, you used %d\nYou can use the 'r' parameter file '%s' which uses option mle=0 to get other results. ",MAXRESULTLINES,nresult,rfileres); |
| break; |
fprintf(ficlog,"ERROR: Current version of IMaCh limits the number of resultlines to %d, you used %d\nYou can use the 'r' parameter file '%s' which uses option mle=0 to get other results. ",MAXRESULTLINES,nresult,rfileres); |
| } else if (num_filled != 1){ |
goto end; |
| printf("ERROR %d: result line! It should be at minimum 'result: V2=0 V1=1 or result:.' %s\n",num_filled, line); |
} |
| fprintf(ficlog,"ERROR %d: result line! It should be at minimum 'result: V2=0 V1=1 or result:.' %s\n",num_filled, line); |
decoderesult(resultline, nresult); /* Fills TKresult[nresult] combination and Tresult[nresult][k4+1] combination values */ |
| } |
fprintf(ficparo,"result: %s\n",resultline); |
| nresult++; /* Sum of resultlines */ |
fprintf(ficres,"result: %s\n",resultline); |
| printf("Result %d: result=%s\n",nresult, resultline); |
fprintf(ficlog,"result: %s\n",resultline); |
| if(nresult > MAXRESULTLINES){ |
break; |
| printf("ERROR: Current version of IMaCh limits the number of resultlines to %d, you used %d\n",MAXRESULTLINES,nresult); |
case 14: |
| fprintf(ficlog,"ERROR: Current version of IMaCh limits the number of resultlines to %d, you used %d\n",MAXRESULTLINES,nresult); |
printf("Error: Unknown command '%s'\n",line); |
| goto end; |
fprintf(ficlog,"Error: Unknown command '%s'\n",line); |
| } |
if(ncovmodel >=2 && nresult==0 ){ |
| decoderesult(resultline, nresult); /* Fills TKresult[nresult] combination and Tresult[nresult][k4+1] combination values */ |
printf("ERROR: no result lines! It should be at minimum 'result: V2=0 V1=1 or result:.' %s\n",line); |
| fprintf(ficparo,"result: %s\n",resultline); |
fprintf(ficlog,"ERROR: no result lines! It should be at minimum 'result: V2=0 V1=1 or result:.' %s\n",line); |
| fprintf(ficres,"result: %s\n",resultline); |
|
| fprintf(ficlog,"result: %s\n",resultline); |
|
| break; |
|
| case 14: |
|
| if(ncovmodel >2 && nresult==0 ){ |
|
| printf("ERROR: no result lines! It should be at minimum 'result: V2=0 V1=1 or result:.' %s\n",line); |
|
| goto end; |
|
| } |
|
| break; |
|
| default: |
|
| nresult=1; |
|
| decoderesult(".",nresult ); /* No covariate */ |
|
| } |
} |
| |
/* goto end; */ |
| |
break; |
| |
default: |
| |
nresult=1; |
| |
decoderesult(".",nresult ); /* No covariate */ |
| } /* End switch parameterline */ |
} /* End switch parameterline */ |
| |
parameterline=0; |
| }while(endishere==0); /* End do */ |
}while(endishere==0); /* End do */ |
| |
|
| /* freqsummary(fileres, agemin, agemax, s, agev, nlstate, imx,Tvaraff,nbcode, ncodemax,mint,anint); */ |
/* freqsummary(fileres, agemin, agemax, s, agev, nlstate, imx,Tvaraff,nbcode, ncodemax,mint,anint); */ |
|
Line 12322 This is probably because your parameter
|
Line 12531 This is probably because your parameter
|
| Please run with mle=-1 to get a correct covariance matrix.\n",ageminpar,agemaxpar); |
Please run with mle=-1 to get a correct covariance matrix.\n",ageminpar,agemaxpar); |
| }else{ |
}else{ |
| /* printinggnuplot(fileresu, optionfilefiname,ageminpar,agemaxpar,fage, prevfcast, backcast, pathc,p, (int)anproj1-(int)agemin, (int)anback1-(int)agemax+1); */ |
/* printinggnuplot(fileresu, optionfilefiname,ageminpar,agemaxpar,fage, prevfcast, backcast, pathc,p, (int)anproj1-(int)agemin, (int)anback1-(int)agemax+1); */ |
| printinggnuplot(fileresu, optionfilefiname,ageminpar,agemaxpar,bage, fage, prevfcast, backcast, pathc,p, (int)anproj1-bage, (int)anback1-fage); |
/* It seems that anprojd which is computed from the mean year at interview which is known yet because of freqsummary */ |
| |
/* date2dmy(dateintmean,&jintmean,&mintmean,&aintmean); */ /* Done in freqsummary */ |
| |
if(prvforecast==1){ |
| |
dateprojd=(jproj1+12*mproj1+365*anproj1)/365; |
| |
jprojd=jproj1; |
| |
mprojd=mproj1; |
| |
anprojd=anproj1; |
| |
dateprojf=(jproj2+12*mproj2+365*anproj2)/365; |
| |
jprojf=jproj2; |
| |
mprojf=mproj2; |
| |
anprojf=anproj2; |
| |
} else if(prvforecast == 2){ |
| |
dateprojd=dateintmean; |
| |
date2dmy(dateprojd,&jprojd, &mprojd, &anprojd); |
| |
dateprojf=dateintmean+yrfproj; |
| |
date2dmy(dateprojf,&jprojf, &mprojf, &anprojf); |
| |
} |
| |
if(prvbackcast==1){ |
| |
datebackd=(jback1+12*mback1+365*anback1)/365; |
| |
jbackd=jback1; |
| |
mbackd=mback1; |
| |
anbackd=anback1; |
| |
datebackf=(jback2+12*mback2+365*anback2)/365; |
| |
jbackf=jback2; |
| |
mbackf=mback2; |
| |
anbackf=anback2; |
| |
} else if(prvbackcast == 2){ |
| |
datebackd=dateintmean; |
| |
date2dmy(datebackd,&jbackd, &mbackd, &anbackd); |
| |
datebackf=dateintmean-yrbproj; |
| |
date2dmy(datebackf,&jbackf, &mbackf, &anbackf); |
| |
} |
| |
|
| |
printinggnuplot(fileresu, optionfilefiname,ageminpar,agemaxpar,bage, fage, prevfcast, prevbcast, pathc,p, (int)anprojd-bage, (int)anbackd-fage); |
| } |
} |
| printinghtml(fileresu,title,datafile, firstpass, lastpass, stepm, weightopt, \ |
printinghtml(fileresu,title,datafile, firstpass, lastpass, stepm, weightopt, \ |
| model,imx,jmin,jmax,jmean,rfileres,popforecast,mobilav,prevfcast,mobilavproj,backcast, estepm, \ |
model,imx,jmin,jmax,jmean,rfileres,popforecast,mobilav,prevfcast,mobilavproj,prevbcast, estepm, \ |
| jprev1,mprev1,anprev1,dateprev1, dateproj1, dateback1,jprev2,mprev2,anprev2,dateprev2,dateproj2, dateback2); |
jprev1,mprev1,anprev1,dateprev1, dateprojd, datebackd,jprev2,mprev2,anprev2,dateprev2,dateprojf, datebackf); |
| |
|
| /*------------ free_vector -------------*/ |
/*------------ free_vector -------------*/ |
| /* chdir(path); */ |
/* chdir(path); */ |
|
Line 12335 Please run with mle=-1 to get a correct
|
Line 12577 Please run with mle=-1 to get a correct
|
| /* free_imatrix(dh,1,lastpass-firstpass+2,1,imx); */ |
/* free_imatrix(dh,1,lastpass-firstpass+2,1,imx); */ |
| /* free_imatrix(bh,1,lastpass-firstpass+2,1,imx); */ |
/* free_imatrix(bh,1,lastpass-firstpass+2,1,imx); */ |
| /* free_imatrix(mw,1,lastpass-firstpass+2,1,imx); */ |
/* free_imatrix(mw,1,lastpass-firstpass+2,1,imx); */ |
| free_lvector(num,1,n); |
free_lvector(num,firstobs,lastobs); |
| free_vector(agedc,1,n); |
free_vector(agedc,firstobs,lastobs); |
| /*free_matrix(covar,0,NCOVMAX,1,n);*/ |
/*free_matrix(covar,0,NCOVMAX,1,n);*/ |
| /*free_matrix(covar,1,NCOVMAX,1,n);*/ |
/*free_matrix(covar,1,NCOVMAX,1,n);*/ |
| fclose(ficparo); |
fclose(ficparo); |
|
Line 12399 Please run with mle=-1 to get a correct
|
Line 12641 Please run with mle=-1 to get a correct
|
| }/* end if moving average */ |
}/* end if moving average */ |
| |
|
| /*---------- Forecasting ------------------*/ |
/*---------- Forecasting ------------------*/ |
| if(prevfcast==1){ |
if(prevfcast==1){ |
| /* if(stepm ==1){*/ |
/* /\* if(stepm ==1){*\/ */ |
| prevforecast(fileresu, anproj1, mproj1, jproj1, agemin, agemax, dateprev1, dateprev2, mobilavproj, mobaverage, bage, fage, firstpass, lastpass, anproj2, p, cptcoveff); |
/* /\* anproj1, mproj1, jproj1 either read explicitly or yrfproj *\/ */ |
| |
/*This done previously after freqsummary.*/ |
| |
/* dateprojd=(jproj1+12*mproj1+365*anproj1)/365; */ |
| |
/* dateprojf=(jproj2+12*mproj2+365*anproj2)/365; */ |
| |
|
| |
/* } else if (prvforecast==2){ */ |
| |
/* /\* if(stepm ==1){*\/ */ |
| |
/* /\* anproj1, mproj1, jproj1 either read explicitly or yrfproj *\/ */ |
| |
/* } */ |
| |
/*prevforecast(fileresu, dateintmean, anproj1, mproj1, jproj1, agemin, agemax, dateprev1, dateprev2, mobilavproj, mobaverage, bage, fage, firstpass, lastpass, anproj2, p, cptcoveff);*/ |
| |
prevforecast(fileresu,dateintmean, dateprojd, dateprojf, agemin, agemax, dateprev1, dateprev2, mobilavproj, mobaverage, bage, fage, firstpass, lastpass, p, cptcoveff); |
| } |
} |
| |
|
| /* Backcasting */ |
/* Prevbcasting */ |
| if(backcast==1){ |
if(prevbcast==1){ |
| ddnewms=matrix(1,nlstate+ndeath,1,nlstate+ndeath); |
ddnewms=matrix(1,nlstate+ndeath,1,nlstate+ndeath); |
| ddoldms=matrix(1,nlstate+ndeath,1,nlstate+ndeath); |
ddoldms=matrix(1,nlstate+ndeath,1,nlstate+ndeath); |
| ddsavms=matrix(1,nlstate+ndeath,1,nlstate+ndeath); |
ddsavms=matrix(1,nlstate+ndeath,1,nlstate+ndeath); |
|
Line 12420 Please run with mle=-1 to get a correct
|
Line 12672 Please run with mle=-1 to get a correct
|
| hBijx(p, bage, fage, mobaverage); |
hBijx(p, bage, fage, mobaverage); |
| fclose(ficrespijb); |
fclose(ficrespijb); |
| |
|
| prevbackforecast(fileresu, mobaverage, anback1, mback1, jback1, agemin, agemax, dateprev1, dateprev2, |
/* /\* prevbackforecast(fileresu, mobaverage, anback1, mback1, jback1, agemin, agemax, dateprev1, dateprev2, *\/ */ |
| mobilavproj, bage, fage, firstpass, lastpass, anback2, p, cptcoveff); |
/* /\* mobilavproj, bage, fage, firstpass, lastpass, anback2, p, cptcoveff); *\/ */ |
| |
/* prevbackforecast(fileresu, mobaverage, anback1, mback1, jback1, agemin, agemax, dateprev1, dateprev2, */ |
| |
/* mobilavproj, bage, fage, firstpass, lastpass, anback2, p, cptcoveff); */ |
| |
prevbackforecast(fileresu, mobaverage, dateintmean, dateprojd, dateprojf, agemin, agemax, dateprev1, dateprev2, |
| |
mobilavproj, bage, fage, firstpass, lastpass, p, cptcoveff); |
| |
|
| |
|
| varbprlim(fileresu, nresult, mobaverage, mobilavproj, bage, fage, bprlim, &ncvyear, ftolpl, p, matcov, delti, stepm, cptcoveff); |
varbprlim(fileresu, nresult, mobaverage, mobilavproj, bage, fage, bprlim, &ncvyear, ftolpl, p, matcov, delti, stepm, cptcoveff); |
| |
|
| |
|
|
Line 12429 Please run with mle=-1 to get a correct
|
Line 12687 Please run with mle=-1 to get a correct
|
| free_matrix(ddnewms, 1, nlstate+ndeath, 1, nlstate+ndeath); |
free_matrix(ddnewms, 1, nlstate+ndeath, 1, nlstate+ndeath); |
| free_matrix(ddsavms, 1, nlstate+ndeath, 1, nlstate+ndeath); |
free_matrix(ddsavms, 1, nlstate+ndeath, 1, nlstate+ndeath); |
| free_matrix(ddoldms, 1, nlstate+ndeath, 1, nlstate+ndeath); |
free_matrix(ddoldms, 1, nlstate+ndeath, 1, nlstate+ndeath); |
| } /* end Backcasting */ |
} /* end Prevbcasting */ |
| |
|
| |
|
| /* ------ Other prevalence ratios------------ */ |
/* ------ Other prevalence ratios------------ */ |
|
Line 12650 Please run with mle=-1 to get a correct
|
Line 12908 Please run with mle=-1 to get a correct
|
| varprlim(fileresu, nresult, mobaverage, mobilavproj, bage, fage, prlim, &ncvyear, ftolpl, p, matcov, delti, stepm, cptcoveff); |
varprlim(fileresu, nresult, mobaverage, mobilavproj, bage, fage, prlim, &ncvyear, ftolpl, p, matcov, delti, stepm, cptcoveff); |
| |
|
| |
|
| free_vector(weight,1,n); |
free_vector(weight,firstobs,lastobs); |
| free_imatrix(Tvard,1,NCOVMAX,1,2); |
free_imatrix(Tvard,1,NCOVMAX,1,2); |
| free_imatrix(s,1,maxwav+1,1,n); |
free_imatrix(s,1,maxwav+1,firstobs,lastobs); |
| free_matrix(anint,1,maxwav,1,n); |
free_matrix(anint,1,maxwav,firstobs,lastobs); |
| free_matrix(mint,1,maxwav,1,n); |
free_matrix(mint,1,maxwav,firstobs,lastobs); |
| free_ivector(cod,1,n); |
free_ivector(cod,firstobs,lastobs); |
| free_ivector(tab,1,NCOVMAX); |
free_ivector(tab,1,NCOVMAX); |
| fclose(ficresstdeij); |
fclose(ficresstdeij); |
| fclose(ficrescveij); |
fclose(ficrescveij); |
|
Line 12675 Please run with mle=-1 to get a correct
|
Line 12933 Please run with mle=-1 to get a correct
|
| free_matrix(oldms, 1,nlstate+ndeath,1,nlstate+ndeath); |
free_matrix(oldms, 1,nlstate+ndeath,1,nlstate+ndeath); |
| free_matrix(newms, 1,nlstate+ndeath,1,nlstate+ndeath); |
free_matrix(newms, 1,nlstate+ndeath,1,nlstate+ndeath); |
| free_matrix(savms, 1,nlstate+ndeath,1,nlstate+ndeath); |
free_matrix(savms, 1,nlstate+ndeath,1,nlstate+ndeath); |
| if(ntv+nqtv>=1)free_ma3x(cotvar,1,maxwav,1,ntv+nqtv,1,n); |
if(ntv+nqtv>=1)free_ma3x(cotvar,1,maxwav,1,ntv+nqtv,firstobs,lastobs); |
| if(nqtv>=1)free_ma3x(cotqvar,1,maxwav,1,nqtv,1,n); |
if(nqtv>=1)free_ma3x(cotqvar,1,maxwav,1,nqtv,firstobs,lastobs); |
| if(nqv>=1)free_matrix(coqvar,1,nqv,1,n); |
if(nqv>=1)free_matrix(coqvar,1,nqv,firstobs,lastobs); |
| free_matrix(covar,0,NCOVMAX,1,n); |
free_matrix(covar,0,NCOVMAX,firstobs,lastobs); |
| free_matrix(matcov,1,npar,1,npar); |
free_matrix(matcov,1,npar,1,npar); |
| free_matrix(hess,1,npar,1,npar); |
free_matrix(hess,1,npar,1,npar); |
| /*free_vector(delti,1,npar);*/ |
/*free_vector(delti,1,npar);*/ |
|
Line 12798 Please run with mle=-1 to get a correct
|
Line 13056 Please run with mle=-1 to get a correct
|
| |
|
| sprintf(plotcmd,"%s %s",pplotcmd, optionfilegnuplot); |
sprintf(plotcmd,"%s %s",pplotcmd, optionfilegnuplot); |
| printf("Starting graphs with: '%s'\n",plotcmd);fflush(stdout); |
printf("Starting graphs with: '%s'\n",plotcmd);fflush(stdout); |
| |
strcpy(pplotcmd,plotcmd); |
| |
|
| if((outcmd=system(plotcmd)) != 0){ |
if((outcmd=system(plotcmd)) != 0){ |
| printf("gnuplot command might not be in your path: '%s', err=%d\n", plotcmd, outcmd); |
printf("Error in gnuplot, command might not be in your path: '%s', err=%d\n", plotcmd, outcmd); |
| printf("\n Trying if gnuplot resides on the same directory that IMaCh\n"); |
printf("\n Trying if gnuplot resides on the same directory that IMaCh\n"); |
| sprintf(plotcmd,"%sgnuplot %s", pathimach, optionfilegnuplot); |
sprintf(plotcmd,"%sgnuplot %s", pathimach, optionfilegnuplot); |
| if((outcmd=system(plotcmd)) != 0) |
if((outcmd=system(plotcmd)) != 0){ |
| printf("\n Still a problem with gnuplot command %s, err=%d\n", plotcmd, outcmd); |
printf("\n Still a problem with gnuplot command %s, err=%d\n", plotcmd, outcmd); |
| |
strcpy(plotcmd,pplotcmd); |
| |
} |
| } |
} |
| printf(" Successful, please wait..."); |
printf(" Successful, please wait..."); |
| while (z[0] != 'q') { |
while (z[0] != 'q') { |
![[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