Learn R Programming

EcoState

Package ecostate fits a state-space mass-balance model EcoState intended for aquatic ecosystems, using mass-balance equations matching those from Ecopath and dynamical equations matching Ecosim. Unlike the Ecopath-with-Ecosim (EwE) package, ecostate fits both biological parameters (e.g., equilibrium biomass and predator-prey vulnerability) and measurement parameters (e.g., catchability coefficients) via fit to time-series data. ecostate also estimates additional process errors representing nonstationarity in growth efficiency, ecotrophic efficient, migration, or other unmodeled processes. These process errors allow biomass patterns to closely match available data, so that resulting consumption (and associated productivity and mortality rates) can accurately be conditioned upon any residual patterns.

Installation

ecostate can be installed from GitHub using:

library(remotes)
install_github( "James-Thorson-NOAA/ecostate" )

Or to access vignettes from your R session, please instead use:

remotes::install_github( "James-Thorson-NOAA/ecostate",
                          build_vignettes = TRUE )
browseVignettes("ecostate")

More details

For more background please read:

Thorson, J. Kristensen, K., Aydin, K., Gaichas, S., Kimmel, D.G., McHuron, E.A., Nielsen, J.N., Townsend, H., Whitehouse, G.A. EcoState: Extending Ecopath with Ecosim to estimate biological parameters and process errors using RTMB and time-series data. Pre-print URL: https://doi.org/10.32942/X2QK81

or the updated version:

Thorson, J. Kristensen, K., Aydin, K., Gaichas, S., Kimmel, D.G., McHuron, E.A., Nielsen, J.N., Townsend, H., Whitehouse, G.A. (In press) The benefits of hierarchical ecosystem models: demonstration using a new state-space mass-balance model EcoState. Fish and Fisheries.

Copy Link

Version

Install

install.packages('ecostate')

Monthly Downloads

120

Version

0.2.0

License

GPL-3

Maintainer

James Thorson

Last Published

November 25th, 2024

Functions in ecostate (0.2.0)

ddirmult

Dirichlet-multinomial
abm3pc_sys

Adams-Bashford-Moulton for system of equations
add_equilibrium

Compute equilibrium values
compute_tracer

Calculate tracers, e.g., trophic level
dBdt

Dynamics from EcoSim
ginv

Penrose-Moore pseudoinverse
eastern_bering_sea

eastern Bering Sea ecosystem data
ecostate

fit EcoState model
compute_nll

Compute negative log-likelihood for EcoState model
ecostate_control

Detailed control for ecostate structure
print_ecopars

Print EcoSim parameters
plot_foodweb

Plot foodweb
print.ecostate

Print fitted ecostate object
rk4sys

Classical Runge-Kutta for system of equations
logLik.ecostate

Marginal log-likelihood
ode23

Non-stiff (and stiff) ODE solvers
stanza_settings

Detailed control for stanza structure
whitehouse_2021

Full rpath inputs for eastern Bering Sea