DynState: 2 species quota dynamic state variable model

Description Usage Arguments Details Value Author(s) Examples

View source: R/RDynState.R

Description

The DynState function simulates decision and outcomes in a 2 species mixed

fisheries under Individual Tradeable Quota constraints or effort constraints. The simulation

technique used is called dynamic state variable modelling (Clark and Mangel,2000). The input to the model is

defined by two DynStateInput objects.

Usage

1
	DynState(inputSpp1,inputSpp2, inputEffort, control)

Arguments

inputSpp1

A DynStateInput object holding the expected mean catch rate and Cv for species 1 (FLQuants).

inputSpp2

A DynStateInput object holding the expected mean catch rate and Cv for species 2 (FLQuants).

inputEffort

A DynStateInput object holding the effort needed (array).

control

A control object of class DynState.control

Details

This function returns an object of class DynState with the decision array and the individual

simulations through this array with information on the effort, catch rates and exceeding of quota.

Value

An object of class DynState.

Author(s)

Jan Jaap Poos

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
 library(FLRDynState)



 inputple4 <- new("DynStateInput")

 inputple4@CatchMean  <- FLQuant(rep(c(25,20,8),each=12), dimnames=list(age="all", year="all", unit="all",season=as.character(1:12), area=c("North","Mid","South")))

 inputple4@CatchSigma <- FLQuant(6, dimnames=dimnames(inputple4@CatchMean))



 inputsol4 <- new("DynStateInput")

 inputsol4@CatchMean <-  FLQuant(rep(c(1,4,6),each=12), dimnames=list(age="all", year="all", unit="all",season=as.character(1:12), area=c("North","Mid","South")))

 inputsol4@CatchSigma  <- FLQuant(2.5, dimnames=dimnames(inputsol4@CatchMean))



 effort <- array(c(9,8,7),dim=c(3,1))



 control <- DynState.control( Increments=20, Spp1Uplimit=20, Spp2Uplimit=30,

                             Spp1DiscardSteps=10, Spp2DiscardSteps=5, EffortUplimit=60,

                             EffortPrice=1, Spp1Price=1, Spp2Price=4,

                             ChoiceDist=1, SimNumber=100)



 z100 <- DynState(inputSpp1=inputple4, inputSpp2=inputsol4, effort, control)

flr/FLRDynState documentation built on Sept. 6, 2021, 5:51 a.m.