# context("Simulation")
# test_that("generated length frequency has no recruits",{
# lh <- create_lh_list(vbk=0.21,
# linf=65,
# t0=-0.01,
# lwa=0.0245,
# lwb=2.79,
# S50=20,
# S95=26,
# selex_input="length",
# M50=34,
# M95=NULL,
# maturity_input="length",
# M=0.27,
# binwidth=1,
# CVlen=0.1,
# SigmaR=0.737,
# SigmaF=0.2,
# SigmaC=0.2,
# SigmaI=0.2,
# R0=1,
# qcoef=1e-5,
# start_ages=0,
# rho=0.43,
# nseasons=1)
# ## Demonstrate data generation option
# true <- generate_data(modpath=NULL,
# data_avail="Index_Catch_LC",
# itervec=1,
# Fdynamics="Ramp",
# Rdynamics="AR",
# lh=lh,
# Nyears=20,
# Nyears_comp=10,
# comp_sample=200,
# init_depl=0.4)
# ## Data input components
# years <- true$years ## total years to model, can be 1:20 or 1998:2017
# LF <- true$LF ## length composition data, years along rows and length bin along columns. year names should match 'year' quantity, e.g. 11:20 or 2008:2017 (can be any years within years to model)
# LF_recruits <- LF[,1]
# expect_true(all(LF_recruits==0))
# }
# test_that("data generation creates a list"){
# lh <- create_lh_list(vbk=0.21,
# linf=65,
# t0=-0.01,
# lwa=0.0245,
# lwb=2.79,
# S50=20,
# S95=26,
# selex_input="length",
# M50=34,
# M95=NULL,
# maturity_input="length",
# M=0.27,
# binwidth=1,
# CVlen=0.1,
# SigmaR=0.737,
# SigmaF=0.2,
# SigmaC=0.2,
# SigmaI=0.2,
# R0=1,
# qcoef=1e-5,
# start_ages=0,
# rho=0.43,
# nseasons=1)
# ## Demonstrate data generation option
# true <- generate_data(modpath=NULL,
# data_avail="Index_Catch_LC",
# itervec=1,
# Fdynamics="Ramp",
# Rdynamics="AR",
# lh=lh,
# Nyears=20,
# Nyears_comp=10,
# comp_sample=200,
# init_depl=0.4)
# expect_true(is.list(true))
# expect_false(all(is.na(true)))
# })
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.