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

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

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

if (FALSE) { # \dontrun{
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
} # }