Denny {fmsb} | R Documentation |
生命表のlxについてのDennyのモデル死亡関数を定義し,それを所与の生命表のlxに自動的に当てはめ,最適パラメータa, b, cを得る。l(0)は100000であることに注意。
Implementing Denny's model mortality function of lx and fitting the model to actual lx of given lifetable, then obtaining the parameters a, b and c for the best fit. Note that l(0) is 100000.
Denny(a, b, c, t) fitDenny(initialpar=rep(0.1, 3), data, mode=3, Method="Nelder-Mead", ...)
a |
Dennyの死亡モデルl(t)=1/(1+a*(t/(105-t))^3+b*sqrt(exp(t/(105-t))-1)+c*(1-exp(-2*t)))のパラメータa The parameter a of the Denny model, l(t)=1/(1+a*(t/(105-t))^3+b*sqrt(exp(t/(105-t))-1)+c*(1-exp(-2*t))). |
b |
Dennyの死亡モデルl(t)=1/(1+a*(t/(105-t))^3+b*sqrt(exp(t/(105-t))-1)+c*(1-exp(-2*t)))のパラメータb The parameter b of the Denny model, l(t)=1/(1+a*(t/(105-t))^3+b*sqrt(exp(t/(105-t))-1)+c*(1-exp(-2*t))). |
c |
Dennyの死亡モデルl(t)=1/(1+a*(t/(105-t))^3+b*sqrt(exp(t/(105-t))-1)+c*(1-exp(-2*t)))のパラメータc The parameter c of the Denny model, l(t)=1/(1+a*(t/(105-t))^3+b*sqrt(exp(t/(105-t))-1)+c*(1-exp(-2*t))). |
t |
年齢(ベクトルで与えることができる)を年単位で与える。tの要素は105未満でなくてはならない。そうでないと,Denny()関数の返値は0になる。 Age (vector OK) in years. The t must be less than 105, otherwise the value by Denny() become 0. |
initialpar |
推定されるべき初期パラメータ。長さ3のベクトルとして,a, b, cの初期値を与える。与えられないとrep(0.1, 3),つまりc(0.1, 0.1, 0.1)が使われる。 Initial value for the parameters to be estimated. If not given, rep(0.1, 3) is used. |
data |
データとしてDennyのモデルを当てはめるのに使われる,生命表関数qxのベクトルをデータとして与える。ただし,105歳以上のqxが入っていても(つまり,105番目以降の要素があっても)当てはめでは無視される。 Actual vector of qx in the lifetable to be used to obtain the best-fit parameters of the Denny's model. If the ages for qx are equal or elder than 105 years old, those will be ignored in fitting. |
mode |
当てはめで平均平方誤差の平方根(RMSE)を計算する際に,生命表関数のどれにデータを当てはめるかを番号で指定する。1ではqx,2だとdx,それ以外だとlxへの当てはめを行う。デフォルトは3なので,lxへの当てはめを行う。 Which of lifetable functions should be used to calculate the RMSE: 1 qx, 2 dx, otherwise lx. Default is 3. |
Method |
平均平方誤差の平方根(RMSE)の最小化をするために,関数optim()に渡す,関数の最小化の方法を文字列で指定する。デフォルトは"Nelder-Mead"である。 The method to be used in optim() function. Default is "Nelder-Mead". |
... |
その他オプションはoptim()関数にそのまま渡される。 Other options to be passed to optim(). |
日本の完全生命表データJlife
Jlife for complete life tables in Japan.