R/hp.nqx9.R

In milokmilo/Stranded: Calculate Biological Parameters from Stranded Cetaceans

#'  Helligman Pollard mortality
#'
#' Helligman Pollard mortality. Modified from hp.nqx.R of HPbayes_0.1
#' @param H.out
#' @param age. Default = seq(0, 85, 1)
#' @keywords Heligman Pollard mortality
#' @export
#' @examples
#' hp.nqx9()

hp.nqx9 <- 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, ] <- mod9p(theta = H.out[i, ], x = age)
  }
  ans <- H.new.hat
  return(ans)
}