| Hadwiger {fmsb} | R Documentation |
Hadwigerの年齢別出生率(ASFR)のモデルを実装し，実際のASFRに当てはめる（自動的に，最も当てはまるパラメータを求める）。
Implementing Hadwiger's model of age-specific fertility rates and fitting the model to actual ASFR.
Hadwiger(a=3.4, b=2.5, c=22.2) fitHad(initialpar=c(3.4, 2.5, 22.2), data, Method="Nelder-Mead", ...)
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(). |
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.
Minato Nakazawa minato-nakazawa@people.kobe-u.ac.jp https://minato.sip21c.org/
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
res <- fitHad(,Jfert$ASFR2000)
FLAG <- res[5]
while (FLAG>0) {
res <- fitHad(res[1:3],Jfert$ASFR2000)
FLAG <- res[5]
}
print(res)