Hadwiger {fmsb}R Documentation

Hadwigerのモデル出生力とその当てはめ
Hadwiger's model fertility and its fitting

説明
Description

Hadwigerの年齢別出生率(ASFR)のモデルを実装し,実際のASFRに当てはめる(自動的に,最も当てはまるパラメータを求める)。

Implementing Hadwiger's model of age-specific fertility rates and fitting the model to actual ASFR.

使用法
Usage

 Hadwiger(a=3.4, b=2.5, c=22.2)
 fitHad(initialpar=c(3.4, 2.5, 22.2), data, Method="Nelder-Mead", ...)

引数
Arguments

a

Hadwigerの年齢別出生率のモデルASFR(x)=a*b/c*(c/x)^1.5*exp(-b^2*(c/x+x/c-2)),年齢xは15歳から54歳まで,のパラメータa

The parameter a of the Hadwiger model, ASFR(x) = a*b/c*(c/x)^1.5*exp(-b^2*(c/x+x/c-2)) for age x from 15 to 54.

b

Hadwigerの年齢別出生率のモデルASFR(x)=a*b/c*(c/x)^1.5*exp(-b^2*(c/x+x/c-2)),年齢xは15歳から54歳まで,のパラメータb

The parameter b of the Hadwiger model, ASFR(x) = a*b/c*(c/x)^1.5*exp(-b^2*(c/x+x/c-2)) for age x from 15 to 54.

c

Hadwigerの年齢別出生率のモデルASFR(x)=a*b/c*(c/x)^1.5*exp(-b^2*(c/x+x/c-2)),年齢xは15歳から54歳まで,のパラメータc。概ね初産年齢を示すと考えられている。

The parameter c of the Hadwiger model, ASFR(x) = a*b/c*(c/x)^1.5*exp(-b^2*(c/x+x/c-2)) for age x from 15 to 54. It may correspond to the age of first birth.

initialpar

推定すべきパラメータa, b, cの初期値。要素3個の数値ベクトルとして与える。与えないと,c(3.4, 2.5, 22.2)が使われる。

Initial value for the parameters to be estimated. If not given, c(3.4, 2.5, 22.2) is used.

data

Hadwigerのモデルを当てはめて最適パラメータを計算する対象となる,実際のASFRデータ(15~54歳まで年齢各歳に与えられる必要がある)。

Actual vector of ASFR (which must be given for from ages from 15 to 54 for each age) to be used to obtain the best-fit parameters of the Hadwiger's model.

Method

optim()関数で使われる関数の最小化の方法。デフォルトは"Nelder-Mead"である。

The method to be used in optim() function. Default is "Nelder-Mead".

...

Other options to be passed to optim().

返り値
Value

Hadwiger()関数は,年齢15歳から54歳までの女性の年齢別出生率の数値ベクトルを返す。fitHad()関数は,データに最も良く当てはまるHadwiger()関数のパラメータa, b, cと平均平方誤差の平方根(RMSR)と,収束判定フラグ(0が収束,それ以外は非収束)という5つの要素からなる数値ベクトルを返す。

Hadwiger() returns model ASFR for ages from 15 to 54. fitHad() returns the numeric vector of fitted parameters a, b and c, RMSE for those values, and the flag of convergence.

Author(s)

Minato Nakazawa (minato-nakazawa[atmark]umin.net)

文献
References

Chandola T, Coleman DA, Horns RW (1999) Recent European fertility patterns: fitting curves to 'distorted' distributions. Population Studies, 53(3): 317-329.

URL: http://www.jstor.org/stable/2584702

See Also

Jfert

Examples

 res <- fitHad(,Jfert$ASFR2000)
 FLAG <- res[5]
 while (FLAG>0) {
   res <- fitHad(res[1:3],Jfert$ASFR2000)
   FLAG <- res[5]
 }
 print(res)

[Package fmsb Index]