/* $Id$
$State$
$Log$
+ Revision 1.319 2022/06/02 04:45:11 brouard
+ * imach.c (Module): Adding the Wald tests from the log to the main
+ htm for better display of the maximum likelihood estimators.
+
Revision 1.318 2022/05/24 08:10:59 brouard
* imach.c (Module): Some attempts to find a bug of wrong estimates
of confidencce intervals with product in the equation modelC
/* Tage[cptcovage]=k 5 8 */ /* Position in the model of ith cov*age */
/* Tvard[1][1]@4={4,3,1,2} V4*V3 V1*V2 */ /* Position in model of the ith prod without age */
/* TvarF TvarF[1]=Tvar[6]=2, TvarF[2]=Tvar[7]=7, TvarF[3]=Tvar[9]=1 ID of fixed covariates or product V2, V1*V2, V1 */
-/* TvarFind; /**< TvarFind[1]=6, TvarFind[2]=7, TvarFind[3]=9 *//* Inverse V2(6) is first fixed (single or prod) */
+/* TvarFind; TvarFind[1]=6, TvarFind[2]=7, TvarFind[3]=9 *//* Inverse V2(6) is first fixed (single or prod) */
/* Type */
/* V 1 2 3 4 5 */
/* F F V V V */
/* # V1=sex, V2=raedyrs Quant Fixed, State=livarnb4..livarnb11, V3=iadl4..iald11, V4=adlw4..adlw11, V5=r4bmi..r11bmi */
/* V5+V4+V3+V4*V3+V5*age+V2+V1*V2+V1*age+V1 */
/* TvarF[1]=Tvar[6]=2, TvarF[2]=Tvar[7]=7, TvarF[3]=Tvar[9]=1 ID of fixed covariates or product V2, V1*V2, V1 */
- /* TvarFind; /**< TvarFind[1]=6, TvarFind[2]=7, TvarFind[3]=9 *//* Inverse V2(6) is first fixed (single or prod) */
+ /* TvarFind; TvarFind[1]=6, TvarFind[2]=7, TvarFind[3]=9 *//* Inverse V2(6) is first fixed (single or prod) */
cov[ioffset+TvarFind[k]]=covar[Tvar[TvarFind[k]]][i];/* V5+V4+V3+V4*V3+V5*age+V2+V1*V2+V1*age+V1, only V1 is fixed (TvarFind[1]=6)*/
/* V1*V2 (7) TvarFind[2]=7, TvarFind[3]=9 */
}
varhe[ij][ji][(int)age] += doldm[ij][ji]*hf*hf;
}
}
- if((int)age ==50){
- printf(" age=%d cij=%d nres=%d varhe[%d][%d]=%f ",(int)age, cij, nres, 1,2,varhe[1][2]);
- }
+ /* if((int)age ==50){ */
+ /* printf(" age=%d cij=%d nres=%d varhe[%d][%d]=%f ",(int)age, cij, nres, 1,2,varhe[1][2]); */
+ /* } */
/* Computing expectancies */
hpxij(p3matm,nhstepm,age,hstepm,x,nlstate,stepm,oldm, savm, cij,nres);
for(i=1; i<=nlstate;i++)