R/hcrSBTD.R

In laurieKell/tMSE: Tools for conducting Management Strategy Evaluation

hcrSBTD<-function(yrs,
                  control=c(k1=2.0,k2=3.0,gamma=1),
                  index,
                  catch,...){

  lambda=as.FLQuant(ddply(as.data.frame(index%/%apply(index,6,mean)), 
                          .(iter), with,  data.frame(data=coefficients(lm(data~year))[2])))
  flag  =lambda<0
  lambda=abs(lambda)
  res   =1+ifelse(flag,-control["k1"],control["k2"])*exp(log(lambda)*ifelse(flag,control["gamma"],1))
  res   =res%*%apply(catch,6,mean)

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

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

  return(res)}