calc_steepness: Steepness (h) と関連するパラメータ (SB0,R0,B0)を計算する関数

In ichimomo/frasyr: Fisheries Research Agency (FRA) provides the method for calculating sustainable yield (SY) with R

Steepness (h) と関連するパラメータ (SB0,R0,B0)を計算する関数

Description

Steepness (h) と関連するパラメータ (SB0,R0,B0)を計算する関数

Usage

calc_steepness(
  SR = "HS",
  rec_pars,
  M,
  waa,
  maa,
  plus_group = TRUE,
  faa = NULL,
  Pope = TRUE
)

Arguments

SR	"HS", "BH", "RI"のいずれか
rec_pars	再生産関係のパラメータでrec_pars$a,rec_pars$bで使用する
M	年齢別自然死亡係数 (ベクトルで与えるか、年齢共通の場合M=0.4のようにしてもよい)
waa	（親魚量の）年齢別体重
maa	年齢別親魚量
plus_group	最高齢がプラスグループかどうか

Value

以下の要素からなるデータフレーム

SPR0

F=0のときのSPR(この逆数がreplacement lineの傾き)

SB0

F=0のときの親魚量

R0

F=0のときの加入量

B0

F=0のときの資源量

h

steepness: BHかRIのときは0.2×SB0のときの加入量がh×R0, HSのときはh=1-b/SB0

Examples

## Not run: 
data(res_vpa)
SRdata <- get.SRdata(res_vpa)
resSR <- fit.SR(SRdata, SR = c("HS","BH","RI")[1],
                method = c("L1","L2")[2], AR = 1,
               out.AR = TRUE)
rec_pars = resSR$pars
year <- "2017"
M = res_vpa$input$dat$M[,year]
waa = res_vpa$input$dat$waa[,year]
maa = res_vpa$input$dat$maa[,year]
Res_h = calc_steepness(SR="HS",rec_pars=rec_pars,M=M,waa=waa,maa=maa,plus_group=TRUE)
Res_h

## End(Not run)

ichimomo/frasyr documentation built on Dec. 11, 2023, 12:50 a.m.