samp.postopt: Multivariate Gaussian Sampling for Heligman-Pollard model...

Description Usage Arguments Value Note References

Description

Samples the nine Heligman-Pollard parameters from the mvnorm distribution for each run of optimizer step where the likelihood for that run exceeds the maximum likelihood from the prior.

Usage

1
2
samp.postopt(opt.cov.d, opt.mu.d, d.keep, prior, B = 400, B0 = 8000,
  d = 10)

Arguments

opt.cov.d

An array containing a covariance matrix for each run of optimizer where the likelihood for that run exceeds the maximum likelihood from the prior.

opt.mu.d

A matrix containing the results of the optimizer step.

d.keep

Number of runs of optimizer where the likelihood for that run exceeds the maximum likelihood from the prior.

prior

A matrix containing the prior distribution (see HP.priors).

B

Sample size at the importance sampling stage.

B0

Sample size of the prior. This is equal to (theta.dim * 1000).

d

Number of optimizer iterations.

Value

H.k

The prior plus new samples.

H.new

The new samples from the multivariate normal.

B1

The number of new samples - should be equal to B * d.keep.

Note

For use within the function HP.mod.

References

Heligman, L. and Pollard, J.H. (1980) The Age Pattern of Mortality. Journal of the Institute of Actuaries 107:49-80.

Poole, D and Raftery, A. (2000). Inference for Deterministic Simulation Models: The Bayesian Melding Approach. Journal of the American Statistical Association 95:1244-1255.

Raftery, A and Bao, L. (2009). "Estimating and Projecting Trends in HIV/AIDS Generalized Epidemics Using Incremental Mixture Importance Sampling." Technical Report 560, Department of Statistics, University of Washington.

Sharrow, D.J., Clark, S.J., Collinson, M.A., Kahn, K. and Tollman, S.M. (2013). The Age Pattern of Increases in Mortality Affected by HIV: Bayesian Fit of the Heligman-Pollard Model to Data from the Agincourt HDSS Field Site in Rural Northeast South Africa. Demogr. Res. 29, 1039<e2><80><93>1096.


strandCet documentation built on May 1, 2019, 8:19 p.m.