Learn R Programming

⚠️There's a newer version (6.0.6) of this package.Take me there.

DLMtool

DLMtool is a R package for management strategy evalaution and implementation of data-limited fishery stock assessment methods.

Find more information about DLMtool at http://datalimitedtoolkit.org/.

Installation

The released version of DLMtool can be downloaded from CRAN:

install.packages("DLMtool")

This GitHub repository contains the development version of DLMtool. The development version may be unstable, and it is not recommended for analyses.

User Guide and Documentation

A DLMtool User Guide and documentation for the DLMtool package are now available: https://dlmtool.github.io/DLMtool/.

Bugs and Issues

Please report any bugs or issues through GitHub or the DLMtool website.

Copy Link

Version

Install

install.packages('DLMtool')

Monthly Downloads

389

Version

3.2.3

License

GPL-2

Maintainer

Tom Carruthers

Last Published

April 10th, 2017

Functions in DLMtool (3.2.3)

BlankSelPlot

Internal function to create a blank plot for mapping selectivity at length
Blue_shark

Blue_shark Stock
Bluefin_tuna

Bluefin_tuna Stock
Bluefin_tuna_WAtl

Bluefin_tuna_WAtl Stock
BK_CC

Beddington and Kirkwood life-history method combined with catch curve
Albacore

Albacore Stock
Atlantic_mackerel

Atlantic_mackerel DLM_data
AvC

Average Catch
BK

Beddington and Kirkwood life-history MP (simple version)
CC4

Constant catch management procedure of Geromont and Butterworth (2014)
CSRA

Catch at size reduction analysis
CSRAfunc

Optimization function for CSRA
Can

What data-limited methods can be applied to this DLM_data object?
ChooseSelect

Manually choose the historical selectivity pattern
Cobia

Cobia DLM_data
BK_ML

Beddington and Kirkwood life-history analysis with mean-length estimator of
Butterfish

Butterfish Stock
CC1

Constant catch management procedure of Geromont and Butterworth (2014)
DAAC

Depletion Adjusted Average Catch
China_rockfish

China_rockfish DLM_data
ChooseEffort

Manually map the historical relative fishing effort trajectory.
Converge

Check Convergence
Cplot

Plot the median biomass and yield relative to last historical year
DD4010

Delay - Difference Stock Assessment with UMSY and MSY leading coupled with a
DDe

Effort control version of DD - Delay - Difference Stock Assessment with UMSY
DLM_fease-class

Class
DBSRA

Depletion-Based Stock Reduction Analysis
DBSRA_ML

Depletion-Based Stock Reduction Analysis using mean length estimator of
DCAC

Depletion Corrected Average Catch
Cant

What methods can't be applied to this DLM data object
CheckConverg

MSE convergence diagnostic
CompSRA

Age-composition-based estimate of current stock depletion given constant Z
CompSRA4010

Age-composition-based estimate of current stock depletion given constant Z
DDe75

Effort control version of DD - Delay - Difference Stock Assessment with UMSY
DBSRA4010

Depletion-Based Stock Reduction Analysis paired with 40-10 harvest control
DBSRA_40

Depletion-Based Stock Reduction Analysis assuming 40 per cent stock
DCAC_ML

Depletion-Based Stock Reduction Analysis using mean-length estimator of
DOM

How dominant is an MP?
DTe40

Effort searching MP aiming for 40 per cent stock depletion
FMSYopt

Internal function FMSY and related metrics
FMSYref

A reference FMSY method (uses perfect information about FMSY)
FlatE_Dom

FlatE_Dom Fleet
DLMtool-package

Data-Limited Methods Toolkit
FMSYref50

A reference FMSY method that fishes at half of FMSY (uses perfect
FMSYref75

A reference FMSY method that fishes at three quarters of FMSY (uses perfect
Fdem_CC

Demographic FMSY method using catch-curve analysis to estimate recent Z
DD

Delay - Difference Stock Assessment with UMSY and MSY leading
DecE_HDom

DecE_HDom Fleet
DecE_NDom

DecE_NDom Fleet
FlatE_NDom

FlatE_NDom Fleet
Fleet-class

Class
DDes

Effort searching version of DD - Delay - Difference Stock Assessment with
DTe50

Effort searching MP aiming for 50 per cent stock depletion
DecE_Dom

DecE_Dom Fleet
DepF

Depletion Corrected Fratio
DynF

Dynamic Fratio MP
Fease

MP feasibility diagnostic
Fease_xl

Read in feasibility parameters from Excel spreadsheet
Generic_DecE

Generic_DecE Fleet
Generic_FlatE

Generic_FlatE Fleet
FlatE_HDom

FlatE_HDom Fleet
Fratio

An FMSY/M ratio method
Fratio4010

An FMSY/M ratio method paired with the 40-10 rule
DLMDataDir

Directory of the installed package on your computer
DLM_data-class

Class
Fadapt

An adaptive MP that uses trajectory in inferred suplus production and
Fdem

Demographic FMSY method
GB_CC

Geromont and Butterworth Constant Catch Harvest Control Rule
GB_slope

Geromont and Butterworth index slope Harvest Control Rule
ITM

Index Target based on natural mortality rate
HDAAC

Hybrid Depletion Adjusted Average Catch
Herring

Herring Stock
ITe5

Index Target Effort-Based 5
Imprecise_Biased

Imprecise_Biased Observation
IT10

Index Target 10
IT5

Index Target 5
Itarget4

A management procedure that incrementally adjusts the TAC (starting from
ItargetE1

A management procedure that incrementally adjusts the effort to reach a
LstepCC1

A management procedure that incrementally adjusts the TAC according to the
DCAC4010

Depletion Corrected Average Catch paired with the 40-10 rule
DCAC_40

Depletion Corrected Average Catch assuming 40 per cent stock depletion
EtargetLopt

Effort MP: adjust effort up/down if mean length above/below Ltarget
Example_datafile

Example_datafile DLM_data
Feasibility

Feasibility DLM_fease
Feasibility2

Feasibility2 DLM_fease
Generic_obs

Generic_obs Observation
Gulf_blue_tilefish

Gulf_blue_tilefish DLM_data
Imprecise_Unbiased

Imprecise_Unbiased Observation
ITe10

Index Target Effort-Based 10
IncE_HDom

IncE_HDom Fleet
IncE_NDom

IncE_NDom Fleet
LBSPR_ItEff

Length-based SPR model with HCR that iteratively adjusts Effort
Fdem_ML

Demographic FMSY method that uses mean length data to estimate recent Z
Fratio_CC

A data-limited method that uses FMSY/M ratio and a naive catch-curve
Fratio_ML

An FMSY/M ratio MP that uses a mean length estimator of recent Z
LBSPR_ItSel

Length-based SPR model with HCR that iteratively adjusts Selectivity
LBSPR_ItTAC

Length-based SPR model with HCR that iteratively adjusts TAC
Low_Effort_Non_Target

Low_Effort_Non_Target Fleet
MCD

Mean Catch Depletion
Input

Function to run a set of input control methods
Islope1

A management procedure that incrementally adjusts the TAC to maintain a
Ltarget1

A management procedure that incrementally adjusts the TAC to reach a target
Ltarget4

A management procedure that incrementally adjusts the TAC to reach a target
MCD4010

Mean Catch Depletion
Perfect_Info

Perfect_Info Observation
Porgy

Porgy Stock
Required

What methods need what data
Rockfish

Rockfish Stock
Kplot

KOBE plot: a projection by projection plot of F/FMSY and B/BMSY
L2A

Length to age conversion
MSE-class

Class
Mackerel

Mackerel Stock
GB_target

Geromont and Butterworth target CPUE and catch MP
Gcontrol

G-control MP
Generic_IncE

Generic_IncE Fleet
Generic_fleet

Generic_fleet Fleet
Islope4

A management procedure that incrementally adjusts the TAC to maintain a
ML2D

Depletion and F estimation from mean length of catches
MPStats

Calculate Statistics for MP Performance
Pplot

A projection by projection plot of F/FMSY and B/BMSY
Pplot2

A projection by projection plot of F/FMSY, B/BMSY, B/B0, and yield
IncE_Dom

IncE_Dom Fleet
ItargetE4

A management procedure that incrementally adjusts the Effort to reach a
Jplot

Joint probability plot
LtargetE1

A management procedure that incrementally adjusts the Effort to reach a
SBT1

SBT simple MP
SBT2

SBT complex MP
SPslope

Slope in surplus production MP
LtargetE4

A management procedure that incrementally adjusts the Effort to reach a
MRnoreal

An marine reserve in area 1 with no spatial reallocation of fishing effort
MRreal

An marine reserve in area 1 with full reallocation of fishing effort
NOAA_plot

National Oceanographic and Atmospheric Administration default plot 1
Needed

Data needed to get MPs running
Precise_Biased

Precise_Biased Observation
Precise_Unbiased

Precise_Unbiased Observation
SelectFun

Selectivity at length function
Sense

Sensitivity analysis
Targeting_Small_Fish

Targeting_Small_Fish Fleet
Toothfish

Toothfish Stock
Itarget1

A management procedure that incrementally adjusts the TAC (starting from
L95target

A management procedure that adjusts the TAC up/down from reference (target)
LBSPR

Apply the Length-based SPR model to DLMtool Data Object
LstepCE1

A management procedure that incrementally adjusts the TAC according to the
Range

Standardize values
Rcontrol

Harvest Control Rule using prior for intrinsic rate of increase
Stock-class

Class
Target_All_Fish

Target_All_Fish Fleet
Tplot

A trade-off plot for an MSE object
condmet

Condition met?
curE

Fishing at current effort levels
genLenComp

Generate length composition of catch
Sam

Conduct stock assessment
SketchFun

Manually map the historical relative fishing effort trajectory.
Snapper

Snapper Stock
LstepCE2

A management procedure that incrementally adjusts the Effort according to
Observation-class

Class
OneRep

Make stochastic variables certain for only one rep
SetRecruitCycle

Function to calculate cyclic recruitment pattern given user-specified values
Simulation_1

Simulation_1 DLM_data
YPR

Yield Per Recruit analysis to get FMSY proxy F01
YPR_CC

Yield Per Recruit analysis to get FMSY proxy F01 paired to a naive catch
betaconv

Calculate beta parameter for beta distribution from mean and standard deviation
boxplot.DLM_data

Boxplot of TAC recommendations
doprojPI

Internal Projection Function
doprojPI_cpp

Rcpp version of the Projection function for calculating Reference Yield
mconv

Get log normal mean from transformed space mean and standard deviation
minlenLopt1

This input control sets the minimum length of fish caught to a fraction of
slotlim

An data-limited method which sets a slot limit
summary,DLM_data-method

Summary of DLM_data object
getAFC

Calculate age at first capture from length at first capture and growth
getSlope2

Internal function to calculate slope
getSlopeFun

Internal function to calculate slope
Sub

Subset MSE object by management procedure (MP) or simulation.
Tplot2

A shorter version of the plot method for MSEs that just shows the overall
TradePlot

Generic Trade-off Plot
boxplot.MSE

Boxplot of MP performance from MSE object
LstepCC4

A management procedure that incrementally adjusts the TAC according to the
NAor0

Is a value NA or zero.
NFref

No Fishing Reference MP
VOI2

Calculate Value Of Information 2
VOIplot

Yet another Value of Information Plot
YPR_ML

Yield Per Recruit analysis to get FMSY proxy F01 paired with a mean-length
alphaconv

Calculate alpha parameter for beta distribution from mean and standard deviation
SPMSY

Catch trend Surplus Production MSY MP
SPSRA

Surplus Production Stock Reduction Analysis
TAC

Calculate TAC recommendations for more than one MP
TACfilter

TAC Filter
TwoSidedFun

Internal Two sided selectivity curve
makePerf

Convert a OM object to one without observation or process error
makeTransparent

Make colors transparent
ourReefFish

Example data object
plot,DLM_data-method

Plot DLM_data object
VOI

Calculate Value Of Information
avail

What objects of this class are available
barplot.MSE

Plot a barplot of MSE results
curE75

Fishing at 75 per cent of current effort levels
cv

Calculate CV from vector of values
getmov2

Optimization function to find a movement model that matches user specified
getq

Optimization function that find the catchability (q where F=qE) value
movfit

Optimization function that returns the squared difference between user
movfit_Rcpp

Rcpp version of the Optimization function that returns the squared difference between user
summary,MSE-method

Summary of MSE object
tdlnorm

Calculate density of log-normally distributed random numbers
OM-class

Class
OM_xl

Read in operating model parameters from Excel spreadsheet
Rcontrol2

MP using prior for intrinsic rate of increase with a quadratic approximation
runMSErobust

Run a Management Strategy Evaluation
sampy

Sample vector
trlnorm

Generate log-normally distributed random numbers
updateMSE

Update an MSE object with new slots
checkMSE

Check that MSE object includes all slots
getFref2

Internal Get Reference F using Rcpp
getSlope1

Internal function to calculate ascending slope of selectivity curve
wormplot

Biomass wormplot
writeCSV

Internal function to write CSVs for objects
getEffhist

Calculate historical fishing mortality
getFMSY

Calculate FMSY and related metrics
getclass

get object class
Red_snapper

Red_snapper DLM_data
SPSRA_ML

Surplus Production Stock Reduction Analysis using a mean-length estimate of
SPmod

Surplus production based catch-limit modifier
Sole

Sole Stock
joinMSE

Join multiple MSE objects together
lmmodel-class

Class
plotOFL

A generic OFL plot for NOAA use
projOpt_cpp

Rcpp version of the Projection Optimizer
Splot

Scatter plot of B/BMSY or B/B0 and F/FMSY for lastYrs
getFMSY2

Calculate FMSY and related metrics using Rcpp code
getFref

Internal Get Reference F
getq2

Optimization function that find the catchability (q where F=qE) value
getmov

Optimization function to find a movement model that matches user specified
matlenlim

A data-limited method in which fishing vulnerability is set according to the
matlenlim2

A data-limited method in which fishing vulnerability is set slightly higher
mypars

mypars DLM_general
sdconv

Get log normal standard deviation from transformed space mean and standard deviation
setup

Setup parallel processing
getroot

Calculate slope from ageM and age95
gettempvar

Creates a time series per simulation that has gradient grad and random normal walk with sigma
initialize-methods

~~ Methods for Function
plot,MSE-method

Plot MSE object
plotFun

Print out plotting functions
runInMP

Runs input control MPs on a DLM_data object.
runMSE

Run a Management Strategy Evaluation
optQ_cpp

Rcpp version of the q Optimizer
qopt

Internal optimization function that find the catchability (q where F=qE)
replic8

Enlarge (replicate) a DLM data object to create an additional dimension for