R/hp.nqx.R

In HPbayes: Heligman Pollard mortality model parameter estimation using Bayesian Melding with Incremental Mixture Importance Sampling

Defines functions hp.nqx

Documented in hp.nqx

hp.nqx <-
function(H.out, age = seq(0, 85, 1)) 
{
    H.new.hat <- matrix(NA, nrow = nrow(H.out), ncol = length(age))
    for (i in 1:nrow(H.out)) {
        H.new.hat[i, ] <- mod8p(theta = H.out[i, ], x = age)
    }
    ans <- H.new.hat
    return(ans)
}

Try the HPbayes package in your browser