createRAVfile: Create rav file.

Description Usage Arguments Value

Description

This function takes the output from the DM model and creates a rav file that can be used forrunning the VRAP simulations. The parameter simInd indicates which MCMC sim to use. This takes the information from the input file in the DM spreadsheet (read in by readDMData) and the results of the MCMC simulations. Current version has no uncertainty in management error and base exploitation rate must be included as a parameter. This could be calculated with the SR params, maturation rates, natMort, and harvest.

Usage

1
2
createRAVfile(bdat, input, tDat, dat, filename = "temp_rav.rav",
  estType = "median", sim = 1, rav.options = list())

Arguments

bdat

data for the bayesian specification of a DM run

input

a list with the other values needed for a DM run. The following are examples naturalMort, analysisType = "DM", SRfunction, includeMarineSurvival, includeFlow

tDat

a data frame of the transformed posteriors and calculated process error and Smsy for each draw. From calculatePEandSmsy(). This is specifically for VRAP input. The productivity (a) and capacity (b) posteriors are not adjusted (centered) if the covariates do not have 0 mean. VRAP needs the SR parameters in this raw form using non-centered (de-meaned) covariates.

dat

data from the A & P file

filename

name to give the rav file

estType

("median")

sim

indicates which MCMC sim from the posteriors to use if estType is not median

rav.options

list of the rav file options to use that do not come from the posteriors.

randomSeed

(0)

numRuns

(1000) The number of simulations to run for computing probabilities in VRAP.

numYears

(25) The number of years to project forward in the simulations.

minAge

(2)

maxAge

(5)

convergeCrit

(0.001)

centerCovFlag

("NO") "NO"/"YES"

marineSurvType

("Autoc") "Autoc"/"Cycle"/"Trend"

flowType

("Autoc") "Autoc"/"Cycle"/"Trend"

modelDepensation

("NO") "NO"/"YES"

depensation

(300)

QETcritical

(63)

depensPar3

(1)

recruitsFromAdultSpawners

("YES") "NO"/"YES"

SRvariation

("YES") "NO"/"YES"

smoltToAdultVar

("NO") "NO"/"YES"

baseExploitationRate

(0.67)

includeManagementError

("YES") "NO"/"YES"

manageErrorA

(65.3946) derived from estimates from Puget Sound for all except STL and WRS

manageErrorB

(0.0158) derived from estimates from Puget Sound for all except STL and WRS

lowerEscThreshold

(200)

numYearsToAvg

(5)

runType

("ER") "ER"/"Pop"

bufferMin

(0)

Value

nothing is returned but the rav file is written.


eeholmes/DM documentation built on May 26, 2019, 3:36 p.m.