(Module): sump fixed, loop imx fixed, and simplifications.
(Module): If the status is missing at the last wave but we know
that the person is alive, then we can code his/her status as -2
(instead of missing=-1 in earlier versions) and his/her
contributions to the likelihood is 1 - Prob of dying from last
health status (= 1-p13= p11+p12 in the easiest case of somebody in
the healthy state at last known wave). Version is 0.98
N. Brouard [Wed, 4 May 2005 17:24:51 +0000 (17:24 +0000)]
Output of make setup.iss using setup.iss.in as a template.
setup.iss can be compile with Inno Setup compiler to get the IMaCh setup.
c:\Program Files\Inno Setup 4\Compile32.exe /cc setup.iss
In fact we cross compile from Linux to get imach.exe
and are using wine under Linux to run Inno Setup :
wine "C:////Program Files////Inno Setup 4////Compil32.exe" /cc setup.iss
N. Brouard [Mon, 24 May 2004 18:32:25 +0000 (18:32 +0000)]
Agnes added a direct estimation of mortality (without the need of
computing period prevalence and differential mortality). Thus here
is version 0.97a which has been distributed to some people at
REVES 16 in Brugge using an Inno setup.exe for PCs. Estimates of
mortality using covariates is not done today. Estimating direct
mortality is a very different process because it doesn't need
interpolation because it is easy to get the lx from the force of
the mortality mux in the simplest case as for a Gompertz (log mux
= a + b*x . But we have been able to incorporate the new code
within former imach program (0.96d) without deteriorating too much
the understanding of the program.
Gnuplot is now installed in the same directory (on Windows) as imach.
Thus the full path of gnuplot is executed in order to access the
current version.
N. Brouard [Sun, 16 May 2004 15:05:56 +0000 (15:05 +0000)]
New version 0.97 . First attempt to estimate force of mortality
directly from the data i.e. without the need of knowing the health
state at each age, but using a Gompertz model: log u =a + b*age .
This is the basic analysis of mortality and should be done before any
other analysis, in order to test if the mortality estimated from the
cross-longitudinal survey is different from the mortality estimated
from other sources like vital statistic data.
The same imach parameter file can be used but the option for mle should be -3.
Agn�s, who wrote this part of the code, tried to keep most of the
former routines in order to include the new code within the former code.
The output is very simple: only an estimate of the intercept and of
the slope with 95% confident intervals.
Current limitations:
A) Even if you enter covariates, i.e. with the
model= V1+V2 equation for example, the programm does only estimate a unique global model without covariates.
B) There is no computation of Life Expectancy nor Life Table.
N. Brouard [Wed, 25 Jun 2003 15:30:29 +0000 (15:30 +0000)]
* imach.c (Repository): Duplicated warning errors corrected.
(Repository): Elapsed time after each iteration is now output. It
helps to forecast when convergence will be reached. Elapsed time
is stamped in powell. We created a new html file for the graphs
concerning matrix of covariance. It has extension -cov.htm.
N. Brouard [Tue, 24 Jun 2003 12:34:15 +0000 (12:34 +0000)]
(Module): Some bugs corrected for windows. Also, when
mle=-1 a template is output in file "or"mypar.txt with the design
of the covariance matrix to be input.
N. Brouard [Tue, 24 Jun 2003 12:30:52 +0000 (12:30 +0000)]
(Module): Some bugs corrected for windows. Also, when
mle=-1 a template is output in file "or"mypar.txt with the design
of the covariance matrix to be input.
N. Brouard [Mon, 23 Jun 2003 17:54:56 +0000 (17:54 +0000)]
* 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.
N. Brouard [Tue, 17 Jun 2003 13:12:43 +0000 (13:12 +0000)]
* imach.c (Repository): Check when date of death was earlier that
current date of interview. It may happen when the death was just
prior to the death. In this case, dh was negative and likelihood
was wrong (infinity). We still send an "Error" but patch by
assuming that the date of death was just one stepm after the
interview.
(Repository): Because some people have very long ID (first column)
we changed int to long in num[] and we added a new lvector for
memory allocation. But we also truncated to 8 characters (left
truncation)
(Repository): No more line truncation errors.
N. Brouard [Fri, 13 Jun 2003 21:44:43 +0000 (21:44 +0000)]
* imach.c (Repository): Replace "freqsummary" at a correct
place. It differs from routine "prevalence" which may be called
many times. Probs is memory consuming and must be used with
parcimony.
Version 0.95a2 (should output exactly the same maximization than 0.8a2)
git://henry.ined.fr / .git / log