Siler {fmsb}R Documentation

Silerの死亡モデルを生命表関数q(x)に当てはめる
Siler's model mortality for qx and its fitting

概要 (Description)

生命表の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.

使用法 (Usage)

 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", ...)

引数 (Arguments)

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().

返り値 (Value)

Siler()
Siler()関数は,tと同じ長さをもつqxのベクトルを返す。
Siler() returns model qx for the same length with t.
fitSiler()
fitSiler()は,optim()関数によって得られた,データに最も当てはまるSiler()のパラメータa1, b1, a2, a3, b3と,そのときのRMSE(平均平方誤差の平方根)と,収束フラグ(0が収束したことを示し,それ以外は収束しなかったことを示す)
fitSiler() returns the numeric vector of fitted parameters a1, b1, a2, a3 and b3, RMSE for those values, and the flag of convergence.

参照 (See Also)

日本の完全生命表データJlife

Jlife for complete life tables in Japan.

用例 (Examples)

qx of 1950 and 2000 males' complete life tables in Japan with Siler's model
[Package fmsb Index]