ssdeltaMVLN: SSdeltaMVLN()

SSdeltaMVLNR Documentation

SSdeltaMVLN()

Description

function to generatbe kobe pdfs from a Multivariate Log-Normal Distribution including plotting option

Usage

SSdeltaMVLN(
  ss3rep,
  Fref = NULL,
  years = NULL,
  virgin = FALSE,
  catch.type = c("Exp", "Obs", "sel_bio", "ret_bio"),
  mc = 5000,
  weight = 1,
  run = "MVLN",
  plot = TRUE,
  addprj = FALSE,
  ymax = NULL,
  xmax = NULL,
  legendcex = 1,
  summary = TRUE,
  verbose = TRUE,
  quantiles = c(0.025, 0.975),
  seed = 123
)

Arguments

ss3rep

from r4ss::SS_output

Fref

Choice of Fratio c("MSY","Btgt","SPR","F01"), correponding to F_MSY and F_Btgt

years

single year or vector of years for mvln

virgin

if FALSE (default) the B0 base for Bratio is SSB_unfished

catch.type

reported catch c("Exp","Obs","sel_bio","ret_bio")

mc

number of monte-carlo simulations

weight

weighting option for model ensembles weight*mc

run

qualifier for model run

plot

option to show plot

addprj

include forecast years

ymax

ylim maximum

xmax

xlim maximum

summary

if TRUE trajectory summary table with MLEs and CIs is added

verbose

Report progress to R GUI?

quantiles

for computing CIs, default c(0.025,0.975)

seed

retains interannual correlation structure like MCMC

legendcex=1

Allows to adjust legend cex

Value

output list of quant posteriors and mle's

Author(s)

Henning Winker (JRC-EC)

Examples

mvn = SSdeltaMVLN(ss3sma,plot=TRUE) 

jabbamodel/ss3diags documentation built on Oct. 8, 2024, 11:49 p.m.