R/hcr_implementation_model_1.R

In einarhjorleifsson/mac: Miscalaneous set of functions related fish stock assessment work

#' @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)
}