Description Usage Arguments Value Note References See Also Examples
View source: R/prior.likewts.R
Calculates the log-likelihood and importance weight for each set (i.e. each row) of Heligman-Pollard parameters in the prior
1 2 | prior.likewts(prior, nrisk, ndeath, theta.dim = 8,
age = c(1e-05, 1, seq(5, 100, 5)))
|
prior |
A ( |
nrisk |
A vector containing the number of persons at risk in each age group |
ndeath |
A vector containing the number of deaths in each age group |
theta.dim |
Number of columns of the prior matrix. Defaults to 8 |
age |
A vector containing the ages at which each age interval begins |
wts.0 |
A vector containing an importance weight for each set of parameters from the prior |
log.ilke.0 |
A vector containing a log likelihood for each set of parameters from the prior |
Used in the loop.optim
function
Heligman, Larry and John H. Pollard. 1980 "The Age Pattern of Mortality." Journal of the Institute of Actuaries 107:49–80.
Poole, David and Adrian Raftery. 2000. "Inference for Deterministic Simulation Models: The Bayesian Melding Approach." Journal of the American Statistical Association 95:1244–1255.
Raftery, Adrian and Le Bao. 2009. "Estimating and Projecting Trends in HIV/AIDS Gen- eralized Epidemics Using Incremental Mixture Importance Sampling." Technical Report 560, Department of Statistics, University of Washington.
1 2 3 4 5 6 7 8 9 | lx <- c(1974, 1906, 1860, 1844, 1834, 1823, 1793, 1700, 1549, 1361,
1181, 1025, 870, 721, 571, 450, 344, 256, 142, 79, 41, 8)
dx <- c(68, 47, 16, 10, 13, 29, 92, 151, 188, 179, 156, 155, 147, 150,
122, 106, 88, 113, 63, 38, 32, 8)
age <- c(1e-05, 1, seq(5, 85, 5))
q0 <- prior.form()
result <- prior.likewts(prior=q0, nrisk=lx, ndeath=dx, age=age)
summary(result$wts.0)
summary(result$log.like.0)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.