CT {fmsb}R Documentation

Coale and Trussellの有配偶出生モデルと,その当てはめ
Coale and Trussell's model marital fertility and its fitting

説明
Description

Coale and Trussell (1978)の年齢別有配偶出生率のモデルを実装し,年齢別有配偶出生率の実データへの当てはめを行う。

Implementing Coale and Trussell's model of age-specific marital fertility rates and fitting the model to actual ASMFR.

使用法
Usage

 CT(M=1, m=0)
 fitCT(initialpar=c(1.0, 1.0), data, Method="Nelder-Mead", ...)

引数
Arguments

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

返り値
Value

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.

Author(s)

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

参考文献
References

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

See Also

Jfert

Examples

 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)

[Package fmsb Index]