--- imach/html/doc/rbiaspar.imach 2005/10/25 11:53:26 1.2 +++ imach/html/doc/rbiaspar.imach 2006/02/09 15:20:33 1.3 @@ -1,41 +1,41 @@ -#Imach version 0.98, September 2005, INED-EUROREVES -#Number of iterations = 34, -2 Log likelihood = 46542.378402526519 -title=1st_example datafile=data1.txt lastobs=8600 firstpass=1 lastpass=4 -ftol=1.000000e-08 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.245504 0.092361 -13 -10.670953 0.060958 -21 -2.645034 -0.022329 -23 -4.775501 0.007884 -# Scales (for hessian or gradient estimation) -12 1.00000e-04 1.00000e-06 -13 1.00000e-04 1.00000e-06 -21 1.00000e-03 1.00000e-06 -23 1.00000e-04 1.00000e-05 -# 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.13401e-01 -122 -1.36090e-03 1.64264e-05 -131 -6.89393e-02 8.35716e-04 3.31264e-01 -132 8.21021e-04 -1.00329e-05 -4.07026e-03 5.03323e-05 -211 7.49121e-02 -8.92959e-04 1.95088e-02 -2.56008e-04 4.75336e-01 -212 -8.93092e-04 1.06987e-05 -2.49893e-04 3.28027e-06 -5.63203e-03 6.70896e-05 -231 5.58947e-02 -6.39542e-04 -1.16708e-01 1.34217e-03 -9.40960e-04 1.38147e-05 2.47512e-01 -232 -6.38539e-04 7.35309e-06 1.37680e-03 -1.59888e-05 2.13084e-05 -2.92279e-07 -2.82744e-03 3.24870e-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 +#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 +# 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