Summary: 0.99r12

diff --git a/html/doc/biaspar.imach b/html/doc/biaspar.imach
index 792d378019f049fef03483aaac5d1d9deee226b2..5367613b0c9c9817e693106159d18ba5dde7fdf8 100644 (file)
--- a/html/doc/biaspar.imach
+++ b/html/doc/biaspar.imach
@@ -1,30 +1,44 @@
-# Imach version 0.98r3, December 2015, INED-EUROREVES \r
-title=1st_example datafile=data1.txt lastobs=8600 firstpass=1 lastpass=4\r
-ftol=1.000000e-008 stepm=1 ncovcol=2 nlstate=2 ndeath=1 maxwav=4 mle=1 weight=0\r
-model=1+age.\r
-12 -12.245240 0.092358 \r
-13 -10.671890 0.060969 \r
-21 -2.645345 -0.022325 \r
-23 -4.773317 0.007859\r
-# Scales\r
-12 0. 0. \r
-13 0. 0. \r
-21 0. 0. \r
-23 0. 0. \r
-#covariance matrix#\r
-121 0.\r
-122 0. 0.\r
-131 0. 0. 0. \r
-132 0. 0. 0. 0. \r
-211 0. 0. 0. 0. 0. \r
-212 0. 0. 0. 0. 0. 0. \r
-231 0. 0. 0. 0. 0. 0. 0. \r
-232 0. 0. 0. 0. 0. 0. 0. 0.\r
-# agemin agemax for life expectancy, bage fage (if mle==0 ie no data nor Max likelihood).\r
-agemin=70 agemax=95 bage=65 fage=95 estepm=1 ftolpl=6e-4\r
-# Observed prevalence period\r
-begin-prev-date=1/1/1984 end-prev-date=1/6/1988 \r
-# Health expectancies computed from stationary prevalence (pop_based=0) or population based (1)\r
-pop_based=1\r
-# Prevalence forecasting\r
-prevforecast=1 starting-proj-date=1/1/1989 final-proj-date=1/1/1990 mov_average=0\r
+#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+.
+# Parameters nlstate*nlstate*ncov a12*1 + b12 * age + ...
+12  -12.3175753    0.0931122 
+13  -10.1631670    0.0548384 
+21   -1.6685685   -0.0339056 
+23   -5.6869201    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.0846319e-01
+122 -1.30062e-03 1.5689741e-05
+131 -5.39098e-02 6.59806e-04 3.0197746e-01
+132 6.43147e-04 -7.94760e-06 -3.73264e-03 4.6435109e-05
+211 7.75092e-02 -9.24613e-04 1.91004e-02 -2.45003e-04 4.7111318e-01
+212 -9.24119e-04 1.10769e-05 -2.40929e-04 3.09854e-06 -5.60404e-03 6.7019135e-05
+231 4.86258e-02 -5.50028e-04 -9.42536e-02 1.06700e-03 -2.57577e-03 2.94283e-05 2.5010261e-01
+232 -5.52529e-04 6.29410e-06 1.11646e-03 -1.27952e-05 3.55428e-05 -4.16204e-07 -2.83674e-03 3.2358065e-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 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.\r
+#result: .\r

diff --git a/html/doc/mypar.imach b/html/doc/mypar.imach
index 06f941bef9b0c3e0da02d6f32d993feafd82d8ff..0fc4e5b7cf0122df38adf5cb208cf13157ced17c 100755 (executable)
--- a/html/doc/mypar.imach
+++ b/html/doc/mypar.imach
@@ -1,7 +1,7 @@
 # MLE 1984-1990\r
 #\r
 title=MLE datafile=../doc/data/mydata.txt lastobs=3000 firstpass=1 lastpass=4\r
-ftol=1e-8 stepm=24 ncovcol=2 nlstate=2 ndeath=1 maxwav=4 mle=1 weight=0\r
+ftol=1e-8 stepm=24 ncovcol=2 nqv=0 ntv=0 nqtv=0 nlstate=2 ndeath=1 maxwav=4 mle=1 weight=0\r
 model=1+age.\r
 #q Parameters nlstate*nlstate*ncov a12*1 + b12 * age + ...\r
 12 0. 0.\r