R/mod.nM.R

In milokmilo/Stranded: Calculate Biological Parameters from Stranded Cetaceans

#' Helligman Pollard mortality function 
#'
#' Helligman Pollard mortality function. Modified from HPbayes_0.1
#' @param age
#' @param hpp
#' @keywords Heligman Pollard mortality
#' @export
#' @examples
#' mod.nM()

mod.nM <- function (theta, x) {
  A <- theta[1]
  B <- theta[2]
  C <- theta[3]
  G <- theta[4]
  H <- theta[5]
  f.x <- A^((x + B)^C) +
    (G * (H^x))/(1 + G * (H^x))
  f.x2 <- f.x
  f.x2[1e-06 > f.x] <- 1e-06
  f.x2[0.999999 < f.x] <- 0.999999
  return(f.x2)
}