Annotation of imach/src/imachrun.log, revision 1.1.1.1
1.1 brouard 1: C:\usr\biais\imach\mytry>c:\usr\biais\imach\bin\imach.exe mypar.txt
2: # MLE 1984-1990
3:
4: #
5:
6: title=MLE datafile=..\data\mydata.txt lastobs=3000 firstpass=1 lastpass=3
7: ftol=1.000000e-008 stepm=24 ncov=2 nlstate=2 ndeath=1 maxwav=4 mle=1 weight=0
8: # Parameters nlstate*nlstate*ncov a12*1 + b12 * age + ...
9:
10: 11 0.000000 0.000000
11: 12 0.000000 0.000000
12: 21 0.000000 0.000000
13: 22 0.000000 0.000000
14: # Scales (for hessian or gradient estimation)
15:
16: 11 0.000000e+000 0.000000e+000
17: 12 0.000000e+000 0.000000e+000
18: 21 0.000000e+000 0.000000e+000
19: 22 0.000000e+000 0.000000e+000
20: #covariance matrix#
21:
22: 121 0.00000e+000
23: 122 0.00000e+000 0.00000e+000
24: 131 0.00000e+000 0.00000e+000 0.00000e+000
25: 132 0.00000e+000 0.00000e+000 0.00000e+000 0.00000e+000
26: 211 0.00000e+000 0.00000e+000 0.00000e+000 0.00000e+000 0.00000e+000
27: 212 0.00000e+000 0.00000e+000 0.00000e+000 0.00000e+000 0.00000e+000 0.00000e+000
28: 231 0.00000e+000 0.00000e+000 0.00000e+000 0.00000e+000 0.00000e+000 0.00000e+000 0.00000e+000
29: 232 0.00000e+000 0.00000e+000 0.00000e+000 0.00000e+000 0.00000e+000 0.00000e+000 0.00000e+000 0.00000e+000
30:
31: Total number of individuals= 2965, Agemin = 70.00, Agemax= 100.92
32:
33: Warning, no any valid information for:126 line=126
34: Warning, no any valid information for:2307 line=2307
35: Delay (in months) between two waves Min=21 Max=51 Mean=24.495826
36:
37: Age 70 1.=230 loss[1]=3.5% 2.=16 loss[2]=12.5% 1.=222 prev[1]=94.1% 2.=14 prev[2]=5.9% 1-1=8 11=200 12=7 13=15 2-1=2 21=6 22=7 23=1
38: Age 71 1.=235 loss[1]=4.7% 2.=7 loss[2]=14.3% 1.=224 prev[1]=97.4% 2.=6 prev[2]=2.6% 1-1=11 11=196 12=6 13=22 2-1=1 21=1 22=4 23=1 33=6
39: Age 72 1.=417 loss[1]=5.8% 2.=24 loss[2]=4.2% 1.=393 prev[1]=94.5% 2.=23 prev[2]=5.5% 1-1=24 11=347 12=18 13=28 2-1=1 21=6 22=14 23=3 33=20
40: Age 73 1.=395 loss[1]=3.3% 2.=21 loss[2]=4.8% 1.=382 prev[1]=95.0% 2.=20 prev[2]=5.0% 1-1=13 11=327 12=15 13=40 2-1=1 21=4 22=12 23=4 33=17
41: Age 74 1.=364 loss[1]=6.3% 2.=37 loss[2]=5.4% 1.=341 prev[1]=90.7% 2.=35 prev[2]=9.3% 1-1=23 11=307 12=11 13=23 2-1=2 21=6 22=23 23=6 3-1=2 33=20
42: Age 75 1.=357 loss[1]=6.4% 2.=23 loss[2]=4.3% 1.=334 prev[1]=93.8% 2.=22 prev[2]=6.2% 1-1=23 11=286 12=18 13=30 2-1=1 21=4 22=13 23=5 3-1=1 33=21
43: Age 76 1.=315 loss[1]=4.4% 2.=29 loss[2]=3.4% 1.=301 prev[1]=91.5% 2.=28 prev[2]=8.5% 1-1=14 11=260 12=20 13=21 2-1=1 21=5 22=18 23=5 3-1=2 33=15
44: Age 77 1.=285 loss[1]=3.5% 2.=46 loss[2]=0.0% 1.=275 prev[1]=85.7% 2.=46 prev[2]=14.3% 1-1=10 11=224 12=19 13=32 21=6 22=30 23=10 3-1=1 33=10
45: Age 78 1.=277 loss[1]=3.2% 2.=41 loss[2]=2.4% 1.=268 prev[1]=87.0% 2.=40 prev[2]=13.0% 1-1=9 11=216 12=28 13=24 2-1=1 21=4 22=28 23=8 3-1=1 33=22
46: Age 79 1.=235 loss[1]=4.7% 2.=47 loss[2]=8.5% 1.=224 prev[1]=83.9% 2.=43 prev[2]=16.1% 1-1=11 11=184 12=20 13=20 2-1=4 21=5 22=28 23=10 33=12
47: Age 80 1.=226 loss[1]=5.8% 2.=40 loss[2]=5.0% 1.=213 prev[1]=84.9% 2.=38 prev[2]=15.1% 1-1=13 11=182 12=16 13=15 2-1=2 21=7 22=25 23=6 33=12
48: Age 81 1.=204 loss[1]=7.4% 2.=40 loss[2]=7.5% 1.=189 prev[1]=83.6% 2.=37 prev[2]=16.4% 1-1=15 11=144 12=25 13=20 2-1=3 21=4 22=16 23=17 3-1=2 33=11
49: Age 82 1.=183 loss[1]=2.7% 2.=40 loss[2]=2.5% 1.=178 prev[1]=82.0% 2.=39 prev[2]=18.0% 1-1=5 11=139 12=21 13=18 2-1=1 21=7 22=19 23=13 3-1=1 33=15
50: Age 83 1.=149 loss[1]=3.4% 2.=35 loss[2]=2.9% 1.=144 prev[1]=80.9% 2.=34 prev[2]=19.1% 1-1=5 11=109 12=15 13=20 2-1=1 21=6 22=20 23=8 3-1=1 33=12
51: Age 84 1.=133 loss[1]=3.0% 2.=53 loss[2]=9.4% 1.=129 prev[1]=72.9% 2.=48 prev[2]=27.1% 1-1=4 11=97 12=18 13=14 2-1=5 21=10 22=23 23=15 3-1=2 33=14
52: Age 85 1.=109 loss[1]=2.8% 2.=34 loss[2]=0.0% 1.=106 prev[1]=75.7% 2.=34 prev[2]=24.3% 1-1=3 11=69 12=17 13=20 21=7 22=19 23=8 3-1=2 33=15
53: Age 86 1.=95 loss[1]=1.1% 2.=43 loss[2]=2.3% 1.=94 prev[1]=69.1% 2.=42 prev[2]=30.9% 1-1=1 11=59 12=17 13=18 2-1=1 21=5 22=22 23=15 3-1=2 33=9
54: Age 87 1.=78 loss[1]=2.6% 2.=44 loss[2]=6.8% 1.=76 prev[1]=65.0% 2.=41 prev[2]=35.0% 1-1=2 11=52 12=14 13=10 2-1=3 21=9 22=20 23=12 33=15
55: Age 88 1.=67 loss[1]=0.0% 2.=32 loss[2]=3.1% 1.=67 prev[1]=68.4% 2.=31 prev[2]=31.6% 11=37 12=19 13=11 2-1=1 21=5 22=20 23=6 33=9
56: Age 89 1.=55 loss[1]=5.5% 2.=31 loss[2]=3.2% 1.=52 prev[1]=63.4% 2.=30 prev[2]=36.6% 1-1=3 11=26 12=14 13=12 2-1=1 21=8 22=17 23=5 3-1=1 33=9
57: Age 90 1.=44 loss[1]=11.4% 2.=27 loss[2]=3.7% 1.=39 prev[1]=60.0% 2.=26 prev[2]=40.0% 1-1=5 11=22 12=9 13=8 2-1=1 21=2 22=14 23=10 3-1=1 33=9
58: Age 91 1.=31 loss[1]=3.2% 2.=23 loss[2]=13.0% 1.=30 prev[1]=60.0% 2.=20 prev[2]=40.0% 1-1=1 11=17 12=7 13=6 2-1=3 21=1 22=11 23=8 33=3
59: Age 92 1.=23 loss[1]=0.0% 2.=21 loss[2]=4.8% 1.=23 prev[1]=53.5% 2.=20 prev[2]=46.5% 11=9 12=8 13=6 2-1=1 21=1 22=13 23=6 33=8
60: Age 93 1.=8 loss[1]=12.5% 2.=19 loss[2]=5.3% 1.=7 prev[1]=28.0% 2.=18 prev[2]=72.0% 1-1=1 11=4 12=3 2-1=1 21=1 22=11 23=6 3-1=1 33=9
61: Age 94 1.=11 loss[1]=9.1% 2.=14 loss[2]=0.0% 1.=10 prev[1]=41.7% 2.=14 prev[2]=58.3% 1-1=1 11=3 12=6 13=1 21=2 22=7 23=5 33=5
62: Age 95 1.=4 loss[1]=0.0% 2.=6 loss[2]=16.7% 1.=4 prev[1]=44.4% 2.=5 prev[2]=55.6% 11=2 12=1 13=1 2-1=1 21=2 23=3 33=2
63: Age 96 1.=4 loss[1]=25.0% 2.=6 loss[2]=0.0% 1.=3 prev[1]=33.3% 2.=6 prev[2]=66.7% 1-1=1 11=1 12=1 13=1 21=2 22=2 23=2 33=2
64: Age 97 1.=0 loss[1]=NaNQ% 2.=1 loss[2]=0.0% 1.=0 prev[1]=0.0% 2.=1 prev[2]=100.0% 23=1
65: Age 98 1.=0 loss[1]=NaNQ% 2.=0 loss[2]=NaNQ% 1.=0 prev[1]=NaNQ% 2.=0 prev[2]=NaNQ%
66: Age 99 1.=0 loss[1]=NaNQ% 2.=0 loss[2]=NaNQ% 1.=0 prev[1]=NaNQ% 2.=0 prev[2]=NaNQ%
67: Age 100 1.=0 loss[1]=NaNQ% 2.=0 loss[2]=NaNQ% 1.=0 prev[1]=NaNQ% 2.=0 prev[2]=NaNQ%
68: Age 101 1.=0 loss[1]=NaNQ% 2.=0 loss[2]=NaNQ% 1.=0 prev[1]=NaNQ% 2.=0 prev[2]=NaNQ%
69: Age 102 1.=0 loss[1]=NaNQ% 2.=0 loss[2]=NaNQ% 1.=0 prev[1]=NaNQ% 2.=0 prev[2]=NaNQ% -1-1=18 -11=188 -12=30 -13=38
70: Total 1.=4534 loss[1]=4.5% 2.=800 loss[2]=4.9% 1.=4328 prev[1]=85.0% 2.=761 prev[2]=15.0% -1-1=18 -11=188 -12=30 -13=38 1-1=206 11=3519 12=373 13=436 2-1=39 21=126 22=436 23=199 3-1=20 33=302
71: Powell
72:
73: Powell iter=1 -2*LL=11531.405658264877 1 0.000000000000 2 0.000000000000 3 0.000000000000 4 0.000000000000 5 0.000000000000 6 0.000000000000 7 0.000000000000 8 0.000000000000
74: 1..........2.................3..........4.................5.........6................7........8...............
75: Powell iter=2 -2*LL=7189.342479571923 1 -1.618034000000 2 0.000243057112 3 -1.618034000000 4 -0.004756278029 5 -0.914196809953 6 0.000204141812 7 -0.703784678359 8 0.000288223851
76: 1........2.................3.........4..................5.........6................7.........8........................
77: Powell iter=3 -2*LL=6999.869835328977 1 -2.331359020254 2 0.000958881180 3 -1.729014288043 4 -0.004451102909 5 -1.234885928233 6 0.000031918936 7 -0.813603344320 8 0.000604615223
78: 1...............2.................3.........4..................5.........6................7.........8..............
79: Powell iter=4 -2*LL=6996.402685207441 1 -2.478603679819 2 0.002320473825 3 -1.774318202733 4 -0.003932303071 5 -1.267360655852 6 -0.000258419142 7 -0.872906870738 8 0.000688621829
80: 1...............2.......3.........4..................5.........6................7.........8..............
81: Powell iter=5 -2*LL=6991.472623520970 1 -2.754862801334 2 0.004818289656 3 -1.909104599415 4 -0.002278897004 5 -1.199219023235 6 -0.001178497892 7 -1.022698614160 8 0.002118851608
82: 1...............2......3.........4..................5.........6................7.........8.............
83: Powell iter=6 -2*LL=6962.896830060558 1 -4.241651416576 2 0.021100764612 3 -2.782715905961 4 0.008746142517 5 -0.661703877114 6 -0.007293414952 7 -2.002149669075 8 0.013161578084
84: 1...............2........3.........4.................5.........6................7.........8..............
85: Powell iter=7 -2*LL=6856.515462890279 1 -9.104604463847 2 0.080963849064 3 -5.984275969066 4 0.049613285500 5 1.350545085198 6 -0.029870618347 7 -5.621492754966 8 0.056727766830
86: 1...............2......3.........4.................5.........6................7.........8.......
87: Powell iter=8 -2*LL=6807.299976649243 1 -9.679671626302 2 0.090172151940 3 -6.501900056622 4 0.056165638832 5 1.617555411708 6 -0.033345210651 7 -6.285041940821 8 0.065876908728
88: 1...............2.......3.........4.................5.........6................7.........8..............
89: Powell iter=9 -2*LL=6765.223979608776 1 -10.944967760181 2 0.110078352295 3 -7.599918978249 4 0.070265582553 5 2.119367547003 6 -0.040632636562 7 -7.697984528675 8 0.082717269362
90: 1...............2........3...........4.................5.........6................7..........8...........
91: Powell iter=10 -2*LL=6764.942925283532 1 -10.995045251042 2 0.110672992239 3 -7.631384326082 4 0.070639289631 5 2.107918420012 6 -0.040631357619 7 -7.717680431349 8 0.082842084637
92: 1...............2........3.........4.................5.........6................7..........8..................
93: Powell iter=11 -2*LL=6758.126928754780 1 -12.649915944424 2 0.130679851428 3 -8.532294367391 4 0.081933260496 5 1.269832277292 6 -0.031153280246 7 -7.541104999513 8 0.080670456845
94: 1........2........3.........4.................5.........6................7.........8.......
95: Powell iter=12 -2*LL=6751.739055814885 1 -12.595143083360 2 0.130107012348 3 -8.389513656808 4 0.080162752729 5 0.471603387369 6 -0.021009924396 7 -6.696076157144 8 0.070744440420
96: 1........2........3.........4.................5.........6................7.........8.............
97: Powell iter=13 -2*LL=6746.696865848237 1 -12.517990378064 2 0.129617715142 3 -8.146403583528 4 0.077181406366 5 -1.054592628099 6 -0.001907385381 7 -5.139371065631 8 0.052382230359
98: 1........2...........3.........4.................5.........6................7.........8........
99: Powell iter=14 -2*LL=6746.679040245447 1 -12.521586079436 2 0.129677462699 3 -8.144304958850 4 0.077159986974 5 -1.065387233819 6 -0.001856671774 7 -5.138397462232 8 0.052326966712
100: 1.................2........3.........4.................5.........6................7........8..............
101: Powell iter=15 -2*LL=6746.584135522181 1 -12.394946753976 2 0.128077576911 3 -8.009745762465 4 0.075460831413 5 -1.045986257644 6 -0.002324513261 7 -5.086087706239 8 0.051646596279
102: 1........2........3.........4.........5..........6................7.........8............
103: Powell iter=16 -2*LL=6746.409688111119 1 -12.337298733327 2 0.127349601718 3 -7.867238995448 4 0.073740146981 5 -1.028502684303 6 -0.002611889339 7 -5.005904152361 8 0.050670578982
104: 1........2........3.........4........5.......6................7.........8..........
105: Powell iter=17 -2*LL=6745.813273626674 1 -12.426249555087 2 0.128447083617 3 -7.474039244737 4 0.068916832631 5 -0.999114433414 6 -0.003020647652 7 -4.804607433487 8 0.048253597224
106: 1........2........3.........4........5.....6................7.........8......
107: Powell iter=18 -2*LL=6745.117193568432 1 -12.807569088249 2 0.133202806969 3 -7.467620316340 4 0.068773114112 5 -0.898426027976 6 -0.004024137111 7 -4.893865605332 8 0.049367546860
108: 1........2........3.........4........5......6................7.........8.......
109: Powell iter=19 -2*LL=6744.983460451949 1 -12.894722466826 2 0.134251952454 3 -7.428738535044 4 0.068222705074 5 -0.767927408730 6 -0.005529901364 7 -4.985241944126 8 0.050472191885
110: 1........2........3.........4........5.......6................7.........8........
111: Powell iter=20 -2*LL=6744.959681724528 1 -12.935337614690 2 0.134745057934 3 -7.413886221380 4 0.068011967304 5 -0.715509808093 6 -0.006132925395 7 -5.022473778604 8 0.050922685102
112: 1............2.........3........4........5.......6................7..........8........
113: Powell iter=21 -2*LL=6744.955290044597 1 -12.952671908339 2 0.134955278220 3 -7.407513503114 4 0.067921460677 5 -0.693105554462 6 -0.006390806673 7 -5.038208232247 8 0.051113073984
114: 1...............2..............3.............4...........5.......6................7............8........
115: Powell iter=22 -2*LL=6744.954431665786 1 -12.960110073059 2 0.135045475315 3 -7.404811455032 4 0.067883051952 5 -0.683487247450 6 -0.006501558124 7 -5.044834063352 8 0.051193262301
116: 1...............2............3...........4...............5.......6.........................7.............8.........................
117: Powell iter=23 -2*LL=6744.954108371555 1 -12.967632334283 2 0.135136681033 3 -7.402109728262 4 0.067844593326 5 -0.673601538129 6 -0.006615504377 7 -5.051341616718 8 0.051272038506
118: 1..............2...........3..............4...........5..........6................7...........8.........
119: #Number of iterations = 23, -2 Log likelihood = 6744.954042573691
120: # Parameters
121: 12 -12.966061 0.135117
122: 13 -7.401109 0.067831
123: 21 -0.672648 -0.006627
124: 23 -5.051297 0.051271
125:
126: Calculation of the hessian matrix. Wait...
127: 12345678.12.13.14.15.16.17.18.23.24.25.26.27.28.34.35.36.37.38.45.46.47.48.56.57.58.67.68.78
128:
129: Inverting the hessian to get the covariance matrix. Wait...
130:
131: #Hessian matrix#
132: 3.344e+002 2.708e+004 -4.586e+001 -3.806e+003 -1.577e+000 -1.313e+002 3.914e-001 3.166e+001
133: 2.708e+004 2.204e+006 -3.805e+003 -3.174e+005 -1.303e+002 -1.091e+004 2.967e+001 2.399e+003
134: -4.586e+001 -3.805e+003 4.044e+002 3.197e+004 2.431e-002 1.995e+000 1.783e-001 1.486e+001
135: -3.806e+003 -3.174e+005 3.197e+004 2.541e+006 2.436e+000 2.051e+002 1.483e+001 1.244e+003
136: -1.577e+000 -1.303e+002 2.431e-002 2.436e+000 1.093e+002 8.979e+003 -3.402e+001 -2.843e+003
137: -1.313e+002 -1.091e+004 1.995e+000 2.051e+002 8.979e+003 7.420e+005 -2.842e+003 -2.388e+005
138: 3.914e-001 2.967e+001 1.783e-001 1.483e+001 -3.402e+001 -2.842e+003 1.494e+002 1.251e+004
139: 3.166e+001 2.399e+003 1.486e+001 1.244e+003 -2.843e+003 -2.388e+005 1.251e+004 1.053e+006
140: # Scales
141: 12 1.00000e-004 1.00000e-006
142: 13 1.00000e-004 1.00000e-006
143: 21 1.00000e-003 1.00000e-005
144: 23 1.00000e-004 1.00000e-005
145: # Covariance
146: 1 5.90661e-001
147: 2 -7.26732e-003 8.98810e-005
148: 3 8.80177e-002 -1.12706e-003 5.15824e-001
149: 4 -1.13082e-003 1.45267e-005 -6.50070e-003 8.23270e-005
150: 5 9.31265e-003 -1.16106e-004 6.00210e-004 -8.04151e-006 1.75753e+000
151: 6 -1.15664e-004 1.44850e-006 -7.79995e-006 1.04770e-007 -2.12929e-002 2.59422e-004
152: 7 1.35103e-003 -1.75392e-005 -6.38237e-004 7.85424e-006 4.02601e-001 -4.86776e-003 1.32682e+000
153: 8 -1.82421e-005 2.35811e-007 7.75503e-006 -9.58687e-008 -4.86589e-003 5.91641e-005 -1.57767e-002 1.88622e-004
154: # agemin agemax for lifexpectancy, bage fage (if mle==0 ie no data nor Max likelihood).
155:
156: # Fin de la lecture
157:
158: agemin=70 agemax=100 bage=50 fage=100
159: Computing prevalence limit: result on file 'plrmypar.txt'
160: Computing pij: result on file 'pijrmypar.txt'
161: Computing Health Expectancies: result on file 'ermypar.txt'
162: Computing Variance-covariance of DFLEs: file 'vrmypar.txt'
163: Computing Total LEs with variances: file 'trmypar.txt'
164: Computing Variance-covariance of Prevalence limit: file 'vplrmypar.txt'
165: End of Imach
166:
167: C:\usr\biais\imach\mytry>pause
168: Appuyez sur une touche pour continuer . . .
169:
170:
171: C:\usr\biais\imach\mytry>ECHO Press a key to get graphs
172: Press a key to get graphs
173:
174: C:\usr\biais\imach\mytry>c:\usr\biais\imach\bin\wgnuplot.exe.lnk graph.gp
175:
176: C:\usr\biais\imach\mytry>pause
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>
![[BACK]](/icons/small/back.gif){kind=icon}
Return to imachrun.log CVS log ![[TXT]](/icons/small/text.gif){kind=icon}
![[DIR]](/icons/small/dir.gif){kind=icon}
Up to [Local Repository] / imach / src