SSMethod.TA1.8: Apply Francis composition weighting method TA1.8

Description Usage Arguments Author(s) References See Also Examples

View source: R/SSMethod.TA1.8.R

Description

Uses method TA1.8 (described in Appendix A of Francis 2011) to do stage-2 weighting of composition data from a Stock Synthesis model. Outputs a mutiplier, w (with bootstrap 95 so that N2y = w x N1y, where N1y and N2y are the stage-1 and stage-2 multinomial sample sizes for the data set in year y. Optionally makes a plot of observed (with confidence limits, based on N1y) and expected mean lengths (or ages).

CAUTIONARY/EXPLANATORY NOTE. The large number of options available in SS makes it very difficult to be sure that what this function does is appropriate for all combinations of options. The following notes might help anyone wanting to check or correct the code.

  1. The code first takes the appropriate database (lendbase, sizedbase, agedbase, or condbase) and removes un-needed rows.

  2. The remaining rows of the database are grouped into individual comps (indexed by vector indx) and relevant statistics (e.g., observed and expected mean length or age), and ancillary data, are calculated for each comp (these are stored in pldat - one row per comp). If the data are to be plotted, the comps are grouped, with each group corresponding to a panel in the plot, and groups are indexed by plindx.

  3. A single multiplier is calculated to apply to all the comps.

Usage

1
2
SSMethod.TA1.8(fit, type, fleet, part = 0:2, pick.gender = 0:3,
  seas = NULL, method = NULL, plotit = TRUE, maxpanel = 1000)

Arguments

fit

Stock Synthesis output as read by r4SS function SS_output

type

either 'len' (for length composition data), 'size' (for generalized size composition data), 'age' (for age composition data), or 'con' (for conditional age at length data)

fleet

vector of one or more fleet numbers whose data are to be analysed simultaneously (the output N multiplier applies to all fleets combined)

part

vector of one or more partition values; analysis is restricted to composition data with one of these partition values. Default is to include all partition values (0, 1, 2).

pick.gender

vector of one or more values for Pick_gender; analysis is restricted to composition data with one of these Pick_gender values. Ignored if type=='con'

seas

string indicating how to treat data from multiple seasons 'comb' - combine seasonal data for each year and plot against Yr 'sep' - treat seasons separately, plotting against Yr.S If is.null(seas) it is assumed that there is only one season in the selected data (a warning is output if this is not true) and option 'comb' is used.

method

a vector of one or more size-frequency method numbers (ignored unless type = 'size'). If !is.null(method), analysis is restricted to size-frequency methods in this vector. NB comps are separated by method

plotit

if TRUE, make an illustrative plot like one or more panels of Fig. 4 in Francis (2011).

maxpanel

maximum number of panels within a plot

Author(s)

Chris Francis, Andre Punt, Ian Taylor

References

Francis, R.I.C.C. (2011). Data weighting in statistical fisheries stock assessment models. Canadian Journal of Fisheries and Aquatic Sciences 68: 1124-1138.

See Also

SSMethod.Cond.TA1.8

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
## Not run: 
Nfleet <- length(myreplist$FleetNames)
for (Ifleet in 1:Nfleet)
  SSMethod.TA1.8(myreplist,"len",fleet=Ifleet,maxpanel=maxpanel)
for (Ifleet in 1:Nfleet)
  SSMethod.TA1.8(myreplist,"age",fleet=Ifleet,maxpanel=maxpanel)
for (Ifleet in 1:Nfleet)
  SSMethod.TA1.8(myreplist,"size",fleet=Ifleet,maxpanel=maxpanel)
for (Ifleet in 1:Nfleet)
  SSMethod.TA1.8(myreplist,"con",fleet=Ifleet,maxpanel=maxpanel)
for (Ifleet in 1:Nfleet)
  SSMethod.Cond.TA1.8(myreplist,fleet=Ifleet,maxpanel=maxpanel)

## End(Not run)

r4ss documentation built on May 29, 2017, 2:24 p.m.