複数のVPAの結果を重ね書きする

plot_vpa(
  vpalist,
  vpatibble = NULL,
  what.plot = NULL,
  plot_year = NULL,
  legend.position = "top",
  vpaname = NULL,
  ncol = 2,
  scale_value = NULL,
  is_minor_ticks = TRUE
)

Arguments

vpalist

vpaの返り値をリストにしたもの; 単独でも可

vpatibble

tibble形式のVPA結果も可。この場合、convert_vpa_tibble関数の出力に準じる。複数のVPA結果がrbindされている場合は、列名"id"で区別する。

what.plot

何の値をプロットするか. NULLの場合、全て(SSB, biomass, U, catch, Recruitment, fish_number, fishing_mortality, weight, maturity, catch_number)をプロットする。what.plot=c("SSB","Recruitment")とすると一部のみプロット。ここで指定された順番にプロットされる。

legend.position

凡例の位置。"top" (上部), "bottom" (下部), "left" (左), "right" (右), "none" (なし)。

vpaname

凡例につけるVPAのケースの名前。vpalistと同じ長さにしないとエラーとなる

ncol

図の列数

is_minor_ticks

補助目盛りをつける場合はTRUE。デフォルトはTRUE

Examples

if (FALSE) { # \dontrun{
data(res_vpa)
res_vpa2 <- res_vpa
res_vpa2$naa <- res_vpa2$naa*1.2

plot_vpa(list(res_vpa, res_vpa2), vpaname=c("True","Dummy"))
plot_vpa(list(res_vpa, res_vpa2), vpaname=c("True","Dummy"),
                 what.plot=c("SSB","fish_number"))

} # }