群馬大学 | 医学部 | サイトトップ | 人口学

Getting MLEs of a Denny's model survival curve

Last updated on 6 January 2012 (Fri)

References

  1. Denny, C. (1997) A model of the probability of survival from birth. Mathematical and Computer Modelling, 26: 69-78.
  2. Kobayashi, K. and Z. Nanjo (1988) Generation life tables for Japan based on period life tables covering the years 1891-1986. Nihon University Population Research Institute, Tokyo.
  3. Nelder, J.A. and R. Mead (1965) A simplex method for function minimization. The Computer Journal, 7(4): 308-313.

ダウンロード

dennys.exe (Win32自己解凍ファイル,27966 bytes),または,dennys.tgz (tar+gzip圧縮ファイル,9055 bytes)のどちらかをダウンロードしてください。

解凍すると,denny.txt(この説明とほぼ同じ),denny.c(Cで書かれたソースプログラム。gccでコンパイル可能),denny.exe(Mingw32でコンパイルしたWin32実行形式),jpcohortm.dat(小林, 南條(1988)による日本人男性のコホート生命表から作成したサンプルデータ)が得られます。

概要説明

Dennyのモデル(Denny, 1997)の3つのパラメータを推計するプログラム。入力ファイルとして,年齢x,タブ記号,生残率l(x)を1行として0歳から年齢各歳別に縦に並べたデータ(サンプルデータ参照)をコマンドラインパラメータとして与えると,データの生残率とDennyモデルでの生残率の差の二乗和の平均の平方根を算出し,それをNelder-Meadの滑降シンプレックス法(Nelder and Mead, 1965)を使って最小化するパラメータを求め,標準出力する。データによっては,局所最小に落ち込んでしまう可能性があるので,初期値を変更する必要があるかもしれない(初期値によって結果が変わる場合は,RMSEを比べて真の最小を探ればよい)。初期値の与え方は,データファイル名に続けて,コマンドラインパラメータとしてスペースで区切って与えればよい。

著作権等

revision 1.0 All right reserved by Minato Nakazawa <minato-nakazawa[atmark]umin.net>.

配布条件

基本的に無料です。配布条件は基本的にGNUのGeneral Public Licenseに従うものとお考えください。

Download

Please download either dennys.exe (Win32 self-executable archive, 27966 bytes), or dennys.tgz (tar+gzip archive, 9055 bytes).

When extracting, you may get denny.txt (explanation like this page), denny.c (the source file written in C, can be compiled by gcc), denny.exe (Win32 executable file, compiled by Mingw32), jpcohortm.dat (sample data based on the cohort life table of Japanese males, Kobayashi and Nanjo, 1988).

Description

The program estimates 3 parameters of Denny's model for survival curve (Denny, 1997), with best fit to any given data. Input data file (survival curve) must have age x, <tab>, survival rate l(x), which are aligned in each line and the lines compose whole survival curve from age 0 to eldest age (see, sample data file, jpcohortm.dat).

When running the program, if the data file is given as 1st command line parameter, then the best fit parameters will be shown in standard output, using so-called downhill simplex method for function minimization to achieve minimum root MSE (Nelder and Mead, 1965). If obtained RMSE is not small probably due to local minimum pocket, please change initial parameter values by giving 2nd, 3rd, and fourth command line parameters.

Copyright

revision 1.0 All right reserved by Minato Nakazawa <minato-nakazawa[atmark]umin.net>.

Distribution

Basically free of charge. This program basically adopts GNU GPL.


リンクと引用について