The function creates two files containing projection summaries, such as the median, the lower and upper bound of the 80 and 90% probability intervals, respectively, the +/- 0.5 child variant and the constant variant. One file is in a user-friendly format, whereas the other is in a UN-specific format with internal coding of the time and the variants. In addition, a file containing some of the model parameters is created.
write.projection.summary(dir = file.path(getwd(), "bayesTFR.output"),
output.dir = NULL, revision = NULL, adjusted = FALSE, est.uncertainty = FALSE, ...)
Directory containing the prediction object. It should correspond to the output.dir argument of the tfr.predict function.
Directory in which the resulting file will be stored. If NULL the same directory is used as for the prediction.
UN WPP revision number. If NULL it is determined from the corresponding WPP year: WPP 2008 corresponds to revision 13, every subsequent WPP increases the revision number by one. Used as a constant in the second file only.
Logical. By default the function writes summary using the original BHM projections. If the projection medians are adjusted (using e.g. tfr.median.set), setting this argument to TRUE causes writing the adjusted projections.
Logical. If TRUE and the simulation was generated with uncertainty around estimation, that uncertainty info is included in the summaries.
Additional arguments passed to the underlying functions. Here, argument precision can be set to determine the number of significant digits (default is 4).
Hana Sevcikova
The first file that the function creates is called projection_summary_user_friendly.csv (or projection_summary_user_friendly_adjusted.csv if adjusted=TRUE), it is a comma-separated table with the following columns:
country code
name of the variant, such as “median”, “lower 80”, “upper 80”, “lower 95”, “upper 95”, “-0.5child”, “+0.5child”, “constant”
e.g. “2005-2010”: TFR for the first time period. If est.uncertainty is TRUE, the first time period is the first observed time period. Otherwise it is the last observed one.
e.g. “2010-2015”: TFR for the second time period
further columns with TFR projections
The second file, called projection_summary.csv (or projection_summary_adjusted.csv if adjusted=TRUE), also comma-separated table, contains the same information as above in a UN-specific format:
variant identifier, extracted from the UN_variants dataset
country code
time identifier, extracted from the UN_time dataset
the total fertility rate for this variant, location and time period
The third comma-separated file, called projection_summary_parameters.csv contains the following columns:
country code
start period of TFR decline
TFR at the onset of the fertitility transition (median of the \(U_c\) parameter)
maximum decrement of TFR decline (median of the \(d_c\) parameter)
median of the end level of the fertility transition (\(\Delta_{c4}\) parameter)
2.5 percentile of the \(\Delta_{c4}\) distribution
97.5 percentile of the \(\Delta_{c4}\) distribution
period of the end decline, measured using the prediction median
Note that this file is not created if adjusted=TRUE.
convert.tfr.trajectories, tfr.predict