Learn R Programming

 

Mizer is an R package to run dynamic multi-species size-spectrum models of fish communities. The package has been developed to model marine ecosystems that are subject to fishing. However, it may also be appropriate for other aquatic ecosystems. By providing a framework for multi-species fisheries modelling as an R package, mizer enhances the accessibility, usability, and reproducibility of models, and thus aims to facilitate collaboration and innovation.

The package contains functions that allow you to set up an ecosystem model and then project it through time under different fishing strategies. Methods are included to explore the results, including plots and calculations of community indicators such as the slope of the size spectrum. Size-based models can be complicated, so mizer contains many default options that you can however change when needed.

Mizer has been supporting research in marine ecology and fisheries science since 2014 (see publications). Mizer is still under active development. Version 2.0 has increased the user-friendliness and the flexibility of the framework. Contributions from the user community are very welcome. There is a sister package called mizerExperimental where user contributions can be checked out and receive feedback from the community. Example mizer models can be contributed to mizerExamples. Follow us on twitter and read our blog to stay up-to-date with new developments. Does your project or publication use mizer? If so, we would love to know.

Recent work on mizer was funded by the European Commission Horizon 2020 Research and Innovation Programme under Grant Agreement No 634495 for the project MINOUW (http://minouw-project.eu/) and the Australian Research Council Discovery Project Rewiring Marine Food Webs.

Installation

The package is on CRAN and therefore available from R’s built-in package manager.

# Install latest released version from CRAN
install.packages("mizer")

# Alternatively, install the development version from GitHub
remotes::install_github("sizespectrum/mizer")

Example

The following code loads the mizer package, loads some information about species in the North Sea that comes as an example with the package, sets up the parameters for the mizer model, and runs a simulation for 10 years.

library(mizer)
params <- newMultispeciesParams(NS_species_params, NS_interaction)
sim <- project(params, t_max = 10, effort = 0)

The results of the simulation can then be analysed, for example via plots:

plot(sim)

See the accompanying Get started page for more details on how the package works, including detailed examples.

Dynamic multi-species size-spectrum model

Size-based multi-species models are important for fisheries science because they provide a more realistic and accurate representation of the dynamics of fish populations and the ecosystems in which they live. In contrast to traditional single-species models, which consider a single fish stock as an isolated unit, size-based multi-species models account for the fact that fish populations are part of a larger ecosystem and interact with other species through predation, competition, and other ecological processes.

One of the key advantages of size-based multi-species models is that they provide a more comprehensive understanding of the impacts of fishing on fish populations and ecosystems. By considering the size distribution of different fish species, these models can capture the effects of fishing on both target and non-target species, and on different life stages of a species. This is particularly important for species that are caught as bycatch or that are indirectly affected by fishing through changes in their food web.

Another advantage of size-based multi-species models is that they can be used to investigate the effects of environmental changes and other perturbations on fish populations and ecosystems. For example, these models can be used to explore the impacts of climate change on the distribution and abundance of fish populations, or the effects of habitat loss or pollution on fish communities. Because mizer is a mechanistic model, it can deduce the complex population-level changes that we are interested in from the simpler changes in the physiological rates and feeding interactions of individual fish species.

Overall, size-based multi-species models provide a more comprehensive and realistic framework for understanding the dynamics of fish populations and ecosystems, and for developing effective fisheries management strategies that account for the complex interactions among species and their environment.

A mizer model captures the interactions between species. The growth rates of fish are determined by the availability of prey and the death rates are influenced by the abundance of predators, as well as fishing. The model starts with the individual-level physiological rates for each species, as well as the predation preferences, and deduces the population-level dynamics from these. Thus quantities like fish diets and fisheries yields emerge dynamically and can be projected into the future.

Because a mizer model tracks the size of individuals as they grow up over several orders of magnitude from their egg size to their maximum size, it correctly tracks the ontogenetic diet shifts. An individual typically moves through several trophic levels during its life time. This is often not correctly captured in other multi-species models.

A mizer model can be set up with only a small amount of information because it uses allometric scaling relations and size-based feeding rules to choose sensible defaults for unknown parameters.

Setting up a new multi-species mizer model is a two-step process, similar to what may be familiar from Ecopath with Ecosim: First one calibrates the model to describe a steady state that is in agreement with current observations (as in Ecopath), then one chooses the additional parameters that determine the dynamics away from the steady state (as in Ecosim). This model can then be used to investigate future effects of changes in fishing policy or of environmental stressors.

A strong theoretical basis

One big advantage of a mizer model is that it is based on a strong mathematical foundation. This allows a degree of a priori understanding of the behaviour of the model that is absent in many other multi-species models. This theoretical foundation is well presented in the book “Fish Ecology, Evolution, and Exploitation” by Ken Andersen.

It is interesting to think of the marine ecosystem as a transport system that moves biomass from the size of primary producers (mostly unicellular plankton) up to the sizes of fish that humans like to consume. Each fish that grows up from egg size to maturity by eating smaller individuals is like a car on this biomass highway. The yield of our fisheries depend on this traffic flowing smoothly and without traffic jams.

An analogy with road traffic may be helpful:

In road traffic, if traffic density gets too high in a section of the highway, drivers slow down, which leads to a pile-up producing even higher traffic density, leading to further slow-down in a potentially vicious cycle known as a traffic jam. Traffic management that ignores how the traffic density affects traffic speed fails. Luckily our mathematical understanding of transport equations has made practical contributions to managing traffic in ways that produce smoother traffic flow and hence higher throughput.

Mizer implements the transport equations for marine ecosystems. The potential for traffic jams is the same: if for example there is a high density of predators of a particular size, which all have preference for prey of a particular smaller size, then due to competition for that prey the growth of those predators slows down, leading to a pile-up which leads to further depletion of prey, leading to further slow-down, in a potentially vicious cycle.

Luckily, the natural ecosystem has evolved to facilitate very smooth traffic on this biomass highway, with resultant high productivity. This state is characterised by an approximate power-law shape of the biomass size spectrum. The purpose of mizer is to allow us to understand how various stressors, like fishing or climate change, affect the size spectrum and hence the flow of biomass and the productivity and resilience of the marine ecosystem. Mizer allows us to investigate how size-based fisheries management strategies can be used to keep the ecosystem close to its natural productive state.

Copy Link

Version

Install

install.packages('mizer')

Monthly Downloads

481

Version

2.5.3

License

GPL-3

Issues

Pull Requests

Stars

Forks

Maintainer

Gustav Delius

Last Published

October 17th, 2024

Functions in mizer (2.5.3)

animateSpectra

Animation of the abundance spectra
age_mat

Calculate age at maturity
calc_selectivity

Calculate selectivity from gear parameters
SheperdRDD

Sheperd function to calculate density-dependent reproduction rate
NS_params

Example MizerParams object for the North Sea example
box_pred_kernel

Box predation kernel
NS_species_params_gears

Example species parameter set based on the North Sea with different gears
NS_species_params

Example species parameter set based on the North Sea
customFunction

Replace a mizer function with a custom version
calibrateYield

Calibrate the model scale to match total observed yield
compareParams

Compare two MizerParams objects and print out differences
constantEggRDI

Choose egg production to keep egg density constant
calibrateNumber

Calibrate the model scale to match total observed number
completeSpeciesParams

Alias for validSpeciesParams()
constant_other

Helper function to keep other components constant
constantRDD

Give constant reproduction rate
calibrateBiomass

Calibrate the model scale to match total observed biomass
default_pred_kernel_params

Set defaults for predation kernel parameters
defaults_edition

Default editions
gear_params

Gear parameters
different

Check whether two objects are different
getBiomass

Calculate the total biomass of each species within a size range at each time step.
distanceMaxRelRDI

Measure distance between current and previous state in terms of RDI
emptyParams

Create empty MizerParams object of the right size
double_sigmoid_length

Length based double-sigmoid selectivity function
finalN

Size spectra at end of simulation
finalNOther

Values of other ecosystem components at end of simulation
distanceSSLogN

Measure distance between current and previous state in terms of fish abundances
getComponent

Get information about other ecosystem components
getDiet

Get diet of predator at size, resolved by prey species
getCommunitySlope

Calculate the slope of the community abundance
getEGrowth

Get energy rate available for growth
getEffort

Fishing effort used in simulation
getEncounter

Get encounter rate
getERepro

Get energy rate available for reproduction
getCriticalFeedingLevel

Get critical feeding level
getESpawning

Alias for getERepro()
getEReproAndGrowth

Get energy rate available for reproduction and growth
getFMort

Get the total fishing mortality rate from all fishing gears by time, species and size.
getInteraction

Deprecated function to get interaction matrix
getM2Background

Alias for getResourceMort()
getM2

Alias for getPredMort()
getMeanMaxWeight

Calculate the mean maximum weight of the community
getFMortGear

Get the fishing mortality by time, gear, species and size
getMort

Get total mortality rate
getMeanWeight

Calculate the mean weight of the community
getFeedingLevel

Get feeding level
getGrowthCurves

Get growth curves giving weight as a function of age
getPredRate

Get predation rate
getParams

Extract the parameter object underlying a simulation
getRates

Get all rates
getRDD

Get density dependent reproduction rate
getRDI

Get density independent rate of egg production
getPredMort

Get total predation mortality rate
getPhiPrey

Get available energy
getProportionOfLargeFish

Calculate the proportion of large fish
getN

Calculate the number of individuals within a size range
getReproductionLevel

Get reproduction level
get_h_default

Get default value for h
getYieldGear

Calculate the rate at which biomass of each species is fished by each gear
getSSB

Calculate the SSB of species
getTimes

Times for which simulation results are available
get_gamma_default

Get default value for gamma
getYield

Calculate the rate at which biomass of each species is fished
getZ

Alias for getMort()
get_f0_default

Get default value for f0
initialNOther<-

Initial values for other ecosystem components
get_phi

Get values from feeding kernel function
indicator_functions

Description of indicator functions
get_initial_n

Calculate initial population abundances
get_size_range_array

Get size range array
getResourceDynamics

Deprecated functions for getting resource parameters
get_required_reproduction

Determine reproduction rate needed for initial egg abundance
getResourceMort

Get predation mortality rate for resource
initialNResource<-

Initial value for resource spectrum
get_ks_default

Get default value for ks
get_time_elements

Get_time_elements
log_breaks

Helper function to produce nice breaks on logarithmic axes
lognormal_pred_kernel

Lognormal predation kernel
matchBiomasses

Match biomasses to observations
matchYields

Match yields to observations
matchNumbers

Match numbers to observations
matchGrowth

Adjust model to produce observed growth
initial_effort

Initial fishing effort
inter

Alias for NS_interaction
mizer-package

mizer: Multi-species size-based modelling in R
mizerPredMort

Get total predation mortality rate needed to project standard mizer model
mizerEncounter

Get encounter rate needed to project standard mizer model
mizerEReproAndGrowth

Get energy rate available for reproduction and growth needed to project standard mizer model
mizerERepro

Get energy rate available for reproduction needed to project standard mizer model
mizerEGrowth

Get energy rate available for growth needed to project standard mizer model
mizerFMortGear

Get the fishing mortality needed to project standard mizer model
knife_edge

Weight based knife-edge selectivity function
mizerMort

Get total mortality rate needed to project standard mizer model
initialN<-

Initial values for fish spectra
mizerFMort

Get the total fishing mortality rate from all fishing gears
l2w

Length-weight conversion
mizerFeedingLevel

Get feeding level needed to project standard mizer model
newSingleSpeciesParams

Set up parameters for a single species in a power-law background
needs_upgrading

Determine whether a MizerParams or MizerSim object needs to be upgraded
mizerRates

Get all rates needed to project standard mizer model
mizerResourceMort

Get predation mortality rate for resource needed to project standard mizer model
newCommunityParams

Set up parameters for a community-type model
newTraitParams

Set up parameters for a trait-based multispecies model
mizerPredRate

Get predation rate needed to project standard mizer model
noRDD

Give density-independent reproduction rate
mizerRDI

Get density-independent rate of reproduction needed to project standard mizer model
newMultispeciesParams

Set up parameters for a general multispecies model
plotGrowthCurves

Plot growth curves
plot,MizerParams,missing-method

Summary plot for MizerParams objects
plot,MizerSim,missing-method

Summary plot for MizerSim objects
plotFMort

Plot total fishing mortality of each species by size
plotDiet

Plot diet, resolved by prey species, as function of predator at size.
plotBiomass

Plot the biomass of species through time
plotDataFrame

Make a plot from a data frame
plotM2

Alias for plotPredMort()
plotFeedingLevel

Plot the feeding level of species by size
plotBiomassObservedVsModel

Plotting observed vs. model biomass data
project_simple

Project abundances by a given number of time steps into the future
plotPredMort

Plot predation mortality rate of each species against size
plotYieldObservedVsModel

Plotting observed vs. model yields
project

Project size spectrum forward in time
plotYield

Plot the total yield of species through time
plotting_functions

Description of the plotting functions
power_law_pred_kernel

Power-law predation kernel
projectToSteady

Project to steady state
plotSpectra

Plot the abundance spectra
plotYieldGear

Plot the total yield of each species by gear through time
scaleModel

Change scale of the model
reexports

Objects exported from other packages
setBevertonHolt

Set Beverton-Holt reproduction without changing the steady state
resource_semichemostat

Project resource using semichemostat model
resource_logistic

Project resource using logistic model
resource_params

Resource parameters
removeSpecies

Remove species
saveParams

Save a MizerParams object to file, and restore it
resource_constant

Keep resource abundance constant
renameSpecies

Rename species
setMetabolicRate

Set metabolic rate
setInteraction

Set species interaction matrix
setMetadata

Set metadata for a model
setFishing

Set fishing parameters
setExtMort

Set external mortality rate
setColours

Set line colours and line types to be used in mizer plots
setExtEncounter

Set external encounter rate
setInitialValues

Set initial values to values from a simulation
setMaxIntakeRate

Set maximum intake rate
setComponent

Add a dynamical ecosystem component
setRmax

Alias for setBevertonHolt()
setSearchVolume

Set search volume
setRateFunction

Set own rate function to replace mizer rate function
set_multispecies_model

Deprecated obsolete function for setting up multispecies parameters
setReproduction

Set reproduction parameters
setPredKernel

Set predation kernel
setParams

Set or change any model parameters
set_community_model

Deprecated function for setting up parameters for a community-type model
set_species_param_default

Set a species parameter to a default value
setResource

Set resource dynamics
steady

Set initial values to a steady state for the model
truncated_lognormal_pred_kernel

Truncated lognormal predation kernel
set_trait_model

Deprecated function for setting up parameters for a trait-based model
summary,MizerParams-method

Summarize MizerParams object
species_params

Species parameters
sigmoid_weight

Weight based sigmoidal selectivity function
summary_functions

Description of summary functions
steadySingleSpecies

Set initial abundances to single-species steady state abundances
summary,MizerSim-method

Summarize MizerSim object
sigmoid_length

Length based sigmoid selectivity function
validSim

Validate MizerSim object and upgrade if necessary
validSpeciesParams

Validate species parameter data frame
upgradeParams

Upgrade MizerParams object from earlier mizer versions
valid_gears_arg

Helper function to assure validity of gears argument
validGearParams

Check validity of gear parameters and set defaults
valid_species_arg

Helper function to assure validity of species argument
validParams

Validate MizerParams object and upgrade if necessary
w

Size bins
upgradeSim

Upgrade MizerSim object from earlier mizer versions
MizerSim-class

A class to hold the results of a simulation
NS_interaction

Example interaction matrix for the North Sea example
BevertonHoltRDD

Beverton Holt function to calculate density-dependent reproduction rate
N

Time series of size spectra
MizerParams

Alias for set_multispecies_model()
NOther

Time series of other components
MizerSim

Constructor for the MizerSim class
MizerParams-class

A class to hold the parameters for a size based model.
RickerRDD

Ricker function to calculate density-dependent reproduction rate
addSpecies

Add new species
NS_sim

Example MizerSim object for the North Sea example
age_mat_vB

Calculate age at maturity from von Bertalanffy growth parameters