Annotation of imach/html/doc/biaspar.log, revision 1.2
1.2 ! brouard 1: Log filename:biaspar.log
! 2:
! 3: Imach version 0.98, September 2005, INED-EUROREVES
! 4: $Revision: 1.104 $ $Date: 2005/09/30 16:11:43 $
! 5: Enter the parameter file name:
! 6: pathimach=../../bin/
! 7: pathtot=biaspar.imach
! 8: path=/Users/brouard/bin/html/doc
! 9: optionfile=biaspar.imach
! 10: optionfilext=imach
! 11: optionfilefiname=biaspar
! 12: Local time (at start): Mon Oct 24 14:51:36 2005
1.1 brouard 13: # Imach version 0.97b, June 2004, INED-EUROREVES
1.2 ! brouard 14: title=1st_example datafile=data1.txt lastobs=8600 firstpass=1 lastpass=4
! 15: ftol=1.000000e-08 stepm=1 ncovcol=2 nlstate=2 ndeath=1 maxwav=4 mle=1 weight=0
! 16: model=.
1.1 brouard 17: # Parameters
1.2 ! brouard 18: 11 0.000000 0.000000
! 19: 12 0.000000 0.000000
! 20: 21 0.000000 0.000000
! 21: 22 0.000000 0.000000
1.1 brouard 22: # Scales
1.2 ! brouard 23: 12 0.000000e+00 0.000000e+00
! 24: 13 0.000000e+00 0.000000e+00
! 25: 21 0.000000e+00 0.000000e+00
! 26: 23 0.000000e+00 0.000000e+00
1.1 brouard 27: #covariance matrix#
1.2 ! brouard 28: 121 0.00000e+00
! 29: 122 0.00000e+00 0.00000e+00
! 30: 131 0.00000e+00 0.00000e+00 0.00000e+00
! 31: 132 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
! 32: 211 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
! 33: 212 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
! 34: 231 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
! 35: 232 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
! 36:
! 37: Total number of individuals= 8270, Agemin = 70.00, Agemax= 104.17
! 38:
! 39: Delay (in months) between two waves Min=1 Max=74 Mean=24.040080
! 40:
! 41: Age 70 1.=665 loss[1]=5.1% 2.=0 loss[2]=NaNQ% 1.=631 prev[1]=100.0% 2.=0 prev[2]=0.0% 1-1=34 11=563 12=33 13=35 33=6
! 42: Age 71 1.=663 loss[1]=5.7% 2.=2 loss[2]=0.0% 1.=625 prev[1]=99.7% 2.=2 prev[2]=0.3% 1-1=38 11=530 12=47 13=48 21=1 22=1 3-1=2 33=31
! 43: Age 72 1.=1181 loss[1]=5.6% 2.=35 loss[2]=5.7% 1.=1115 prev[1]=97.1% 2.=33 prev[2]=2.9% 1-1=66 11=969 12=70 13=76 2-1=2 21=9 22=18 23=6 3-1=2 33=87
! 44: Age 73 1.=1119 loss[1]=5.1% 2.=48 loss[2]=8.3% 1.=1062 prev[1]=96.0% 2.=44 prev[2]=4.0% 1-1=57 11=912 12=71 13=79 2-1=4 21=12 22=26 23=6 3-1=12 33=94
! 45: Age 74 1.=1563 loss[1]=4.2% 2.=90 loss[2]=4.4% 1.=1497 prev[1]=94.6% 2.=86 prev[2]=5.4% 1-1=66 11=1277 12=104 13=116 2-1=4 21=30 22=41 23=15 3-1=4 33=125
! 46: Age 75 1.=1463 loss[1]=6.2% 2.=98 loss[2]=7.1% 1.=1372 prev[1]=93.8% 2.=91 prev[2]=6.2% 1-1=91 11=1155 12=98 13=119 2-1=7 21=15 22=55 23=21 3-1=4 33=124
! 47: Age 76 1.=1390 loss[1]=4.5% 2.=112 loss[2]=3.6% 1.=1328 prev[1]=92.5% 2.=108 prev[2]=7.5% 1-1=62 11=1090 12=111 13=127 2-1=4 21=25 22=54 23=29 3-1=4 33=142
! 48: Age 77 1.=1271 loss[1]=5.0% 2.=101 loss[2]=5.9% 1.=1208 prev[1]=92.7% 2.=95 prev[2]=7.3% 1-1=63 11=962 12=131 13=115 2-1=6 21=24 22=52 23=19 3-1=7 33=145
! 49: Age 78 1.=1166 loss[1]=4.9% 2.=122 loss[2]=4.1% 1.=1109 prev[1]=90.5% 2.=117 prev[2]=9.5% 1-1=57 11=894 12=111 13=104 2-1=5 21=21 22=75 23=21 3-1=4 33=146
! 50: Age 79 1.=1045 loss[1]=5.7% 2.=133 loss[2]=5.3% 1.=985 prev[1]=88.7% 2.=126 prev[2]=11.3% 1-1=60 11=773 12=109 13=103 2-1=7 21=31 22=68 23=27 3-1=8 33=136
! 51: Age 80 1.=986 loss[1]=4.1% 2.=128 loss[2]=3.1% 1.=946 prev[1]=88.4% 2.=124 prev[2]=11.6% 1-1=40 11=718 12=140 13=88 2-1=4 21=21 22=71 23=32 3-1=6 33=121
! 52: Age 81 1.=872 loss[1]=4.9% 2.=130 loss[2]=3.8% 1.=829 prev[1]=86.9% 2.=125 prev[2]=13.1% 1-1=43 11=629 12=102 13=98 2-1=5 21=24 22=74 23=27 3-1=5 33=126
! 53: Age 82 1.=803 loss[1]=4.4% 2.=146 loss[2]=4.1% 1.=768 prev[1]=84.6% 2.=140 prev[2]=15.4% 1-1=35 11=546 12=131 13=91 2-1=6 21=15 22=77 23=48 3-1=11 33=100
! 54: Age 83 1.=726 loss[1]=6.5% 2.=116 loss[2]=0.9% 1.=679 prev[1]=85.5% 2.=115 prev[2]=14.5% 1-1=47 11=456 12=124 13=99 2-1=1 21=19 22=71 23=25 3-1=8 33=133
! 55: Age 84 1.=622 loss[1]=5.6% 2.=147 loss[2]=4.1% 1.=587 prev[1]=80.6% 2.=141 prev[2]=19.4% 1-1=35 11=394 12=112 13=81 2-1=6 21=18 22=82 23=41 3-1=5 33=114
! 56: Age 85 1.=526 loss[1]=4.4% 2.=127 loss[2]=6.3% 1.=503 prev[1]=80.9% 2.=119 prev[2]=19.1% 1-1=23 11=315 12=122 13=66 2-1=8 21=13 22=67 23=39 3-1=4 33=137
! 57: Age 86 1.=464 loss[1]=5.2% 2.=130 loss[2]=6.2% 1.=440 prev[1]=78.3% 2.=122 prev[2]=21.7% 1-1=24 11=271 12=103 13=66 2-1=8 21=18 22=70 23=34 3-1=3 33=112
! 58: Age 87 1.=372 loss[1]=4.8% 2.=143 loss[2]=2.8% 1.=354 prev[1]=71.8% 2.=139 prev[2]=28.2% 1-1=18 11=219 12=79 13=56 2-1=4 21=15 22=71 23=53 3-1=1 33=81
! 59: Age 88 1.=283 loss[1]=3.9% 2.=164 loss[2]=6.7% 1.=272 prev[1]=64.0% 2.=153 prev[2]=36.0% 1-1=11 11=172 12=63 13=37 2-1=11 21=25 22=84 23=44 3-1=8 33=102
! 60: Age 89 1.=227 loss[1]=3.5% 2.=152 loss[2]=4.6% 1.=219 prev[1]=60.2% 2.=145 prev[2]=39.8% 1-1=8 11=135 12=45 13=39 2-1=7 21=21 22=73 23=51 3-1=2 33=77
! 61: Age 90 1.=181 loss[1]=3.9% 2.=143 loss[2]=4.9% 1.=174 prev[1]=56.1% 2.=136 prev[2]=43.9% 1-1=7 11=104 12=41 13=29 2-1=7 21=18 22=82 23=36 3-1=1 33=75
! 62: Age 91 1.=146 loss[1]=4.1% 2.=109 loss[2]=8.3% 1.=140 prev[1]=58.3% 2.=100 prev[2]=41.7% 1-1=6 11=78 12=36 13=26 2-1=9 21=11 22=51 23=38 3-1=3 33=77
! 63: Age 92 1.=114 loss[1]=6.1% 2.=101 loss[2]=2.0% 1.=107 prev[1]=51.9% 2.=99 prev[2]=48.1% 1-1=7 11=47 12=43 13=17 2-1=2 21=11 22=52 23=36 3-1=2 33=61
! 64: Age 93 1.=78 loss[1]=3.8% 2.=89 loss[2]=12.4% 1.=75 prev[1]=49.0% 2.=78 prev[2]=51.0% 1-1=3 11=36 12=25 13=14 2-1=11 21=10 22=46 23=22 3-1=2 33=62
! 65: Age 94 1.=51 loss[1]=3.9% 2.=71 loss[2]=1.4% 1.=49 prev[1]=41.2% 2.=70 prev[2]=58.8% 1-1=2 11=24 12=8 13=17 2-1=1 21=4 22=39 23=27 3-1=1 33=37
! 66: Age 95 1.=38 loss[1]=5.3% 2.=57 loss[2]=5.3% 1.=36 prev[1]=40.0% 2.=54 prev[2]=60.0% 1-1=2 11=16 12=13 13=7 2-1=3 21=4 22=31 23=19 3-1=1 33=39
! 67: Age 96 1.=25 loss[1]=4.0% 2.=43 loss[2]=4.7% 1.=24 prev[1]=36.9% 2.=41 prev[2]=63.1% 1-1=1 11=11 12=9 13=4 2-1=2 21=2 22=22 23=17 3-1=2 33=25
! 68: Age 97 1.=11 loss[1]=0.0% 2.=36 loss[2]=8.3% 1.=11 prev[1]=25.0% 2.=33 prev[2]=75.0% 11=4 12=7 2-1=3 21=2 22=22 23=9 33=23
! 69: Age 98 1.=5 loss[1]=0.0% 2.=21 loss[2]=0.0% 1.=5 prev[1]=19.2% 2.=21 prev[2]=80.8% 11=1 12=4 21=1 22=13 23=7 3-1=1 33=17
! 70: Age 99 1.=6 loss[1]=0.0% 2.=14 loss[2]=0.0% 1.=6 prev[1]=30.0% 2.=14 prev[2]=70.0% 11=3 12=1 13=2 21=2 22=2 23=10 33=6
! 71: Age 100 1.=2 loss[1]=0.0% 2.=9 loss[2]=0.0% 1.=2 prev[1]=18.2% 2.=9 prev[2]=81.8% 11=1 12=1 21=2 22=5 23=2 33=2
! 72: Age 101 1.=2 loss[1]=50.0% 2.=0 loss[2]=NaNQ% 1.=1 prev[1]=100.0% 2.=0 prev[2]=0.0% 1-1=1 12=1 3-1=1 33=3
! 73: Age 102 1.=2 loss[1]=0.0% 2.=1 loss[2]=0.0% 1.=2 prev[1]=66.7% 2.=1 prev[2]=33.3% 12=2 22=1 33=1
! 74: Age 103 1.=0 loss[1]=NaNQ% 2.=0 loss[2]=NaNQ% 1.=0 prev[1]=NaNQ% 2.=0 prev[2]=NaNQ%
! 75: Age 104 1.=0 loss[1]=NaNQ% 2.=0 loss[2]=NaNQ% 1.=0 prev[1]=NaNQ% 2.=0 prev[2]=NaNQ%
! 76: Age 105 1.=0 loss[1]=NaNQ% 2.=0 loss[2]=NaNQ% 1.=0 prev[1]=NaNQ% 2.=0 prev[2]=NaNQ%
! 77: Age 106 1.=0 loss[1]=NaNQ% 2.=0 loss[2]=NaNQ% 1.=0 prev[1]=NaNQ% 2.=0 prev[2]=NaNQ% -1-1=69 -11=718 -12=195 -13=261
! 78: Total 1.=18068 loss[1]=5.0% 2.=2818 loss[2]=4.9% 1.=17161 prev[1]=86.5% 2.=2681 prev[2]=13.5% -1-1=69 -11=718 -12=195 -13=261 1-1=907 11=13305 12=2097 13=1759 2-1=137 21=424 22=1496 23=761 3-1=114 33=2567
! 79: Powell
! 80: 1 0.000000000000 2 0.000000000000 3 0.000000000000 4 0.000000000000 5 0.000000000000 6 0.000000000000 7 0.000000000000 8 0.000000000000
! 81:
! 82: Considering the time needed for this last iteration #1: 5 seconds,
! 83: - if your program needs 10 iterations to converge, convergence will be
! 84: reached in 0 day(s) 0 hour(s) 0 minute(s) 45 second(s) i.e.
! 85: on Mon Oct 24 14:52:26 2005 (current time is Mon Oct 24 14:51:41 2005);
! 86: - if your program needs 20 iterations to converge, convergence will be
! 87: reached in 0 day(s) 0 hour(s) 1 minute(s) 35 second(s) i.e.
! 88: on Mon Oct 24 14:53:16 2005 (current time is Mon Oct 24 14:51:41 2005);
! 89: - if your program needs 30 iterations to converge, convergence will be
! 90: reached in 0 day(s) 0 hour(s) 2 minute(s) 25 second(s) i.e.
! 91: on Mon Oct 24 14:54:06 2005 (current time is Mon Oct 24 14:51:41 2005);
! 92: 1........2..............3.........4.......................5.....6............7.........8.............................. 1 8.981830530002 2 0.016741746069 3 -0.434889243858 4 0.000506810742 5 3.199437430882 6 -0.000133684844 7 -0.817356206705 8 0.000414700052
! 93:
! 94: Considering the time needed for this last iteration #2: 167 seconds,
! 95: - if your program needs 10 iterations to converge, convergence will be
! 96: reached in 0 day(s) 0 hour(s) 22 minute(s) 16 second(s) i.e.
! 97: on Mon Oct 24 15:16:44 2005 (current time is Mon Oct 24 14:54:28 2005);
! 98: - if your program needs 20 iterations to converge, convergence will be
! 99: reached in 0 day(s) 0 hour(s) 50 minute(s) 6 second(s) i.e.
! 100: on Mon Oct 24 15:44:34 2005 (current time is Mon Oct 24 14:54:28 2005);
! 101: - if your program needs 30 iterations to converge, convergence will be
! 102: reached in 0 day(s) 1 hour(s) 17 minute(s) 56 second(s) i.e.
! 103: on Mon Oct 24 16:12:24 2005 (current time is Mon Oct 24 14:54:28 2005);
! 104: 1...........2.....................3..........4..................5.................6................7..........8................. 1 -2.281804650756 2 0.011363116543 3 -3.509368523070 4 -0.009363669124 5 2.547946952776 6 -0.006890085832 7 -1.246147647315 8 -0.012650266592
! 105:
! 106: Considering the time needed for this last iteration #3: 168 seconds,
! 107: - if your program needs 10 iterations to converge, convergence will be
! 108: reached in 0 day(s) 0 hour(s) 19 minute(s) 36 second(s) i.e.
! 109: on Mon Oct 24 15:16:52 2005 (current time is Mon Oct 24 14:57:16 2005);
! 110: - if your program needs 20 iterations to converge, convergence will be
! 111: reached in 0 day(s) 0 hour(s) 47 minute(s) 36 second(s) i.e.
! 112: on Mon Oct 24 15:44:52 2005 (current time is Mon Oct 24 14:57:16 2005);
! 113: - if your program needs 30 iterations to converge, convergence will be
! 114: reached in 0 day(s) 1 hour(s) 15 minute(s) 36 second(s) i.e.
! 115: on Mon Oct 24 16:12:52 2005 (current time is Mon Oct 24 14:57:16 2005);
! 116: 1........2.................3......4....................5........6....................7..........8........... 1 -2.181985096129 2 0.011465880658 3 -4.732211491614 4 -0.008150559262 5 2.446857218233 6 -0.008350202510 7 -1.338377277831 8 -0.012700875334
! 117: 1........2................3..........4...................5.........6...................7..........8....................... 1 -3.070133451540 2 0.018029818494 3 -6.179866175818 4 0.001935898306 5 2.212009143103 6 -0.027638638041 7 -2.308991353528 8 -0.013146236293
! 118: 1........2...............3..........4....................5........6..........7..........8........... 1 -3.626502091197 2 0.020253035211 3 -6.684342507080 4 0.006057014157 5 2.123782483735 6 -0.032645573113 7 -2.451363857789 8 -0.013247681665
! 119: 1.........2................3..........4....................5........6..........7..........8................ 1 -7.192534317905 2 0.034624192693 3 -10.542205375135 4 0.036666250158 5 1.541870117892 6 -0.063020402356 7 -2.948321308915 8 -0.013800448747
! 120: 1.........2...........3.........4...................5.........6............7..........8......... 1 -7.695451703104 2 0.037194072253 3 -9.230251380079 4 0.039357187891 5 1.401560774610 6 -0.065175085233 7 -2.866616489626 8 -0.013803976993
! 121: 1.........2.............3..........4...................5.........6...........7...........8......... 1 -7.940212354544 2 0.039881323201 3 -9.188932954744 4 0.041398021268 5 1.337148135578 6 -0.067668349287 7 -2.901457910987 8 -0.013846848839
! 122: 1.........2...............3...........4...................5.........6...........7...........8.................. 1 -8.176431117205 2 0.042833923530 3 -9.304312328468 4 0.043720071441 5 1.276877306765 6 -0.070125167890 7 -2.891118730399 8 -0.013891797083
! 123: 1.........2.........3..........4...................5..........6...........7..........8........ 1 -8.179934732156 2 0.042989046121 3 -9.366297995893 4 0.044103933684 5 1.279115322988 6 -0.070236292085 7 -2.878917655261 8 -0.013898162887
! 124: 1..........2.........3..........4...................5..........6...........7..........8................ 1 -8.254585369807 2 0.044560894850 3 -9.868936732048 4 0.048207926200 5 1.293598356293 6 -0.070738582582 7 -2.797523985671 8 -0.013960095547
! 125: 1..........2.........3........4....................5..........6...........7.........8........ 1 -8.301539191723 2 0.045052305123 3 -10.022804846924 4 0.050130364363 5 1.294397184532 6 -0.070701522071 7 -2.801520599321 8 -0.013984314301
! 126: 1..........2.........3........4....................5..........6...........7.........8................ 1 -9.366193466838 2 0.057625368748 3 -13.302743404812 4 0.092879588144 5 1.298194338761 6 -0.069458719670 7 -2.845793360462 8 -0.014517392773
! 127: 1.........2.........3........4...................5..........6....................7........8......... 1 -9.460968984193 2 0.058999781798 3 -13.723675221324 4 0.097150765005 5 1.308211839961 6 -0.069448653825 7 -2.836187353406 8 -0.014578542733
! 128: 1..........2.........3........4...................5..................6...........7.........8.......... 1 -9.558132338168 2 0.060212099576 3 -13.942989603808 4 0.099311940480 5 1.302080125134 6 -0.069569168801 7 -2.798255305774 8 -0.014606029284
! 129: 1......................2.........3........4...................5...........6......................7.........8....... 1 -9.572161154893 2 0.060361245477 3 -13.929855173962 4 0.099310593326 5 1.298209054802 6 -0.069455399052 7 -2.795021923613 8 -0.014604875644
! 130: 1.........2.........3........4...................5..............6...........7.........8............. 1 -9.521658848729 2 0.059772467720 3 -13.378117941017 4 0.092719757070 5 1.273997949530 6 -0.068987502478 7 -2.776883693013 8 -0.014526513051
! 131: 1.........2.........3........4...................5..........6...........7........8....... 1 -9.561633452062 2 0.060333173985 3 -13.313444344667 4 0.091835856342 5 1.263732018431 6 -0.068847034820 7 -2.776256237312 8 -0.014524728372
! 132: 1.........2.........3........4...................5..........6..........7........8............... 1 -11.211177020816 2 0.080500005299 3 -11.537258013108 4 0.070983535362 5 0.867292040233 6 -0.064099822721 7 -2.823227325066 8 -0.014537371296
! 133: 1.........2.........3........4...................5..........6...........7........8........ 1 -11.325506728779 2 0.081619385497 3 -11.381316058273 4 0.068836396199 5 0.837536179182 6 -0.063930299171 7 -2.827086277323 8 -0.014528843884
! 134: 1.........2.........3........4......................5..........6...........7........8........... 1 -11.341676688719 2 0.081802878533 3 -11.416589557874 4 0.069266621366 5 0.836897464135 6 -0.063871223031 7 -2.826483099116 8 -0.014533771831
! 135: 1.....................2............3.........4....................5.................6...........7.........8........................ 1 -11.341583321858 2 0.081813100925 3 -11.417800327848 4 0.069285874759 5 0.836665337138 6 -0.063830701725 7 -2.824804546159 8 -0.014530472468
! 136: 1.................2.................3...........4........................5.............6..........7........8...... 1 -11.343601882043 2 0.081830552688 3 -11.402331269733 4 0.069097487762 5 0.835136563121 6 -0.063812998546 7 -2.824904482410 8 -0.014527159992
! 137: 1..............2.........3.........4...................5............6....................7..........8............. 1 -11.344108280017 2 0.081798969097 3 -11.289229702040 4 0.067698422410 5 0.821725110822 6 -0.063651576307 7 -2.828095083276 8 -0.014467808496
! 138: 1..........2.........3........4...................5..........6........7..........8............ 1 -11.346884442242 2 0.081798944160 3 -11.190636099860 4 0.066481588797 5 0.799901616320 6 -0.063393443465 7 -2.836778391487 8 -0.014350132999
! 139: 1.........2.........3........4...................5..........6........7...........8............ 1 -11.374550257257 2 0.082018917627 3 -10.813465698425 4 0.061852349060 5 0.649704022353 6 -0.061545565813 7 -2.904867373668 8 -0.013465068949
! 140: 1.........2.........3........4...................5..........6........7..........8............. 1 -11.497299557475 2 0.083382209466 3 -10.214520395463 4 0.054601478423 5 0.106321936126 6 -0.054800998817 7 -3.176294913042 8 -0.010075948370
! 141: 1.........2.........3........4...................5..........6........7........8............... 1 -11.841632967445 2 0.087600665563 3 -9.764976982488 4 0.049419193834 5 -1.249466269151 6 -0.038183184473 7 -3.898766321772 8 -0.001377882469
! 142: 1.........2.........3.......4....................5..........6........7.......8........ 1 -11.994458882977 2 0.089578157478 3 -10.199298537104 4 0.054876691393 5 -1.732627878219 6 -0.032656820997 7 -4.194520508882 8 0.001844755037
! 143: 1.........2.........3........4...................5..........6........7.......8......... 1 -12.258447292779 2 0.092595784728 3 -10.661249514758 4 0.060826936582 5 -2.635282790368 6 -0.022513797598 7 -4.762328938804 8 0.007796187842
! 144: 1..........2....................3........4...................5...........6...........7..........8........ 1 -12.259301101757 2 0.092547594732 3 -10.639466323537 4 0.060551352623 5 -2.684448479366 6 -0.021830603905 7 -4.795109470365 8 0.008128524946
! 145: 1..........2........3..........4....................5............6...........7............8.................. 1 -12.250033234263 2 0.092418316868 3 -10.671177063275 4 0.060960172903 5 -2.649386951711 6 -0.022262295734 7 -4.777729640313 8 0.007908146901
! 146: 1...........2..........3..........4....................5............6............7..............8................. 1 -12.245968180439 2 0.092366839432 3 -10.671951574439 4 0.060969690741 5 -2.644234071218 6 -0.022336213263 7 -4.775102360558 8 0.007877895201
! 147: 1...........2..........................3...................4.........................5......................6...................7............8..................
! 148: #Number of iterations = 34, -2 Log likelihood = 46542.378402526519
! 149: # Parameters nlstate*nlstate*ncov a12*1 + b12 * age + ...
! 150: 12 -12.245504 0.092361
! 151: 13 -10.670953 0.060958
! 152: 21 -2.645034 -0.022329
! 153: 23 -4.775501 0.007884
! 154:
! 155: Calculation of the hessian matrix. Wait...
! 156: 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
! 157:
! 158: Inverting the hessian to get the covariance matrix. Wait...
! 159:
! 160: #Hessian matrix#
! 161: 2.314e+03 1.908e+05 4.106e+02 3.377e+04 -3.590e+02 -2.993e+04 -3.182e+02 -2.620e+04
! 162: 1.908e+05 1.581e+07 3.349e+04 2.769e+06 -2.988e+04 -2.505e+06 -2.611e+04 -2.161e+06
! 163: 4.106e+02 3.349e+04 8.126e+02 6.556e+04 -8.087e+01 -6.783e+03 3.862e+02 3.192e+04
! 164: 3.377e+04 2.769e+06 6.556e+04 5.317e+06 -6.823e+03 -5.756e+05 3.190e+04 2.651e+06
! 165: -3.590e+02 -2.988e+04 -8.087e+01 -6.823e+03 4.524e+02 3.800e+04 4.993e+01 4.166e+03
! 166: -2.993e+04 -2.505e+06 -6.783e+03 -5.756e+05 3.800e+04 3.210e+06 4.173e+03 3.500e+05
! 167: -3.182e+02 -2.611e+04 3.862e+02 3.190e+04 4.993e+01 4.173e+03 1.039e+03 8.943e+04
! 168: -2.620e+04 -2.161e+06 3.192e+04 2.651e+06 4.166e+03 3.500e+05 8.943e+04 7.740e+06
! 169: # Scales (for hessian or gradient estimation)
! 170: 12 1.00000e-04 1.00000e-06
! 171: 13 1.00000e-04 1.00000e-06
! 172: 21 1.00000e-03 1.00000e-06
! 173: 23 1.00000e-04 1.00000e-05
! 174: # Covariance matrix
! 175: # 121 Var(a12)
! 176: # 122 Cov(b12,a12) Var(b12)
! 177: # ...
! 178: # 232 Cov(b23,a12) Cov(b23,b12) ... Var (b23)
! 179: #121 Var(a12)
! 180: #122 Cov(b12,a12) Var(b12)
! 181: #131 Cov(a13,a12) Cov(a13,b12) Var(a13)
! 182: #132 Cov(b13,a12) Cov(b13,b12) Cov(b13,a13) Var(b13)
! 183: #211 Cov(a21,a12) Cov(a21,b12) Cov(a21,a13) Cov(a21,b13) Var(a21)
! 184: #212 Cov(b21,a12) Cov(b21,b12) Cov(b21,a13) Cov(b21,b13) Cov(b21,a21) Var(b21)
! 185: #231 Cov(a23,a12) Cov(a23,b12) Cov(a23,a13) Cov(a23,b13) Cov(a23,a21) Cov(a23,b21) Var(a23)
! 186: #232 Cov(b23,a12) Cov(b23,b12) Cov(b23,a13) Cov(b23,b13) Cov(b23,a21) Cov(b23,b21) Cov(b23,a23) Var(b23)
! 187: 121 1.13401e-01
! 188: 122 -1.36090e-03 1.64264e-05
! 189: 131 -6.89393e-02 8.35716e-04 3.31264e-01
! 190: 132 8.21021e-04 -1.00329e-05 -4.07026e-03 5.03323e-05
! 191: 211 7.49121e-02 -8.92959e-04 1.95088e-02 -2.56008e-04 4.75336e-01
! 192: 212 -8.93092e-04 1.06987e-05 -2.49893e-04 3.28027e-06 -5.63203e-03 6.70896e-05
! 193: 231 5.58947e-02 -6.39542e-04 -1.16708e-01 1.34217e-03 -9.40960e-04 1.38147e-05 2.47512e-01
! 194: 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
! 195: begin-prev-date=1/1/1984 end-prev-date=1/6/1988 mov_average=0
! 196: prevforecast=0 starting-proj-date=1/1/2000 final-proj-date=1/1/2000 mobil_average=0
! 197: Computing stable prevalence: result on file 'plrbiaspar.txt'
! 198:
! 199: #************
! 200: Computing pij: result on file 'pijrbiaspar.txt'
! 201: Computing standard deviation of one-step probabilities: result on file 'probrbiaspar.txt'
! 202: Computing matrix of variance covariance of one-step probabilities: result on file 'probcovrbiaspar.txt'
! 203: and correlation matrix of one-step probabilities: result on file 'probcorrbiaspar.txt'
! 204: 65 13-12 mu 1.4597e-02 2.3277e-02 Var 3.1453e-06 3.1816e-06 cor -0.393 cov -1.2429e-06 Eig 4.407e-06 1.920e-06 1stv 0.702 -0.712 tan -1.015
! 205: 70 13-12 mu 1.9767e-02 3.6881e-02 Var 3.1424e-06 4.5584e-06 cor -0.413 cov -1.5617e-06 Eig 5.565e-06 2.136e-06 1stv 0.542 -0.841 tan -1.551
! 206: 75 13-12 mu 2.6747e-02 5.8389e-02 Var 2.7455e-06 5.6383e-06 cor -0.448 cov -1.7608e-06 Eig 6.471e-06 1.913e-06 1stv 0.427 -0.904 tan -2.116
! 207: 80 13-12 mu 3.6146e-02 9.2323e-02 Var 2.8030e-06 6.6320e-06 cor -0.477 cov -2.0574e-06 Eig 7.528e-06 1.907e-06 1stv 0.399 -0.917 tan -2.297
! 208: 85 13-12 mu 4.8754e-02 1.4570e-01 Var 7.0382e-06 1.5211e-05 cor -0.422 cov -4.3687e-06 Eig 1.711e-05 5.143e-06 1stv 0.398 -0.917 tan -2.305
! 209: 90 13-12 mu 6.5564e-02 2.2924e-01 Var 2.7005e-05 7.6142e-05 cor -0.376 cov -1.7053e-05 Eig 8.148e-05 2.167e-05 1stv 0.299 -0.954 tan -3.195
! 210: 95 13-12 mu 8.7773e-02 3.5907e-01 Var 9.3321e-05 3.7748e-04 cor -0.369 cov -6.9221e-05 Eig 3.934e-04 7.736e-05 1stv 0.225 -0.974 tan -4.336
! 211: 65 21-12 mu 1.9364e-01 2.3277e-02 Var 9.6637e-04 3.1816e-06 cor 0.317 cov 1.7554e-05 Eig 9.667e-04 2.862e-06 1stv 1.000 0.018 tan 0.018
! 212: 70 21-12 mu 1.7338e-01 3.6881e-02 Var 4.5522e-04 4.5584e-06 cor 0.314 cov 1.4302e-05 Eig 4.557e-04 4.105e-06 1stv 0.999 0.032 tan 0.032
! 213: 75 21-12 mu 1.5521e-01 5.8389e-02 Var 1.8567e-04 5.6383e-06 cor 0.307 cov 9.9468e-06 Eig 1.862e-04 5.090e-06 1stv 0.998 0.055 tan 0.055
! 214: 80 21-12 mu 1.3892e-01 9.2323e-02 Var 6.7602e-05 6.6320e-06 cor 0.290 cov 6.1458e-06 Eig 6.822e-05 6.019e-06 1stv 0.995 0.099 tan 0.100
! 215: 85 21-12 mu 1.2432e-01 1.4570e-01 Var 3.9595e-05 1.5211e-05 cor 0.281 cov 6.8846e-06 Eig 4.140e-05 1.340e-05 1stv 0.967 0.254 tan 0.263
! 216: 90 21-12 mu 1.1124e-01 2.2924e-01 Var 6.0749e-05 7.6142e-05 cor 0.298 cov 2.0238e-05 Eig 9.010e-05 4.679e-05 1stv 0.568 0.823 tan 1.450
! 217: 95 21-12 mu 9.9518e-02 3.5907e-01 Var 1.0464e-04 3.7748e-04 cor 0.306 cov 6.0852e-05 Eig 3.904e-04 9.169e-05 1stv 0.208 0.978 tan 4.697
! 218: 65 23-12 mu 1.6392e-01 2.3277e-02 Var 4.4973e-04 3.1816e-06 cor 0.380 cov 1.4384e-05 Eig 4.502e-04 2.719e-06 1stv 0.999 0.032 tan 0.032
! 219: 70 23-12 mu 1.7070e-01 3.6881e-02 Var 3.0743e-04 4.5584e-06 cor 0.402 cov 1.5054e-05 Eig 3.082e-04 3.812e-06 1stv 0.999 0.050 tan 0.050
! 220: 75 23-12 mu 1.7774e-01 5.8389e-02 Var 1.8810e-04 5.6383e-06 cor 0.435 cov 1.4152e-05 Eig 1.892e-04 4.547e-06 1stv 0.997 0.077 tan 0.077
! 221: 80 23-12 mu 1.8502e-01 9.2323e-02 Var 1.0082e-04 6.6320e-06 cor 0.456 cov 1.1795e-05 Eig 1.023e-04 5.178e-06 1stv 0.992 0.122 tan 0.123
! 222: 85 23-12 mu 1.9258e-01 1.4570e-01 Var 5.6251e-05 1.5211e-05 cor 0.358 cov 1.0482e-05 Eig 5.877e-05 1.269e-05 1stv 0.972 0.234 tan 0.241
! 223: 90 23-12 mu 2.0041e-01 2.2924e-01 Var 6.6788e-05 7.6142e-05 cor 0.266 cov 1.8971e-05 Eig 9.100e-05 5.193e-05 1stv 0.617 0.787 tan 1.276
! 224: 95 23-12 mu 2.0853e-01 3.5907e-01 Var 1.4685e-04 3.7748e-04 cor 0.256 cov 6.0241e-05 Eig 3.923e-04 1.321e-04 1stv 0.238 0.971 tan 4.074
! 225: 65 21-13 mu 1.9364e-01 1.4597e-02 Var 9.6637e-04 3.1453e-06 cor 0.029 cov 1.6142e-06 Eig 9.664e-04 3.143e-06 1stv 1.000 0.002 tan 0.002
! 226: 70 21-13 mu 1.7338e-01 1.9767e-02 Var 4.5522e-04 3.1424e-06 cor 0.021 cov 7.7634e-07 Eig 4.552e-04 3.141e-06 1stv 1.000 0.002 tan 0.002
! 227: 75 21-13 mu 1.5521e-01 2.6747e-02 Var 1.8567e-04 2.7455e-06 cor 0.010 cov 2.1716e-07 Eig 1.857e-04 2.745e-06 1stv 1.000 0.001 tan 0.001
! 228: 80 21-13 mu 1.3892e-01 3.6146e-02 Var 6.7602e-05 2.8030e-06 cor 0.018 cov 2.4908e-07 Eig 6.760e-05 2.802e-06 1stv 1.000 0.004 tan 0.004
! 229: 85 21-13 mu 1.2432e-01 4.8754e-02 Var 3.9595e-05 7.0382e-06 cor 0.079 cov 1.3170e-06 Eig 3.965e-05 6.985e-06 1stv 0.999 0.040 tan 0.040
! 230: 90 21-13 mu 1.1124e-01 6.5564e-02 Var 6.0749e-05 2.7005e-05 cor 0.099 cov 4.0068e-06 Eig 6.122e-05 2.654e-05 1stv 0.993 0.116 tan 0.117
! 231: 95 21-13 mu 9.9518e-02 8.7773e-02 Var 1.0464e-04 9.3321e-05 cor 0.091 cov 8.9980e-06 Eig 1.096e-04 8.835e-05 1stv 0.875 0.483 tan 0.552
! 232: 65 23-13 mu 1.6392e-01 1.4597e-02 Var 4.4973e-04 3.1453e-06 cor -0.473 cov -1.7779e-05 Eig 4.504e-04 2.439e-06 1stv 0.999 -0.040 tan -0.040
! 233: 70 23-13 mu 1.7070e-01 1.9767e-02 Var 3.0743e-04 3.1424e-06 cor -0.506 cov -1.5725e-05 Eig 3.082e-04 2.332e-06 1stv 0.999 -0.051 tan -0.052
! 234: 75 23-13 mu 1.7774e-01 2.6747e-02 Var 1.8810e-04 2.7455e-06 cor -0.561 cov -1.2754e-05 Eig 1.890e-04 1.872e-06 1stv 0.998 -0.068 tan -0.068
! 235: 80 23-13 mu 1.8502e-01 3.6146e-02 Var 1.0082e-04 2.8030e-06 cor -0.595 cov -1.0008e-05 Eig 1.018e-04 1.792e-06 1stv 0.995 -0.101 tan -0.101
! 236: 85 23-13 mu 1.9258e-01 4.8754e-02 Var 5.6251e-05 7.0382e-06 cor -0.519 cov -1.0318e-05 Eig 5.833e-05 4.963e-06 1stv 0.980 -0.197 tan -0.201
! 237: 90 23-13 mu 2.0041e-01 6.5564e-02 Var 6.6788e-05 2.7005e-05 cor -0.461 cov -1.9572e-05 Eig 7.480e-05 1.899e-05 1stv 0.925 -0.379 tan -0.409
! 238: 95 23-13 mu 2.0853e-01 8.7773e-02 Var 1.4685e-04 9.3321e-05 cor -0.418 cov -4.8891e-05 Eig 1.758e-04 6.435e-05 1stv 0.860 -0.510 tan -0.593
! 239: 65 23-21 mu 1.6392e-01 1.9364e-01 Var 4.4973e-04 9.6637e-04 cor -0.027 cov -1.7471e-05 Eig 9.670e-04 4.491e-04 1stv 0.034 -0.999 tan -29.605
! 240: 70 23-21 mu 1.7070e-01 1.7338e-01 Var 3.0743e-04 4.5522e-04 cor -0.023 cov -8.6175e-06 Eig 4.557e-04 3.069e-04 1stv 0.058 -0.998 tan -17.208
! 241: 75 23-21 mu 1.7774e-01 1.5521e-01 Var 1.8810e-04 1.8567e-04 cor -0.021 cov -3.9354e-06 Eig 1.910e-04 1.828e-04 1stv 0.805 -0.594 tan -0.738
! 242: 80 23-21 mu 1.8502e-01 1.3892e-01 Var 1.0082e-04 6.7602e-05 cor -0.028 cov -2.2822e-06 Eig 1.010e-04 6.745e-05 1stv 0.998 -0.068 tan -0.068
! 243: 85 23-21 mu 1.9258e-01 1.2432e-01 Var 5.6251e-05 3.9595e-05 cor -0.059 cov -2.8004e-06 Eig 5.671e-05 3.914e-05 1stv 0.987 -0.161 tan -0.164
! 244: 90 23-21 mu 2.0041e-01 1.1124e-01 Var 6.6788e-05 6.0749e-05 cor -0.076 cov -4.8520e-06 Eig 6.948e-05 5.805e-05 1stv 0.874 -0.486 tan -0.556
! 245: 95 23-21 mu 2.0853e-01 9.9518e-02 Var 1.4685e-04 1.0464e-04 cor -0.064 cov -7.9663e-06 Eig 1.483e-04 1.032e-04 1stv 0.984 -0.179 tan -0.182
! 246: Computing Total LEs with variances: file 'trbiaspar.txt'
! 247: Computing Health Expectancies: result on file 'erbiaspar.txt'
! 248: Computing Variance-covariance of DFLEs: file 'vrbiaspar.txt'
! 249: 65|66|67|68|69|70|71|72|73|74|75|76|77|78|79|80|81|82|83|84|85|86|87|88|89|90|91|92|93|94|95|
! 250: Computing total mortality p.j=w1*p1j+w2*p2j+..: result on file 'prmorprev1-stablbased-rbiaspar.txt'
! 251: End of Imach
! 252: Local time at start Mon Oct 24 14:51:36 2005
! 253:
! 254: Local time at end Mon Oct 24 16:19:18 2005
! 255:
! 256: Total time used 0 day(s) 1 hour(s) 27 minute(s) 42 second(s)
! 257: Total time was 5262 Sec.
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>