Management Procedures return S4 objects of class Rec containing the management recommendations returned by an MP.
Rec objects have the following slots:
slotNames('Rec')## [1] "TAC" "Effort" "Spatial" "Allocate" "LR5" "LFR"
## [7] "HS" "Rmaxlen" "L5" "LFS" "Vmaxlen" "Fdisc"
## [13] "DR" "Misc"Information on the contents the Rec object is available in the help documentation or by typing class?Rec in the R console.
We demonstrate the Rec object by applying several management procedures of different types to the Cobia example Data object.
First, we use the Can function to determine which MPs can be run on the Cobia data:
MPs <- Can(Cobia, silent=TRUE)
MPtype(MPs)## MP Type Recs
## 1 curE Input TAE
## 2 curE75 Input TAE
## 3 ItargetE1 Input TAE
## 4 ItargetE2 Input TAE
## 5 ItargetE3 Input TAE
## 6 ItargetE4 Input TAE
## 7 matlenlim Input Retention
## 8 matlenlim2 Input Retention
## 9 MRnoreal Input Spatial
## 10 MRreal Input Spatial
## 11 slotlim Input Retention
## 12 AvC_MLL Mixed TAC, Retention
## 13 Itarget1_MPA Mixed TAC, Spatial
## 14 AvC Output TAC
## 15 BK_CC Output TAC
## 16 CC1 Output TAC
## 17 CC2 Output TAC
## 18 CC3 Output TAC
## 19 CC4 Output TAC
## 20 CC5 Output TAC
## 21 CurC Output TAC
## 22 GB_slope Output TAC
## 23 ICI Output TAC
## 24 ICI2 Output TAC
## 25 Iratio Output TAC
## 26 Islope1 Output TAC
## 27 Islope2 Output TAC
## 28 Islope3 Output TAC
## 29 Islope4 Output TAC
## 30 Itarget1 Output TAC
## 31 Itarget2 Output TAC
## 32 Itarget3 Output TAC
## 33 Itarget4 Output TAC
## 34 SBT1 Output TAC
## 35 SPmod Output TAC
## 36 SPMSY Output TAC
## 37 curEref Reference TAE
## 38 NFref Reference TACNext we apply several MPs in turn to the Cobia data.
Effort Control MP
Apply ItargetE1 MP to Cobia:
Rec <- ItargetE1(1, Cobia)
class(Rec)## [1] "Rec"
## attr(,"package")
## [1] "MSEtool"Rec## Effort
## 0.85The MP returns an object of class Rec. By default, printing the Rec object to the console will only show the slots that are populated, in this case the Effort slot is the only slot that has been populated by the MP:
Rec@Effort## [1] 0.85Rec@TAC## numeric(0)Size Limit MP
Apply matlenlim MP which sets the length-at-retention to the length-at-maturity:
Rec <- matlenlim(1, Cobia)
Rec## LR5 LFR
## 611.8 644.0Catch Limit MP
Apply AvC MP which sets the TAC to the average historical catch:
Rec <- AvC(1, Cobia, reps=1000)
Rec## TAC (median)
## 461.5988Note that the TAC slot contains 1000 reps of the average catch (assuming a CV of 0.2, see ?AvC for details):
length(Rec@TAC)## [1] 1000hist(Rec@TAC)
Printing the Rec object returns the median of the 1000 reps:
Rec## TAC (median)
## 461.5988Mixed MP (TAC and Spatial)
Apply Itarget1_MPA MP which sets a TAC and closes Area 1:
Rec <- Itarget1_MPA(1, Cobia, reps=1000)
Rec## TAC (median) Area 1 Area 2
## 572.0487 0.0000 1.0000