GoF.dl: Goodness of Fit of the Double Logistic Function

In jasonacollins/bayesTFR1: Bayesian Fertility Projection

Description

The function computes coverage, i.e. the ratio of observed data fitted within the given probability intervals of the predictive posterior distribution of the double logistic function, as well as the root mean square error and mean absolute error of the simulation.

Usage

tfr.dl.coverage(sim.dir, pi = c(80, 90, 95), burnin = 2000, verbose = TRUE)

Arguments

sim.dir	Directory with the MCMC simulation results. If a prediction and its corresponding thinned MCMCs are available in the simulation directory, those are taken for assessing the goodness of fit.
pi	Probability interval. It can be a single number or an array.
burnin	Burnin. Only relevant if sim.dir does not contain thinned chains.
verbose	Logical switching log messages on and off.

Value

List with the following components:

total.coverage	Vector of the coverage, one element per probability interval. For each pi, it is the ratio of the number of observed data points that fall within the probability interval of the posterior distribution over the total number of data points, i.e. TFR for all countries and historical time periods.
time.coverage	Matrix corresponding to the coverage computed per time period. (Rows correspond to probability intervals, columns correspond to time.) It is derived like total.coverage except that both, the nominator and denominator, contain only data points belonging to the corresponding time period.
country.coverage	Matrix corresponding to the coverage computed per country. (Rows correspond to probability intervals, columns correspond to countries.) It is derived like total.coverage except that both, the nominator and denominator, contain only data points belonging to the corresponding country.
total.rmse	Root mean square error as √{(1/n∑(x-m)^2)} where x are observed data points, m is the mean of the posterior distribution and n is the number of data points. Here the sum is taken over all countries and historical time periods.
time.rmse	Like total.rmse except that each time period is considered separately.
country.rmse	Like total.rmse except that each country is considered separately.
total.mae	Mean absolute error as 1/n∑|x-m| where x are observed data points, m is the median of the posterior distribution and n is the number of data points. Here the sum is taken over all countries and historical time periods.
time.mae	Like total.mae except that each time period is considered separately.
country.mae	Like total.mae except that each country is considered separately.
pred.cdf	TxC matrix (with T being the number of time periods and C being the number of countries), containing the predictive CDF of the observation, i.e. the quantile of each data point within the predictive posterior distribution.
n	0-1 TxC matrix indicating if the corresponding data point was included in the goodness of fit computation. Zeros indicate missing historical values.

Note

To see the fit visually per country, use DLcurve.plot(..., predictive.distr=TRUE,...).

Author(s)

Hana Sevcikova

See Also

DLcurve.plot

Examples

## Not run: 
sim.dir <- file.path(find.package("bayesTFR"), "ex-data", "bayesTFR.output")
tfr <- get.tfr.mcmc(sim.dir)
# Note that this simulation is a toy example and thus has not converged.
gof <- tfr.dl.coverage(sim.dir)
gof$time.coverage
DLcurve.plot(tfr, country=608, predictive.distr=TRUE, pi=c(80, 90, 95))

## End(Not run)

jasonacollins/bayesTFR1 documentation built on May 23, 2019, 7:34 a.m.