Produces an 8 panel plot of both the prior and posterior distribution for the eight parameters of the Heligman-Pollard model
● Data Source:
CranContrib
● Keywords: misc
● Alias: postpri.plot
●
2 images
|
Forms a prior distribution (drawn from a uniform distribution) for each of the eight Heligman-Pollard parameters. First, using optim() , mle estimates of the parameters are fitted to the deaths and persons at risk supplied by the user. Once these estimates (returned as mle ) and their standard errors (se ) are obtained, 8000 (See documentation for prior.form ) draws from a uniform distribution with bounds mle[i] +/- senum*se are taken to form a prior distribution for each parameter.
● Data Source:
CranContrib
● Keywords: misc
● Alias: pri.mle
●
0 images
|
Calculates the log-likelihood and importance weight for each set (i.e. each row) of Heligman-Pollard parameters in the prior
● Data Source:
CranContrib
● Keywords: misc
● Alias: prior.likewts
●
0 images
|
Performs the optimizer step in the IMIS procedure for the eight Heligman-Pollard parameters
● Data Source:
CranContrib
● Keywords: misc
● Alias: loop.optim
●
0 images
|
Performs the final re-sampling step in the Bayesian Melding with IMIS procedure for the eight Heligman-Pollard parameters.
● Data Source:
CranContrib
● Keywords: misc
● Alias: final.resamp
●
0 images
|
samp.postopt
(Package: HPbayes) :
Multivariate Gaussian Sampling for Heligman-Pollard model estimated via Bayesian Melding
Samples the eight 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
● Data Source:
CranContrib
● Keywords: misc
● Alias: samp.postopt
●
0 images
|
This function converts the posterior Heligman-Pollard parameter distribution resulting from a Bayesian Melding procedure to probabilities of death over a specified age range. In addition, this function also calculates and plots a CI over the specified age range. The user can also elect to plot the data (ndeath /nrisk ) from which the model was estimated.
● Data Source:
CranContrib
● Keywords: misc
● Alias: hpbayes.plot
●
2 images
|
Converts a set of Heligman-Pollard mortality model parameters into age-specific probabilities of death
● Data Source:
CranContrib
● Keywords: misc
● Alias: hp.nqx
●
0 images
|
Calculates a log likelihood from a binomial distribution
● Data Source:
CranContrib
● Keywords: misc
● Alias: ll.binom
●
0 images
|
Calculates the log-likelihood based on a binomial distribution of observing a set of age-specific probabilities of death resulting from a set of eight parameters of the Heligman-Pollard Model given a certain number of age-specific death counts and persons at risk
● Data Source:
CranContrib
● Keywords: misc
● Alias: mp8.ll
●
0 images
|