Siler {fmsb} | R Documentation |
生命表のq(x)についてのSilerのモデル死亡関数を実装し,それを所与のq(x)に自動的に当てはめ,最適パラメータa1, b1, a2, a3, b3を得る。
Implementing Siler's model mortality function of qx and fitting the model to actual qx of given lifetable.
Siler(a1, b1, a2, a3, b3, t) fitSiler(initialpar=c(0.01,3,1e-4,1e-5,0.1), data, mode=1, Method="Nelder-Mead", ...)
a1 |
Silerの死亡モデルq(t)=a1*exp(-b1*t)+a2+a3*exp(b3*t)のパラメータa1。 The parameter a1 of the Siler model, q(t)=a1*exp(-b1*t)+a2+a3*exp(b3*t). |
b1 |
Silerの死亡モデルq(t)=a1*exp(-b1*t)+a2+a3*exp(b3*t)のパラメータb1。 The parameter b1 of the Siler model, q(t)=a1*exp(-b1*t)+a2+a3*exp(b3*t). |
a2 |
Silerの死亡モデルq(t)=a1*exp(-b1*t)+a2+a3*exp(b3*t)のパラメータa2。 The parameter a2 of the Siler model, q(t)=a1*exp(-b1*t)+a2+a3*exp(b3*t). |
a3 |
Silerの死亡モデルq(t)=a1*exp(-b1*t)+a2+a3*exp(b3*t)のパラメータa3。 The parameter a3 of the Siler model, q(t)=a1*exp(-b1*t)+a2+a3*exp(b3*t). |
b3 |
Silerの死亡モデルq(t)=a1*exp(-b1*t)+a2+a3*exp(b3*t)のパラメータb3。 The parameter b3 of the Siler model, q(t)=a1*exp(-b1*t)+a2+a3*exp(b3*t). |
t |
年齢(ベクトルでもOK)を年単位で。 Age (vector OK) in years |
initialpar |
推定されるべきパラメータa1, b1, a2, a3, b3の初期値を長さ5のベクトルとして与える。省略すると,c(0.01, 3, 1e-04, 1e-05, 0.1)が使われる。 Initial value for the parameters to be estimated. If not given, c(0.01, 3, 1e-04, 1e-05, 0.1) is used. |
data |
Silerの死亡モデルを当てはめるべき生命表のqx関数データ。ベクトルとして与える。 Actual vector of qx in the lifetable to be used to obtain the best-fit parameters of the Siler's model. |
mode |
モデルの当てはめで最小化すべきRMSE(平均平方誤差の平方根)を,生命表関数の何を使って計算するかを指定する整数値。1ならqx,2ならdx,それ以外ならlxを用いる。デフォルトは1。 Which of lifetable functions should be used to calculate the RMSE: 1 qx, 2 dx, otherwise lx. Default is 1. |
Method |
optim()関数への引数として渡される,モデルの当てはめで関数の最小値を求める方法。デフォルトはNelder-Mead法。他に"SANN"や"BFGS"が使える。 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.