--- imach/html/doc/rbiaspar.imach 2006/02/09 15:20:33 1.3
+++ imach/html/doc/rbiaspar.imach 2019/05/22 23:12:10 1.5
@@ -1,41 +1,47 @@
-#Imach version 0.98b, January 2006, INED-EUROREVES
-#Number of iterations = 34, -2 Log likelihood = 46542.395030424377
-title=1st_example datafile=data1.txt lastobs=8600 firstpass=1 lastpass=4
-ftol=1.000000e-008 stepm=1 ncovcol=2 nlstate=2 ndeath=1 maxwav=4 mle= 0 weight=0
-model=.
-# Parameters nlstate*nlstate*ncov a12*1 + b12 * age + ...
-12 -12.245240 0.092358
-13 -10.671890 0.060969
-21 -2.645345 -0.022325
-23 -4.773317 0.007859
-# Scales (for hessian or gradient estimation)
-12 1.00000e-004 1.00000e-006
-13 1.00000e-004 1.00000e-006
-21 1.00000e-003 1.00000e-006
-23 1.00000e-004 1.00000e-005
-# Covariance matrix
-# 121 Var(a12)
-# 122 Cov(b12,a12) Var(b12)
-# ...
-# 232 Cov(b23,a12) Cov(b23,b12) ... Var (b23)
-#121 Var(a12)
-#122 Cov(b12,a12) Var(b12)
-#131 Cov(a13,a12) Cov(a13,b12) Var(a13)
-#132 Cov(b13,a12) Cov(b13,b12) Cov(b13,a13) Var(b13)
-#211 Cov(a21,a12) Cov(a21,b12) Cov(a21,a13) Cov(a21,b13) Var(a21)
-#212 Cov(b21,a12) Cov(b21,b12) Cov(b21,a13) Cov(b21,b13) Cov(b21,a21) Var(b21)
-#231 Cov(a23,a12) Cov(a23,b12) Cov(a23,a13) Cov(a23,b13) Cov(a23,a21) Cov(a23,b21) Var(a23)
-#232 Cov(b23,a12) Cov(b23,b12) Cov(b23,a13) Cov(b23,b13) Cov(b23,a21) Cov(b23,b21) Cov(b23,a23) Var(b23)
-121 1.11832e-001
-122 -1.34196e-003 1.61978e-005
-131 -6.31135e-002 7.65687e-004 3.32706e-001
-132 7.49702e-004 -9.17551e-006 -4.08921e-003 5.05801e-005
-211 7.93255e-002 -9.47003e-004 2.83678e-002 -3.63380e-004 4.85740e-001
-212 -9.45906e-004 1.13454e-005 -3.55777e-004 4.56374e-006 -5.75682e-003 6.85863e-005
-231 5.35302e-002 -6.11201e-004 -1.16519e-001 1.34050e-003 -7.78502e-003 9.53787e-005 2.47205e-001
-232 -6.10827e-004 7.02090e-006 1.37532e-003 -1.59782e-005 1.00595e-004 -1.23725e-006 -2.82423e-003 3.24537e-005
+#IMaCh 0.99r19
+#0.99r10
+#Number of iterations & function calls = 33 & 6007, -2 Log likelihood = 46617.105944307616
+title=1st_example datafile=data1.txt lastobs=8600 firstpass=1 lastpass=4
+ftol=1.000000e-08 stepm=1 ncovcol=2 nqv=0 ntv=0 nqtv=0 nlstate=2 ndeath=1 maxwav=4 mle=0 weight=0
+model=1+age+
+#Number of iterations & function calls = 2 & 425, -2 Log likelihood = 46617.105944307536
+# Parameters nlstate*nlstate*ncov a12*1 + b12 * age + ...
+12 -12.3175784 0.0931122
+13 -10.1631705 0.0548384
+21 -1.6685685 -0.0339056
+23 -5.6869199 0.0181350
+# Scales (for hessian or gradient estimation)
+12 2.00000e-04 2.00000e-06
+13 3.00000e-04 4.00000e-06
+21 1.00000e-04 1.00000e-06
+23 3.00000e-04 3.00000e-06
+# Covariance matrix
+# 121 Var(a12)
+# 122 Cov(b12,a12) Var(b12)
+# ...
+# 232 Cov(b23,a12) Cov(b23,b12) ... Var (b23)
+#121 Var(a12)
+#122 Cov(b12,a12) Var(b12)
+#131 Cov(a13,a12) Cov(a13,b12) Var(a13)
+#132 Cov(b13,a12) Cov(b13,b12) Cov(b13,a13) Var(b13)
+#211 Cov(a21,a12) Cov(a21,b12) Cov(a21,a13) Cov(a21,b13) Var(a21)
+#212 Cov(b21,a12) Cov(b21,b12) Cov(b21,a13) Cov(b21,b13) Cov(b21,a21) Var(b21)
+#231 Cov(a23,a12) Cov(a23,b12) Cov(a23,a13) Cov(a23,b13) Cov(a23,a21) Cov(a23,b21) Var(a23)
+#232 Cov(b23,a12) Cov(b23,b12) Cov(b23,a13) Cov(b23,b13) Cov(b23,a21) Cov(b23,b21) Cov(b23,a23) Var(b23)
+121 1.0942155e-01
+122 -1.31211e-03 1.5827559e-05
+131 -5.42260e-02 6.63610e-04 3.0218559e-01
+132 6.46897e-04 -7.99273e-06 -3.73518e-03 4.6465922e-05
+211 8.05586e-02 -9.61087e-04 1.88983e-02 -2.43008e-04 4.8553004e-01
+212 -9.60482e-04 1.15119e-05 -2.38550e-04 3.07513e-06 -5.77598e-03 6.9069819e-05
+231 4.89528e-02 -5.53986e-04 -9.43826e-02 1.06862e-03 -1.83095e-03 2.05848e-05 2.5001784e-01
+232 -5.56252e-04 6.33917e-06 1.11795e-03 -1.28139e-05 2.71734e-05 -3.16845e-07 -2.83579e-03 3.2347405e-05
# agemin agemax for life expectancy, bage fage (if mle==0 ie no data nor Max likelihood).
-agemin=70 agemax=95 bage=65 fage=95 estepm=1
-begin-prev-date=1/1/1984 end-prev-date=1/6/1988 mov_average=0
-pop_based=0
-prevforecast=0 starting-proj-date=1/1/2000 final-proj-date=1/1/2000 mobil_average=0
+# agemin agemax for life expectancy, bage fage (if mle==0 ie no data nor Max likelihood).
+agemin=70 agemax=95 bage=70 fage=95 estepm=1 ftolpl=6.000000e-04
+begin-prev-date=1/1/1984 end-prev-date=1/6/1988 mov_average=0
+pop_based=1
+prevforecast=1 starting-proj-date=1/1/1989 final-proj-date=1/1/2000 mobil_average=0
+#backcast=1 starting-back-date=1/1/1990 final-back-date=1/1/1970 mobil_average=1
+# Results model=1+age.
+result: .