Chapter 22 Applying Management Procedures

The runMP function can be used to apply a MP to a Data object. For example, to apply the AvC method:

runMP(Atlantic_mackerel, "AvC", reps=1000)
##       TAC
## AvC 24.33

The runMP prints out the MP recommendations to the console. In the case of a TAC, where multiple repititions where (see reps = 1000 above) used the runMP function prints the median TAC recommendation.

Although it only displays a summary, runMP invisibly returns the Data object with the TAC slot populated:

Recs <- runMP(Atlantic_mackerel, "AvC")
##       TAC
## AvC 24.76
hist(Recs@TAC)

runMP can be used to run several MPs:

runMP(Atlantic_mackerel, c("AvC", "AvC_MLL"))
##           TAC   LR5 LFR
## AvC     24.54          
## AvC_MLL 23.88 90.25  95

Or all available MPs:

Atlantic_mackerel <- runMP(Atlantic_mackerel, reps=1000)
## running all available MPs
## Warning in rlnorm(reps, mconv(mu, mu * cv), sdconv(mu, mu * cv)): NAs produced

## Warning in rlnorm(reps, mconv(mu, mu * cv), sdconv(mu, mu * cv)): NAs produced

## Warning in rlnorm(reps, mconv(mu, mu * cv), sdconv(mu, mu * cv)): NAs produced

## Warning in rlnorm(reps, mconv(mu, mu * cv), sdconv(mu, mu * cv)): NAs produced

## Warning in rlnorm(reps, mconv(mu, mu * cv), sdconv(mu, mu * cv)): NAs produced

## Warning in rlnorm(reps, mconv(mu, mu * cv), sdconv(mu, mu * cv)): NAs produced

## Warning in rlnorm(reps, mconv(mu, mu * cv), sdconv(mu, mu * cv)): NAs produced

## Warning in rlnorm(reps, mconv(mu, mu * cv), sdconv(mu, mu * cv)): NAs produced

## Warning in rlnorm(reps, mconv(mu, mu * cv), sdconv(mu, mu * cv)): NAs produced

## Warning in rlnorm(reps, mconv(mu, mu * cv), sdconv(mu, mu * cv)): NAs produced

## Warning in rlnorm(reps, mconv(mu, mu * cv), sdconv(mu, mu * cv)): NAs produced

## Warning in rlnorm(reps, mconv(mu, mu * cv), sdconv(mu, mu * cv)): NAs produced
## Method CC1 produced greater than 50% NA values
## Warning in rlnorm(reps, mconv(mu, mu * cv), sdconv(mu, mu * cv)): NAs produced
## Method CC2 produced greater than 50% NA values
## Warning in rlnorm(reps, mconv(mu, mu * cv), sdconv(mu, mu * cv)): NAs produced
## Method CC3 produced greater than 50% NA values
## Warning in rlnorm(reps, mconv(mu, mu * cv), sdconv(mu, mu * cv)): NAs produced
## Method CC4 produced greater than 50% NA values
## Warning in rlnorm(reps, mconv(mu, mu * cv), sdconv(mu, mu * cv)): NAs produced
## Method CC5 produced greater than 50% NA values
## Warning in rlnorm(reps, mconv(mu, mu * cv), sdconv(mu, mu * cv)): NAs produced
## Method CurC produced greater than 50% NA values
## Warning in rlnorm(reps, mconv(mu, mu * cv), sdconv(mu, mu * cv)): NAs produced
## Method Islope1 produced greater than 50% NA values
## Warning in rlnorm(reps, mconv(mu, mu * cv), sdconv(mu, mu * cv)): NAs produced
## Method Islope2 produced greater than 50% NA values
## Warning in rlnorm(reps, mconv(mu, mu * cv), sdconv(mu, mu * cv)): NAs produced
## Method Islope3 produced greater than 50% NA values
## Warning in rlnorm(reps, mconv(mu, mu * cv), sdconv(mu, mu * cv)): NAs produced
## Method Islope4 produced greater than 50% NA values
## Method Rcontrol produced greater than 50% NA values
## Warning in rlnorm(reps, mconv(mu, mu * cv), sdconv(mu, mu * cv)): NAs produced
## Method SBT1 produced greater than 50% NA values
##                TAC Effort   LR5    LFR     HS Area 1 Area 2
## AvC          24.70                                         
## AvC_MLL      24.44        90.25  95.00                     
## BK           11.13                                         
## CC1          23.74                                         
## CC2          20.48                                         
## CC3          18.82                                         
## CC4          16.50                                         
## CC5          14.09                                         
## CurC         18.28                                         
## curE                 1.00                                  
## curE75               0.75                                  
## DAAC          2.87                                         
## DBSRA         6.22                                         
## DBSRA_40      9.83                                         
## DBSRA4010     3.97                                         
## DCAC          5.70                                         
## DCAC_40       7.16                                         
## DCAC4010      1.48                                         
## DCACs         5.81                                         
## DD            1.99                                         
## DD4010        0.45                                         
## DDe                  0.09                                  
## DDe75                0.07                                  
## DDes                 0.90                                  
## DepF          2.88                                         
## DTe40                0.90                                  
## DTe50                0.90                                  
## DynF          4.92                                         
## Fadapt        8.13                                         
## Fdem          4.70                                         
## Fratio        3.70                                         
## Fratio4010    2.24                                         
## GB_slope     15.14                                         
## Gcontrol      6.53                                         
## HDAAC         2.86                                         
## ICI          13.75                                         
## ICI2         13.98                                         
## Iratio       10.24                                         
## Islope1      16.90                                         
## Islope2      15.22                                         
## Islope3      12.47                                         
## Islope4      13.23                                         
## Itarget1      7.10                                         
## Itarget1_MPA  7.03                                 0      1
## Itarget2      7.08                                         
## Itarget3      7.09                                         
## Itarget4      4.94                                         
## ItargetE1            0.85                                  
## ItargetE2            0.85                                  
## ItargetE3            0.85                                  
## ItargetE4            0.85                                  
## matlenlim                 90.25  95.00                     
## matlenlim2                99.28 104.50                     
## MCD          14.29                                         
## MCD4010       8.88                                         
## minlenLopt1               89.75  99.72                     
## MRnoreal                                           0      1
## MRreal                                             0      1
## NFref         0.01                                         
## Rcontrol      6.53                                         
## Rcontrol2     6.53                                         
## SBT1         14.16                                         
## slotlim                   99.28 104.50 142.55              
## SPmod        15.90                                         
## SPMSY        12.04                                         
## SPslope      13.85                                         
## SPSRA         5.75                                         
## YPR           8.09                                         
## avgMP         8.62                                         
## TCPUE_e              0.95                                  
## THC          42.98

The TAC recommendations from each Output control can be plotted:

boxplot(Atlantic_mackerel)

##              MP Median   SD           Units
## 1        DD4010   0.45 0.28 thousand tonnes
## 2      DCAC4010   1.48 1.01 thousand tonnes
## 3            DD   1.99 0.85 thousand tonnes
## 4    Fratio4010   2.24 2.08 thousand tonnes
## 5         HDAAC   2.86 1.30 thousand tonnes
## 6          DAAC   2.87 1.31 thousand tonnes
## 7          DepF   2.88 2.09 thousand tonnes
## 8        Fratio   3.70 2.27 thousand tonnes
## 9     DBSRA4010   3.97 5.27 thousand tonnes
## 10         Fdem   4.70 1.83 thousand tonnes
## 11         DynF   4.92 2.14 thousand tonnes
## 12     Itarget4   4.94 0.46 thousand tonnes
## 13         DCAC   5.70 1.86 thousand tonnes
## 14        SPSRA   5.75 3.37 thousand tonnes
## 15        DCACs   5.81 1.84 thousand tonnes
## 16        DBSRA   6.22 4.42 thousand tonnes
## 17     Gcontrol   6.53 0.04 thousand tonnes
## 18     Rcontrol   6.53 0.00 thousand tonnes
## 19    Rcontrol2   6.53 0.00 thousand tonnes
## 20 Itarget1_MPA   7.03 0.62 thousand tonnes
## 21     Itarget2   7.08 0.65 thousand tonnes
## 22     Itarget3   7.09 0.63 thousand tonnes
## 23     Itarget1   7.10 0.64 thousand tonnes
## 24      DCAC_40   7.16 2.07 thousand tonnes
## 25          YPR   8.09 3.59 thousand tonnes
## 26       Fadapt   8.13 0.05 thousand tonnes
## 27        avgMP   8.62 1.94 thousand tonnes
## 28      MCD4010   8.88 8.07 thousand tonnes
## 29     DBSRA_40   9.83 3.00 thousand tonnes
## 30       Iratio  10.24 4.14 thousand tonnes
## 31           BK  11.13 4.16 thousand tonnes
## 32        SPMSY  12.04 7.96 thousand tonnes
## 33      Islope3  12.47 1.21 thousand tonnes
## 34      Islope4  13.23 1.15 thousand tonnes
## 35          ICI  13.75 2.79 thousand tonnes
## 36      SPslope  13.85 1.46 thousand tonnes
## 37         ICI2  13.98 2.89 thousand tonnes
## 38          CC5  14.09 1.32 thousand tonnes
## 39         SBT1  14.16 3.29 thousand tonnes
## 40          MCD  14.29 6.23 thousand tonnes
## 41     GB_slope  15.14 1.52 thousand tonnes
## 42      Islope2  15.22 1.45 thousand tonnes
## 43        SPmod  15.90 3.57 thousand tonnes
## 44          CC4  16.50 1.39 thousand tonnes
## 45      Islope1  16.90 1.37 thousand tonnes
## 46         CurC  18.28 3.88 thousand tonnes
## 47          CC3  18.82 1.41 thousand tonnes
## 48          CC2  20.48 1.94 thousand tonnes
## 49          CC1  23.74 2.07 thousand tonnes
## 50      AvC_MLL  24.44 5.01 thousand tonnes
## 51          AvC  24.70 5.10 thousand tonnes
## 52          THC  42.98 4.33 thousand tonnes