R/hcrSBTP.R

In flr/mydas: Tools for conducting Management Strategy Evaluation for data limited stocks

#' hcrSBTP
#' 
#' @title hcrSBTP
#' 
#' @description 
#' @author Laurence Kell, Sea++
#' 
#' @param yrs  blah,blah,blah,...
#' @param control  blah,blah,blah,...
#' @param ref  blah,blah,blah,...
#' @param target  blah,blah,blah,...
#' @param cpue  blah,blah,blah,...
#' @param catch  blah,blah,blah,...
#' @param wt default 0.5  blah,blah,blah,...
#'  
#' @name hcrSBTP
#' 
#' @export hcrSBTP
#' @docType methods
#' 
#' @rdname hcrSBTP
#' 
#' @examples
#' \dontrun{
#' data(pl4)
#' }

hcrSBTP<-function(yrs,
                  control,
                  ref,
                  target,
                  cpue,
                  catch,
                  wt=0.5){

  #SBT  
  coef       =1+control["k1"]
  flag       =c(cpue<ref)
  coef[,flag]=1-control["k2",flag]
  bit        =target%*%((cpue%/%ref)%^%(coef))
  res        =catch*(1-wt)+bit*(wt)

  dmns=dimnames(catch)
  dmns$year=yrs
  dmns$iter=dimnames(cpue)$iter

  res=FLQuant(rep(c(res),each=length(yrs)),dimnames=dmns)

  return(res)}