imach/src/ChangeLog - view

File: [Local Repository] / imach / src / ChangeLog
Revision 1.65: download - view: text, annotated - select for diffs
Mon May 8 17:46:42 2023 UTC (17 months, 4 weeks ago) by brouard
Branches: MAIN
CVS tags: HEAD

*** empty log message ***

1: 2023-04-29 Nicolas Brouard <brouard@ined.fr> 2: 3: * imach.c (Module): Inverting the model equation and fixingg bugs 4: in the drawings 5: "1+age+V7*V4*age+V6*V4*age+V7*V3*age+V6*V3*age+V6*V2*age+V7*age+V6*age+V4*age+V3*age+V2*age+V7*V4+V6*V4+V7*V3+V6*V3+V7*V2+V6*V2+V7+V6+V4+V3+V2\000\0003+V2" 6: 7: 2023-04-24 Nicolas Brouard <brouard@ined.fr> 8: 9: * imach.c (Module): some bugs in printinggnuplot 10: 11: 2022-09-19 Nicolas Brouard <brouard@ined.fr> 12: 13: * imach.c (Module): Version 0.99r42 needed a newer version of 14: Gnuplot. But newer version 0.99r43 should run with the Gnuplot 15: version 5.0 or 5.1 distributed with IMaCh. 16: 17: 2022-09-18 Nicolas Brouard <brouard@ined.fr> 18: 19: * imach.c (Module): Fixing when a state in the data is 0 or higher 20: than lstate+ndeath. Fixing the plot of contribution to Likelihood. 21: 22: 2022-09-16 Nicolas Brouard <brouard@ined.fr> 23: 24: * imach.c (Module): 0.99r41 Was an error when product of timevarying and fixed. Using FixedV[of name] now. Thank you Feinuo 25: 26: 2022-09-14 Nicolas Brouard <brouard@ined.fr> 27: 28: * imach.c (Module): Fixing names of variables in T_ (thanks to Feinuo) 29: 30: 2022-09-14 Nicolas Brouard <brouard@ined.fr> 31: 32: * imach.c (Module): Version 0.99r39 with colored dummy covariates 33: (fixed or time varying), using new last columns of 34: ILK_parameter.txt file. 35: 36: 2022-09-11 Nicolas Brouard <brouard@ined.fr> 37: 38: * imach.c (Module): Adding timevarying products of any kinds, 39: should work before shifting cotvar from ncovcol+nqv columns in 40: order to have a correspondance between the column of cotvar and 41: the id of column. 42: (Module): Some cleaning and adding covariates in ILK.txt 43: 44: 2022-09-09 Nicolas Brouard <brouard@ined.fr> 45: 46: * imach.c (Module): Many improvements for fixing products of fixed 47: timevarying as well as fixed * fixed, and test with quantitative 48: covariate. 49: 50: 2022-09-04 Nicolas Brouard <brouard@ined.fr> 51: 52: * imach.c (Module): Now the easy runs i.e. without result or 53: model=1+age only did not work. The defautl combination should be 1 54: and not 0 because everything hasn't been tranformed yet. 55: 56: 2022-09-02 Nicolas Brouard <brouard@ined.fr> 57: 58: * imach.c: Version 0.99r35 because it outputs same results with 59: 1+age+V1+V1*age for females and 1+age for females only 60: (education=1 noweight) 61: 62: 2022-08-31 Nicolas Brouard <brouard@ined.fr> 63: 64: * imach.c (Module): Some improvments in fichtm and many verifications 0.99r34 65: 66: 2022-08-21 Nicolas Brouard <brouard@ined.fr> 67: 68: * imach.c (Module): Version 0.99r33 A lot of changes in 69: reassigning covariates: my first idea was that people will always 70: use the first covariate V1 into the model but in fact they are 71: producing data with many covariates and can use an equation model 72: with some of the covariate; it means that in a model V2+V3 instead 73: of codtabm(k,Tvaraff[j]) which calculates for combination k, for 74: three covariates (V1, V2, V3) the value of Tvaraff[j], but in fact 75: the equation model is restricted to two variables only (V2, V3) 76: and the combination for V2 should be codtabm(k,1) instead of 77: (codtabm(k,2), and the code should be 78: codtabm(k,TnsdVar[Tvaraff[j]]. Many many changes have been 79: made. All of these should be simplified once a day like we did in 80: hpxij() for example by using precov[nres] which is computed in 81: decoderesult for each nres of each resultline. Loop should be done 82: on the equation model globally by distinguishing only product with 83: age (which are changing with age) and no more on type of 84: covariates, single dummies, single covariates. 85: 86: 2022-08-06 Nicolas Brouard <brouard@ined.fr> 87: 88: * imach.c (Module): Version of imach using partly decoderesult to rebuild xpxij function 89: 90: 2022-08-03 Nicolas Brouard <brouard@ined.fr> 91: 92: * imach.c (Module): Many errors in graphs fixed with Vn*age covariates. 93: 94: 2022-07-25 Brouard Nicolas <brouard@brouard.name> 95: 96: * imach.c (Module): Error cptcovn instead of nsd in bmij (was 97: coredumped, revealed by Feiuno, thank you. 98: 99: 2022-07-23 Nicolas Brouard <brouard@ined.fr> 100: 101: * r29 W and not sqrt(Wald) 102: 103: 2022-07-22 Nicolas Brouard <brouard@ined.fr> 104: 105: * imach.c (Module): Output of Wald test in the htm file and not only in the log. 106: 107: 2022-06-02 Brouard Nicolas <brouard@brouard.name> 108: 109: * imach.c (Module): Adding the Wald tests from the log to the main 110: htm for better display of the maximum likelihood estimators. 111: 112: 2022-05-30 Brouard Nicolas <brouard@brouard.name> 113: 114: * imach.c: With products of covariates (age or dummies or 115: quantitatives), estimates of parameters were good but when 116: estimating variances the positions of newly created covariates 117: were wrongly assigned and results were wrong. Thank to Feinua_Sun! 118: 119: 2022-05-24 Nicolas Brouard <brouard@ined.fr> 120: 121: * imach.c (Module): Some attempts to find a bug of wrong estimates 122: of confidencce intervals with product in the equation modelC 123: 124: 2022-05-15 Nicolas Brouard <brouard@ined.fr> 125: 126: * imach.c (Module): Some minor improvements 127: 128: 2022-04-13 Brouard Nicolas <brouard@brouard.name> 129: 130: * imach.c (Module): Adding link to text data files 131: 132: 2022-04-11 Nicolas Brouard <brouard@ined.fr> 133: 134: * imach.c (Module): Error in rewriting the 'r' file with yearsfproj or yearsbproj fixed 135: 136: 2022-04-05 Brouard Nicolas <brouard@ined.fr> 137: 138: * imach.c (Module): Fixed covariates (dummy or quantitative) 139: with missing values have never been allowed but are ERRORS and 140: program quits. Standard deviations of fixed covariates were 141: wrongly computed. Mean and standard deviations of time varying 142: covariates are still not computed. 143: 144: 2022-03-17 Brouard Nicolas <brouard@ined.fr> 145: 146: * 99r25 Improvements in fixing discrepancies between covariates 147: missing in result line but in model. 148: 149: *Wed Mar 31 2021 Nicolas Brouard <brouard@ined.fr> 150: * imach.c (Module): Still bugs in the result loop. Thank to Holly Benett 151: 152: *2021-03-08 Nicolas Brouard <brouard@ined.fr> 153: 154: * ./ (Module): Fixed bug on result: 155: 156: *2021-02-20 Nicolas Brouard <brouard@ined.fr> 157: 158: * imach.c (Module): Fix bug on quitting after result lines! 159: (Module): Version 0.99r21 160: 161: 2021-02-12 Nicolas Brouard <brouard@ined.fr> 162: 163: * imach.c (Module): The use of a Windows BOM (huge) file is now an error 164: 165: 166: 2021-02-11 Nicolas Brouard <brouard@ined.fr> 167: 168: * (Module): imach.c Someone entered 'results:' instead of 'result:'. Now it is an error which is printed. 169: 170: 2020-02-22 Brouard Nicolas <brouard@ined.fr> 171: 172: * (Module): imach.c Update mle=-3 (for computing Life expectancy 173: and life table from the data without any state) 174: 175: 2019-05-20 Brouard Nicolas <brouard@ined.fr> 176: 177: * version.h: Summary: version 0.99.r19 178: 179: * imach.c: Summary: Projection syntax simplified 180: We can now start projections, forward or backward, from the mean date 181: of inteviews up to or down to a number of years of projection: 182: prevforecast=1 yearsfproj=15.3 mobil_average=0 183: or 184: prevforecast=1 starting-proj-date=1/1/2007 final-proj-date=12/31/2017 mobil_average=0 185: or 186: prevbackcast=1 yearsbproj=12.3 mobil_average=1 187: or 188: prevbackcast=1 starting-back-date=1/10/1999 final-back-date=1/1/1985 mobil_average=1 189: 190: 2019-05-18 Brouard Nicolas <brouard@ined.fr> 191: 192: * imach.c: Summary: doxygen tex bug 193: 194: 2019-05-16 Brouard Nicolas <brouard@ined.fr> 195: 196: * imach.c: Summary: There was some wrong lines added 197: 198: 2019-05-09 Brouard Nicolas <brouard@ined.fr> 199: 200: * Makefile, imach.c: *** empty log message *** 201: 202: * ChangeLog: Summary: C= 203: 204: * imach.c: Summary: Some updates 205: 206: * imach.c: Summary: Before ncovmax 207: 208: * version.h: *** empty log message *** 209: 210: * imach.c: Summary: 0.99r18 unlimited number of individuals 211: The number n which was limited to 20,000 cases is now unlimited, from firstobs to lastobs. If the number is too for the virtual memory, probably an error will occur. 212: 213: 2018-12-13 Brouard Nicolas <brouard@ined.fr> 214: 215: * imach.c: Summary: Bug for young ages (<-30) will be in r17 216: 217: 2018-12-05 Brouard Nicolas <brouard@ined.fr> 218: 219: * Makefile: Summary: 0.99r16 220: 221: 2018-05-02 Brouard Nicolas <brouard@ined.fr> 222: 223: * imach.c: Summary: Some bugs fixed 224: 225: 2018-05-02 Brouard Nicolas <brouard@ined.fr> 226: 227: * imach.c: Summary: Bug for young ages (<-30) will be in r17 228: 229: 2018-05-02 Brouard Nicolas <brouard@ined.fr> 230: 231: * imach.c: Summary: Some bugs fixed 232: 233: 2018-05-01 Brouard Nicolas <brouard@ined.fr> 234: 235: * imach.c: 236: Summary: Bug fixed by providing frequencies only for non missing covariates 237: 238: 2018-04-27 Brouard Nicolas <brouard@ined.fr> 239: 240: * imach.c: Summary: some minor bugs 241: 242: 2018-04-21 Brouard Nicolas <brouard@ined.fr> 243: 244: * README.txt: *** empty log message *** 245: 246: * imach.c: Summary: Some bugs fixed, valgrind tested 247: 248: 2018-04-20 Brouard Nicolas <brouard@ined.fr> 249: 250: * Makefile, README.txt, setup.iss: Summary: imach 0.99r16 251: 252: * imach.c: 253: Summary: Computing mean and stdeviation of fixed quantitative variables 254: 255: 2018-04-19 Brouard Nicolas <brouard@ined.fr> 256: 257: * imach.c: Summary: Some minor bugs fixed 258: 259: * version.h: *** empty log message *** 260: 261: 2018-02-27 Brouard Nicolas <brouard@ined.fr> 262: 263: * imach.c: *** empty log message *** 264: 265: * imach.c: Summary: Adding second argument for quitting 266: 267: 2018-02-21 Brouard Nicolas <brouard@ined.fr> 268: 269: * Makefile, README.txt, imach.c: Summary: 0.99r15 270: New Makefile with recent VirtualBox 5.26. Bug in sqrt negatve in imach.c 271: 272: 2017-07-20 Brouard Nicolas <brouard@ined.fr> 273: 274: * imach.c: Summary: temporary working 275: 276: 2017-07-19 Brouard Nicolas <brouard@ined.fr> 277: 278: * imach.c: Summary: Bug for mobil_average=0 and prevforecast fixed(?) 279: 280: 2017-07-17 Brouard Nicolas <brouard@ined.fr> 281: 282: * imach.c: Summary: BOM files can be read now 283: 284: 2017-06-30 Brouard Nicolas <brouard@ined.fr> 285: 286: * imach.c: Summary: Graphs improvements 287: 288: * imach.c: Summary: Saito's color 289: 290: 2017-06-29 Brouard Nicolas <brouard@ined.fr> 291: 292: * imach.c: Summary: Version 0.99r14 293: 294: 2017-06-27 Brouard Nicolas <brouard@ined.fr> 295: 296: * imach.c: Summary: More documentation on projections 297: 298: * imach.c: Summary: Color of backprojection changed from 6 to 5(yellow) 299: 300: * imach.c: Summary: Some bug with rint 301: 302: 2017-05-24 Brouard Nicolas <brouard@ined.fr> 303: 304: * imach.c: *** empty log message *** 305: 306: 2017-05-23 Brouard Nicolas <brouard@ined.fr> 307: 308: * imach.c: Summary: Code into subroutine, cleanings 309: 310: 2017-05-18 Brouard Nicolas <brouard@ined.fr> 311: 312: * imach.c: 313: Summary: backprojection and confidence intervals of backprevalence 314: 315: 2017-05-13 Brouard Nicolas <brouard@ined.fr> 316: 317: * imach.c: Summary: temporary save for backprojection 318: 319: * imach.c: Summary: Version 0.99r13 (improvements and bugs fixed) 320: 321: 2017-04-26 Brouard Nicolas <brouard@ined.fr> 322: 323: * imach.c: Summary: imach 0.99r13 Some bugs fixed 324: 325: * imach.c: Summary: Labels in graphs 326: 327: 2017-04-24 Brouard Nicolas <brouard@ined.fr> 328: 329: * imach.c: Summary: to save 330: 331: 2017-04-18 Brouard Nicolas <brouard@ined.fr> 332: 333: * imach.c: *** empty log message *** 334: 335: 2017-04-05 Brouard Nicolas <brouard@ined.fr> 336: 337: * imach.c: Summary: Bug in E_ as well as in T_ fixed nres-1 vs k1-1 338: 339: 2017-04-04 Brouard Nicolas <brouard@ined.fr> 340: 341: * imach.c: Summary: Gnuplot indexations fixed (humm) 342: 343: * imach.c: 344: Summary: Some errors to warnings only if date of death is unknown but status is death we could set to pi3 345: 346: 2017-04-03 Brouard Nicolas <brouard@ined.fr> 347: 348: * imach.c: Summary: Version 0.99r12 349: Some cleanings, conformed with updated documentation. 350: 351: 2017-03-30 Brouard Nicolas <brouard@ined.fr> 352: 353: * version.h: *** empty log message *** 354: 355: * version.h: Summary:0.99r11 356: 357: 2017-03-29 Brouard Nicolas <brouard@ined.fr> 358: 359: * imach.c: Summary: Temp 360: 361: 2017-03-27 Brouard Nicolas <brouard@ined.fr> 362: 363: * imach.c: Summary: Temporary 364: 365: 2017-03-08 Brouard Nicolas <brouard@ined.fr> 366: 367: * imach.c: Summary: IMaCh version 0.99r10 bugs in gnuplot fixed 368: 369: * imach.c: Summary: Fixing data parameter line 370: 371: 2016-12-15 Brouard Nicolas <brouard@ined.fr> 372: 373: * imach.c: Summary: 0.99 in progress 374: 375: 2016-09-15 Brouard Nicolas <brouard@ined.fr> 376: 377: * imach.c: *** empty log message *** 378: 379: * imach.c: Summary: not working 380: 381: 2016-09-08 Brouard Nicolas <brouard@ined.fr> 382: 383: * imach.c: Summary: continue 384: 385: 2016-09-07 Brouard Nicolas <brouard@ined.fr> 386: 387: * imach.c: Summary: Starting values from frequencies 388: 389: * imach.c: *** empty log message *** 390: 391: 2016-09-02 Brouard Nicolas <brouard@ined.fr> 392: 393: * imach.c: *** empty log message *** 394: 395: 2016-08-30 Brouard Nicolas <brouard@ined.fr> 396: 397: * imach.c: Summary: Fixing a lots 398: 399: 2016-08-29 Brouard Nicolas <brouard@ined.fr> 400: 401: * imach.c: Summary: gnuplot problem in Back projection to fix 402: 403: * imach.c: Summary: Better 404: 405: 2016-08-26 Brouard Nicolas <brouard@ined.fr> 406: 407: * imach.c: 408: Summary: Improvement in Powell output in order to copy and paste 409: 410: * imach.c: Summary: Starting tests of 0.99 411: 412: * imach.c: Summary: to valgrind 413: 414: 2016-08-25 Brouard Nicolas <brouard@ined.fr> 415: 416: * imach.c: *** empty log message *** 417: 418: 2016-08-23 Brouard Nicolas <brouard@ined.fr> 419: 420: * imach.c: *** empty log message *** 421: 422: * imach.c: Summary: not working 423: 424: 2016-08-22 Brouard Nicolas <brouard@ined.fr> 425: 426: * imach.c: Summary: not working 427: 428: * imach.c: Summary: Not working 429: 430: 2016-07-23 Brouard Nicolas <brouard@ined.fr> 431: 432: * imach.c: Summary: Completing for func too 433: 434: 2016-07-22 Brouard Nicolas <brouard@ined.fr> 435: 436: * imach.c: Summary: Fixing some arrays, still debugging 437: 438: 2016-07-21 Brouard Nicolas <brouard@ined.fr> 439: 440: * imach.c: 441: Summary: 0.99 working (more or less) for Asian Workshop on multitate methods 442: 443: 2016-07-12 Brouard Nicolas <brouard@ined.fr> 444: 445: * imach.c: Summary: temp 446: 447: * imach.c: Summary: saving but not running 448: 449: 2016-07-01 Brouard Nicolas <brouard@ined.fr> 450: 451: * imach.c: Summary: Fixes 452: 453: 2016-02-19 Brouard Nicolas <brouard@ined.fr> 454: 455: * imach-0.98r.c: Summary: Kind of 0.98r? series, starting with r7 456: 457: * setup.iss, Makefile, version.h: *** empty log message *** 458: 459: * imach.c: Summary: temporary 460: 461: 2016-02-17 Brouard Nicolas <brouard@ined.fr> 462: 463: * imach.c: Summary: Probably last 0.98 stable version 0.98r6 464: 465: 2016-02-16 Brouard Nicolas <brouard@ined.fr> 466: 467: * imach.c: Summary: minor bug 468: 469: * imach.c: Summary: 0.99r2 470: 471: 2016-02-15 Brouard Nicolas <brouard@ined.fr> 472: 473: * imach.c: *** empty log message *** 474: 475: 2016-02-12 Brouard Nicolas <brouard@ined.fr> 476: 477: * imach.c: Summary: 0.99 Back projections 478: 479: 2015-12-23 Brouard Nicolas <brouard@ined.fr> 480: 481: * imach.c: Summary: Experimental backcast 482: 483: 2015-12-18 Brouard Nicolas <brouard@ined.fr> 484: 485: * imach.c: Summary: 0.98r4 Warning and status=-2 486: 487: Version 0.98r4 is now: 488: - displaying an error when status is -1, date of interview unknown and date of death known; 489: - permitting a status -2 when the vital status is unknown at a known date of right truncation. 490: Older changes concerning s=-2, dating from 2005 have been supersed. 491: 492: * Makefile: *** empty log message *** 493: 494: 2015-12-16 Brouard Nicolas <brouard@ined.fr> 495: 496: * Makefile: *** empty log message *** 497: 498: * imach.c: Summary: 0.98r4 working 499: 500: * Makefile: Summary: linuxrpm added 501: 502: * Makefile: *** empty log message *** 503: 504: * imach.c: Summary: temporary not working 505: 506: 2015-12-11 Brouard Nicolas <brouard@ined.fr> 507: 508: * imach.c: Summary: 0.98r4 509: 510: 2015-12-08 Brouard Nicolas <brouard@ined.fr> 511: 512: * Makefile: *** empty log message *** 513: 514: 2015-11-21 Brouard Nicolas <brouard@ined.fr> 515: 516: * Makefile: 517: Summary: VBoxManage change from execute --image to run --exe and --wait-exit doesn't exist any more 518: 519: * imach.c: Summary: minor typo 520: 521: * imach.c: Summary: 0.98r3 with some graph of projected cross-sectional 522: 523: Author: Nicolas Brouard 524: 525: 2015-11-18 Brouard Nicolas <brouard@ined.fr> 526: 527: * imach.c: Summary: Start working on projected prevalences 528: 529: 2015-11-17 Brouard Nicolas <brouard@ined.fr> 530: 531: * imach.c: Summary: Adding ftolpl parameter 532: Author: N Brouard 533: 534: We had difficulties to get smoothed confidence intervals. It was due 535: to the period prevalence which wasn't computed accurately. The inner 536: parameter ftolpl is now an outer parameter of the .imach parameter 537: file after estepm. If ftolpl is small 1.e-4 and estepm too, 538: computation are long. 539: 540: * version.h: *** empty log message *** 541: 542: * imach.c: Summary: temporary 543: 544: 2015-10-27 Brouard Nicolas <brouard@ined.fr> 545: 546: * Makefile: *** empty log message *** 547: 548: * setup.iss.in: Summary: some new subdirs 549: 550: * imach.c: *** empty log message *** 551: 552: * setup.iss.in: 553: Summary: Adding gnuplot 5.1 with bin and etc for pango fonts 554: 555: 2015-10-24 Brouard Nicolas <brouard@ined.fr> 556: 557: * imach.c: *** empty log message *** 558: 559: 2015-10-23 Brouard Nicolas <brouard@ined.fr> 560: 561: * imach.c: 562: Summary: 0.98r3 some clarification for graphs on likelihood contributions 563: 564: 2015-10-01 Brouard Nicolas <brouard@ined.fr> 565: 566: * imach.c: Summary: Some new graphs of contribution to likelihood 567: 568: 2015-09-30 Brouard Nicolas <brouard@ined.fr> 569: 570: * Makefile: 571: Summary: Version read from version.h, itself read with cmake 572: 573: * imach.c: Summary: looking at better estimation of the hessian 574: 575: Also a better criteria for convergence to the period prevalence And 576: therefore adding the number of years needed to converge. (The 577: prevalence in any alive state shold sum to one 578: 579: * Makefile: *** empty log message *** 580: 581: 2015-09-22 Brouard Nicolas <brouard@ined.fr> 582: 583: * imach.c: 584: Summary: Adding some overall graph on contribution to likelihood. Might change 585: 586: 2015-09-15 Brouard Nicolas <brouard@ined.fr> 587: 588: * Makefile: *** empty log message *** 589: 590: * imach.c: Summary: 0.98r0 591: 592: - Some new graphs like survival functions 593: - Some bugs fixed like model=1+age+V2. 594: 595: * ChangeLog, Makefile: Summary: 0.98r0 596: 597: * CMakeLists.txt: Summary: Some tests for setup.iss 598: 599: 2015-09-15 Nicolas J Brouard <nbrouard@tugault.ined.fr> 600: 601: * (Module): 0.98r0 Some new graphs, some bugs fixed: model=1+age+V2. etc 602: 603: 2015-08-18 Nicolas J Brouard <nbrouard@tugault.local> 604: 605: * (Module): Adding error when the covariance matrix doesn't contain the exact number of lines required by the model line. 606: 607: 2015-08-03 Nicolas J Brouard <nbrouard@tugault.ined.fr> 608: 609: * (Module): Changing Variance of one-step probabilities into Standard deviation (thanks to Yao-Chi Shih) 610: 611: 2014-12-16 Nicolas J Brouard <nbrouard@tugault.local> 612: 613: * imach.c (Module): Merging 1.61 to 1.162 614: 615: 2014-09-02 Nicolas J Brouard <nbrouard@njbrouard.ined.fr> 616: 617: * prevlim.h (Module): 618: 619: 2014-01-26 Nicolas J Brouard <nbrouard@ledermann-2.local> 620: 621: * imach.c (Module): Trying to merge old staffs together while being at Tokyo. Not tested... 622: (Module): Version 0.98nR Running ok, but output format still only works for three covariates. 623: 624: 2010-04-29 brouard <brouard@ined.fr> 625: 626: * imach.c (Module): Checking covariates for more complex models 627: than V1+V2. A lot of change to be done. Unstable. 628: 629: 2010-04-26 brouard <brouard@ined.fr> 630: 631: * imach.c (Module): merging some libgsl code. Fixing computation 632: of likelione (using inter/intrapolation if mle = 0) in order to 633: get same likelihood as if mle=1. 634: Some cleaning of code and comments added. 635: 636: 2009-10-29 brouard <brouard@ined.fr> 637: 638: * imach.c (Module): Now imach stops if date of birth, at least year of birth, is not given. Some cleaning of the code. 639: 640: 2006-06-30 Brouard Nicolas <brouard@localhost> 641: 642: * imach.c (Module): Clarifications on computing e.j 643: 644: 2006-04-28 Brouard Nicolas <brouard@localhost> 645: 646: * imach.c (Module): Yes the sum of survivors was wrong since 647: imach-114 because nhstepm was no more computed in the age 648: loop. Now we define nhstepma in the age loop. 649: (Module): In order to speed up (in case of numerous covariates) we 650: compute health expectancies (without variances) in a first step 651: and then all the health expectancies with variances or standard 652: deviation (needs data from the Hessian matrices) which slows the 653: computation. 654: In the future we should be able to stop the program is only health 655: expectancies and graph are needed without standard deviations. 656: 657: 2006-04-04 LIEVRE Agnes <lievre@PC_00697.ined.fr> 658: 659: * imach.iss (Module): 660: 661: * imach.c (Module): bugs in evsi + population file added in htm output file 662: 663: 2006-03-22 LIEVRE Agnes <lievre@PC_00697.ined.fr> 664: 665: * imach.c (Module): 666: 667: 2006-03-20 Brouard Nicolas <brouard@localhost> 668: 669: * imach.c (Module): <title> changed, corresponds to .htm file 670: name. <head> headers where missing. 671: 672: * Makefile (Module): 0.98g 673: 674: * imach.c (Module): Weights can have a decimal point as for 675: English (a comma might work with a correct LC_NUMERIC environment, 676: otherwise the weight is truncated). 677: Modification of warning when the covariates values are not 0 or 678: 1. 679: Version 0.98g 680: 681: 2006-03-16 LIEVRE Agnes <lievre@PC_00697.ined.fr> 682: 683: * imach.c (Module): Comments concerning covariates added 684: 685: * noreg-setup.iss (Module): Version 0.98f 686: 687: * Makefile (Module): Version 0.98f make _windows=1 imachcyg-setup 688: 689: * imach.c (Module): refinements in the computation of lli if 690: status=-2 in order to have more reliable computation if stepm is 691: not 1 month. Version 0.98f 692: 693: 2006-03-15 Brouard Nicolas <brouard@localhost> 694: 695: * imach.c (Module): Bug if status = -2, the loglikelihood was 696: computed as likelihood omitting the logarithm. Version O.98e 697: 698: 2006-03-14 Brouard Nicolas <brouard@localhost> 699: 700: * imach.c (Module): varevsij Comments added explaining the second 701: table of variances if popbased=1 . 702: (Module): Covariances of eij, ekl added, graphs fixed, new html link. 703: (Module): Function pstamp added 704: (Module): Version 0.98d 705: 706: 2006-03-06 Brouard Nicolas <brouard@localhost> 707: 708: * imach.c (Module): Variance-covariance wrong links and 709: varian-covariance of ej. is needed (Saito). 710: 711: 2006-02-27 Brouard Nicolas <brouard@localhost> 712: 713: * imach.c (Module): One freematrix added in mlikeli! 714: 715: 2006-02-26 Brouard Nicolas <brouard@localhost> 716: 717: * imach.c (Module): Some improvements in processing parameter 718: filename with strsep. 719: 720: 2006-02-24 Brouard Nicolas <brouard@localhost> 721: 722: * imach.c (Module): Memory leaks checks with valgrind and: 723: datafile was not closed, some imatrix were not freed and on matrix 724: allocation too. 725: 726: 2006-02-09 LIEVRE Agnes <lievre@PC_00697.ined.fr> 727: 728: * strsep.c (Module): Inclusion because not in mingw32!!! 729: 730: * setup.iss.in (Module): IMACHSETUPVERSION added 731: 732: * noreg-setup.iss (Module): added to cvs tree 733: 734: * Makefile (Module): IMACHSETUPVERSION added 735: 736: 2006-01-30 BROUARD Nicolas <brouard@localhost> 737: 738: * Makefile (Module, Module): Latest Makefile 739: 740: * imach.c (Module): Back to gnuplot.exe instead of wgnuplot.exe 741: 742: 2006-01-27 BROUARD Nicolas <brouard@localhost> 743: 744: * noreg-setup.iss (Module): ADDED 745: 746: 2006-01-24 BROUARD Nicolas <brouard@localhost> 747: 748: * imach.c (Module): Comments (lines starting with a #) are allowed in data. 749: 750: 2005-10-25 Nicolas Brouard <brouard@dhcp171.recherche.ined.fr> 751: 752: * Makefile (Module): Added Imach98.pmsp and howto crete a dmg 753: 754: * createdmg.sh (Module): Creates a dmg disk for MaC OS/X 755: 756: 2005-09-30 <lievre@PC_00595> 757: 758: * imach.c (Module): sump fixed, loop imx fixed, and simplifications. 759: (Module): If the status is missing at the last wave but we know 760: that the person is alive, then we can code his/her status as -2 761: (instead of missing=-1 in earlier versions) and his/her 762: contributions to the likelihood is 1 - Prob of dying from last 763: health status (= 1-p13= p11+p12 in the easiest case of somebody in 764: the healthy state at last known wave). Version is 0.98 765: 766: 2004-05-20 Brouard Nicolas <brouard@localhost> 767: 768: * imach.c (Repository): 769: Agnes added a direct estimation of mortality (without the need of 770: computing period prevalence and differential mortality). Thus here 771: is version 0.97a which has been distributed to some people at 772: REVES 16 in Brugge using an Inno setup.exe for PCs. Estimates of 773: mortality using covariates is not done today. Estimating direct 774: mortality is a very different process because it doesn't need 775: interpolation because it is easy to get the lx from the force of 776: the mortality mux in the simplest case as for a Gompertz (log mux 777: = a + b*x . But we have been able to incorporate the new code 778: within former imach program (0.96d) without deteriorating too much 779: the understanding of the program. 780: 781: 2003-06-25 <brouard@BROUARD> 782: 783: * imach.c (Module): On windows (cygwin) function asctime_r doesn't 784: exist so I changed back to asctime which exists. 785: (Module): Version 0.96b 786: 787: 2003-06-24 <brouard@BROUARD> 788: 789: * imach.c (Module): Some bugs corrected for windows. Also, when 790: mle=-1 a template is output in file "or"mypar.txt with the design 791: of the covariance matrix to be input. 792: 793: 2003-06-25 Brouard Nicolas <brouard@homebrou> 794: 795: * imach.c (Repository): Duplicated warning errors corrected. 796: (Repository): Elapsed time after each iteration is now output. It 797: helps to forecast when convergence will be reached. Elapsed time 798: is stamped in powell. We created a new html file for the graphs 799: concerning matrix of covariance. It has extension -cov.htm. 800: 801: 2003-06-23 Brouard Nicolas <brouard@homebrou> 802: 803: * imach.c (Repository): Create a sub-directory where all the secondary files are. Only imach, htm, gp and r(imach) are on the main directory. Correct time and other things. 804: 805: 2003-06-18 Brouard Nicolas <brouard@homebrou> 806: 807: * imach.c (Repository): Suppress abusive calls to append to file html (solution was to put fichtm as a global variable. Version 0.96 808: 809: 2003-06-17 Brouard Nicolas <brouard@homebrou> 810: 811: 812: * imach.c (Repository): Check when date of death was earlier that 813: current date of interview. It may happen when the death was just 814: prior to the death. In this case, dh was negative and likelihood 815: was wrong (infinity). We still send an "Error" but patch by 816: assuming that the date of death was just one stepm after the 817: interview. 818: (Repository): Because some people have very long ID (first column) 819: we changed int to long in num[] and we added a new lvector for 820: memory allocation. But we also truncated to 8 characters (left 821: truncation) 822: 823: (Repository): No more line truncation errors. 824: 825: 2003-06-14 Brouard Nicolas <brouard@homebrou> 826: 827: * imach.c (Repository): Add of a routine likelione (likelihood 828: only once) which prints on a text file 829: (ilk) the contributions to the likelihood for each 830: individual/wave. 831: 832: 2003-06-13 Brouard Nicolas <brouard@homebrou> 833: * imach.c (Repository): Replace "freqsummary" at a correct 834: place. It differs from routine "prevalence" which may be called 835: many times. Probs is memory consuming and must be used with 836: parcimony. 837: 838: 2003-06-17 <brouard@BROUARD> 839: 840: * timeval.h (Module): Added included file to make use of 841: gettimeofday working on win32 with cygwin. 842: 843: 2003-05-16 Brouard Nicolas <brouard@brouard> 844: * imach.c (Module): Insert of a warning if the delay between two 845: waves is negative. 846: 847: * imach.c (Module): 848: 849: 2003-05-03 <brouard@BROUARD> 850: 851: * imach.c: Some cleaning 852: 853: 2003-03-28 <brouard@BROUARD> 854: * imach.c (Module): In version up to 0.92 likelihood was computed 855: as if date of death was unknown. Death was treated as any other 856: health state: the date of the interview describes the actual state 857: and not the date of a change in health state. The former idea was 858: to consider that at each interview the state was recorded 859: (healthy, disable or death) and IMaCh was corrected; but when we 860: introduced the exact date of death then we should have modified 861: the contribution of an exact death to the likelihood. This new 862: contribution is smaller and very dependent of the step unit 863: stepm. It is no more the probability to die between last interview 864: and month of death but the probability to survive from last 865: interview up to one month before death multiplied by the 866: probability to die within a month. Thanks to Chris 867: Jackson for correcting this bug. Former versions increased 868: mortality artificially. The bad side is that we add another loop 869: which slows down the processing. The difference can be up to 10% 870: lower mortality. 871: 872: 2002-05-30 <brouard@BROUARD> 873: 874: * imach.c (Module): Add correlation matrix of one-step 875: probabilities and covariance matrix 876: