--- imach096d/doc/imach.htm 2002/03/11 15:24:05 1.9 +++ imach096d/doc/imach.htm 2002/03/11 22:52:27 1.11 @@ -1,4 +1,4 @@ - + @@ -430,9 +430,40 @@ a large stepm i.e stepm=12 or 24 and the stepm=1 month. If newstepm is the new shorter stepm and stepm can be expressed as a multiple of newstepm, like newstepm=n stepm, then the following approximation holds: -
aij(n stepm) = aij(stepm) +ln(n)
+
aij(stepm) = aij(n . stepm) - ln(n)
 
 and
-
bij(n stepm) = bij(stepm) .

+
bij(stepm) = bij(n . stepm) .

+
+
 For example if you already ran for a 6 months interval and
+got:

+ 
# Parameters
+12 -13.390179  0.126133 
+13  -7.493460  0.048069 
+21   0.575975 -0.041322 
+23  -4.748678  0.030626 
+

+If you now want to get the monthly estimates, you can guess the aij by
+substracting ln(6)= 1,7917
 and running

+
12 -15.18193847  0.126133 
+13 -9.285219469  0.048069
+21 -1.215784469 -0.041322
+23 -6.540437469  0.030626
+

+and get

+
12 -15.029768 0.124347 
+13 -8.472981 0.036599 
+21 -1.472527 -0.038394 
+23 -6.553602 0.029856 
+

+which is closer to the results. The approximation is probably useful
+only for very small intervals and we don't have enough experience to
+know if you will speed up the convergence or not.
+
         -ln(12)= -2.484
+ -ln(6/1)=-ln(6)= -1.791
+ -ln(3/1)=-ln(3)= -1.0986
+-ln(12/6)=-ln(2)= -0.693
+

+
 
Guess values for computing variances

 
 
This is an output if mle=1. But it can be
@@ -446,20 +477,13 @@ matrix of the parameters, that is the in
 matrix, and the variances of health expectancies. Each line
 consists in indices "ij" followed by the initial scales
 (zero to simplify) associated with aij and bij. 

-
-
-
-

-    
# Scales (for hessian or gradient estimation)
+
     
  • If mle=1 you can enter zeros:
    • 
+
    # Scales (for hessian or gradient estimation)
 12 0. 0. 
 13 0. 0. 
 21 0. 0. 
 23 0. 0. 
    
 
    
-
-
      
     
    • If mle=0 you must enter a covariance matrix (usually
         obtained from an earlier run).
      • 
 
    
@@ -469,22 +493,16 @@ consists in indices "ij" follo
 
    This is an output if mle=1. But it can be
 used as an input to get the various output data files (Health
 expectancies, stationary prevalence etc.) and figures without
-rerunning the rather long maximisation phase (mle=0). 
    
-
-
    Each line starts with indices "ijk" followed by the
-covariances between aij and bij: 
    
-
+rerunning the rather long maximisation phase (mle=0). 
    
+Each line starts with indices "ijk" followed by the
+covariances between aij and bij:
    
 
        121 Var(a12) 
    122 Cov(b12,a12)  Var(b12) 
           ...
    232 Cov(b23,a12)  Cov(b23,b12) ... Var (b23) 
    
-
 
      
     
    • If mle=1 you can enter zeros. 
      • 
-
    
-
-
    
     
    # Covariance matrix
 121 0.
 122 0. 0.
@@ -494,12 +512,8 @@ covariances between aij and bij: 
    
 212 0. 0. 0. 0. 0. 0. 
 231 0. 0. 0. 0. 0. 0. 0. 
 232 0. 0. 0. 0. 0. 0. 0. 0.
    
-
    
-
-
      
     
    • If mle=0 you must enter a covariance matrix (usually
-        obtained from an earlier run).
      
-        
      • 
+        obtained from an earlier run). 
 
    
 
 
    Age range for calculation of stationary
@@ -507,20 +521,19 @@ prevalences and health expectanciesagemin=70 agemax=100 bage=50 fage=100

 
-
Once we obtained the estimated parameters, the program is able
+
Once we obtained the estimated parameters, the program is able
 to calculated stationary prevalence, transitions probabilities
 and life expectancies at any age. Choice of age range is useful
 for extrapolation. In our data file, ages varies from age 70 to
 102. It is possible to get extrapolated stationary prevalence by
-age ranging from agemin to agemax. 

+age ranging from agemin to agemax.
 
-
Setting bage=50 (begin age) and fage=100 (final age), makes
+
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!

-
+Program is limited to around 120 for upper age!