calc_H.ty: (SPFI) Calculate the AEQ stratum harvest rate parameters...

In MichaelFolkes/ctctools: Analytical tools For The Pacific Salmon Commission's Chinook Technical Committee.

Description

This is equivalent to equation 6 in the draft SPFI document. In the Visual Basic this is referred to as the AEQScaler.

Usage

calc_H.ty(c.ty.sum, r.ty.sum, hcwt.ty)

Arguments

c.ty.sum	Output from calc_ty.sum.
r.ty.sum	Output from calc_ty.sum.
hcwt.ty	Output from calc_hcwt.ty.

Value

A data frame of the AEQ stratum harvest rate parameter estimates grouped by fishery stratum and year.

Examples

## Not run: 
aeqcwt <- hrj.df[hrj.df$data.type==data.type & hrj.df$fishery.index %in% fishery.subset & hrj.df$Stock.Number %in% stock.subset,]
if(region=="seak") aeqcwt <- adjustAlaska(x = aeqcwt, data.catch = data.catch)
c.ty.sum <- calc_ty.sum(x = aeqcwt, newvar.name = "c.ty.sum")

cwtcatch <- hrj.df[hrj.df$data.type=="NomCat" & hrj.df$fishery.index %in% fishery.subset & hrj.df$Stock.Number %in% stock.subset,]
if(region=="seak") cwtcatch <- adjustAlaska(x = cwtcatch, data.catch = data.catch)
r.ty.sum <- calc_ty.sum(x = cwtcatch, newvar.name = "r.ty.sum")

hcwt.ty <- calc_hcwt.ty(...)

H.ty <- calc_H.ty(c.ty.sum = c.ty.sum, r.ty.sum = r.ty.sum, hcwt.ty = hcwt.ty)

## End(Not run)

MichaelFolkes/ctctools documentation built on May 7, 2019, 4:56 p.m.