sample_survey_numbers: Sample a numbers index of abundance from an atlantis scenario

In r4atlantis/atlantisom: Generates Data Sets From Atlantis Scenarios

Sample a numbers index of abundance from an atlantis scenario

Description

The function takes numbers-at-age data from an Atlantis scenario where the data was read in from Atlantis output using load_nc within run_truth. One does not need to use these functions to create dat, rather you must only ensure that the structure of dat is the same.

Usage

sample_survey_numbers(dat, cv)

Arguments

dat	A data.frame of numbers at age containing the following columns: species agecl polygon layer time atoutput The data.frame is generated from either create_survey or create_fishery_subset.
cv	Coefficient of variation for the entire species specific abundance a matrix with columns: species, cv

Details

This function simply sums numbers-at-age over polygons, and then applies user defined error to the total numbers. The result is a coastwide numbers estimate from the survey Improvements could be to provide polygon specific numbers, but the cv will need to be thought about.

Value

The standard dataframe as specified used in dat. The function sums over layers and makes $layers is NA.

Author(s)

Poseidon

Examples

d <- system.file("extdata", "SETAS_Example", package = "atlantisom")
species <- c("Pisciv_T_Fish","Pisciv_S_Fish")
truth <- run_truth(scenario = "outputs",
  dir = d,
  file_fgs = "Functional_groups.csv",
  file_bgm = "Geography.bgm",
  select_groups = species,
  file_init = "Initial_condition.nc",
  file_biolprm = "Biology.prm",
  file_runprm = "Run_settings.xml")

 boxes <- 1:3
effic <- data.frame(species=c("Pisciv_T_Fish","Pisciv_S_Fish"), efficiency=c(0.3,0.1))
selex <- data.frame(species=c(rep("Pisciv_T_Fish",10),rep("Pisciv_S_Fish",10)),
               agecl=c(1:10,1:10),
                selex=c(0,0,0.1,0.5,0.8,1,1,1,1,1,0,0,0.1,0.3,0.5,0.7,0.9,1,1,1))

tmp <- create_survey(dat=truth$nums, time=seq(10,55,3), species=species, boxes=boxes, effic=effic, selex=selex)


cv <- data.frame(species=species, cv=c(0.2,0.3))

survObsN <- sample_survey_numbers(dat=tmp,cv=cv)

r4atlantis/atlantisom documentation built on Nov. 12, 2023, 2:59 a.m.