Intercept and age are systematically included in the model. -Additional covariates (actually two) can be included with the command:
+Additional covariates can be included with the command:model=list of covariates@@ -368,6 +373,19 @@ Additional covariates (actually two) can the product covariate*age +
In this example, we have two covariates in the data file +(fields 2 and 3). The number of covariates is defined with +statement ncov=2. If now you have 3 covariates in the datafile +(fields 2, 3 and 4), you have to set ncov=3. Then you can run the +programme with a new parametrisation taking into account the +third covariate. For example, model=V1+V3 estimates +a model with the first and third covariates. More complicated +models can be used, but it will takes more time to converge. With +a simple model (no covariates), the programme estimates 8 +parameters. Adding covariates increases the number of parameters +: 12 for model=V1, 16 for model=V1+V1*age +and 20 for model=V1+V2+V3.
+or, to simplify (in most of cases it converges but there is no warranty!):
+or, to simplify (in most of cases it converges but there is no +warranty!):
12 0.0 0.0 @@ -484,15 +503,15 @@ prevalences and health expectancies - +102. It is possible to get extrapolated stationary prevalence by +age ranging from agemin to agemax. -Setting bage=50 (begin age) and fage=100 (final age), makes the program computing -life expectancy from age 'bage' to age 'fage'. As we use a model, we -can interessingly compute life expectancy on a wider age range than the age -range from the data. But the model can be rather wrong on much larger -intervals. Program is limited to around 120 for upper age!
+Setting bage=50 (begin age) and fage=100 (final age), makes +the program computing life expectancy from age 'bage' to age +'fage'. As we use a model, we can interessingly compute life +expectancy on a wider age range than the age range from the data. +But the model can be rather wrong on much larger intervals. +Program is limited to around 120 for upper age!
- agemin= Minimum age for calculation of the @@ -527,47 +546,47 @@ expectancies
pop_based=0-The program computes status-based health expectancies, i.e health -expectancies which depends on your initial health state. If you are -healthy your healthy life expectancy (e11) is higher than if you were -disabled (e21, with e11 > e21).
-To compute a healthy life expectancy independant of the initial status -we have to weight e11 and e21 according to the probability to be in -each state at initial age or, with other word, according to the -proportion of people in each state.
- -We prefer computing a 'pure' period healthy life expectancy based only -on the transtion forces. Then the weights are simply the stationnary -prevalences or 'implied' prevalences at the initial age.
- -Some other people would like to use the cross-sectional prevalences -(the "Sullivan prevalences") observed at the initial age during a -period of time defined just above. +The program computes status-based health expectancies, i.e +health expectancies which depends on your initial health state. +If you are healthy your healthy life expectancy (e11) is higher +than if you were disabled (e21, with e11 > e21).
+To compute a healthy life expectancy independant of the initial +status we have to weight e11 and e21 according to the probability +to be in each state at initial age or, with other word, according +to the proportion of people in each state.
+We prefer computing a 'pure' period healthy life expectancy based +only on the transtion forces. Then the weights are simply the +stationnary prevalences or 'implied' prevalences at the initial +age.
+Some other people would like to use the cross-sectional +prevalences (the "Sullivan prevalences") observed at +the initial age during a period of time defined +just above.-
- -- popbased= 0 Health expectancies are computed - at each age from stationary prevalences 'expected' at this initial age.
-- popbased= 1 Health expectancies are computed - at each age from cross-sectional 'observed' prevalence at this - initial age. As all the population is not observed at the same exact date we - define a short period were the observed prevalence is computed.
+- popbased= 0 Health expectancies are + computed at each age from stationary prevalences + 'expected' at this initial age.
+- popbased= 1 Health expectancies are + computed at each age from cross-sectional 'observed' + prevalence at this initial age. As all the population is + not observed at the same exact date we define a short + period were the observed prevalence is computed.
Prevalence forecasting ( Experimental)
starting-proj-date=1/1/1989 final-proj-date=1/1/1992 mov_average=0-Prevalence and population projections are only available if the -interpolation unit is a month, i.e. stepm=1 and if there are no -covariate. The programme estimates the prevalence in each state at a -precise date expressed in day/month/year. The programme computes one -forecasted prevalence a year from a starting date (1 january of 1989 -in this example) to a final date (1 january 1992). The statement -mov_average allows to compute smoothed forecasted prevalences with a -five-age moving average centered at the mid-age of the five-age -period.
+Prevalence and population projections are only available if +the interpolation unit is a month, i.e. stepm=1 and if there are +no covariate. The programme estimates the prevalence in each +state at a precise date expressed in day/month/year. The +programme computes one forecasted prevalence a year from a +starting date (1 january of 1989 in this example) to a final date +(1 january 1992). The statement mov_average allows to compute +smoothed forecasted prevalences with a five-age moving average +centered at the mid-age of the five-age period.
- starting-proj-date= starting date @@ -1000,7 +1019,7 @@ are in state 2. One year latter, 512892
-Trying an example
+Trying an example
Since you know how to run the program, it is time to test it on your own computer. Try for example on a parameter file named mailto:bro href="mailto:lievre@ined.fr">mailto:lievre@ined.fr .
Latest version (0.71a of March 2002) can be accessed at http://euroreves.ined.fr/imach
+href="http://euroreves.ined.fr/imach">http://euroreves.ined.fr/imach