| CT {fmsb} | R Documentation |
Coale and Trussell (1978)の年齢別有配偶出生率のモデルを実装し，年齢別有配偶出生率の実データへの当てはめを行う。
Implementing Coale and Trussell's model of age-specific marital fertility rates and fitting the model to actual ASMFR.
CT(M=1, m=0) fitCT(initialpar=c(1.0, 1.0), data, Method="Nelder-Mead", ...)
M |
The parameter M of the CT model, the scale (peak height) parameter of fertility |
m |
The parameter m of the CT model, the strength of downward discordance from natural fertility with aging |
initialpar |
Initial value for the parameters to be estimated. If not given, c(1.0, 1.0) is used. |
data |
Actual vector of ASMFR (which must be given for from age 12 to age 49 for each age) to be used to obtain the best-fit parameters of the CT's model. |
Method |
The method to be used in optim() function. Default is "Nelder-Mead". |
... |
Other options to be passed to optim(). |
CT()は12歳から49歳までの年齢各歳について，年齢別有配偶出生率を返す。fitCT()は，データにもっとも当てはまるモデルを返すパラメータMとm，平均二乗誤差の平方根（RMSE），収束したかどうかのフラグという4つの値からなるベクトルを返す。
CT() returns model ASMFR for ages from 12 to 49. fitCT() returns the numeric vector of fitted parameters M and m, RMSE for those values, and the flag of convergence.
Minato Nakazawa minato-nakazawa@people.kobe-u.ac.jp https://minato.sip21c.org/
Coale AJ, Trussell TJ (1978) Technical Note: Finding the Two Parameters That Specify a Model Schedule of Marital Fertility. Population Index, 44(2): 203-213.
URL: http://www.jstor.org/stable/2735537
ASMFR <- c(0,0,0,Jfert$ASMFR2000[1:35]) # Jfert gives ASMFR for ages 15:54, so that it should be rearranged to 12:49
res <- fitCT(,ASMFR)
FLAG <- res[4]
while (FLAG>0) {
res <- fitCT(res[1:2],ASMFR)
FLAG <- res[4]
}
print(res)