Annotation of imach/src/imachrun.log, revision 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>