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)
## Attempting to run 1 MPs:
## AvC
## TAC
## AvC 24.83
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")
## Attempting to run 1 MPs:
## AvC
## TAC
## AvC 24.54
hist(Recs@TAC)
runMP
can be used to run several MPs:
runMP(Atlantic_mackerel, c("AvC", "AvC_MLL"))
## Attempting to run 2 MPs:
## AvC
## AvC_MLL
## TAC LR5 LFR
## AvC 25.1
## AvC_MLL 25.1 90.25 95
Or all available MPs:
Atlantic_mackerel <- runMP(Atlantic_mackerel, reps=1000, chkMPs=TRUE, silent=TRUE)
The TAC recommendations from each Output control can be plotted:
boxplot(Atlantic_mackerel)
## MP Median SD Units
## 1 SP_4010 0.00 0.00 thousand tonnes
## 2 SP_75MSY 0.00 0.00 thousand tonnes
## 3 SP_MSY 0.00 0.00 thousand tonnes
## 4 DDSS_4010 0.45 0.00 thousand tonnes
## 5 DD4010 0.46 0.26 thousand tonnes
## 6 DCAC4010 1.43 1.01 thousand tonnes
## 7 Fratio4010 2.12 2.10 thousand tonnes
## 8 DD 2.12 0.82 thousand tonnes
## 9 DDSS_75MSY 2.60 0.00 thousand tonnes
## 10 DAAC 2.82 1.27 thousand tonnes
## 11 HDAAC 2.86 1.25 thousand tonnes
## 12 DepF 2.93 2.06 thousand tonnes
## 13 DDSS_MSY 3.47 0.00 thousand tonnes
## 14 Fratio 3.74 2.59 thousand tonnes
## 15 DBSRA4010 3.81 4.91 thousand tonnes
## 16 Fdem 4.63 1.98 thousand tonnes
## 17 Itarget4 4.93 0.43 thousand tonnes
## 18 DynF 4.97 2.23 thousand tonnes
## 19 DCACs 5.57 1.88 thousand tonnes
## 20 DCAC 5.66 1.89 thousand tonnes
## 21 DBSRA 5.81 4.56 thousand tonnes
## 22 SPSRA 5.84 3.61 thousand tonnes
## 23 SSS_75MSY 6.43 0.00 thousand tonnes
## 24 Gcontrol 6.53 0.03 thousand tonnes
## 25 Rcontrol 6.53 0.00 thousand tonnes
## 26 Rcontrol2 6.53 0.00 thousand tonnes
## 27 SSS_4010 6.77 0.00 thousand tonnes
## 28 Itarget3 7.08 0.64 thousand tonnes
## 29 Itarget1_MPA 7.09 0.65 thousand tonnes
## 30 Itarget1 7.10 0.61 thousand tonnes
## 31 Itarget2 7.13 0.60 thousand tonnes
## 32 DCAC_40 7.15 2.16 thousand tonnes
## 33 YPR 8.00 3.41 thousand tonnes
## 34 Fadapt 8.13 0.05 thousand tonnes
## 35 SSS_MSY 8.57 0.00 thousand tonnes
## 36 MCD4010 9.24 7.73 thousand tonnes
## 37 DBSRA_40 9.49 2.99 thousand tonnes
## 38 Iratio 10.02 4.08 thousand tonnes
## 39 BK 11.16 4.38 thousand tonnes
## 40 SPMSY 12.47 7.88 thousand tonnes
## 41 Islope3 12.70 1.14 thousand tonnes
## 42 Islope4 13.31 1.17 thousand tonnes
## 43 SPslope 13.86 1.39 thousand tonnes
## 44 MCD 13.93 5.81 thousand tonnes
## 45 ICI 13.94 3.03 thousand tonnes
## 46 ICI2 14.00 2.87 thousand tonnes
## 47 CC5 14.02 1.28 thousand tonnes
## 48 SBT1 14.45 2.94 thousand tonnes
## 49 Islope2 14.67 1.34 thousand tonnes
## 50 GB_slope 15.14 1.58 thousand tonnes
## 51 SPmod 15.86 3.87 thousand tonnes
## 52 CC4 16.34 1.43 thousand tonnes
## 53 Islope1 16.90 1.53 thousand tonnes
## 54 CurC 18.58 3.89 thousand tonnes
## 55 CC3 18.68 1.73 thousand tonnes
## 56 CC2 20.91 1.93 thousand tonnes
## 57 CC1 23.33 2.03 thousand tonnes
## 58 AvC_MLL 24.70 4.95 thousand tonnes
## 59 AvC 24.87 4.76 thousand tonnes