This build method fits fertility rates to Chandola, Coleman and Hiorns Hadwiger mixture model \begin{equation} \begin{array}{lcll} \frate(x)&=&R_1\left(\frac{b_1}{c_1}\right)\left(\frac{c_1}{x}\right)^{3/2} \exp\left[-b_1^2\left (\frac{c_1}{x}+\frac{x}{c_1}-2\right)\right] &\eqnote{non-marital}\\ \\ &&+R_2\left(\frac{b_2}{c_2}\right)\left(\frac{c_2}{x}\right)^{3/2} \exp\left[-b_2^2\left (\frac{c_2}{x}+\frac{x}{c_2}-2\right)\right] & \eqnote{marital} \end{array} \end{equation}
Fitting is done using a nonlinear optimisation routine. Its default behaviour can be changed by appending an extrainfo string to the buildmethod value, HADWIGER_MIX:ExtraInfo. Here ExtraInfo is a comma-separated list of quoted strings of the form "Parameter=Value".
The fitting routine proceeds iteratively from a set of starting values for the parameters. The user can set the starting value in the extrainfo string. For example, "C2=33" sets the starting value of \(c_2\) to 33.
Fit flags can be used to specify whether a parameter is adjusted to fit the fertility probabilities or held fixed at its starting value. For example, "FitC2=0" specifies \(c_2\) to be fixed at its starting value, and "FitC2=1" specifies \(c_2\) to be adjusted. The default behaviour is to fit all parameters. Entire components can be held fixed using the parameters FitNonMarital and FitMarital. For example the extrainfo "FitMarital=0","R2=0" can be used to remove the entire marital component from the fit.
Given a set of observed fertility rates \({_n\ofrate_x}\) the build method searches for parameters that minimise the objective function \begin{equation} \sse=\sum_x \FitError({_n\ofrate_x},{_n\frate_x}) . \end{equation} The form of the function \(\FitError(\ofrate,\frate)\) is controlled by the ErrType extrainfo property \begin{equation} \FitError(\ofrate,\frate)=\left\{% \begin{array}{ll} (\frate - \ofrate)^2 & \eqnote{ABS}\\ \weight(\frate - \ofrate\log\frate) & \eqnote{POISSON}\\ \end{array}\right. \end{equation} The value of the the weight \(\weight\) is controlled by the PCurve extrainfo property. If this is not specified then it defaults to \(\weight=1\). The default error type is "ErrType=ABS".
The error type POISSON is for use when the observed fertilities \(_n\ofrate_x\) display appreciable sample noise. The weight \(\weight\) then corresponds to the population in age group \([x,x+n)\) at risk. An extrainfo string of the form "PCurve=PCurveHandle" specifies that these weights are to be taken from the PCurve object PCurveHandle.