Description Usage Arguments Value Examples
This is equivalent to equation 2 in the draft SPFI document. This function does a single estimate of the harvest rate values. It does not perform iterations to optimize.
1 | calc_hcwt.ty(r.ty.sum, d.tsa, n.ysa)
|
r.ty.sum |
Output from |
d.tsa |
Output from |
n.ysa |
Synonymous with |
A data frame of the harvest rate parameter estimates grouped by fishery stratum and year.
1 2 3 4 5 6 7 8 9 10 11 12 | ## Not run:
#look to \code{\link{buildSPFIscript}} for creating hrj.df
hrj.df <- hrj.df[hrj.df$spfiflag==1,]
cwtpop <- hrj.df[hrj.df$data.type=="Pop" & hrj.df$fishery.index == 1 & hrj.df$Stock.Number %in% stock.subset,]
cwtpop <- subset(cwtpop,select = -fishery.index) #n.ysa
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")
d.tsa <- calc_d.tsa(r.tsa.sum = r.tsa.sum, n.ysa = cwtpop, standardize.bol = TRUE)
hcwt.ty <- calc_hcwt.ty(r.ty.sum=r.ty.sum, d.tsa = d.tsa, n.ysa = cwtpop)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.