|
version 1.10, 2001/05/09 14:25:42
|
version 1.11, 2001/05/17 16:07:14
|
|
Line 8
|
Line 8
|
| Health expectancies are computed from the transistions observed between
|
Health expectancies are computed from the transistions observed between
|
| waves and are computed for each degree of severity of disability (number
|
waves and are computed for each degree of severity of disability (number
|
| of life states). More degrees you consider, more time is necessary to
|
of life states). More degrees you consider, more time is necessary to
|
| reach the Maximum Likelihood of the parameters involved in the model.
|
reach the Maximum Likelihood of the parameters involved in the model.
|
| The simplest model is the multinomial logistic model where pij is
|
The simplest model is the multinomial logistic model where pij is
|
| the probabibility to be observed in state j at the second wave conditional
|
the probabibility to be observed in state j at the second wave conditional
|
| to be observed in state i at the first wave. Therefore the model is:
|
to be observed in state i at the first wave. Therefore the model is:
|
|
Line 1270 void concatwav(int wav[], int **dh, int
|
Line 1270 void concatwav(int wav[], int **dh, int
|
| /* int j, k=0,jk, ju, jl,jmin=1e+5, jmax=-1;
|
/* int j, k=0,jk, ju, jl,jmin=1e+5, jmax=-1;
|
| double sum=0., jmean=0.;*/
|
double sum=0., jmean=0.;*/
|
|
|
|
| int j, k=0,jk, ju, jl;
|
int j, k=0,jk, ju, jl;
|
| double sum=0.;
|
double sum=0.;
|
| jmin=1e+5;
|
jmin=1e+5;
|
| jmax=-1;
|
jmax=-1;
|
| jmean=0.;
|
jmean=0.;
|
| for(i=1; i<=imx; i++){
|
for(i=1; i<=imx; i++){
|
| mi=0;
|
mi=0;
|
| m=firstpass;
|
m=firstpass;
|
|
Line 1309 jmean=0.;
|
Line 1309 jmean=0.;
|
| if(j==0) j=1; /* Survives at least one month after exam */
|
if(j==0) j=1; /* Survives at least one month after exam */
|
| k=k+1;
|
k=k+1;
|
| if (j >= jmax) jmax=j;
|
if (j >= jmax) jmax=j;
|
| else if (j <= jmin)jmin=j;
|
if (j <= jmin) jmin=j;
|
| sum=sum+j;
|
sum=sum+j;
|
| if (j<0) printf("j=%d num=%d ",j,i);
|
if (j<0) printf("j=%d num=%d ",j,i);
|
| }
|
}
|
|
Line 2671 strcpy(fileresvpl,"vpl");
|
Line 2671 strcpy(fileresvpl,"vpl");
|
| chdir(pathcd);
|
chdir(pathcd);
|
| #endif
|
#endif
|
| /*system("wgnuplot graph.plt");*/
|
/*system("wgnuplot graph.plt");*/
|
| system("../gp37mgw/wgnuplot graph.plt");
|
/*system("../gp37mgw/wgnuplot graph.plt");*/
|
| |
/*system("cd ../gp37mgw");*/
|
| |
system("..\\gp37mgw\\wgnuplot graph.plt");
|
|
|
|
| #ifdef windows
|
#ifdef windows
|
| 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