version 1.326, 2022/07/26 17:33:55
|
version 1.328, 2022/07/27 17:40:48
|
Line 1
|
Line 1
|
/* $Id$ |
/* $Id$ |
$State$ |
$State$ |
$Log$ |
$Log$ |
|
Revision 1.328 2022/07/27 17:40:48 brouard |
|
Summary: valgrind bug fixed by initializing to zero DummyV as well as Tage |
|
|
|
Revision 1.327 2022/07/27 14:47:35 brouard |
|
Summary: Still a problem for one-step probabilities in case of quantitative variables |
|
|
Revision 1.326 2022/07/26 17:33:55 brouard |
Revision 1.326 2022/07/26 17:33:55 brouard |
Summary: some test with nres=1 |
Summary: some test with nres=1 |
|
|
Line 6968 To be simple, these graphs help to under
|
Line 6974 To be simple, these graphs help to under
|
cov[2+nagesqr+Tage[k]]=nbcode[Tvar[Tage[k]]][codtabm(j1,k)]*cov[2]; |
cov[2+nagesqr+Tage[k]]=nbcode[Tvar[Tage[k]]][codtabm(j1,k)]*cov[2]; |
/* cov[++k1]=nbcode[Tvar[Tage[k]]][codtabm(ij,k)]*cov[2]; */ |
/* cov[++k1]=nbcode[Tvar[Tage[k]]][codtabm(ij,k)]*cov[2]; */ |
} else if(Dummy[Tage[k]]==3){ /* quantitative with age */ |
} else if(Dummy[Tage[k]]==3){ /* quantitative with age */ |
cov[2+nagesqr+Tage[k]]=Tqresult[nres][k]; |
printf("Internal IMaCh error, don't know which value for quantitative covariate with age, Tage[k]%d, k=%d, Tvar[Tage[k]]=V%d, age=%d\n",Tage[k],k ,Tvar[Tage[k]], (int)cov[2]); |
|
exit(1); |
|
/* cov[2+nagesqr+Tage[k]]=meanq[k]/idq[k]*cov[2];/\* Using the mean of quantitative variable Tvar[Tage[k]] /\* Tqresult[nres][k]; *\/ */ |
/* cov[++k1]=Tqresult[nres][k]; */ |
/* cov[++k1]=Tqresult[nres][k]; */ |
} |
} |
/* cov[2+Tage[k]+nagesqr]=nbcode[Tvar[Tage[k]]][codtabm(ij,k)]*cov[2]; */ |
/* cov[2+Tage[k]+nagesqr]=nbcode[Tvar[Tage[k]]][codtabm(ij,k)]*cov[2]; */ |
Line 9595 int readdata(char datafile[], int firsto
|
Line 9603 int readdata(char datafile[], int firsto
|
|
|
DummyV=ivector(1,NCOVMAX); /* 1 to 3 */ |
DummyV=ivector(1,NCOVMAX); /* 1 to 3 */ |
FixedV=ivector(1,NCOVMAX); /* 1 to 3 */ |
FixedV=ivector(1,NCOVMAX); /* 1 to 3 */ |
|
for(v=1;v<NCOVMAX;v++){ |
|
DummyV[v]=0; |
|
FixedV[v]=0; |
|
} |
|
|
for(v=1; v <=ncovcol;v++){ |
for(v=1; v <=ncovcol;v++){ |
DummyV[v]=0; |
DummyV[v]=0; |
Line 11997 Please run with mle=-1 to get a correct
|
Line 12009 Please run with mle=-1 to get a correct
|
Tage=ivector(1,NCOVMAX); /* Gives the covariate id of covariates associated with age: V2 + V1 + age*V4 + V3*age |
Tage=ivector(1,NCOVMAX); /* Gives the covariate id of covariates associated with age: V2 + V1 + age*V4 + V3*age |
4 covariates (3 plus signs) |
4 covariates (3 plus signs) |
Tage[1=V3*age]= 4; Tage[2=age*V4] = 3 |
Tage[1=V3*age]= 4; Tage[2=age*V4] = 3 |
*/ |
*/ |
|
for(i=1;i<NCOVMAX;i++) |
|
Tage[i]=0; |
Tmodelind=ivector(1,NCOVMAX);/** gives the k model position of an |
Tmodelind=ivector(1,NCOVMAX);/** gives the k model position of an |
* individual dummy, fixed or varying: |
* individual dummy, fixed or varying: |
* Tmodelind[Tvaraff[3]]=9,Tvaraff[1]@9={4, |
* Tmodelind[Tvaraff[3]]=9,Tvaraff[1]@9={4, |