R/hcr_implementation_model_1.R

In einarhjorleifsson/fpress2: A collection functions for hcr evaluation

#' @title Add implementation error to TAC
#'
#' @description Dummy function
#'
#' @export
#'

hcr_implementation_model_1 <- function () #(TACy2, TACy1, beta1, beta2)
{
  #i <- TACy2 < TACy1
  #if (sum(i) > 0)
  #  TACy2[i] <- TACy2[i] * (TACy1[i]/TACy2[i])^(rbeta(1,beta1, beta2))
  #return(TACy2)
}