--- imach/html/doc/rbiaspar.imach 2004/06/16 21:43:53 1.1 +++ imach/html/doc/rbiaspar.imach 2022/09/11 07:59:43 1.6 @@ -1,41 +1,47 @@ -#Imach version 0.07a, May 2004, INED-EUROREVES -#Number of iterations = 34, -2 Log likelihood = 46542.397795763849 -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.245160 0.092357 -13 -10.672078 0.060971 -21 -2.645815 -0.022320 -23 -4.773208 0.007857 -# 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.12641e-001 -122 -1.35168e-003 1.63147e-005 -131 -6.15783e-002 7.47501e-004 3.20023e-001 -132 7.30656e-004 -8.95009e-006 -3.93444e-003 4.86922e-005 -211 8.34954e-002 -9.96445e-004 2.19213e-002 -2.86865e-004 4.89982e-001 -212 -9.95687e-004 1.19356e-005 -2.78360e-004 3.64465e-006 -5.80715e-003 6.91834e-005 -231 5.25828e-002 -5.99555e-004 -1.07947e-001 1.23483e-003 -4.52302e-004 7.66983e-006 2.43192e-001 -232 -5.99669e-004 6.88392e-006 1.27573e-003 -1.47509e-005 1.65541e-005 -2.31945e-007 -2.77700e-003 3.18981e-005 +#IMaCh 0.99r24 +#0.99r24 +#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: .