R/Rev_solow93.R

In Clara-Casabona/ExtinctModels: Models to infer the date of extinction

#' Rev Solow1993.eq2
#'
#' @param dd
#' @param g
#'
#' @return
#' @export
#'
#' @examples

Sol93_col <- function(dd,g){
  names(dd) <-c("yrs", "sights")
  Delta_t <- max(dd$yrs) - min(dd$yrs)
  n <- sum(na.omit(dd$sights))
  Max_year <-  g - min(dd$yrs)
  eq <- (1/(1-(Max_year/Delta_t)))^n
  return(eq)
}