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