All functions |
|
---|---|
VPAの結果に格納されている生物パラメータから世代時間を計算 |
|
Dispacher function for making table |
|
サブ目盛りの追加 |
|
列が年である行列に対して、年を指定するとその年のあいだの平均値(等)を返す関数 |
|
モデルのAICから計算される重みに応じて、複数のモデルからシミュレーション結果を抜き出す |
|
リッジVPAのλ探索の自動化 |
|
再生産関係における残差の時系列の自己相関等についてのプロット |
|
モデル平均的な再生産関係を与える |
|
様々なベータやinput設定をもとに複数の将来予測を繰り返して結果をtibbleで返す関数 |
|
bootstrapを実施する |
|
再生産関係のブートストラップ |
|
再生産関係のブートストラップの結果をggplotで生成する関数する関数 |
|
再生産関係の残差ブートストラップをプロットする関数 |
|
再生産関係の推定結果からbootstrapを実行し、hの分布を計算する関数 |
|
半期VPA専用のバブルプロット作図の関数 |
|
VPAの結果を受けて,バブルプロット(横軸が年,縦軸が年齢,バブルの大きさと色で値を示す)を生成する |
|
年齢別生物パラメータとFと漁獲量を与えると与えた漁獲量と一致するFへの乗数を返す |
|
年齢別生物パラメータとFと漁獲量を与えると与えた漁獲量と一致するFへの乗数を返す (optimを使わないでやろうとしたけどだめだったやつ) |
|
再生産関係の推定における標準化残差を計算する関数 |
|
年齢別パラメータを与えて、1年分前進計算する |
|
再生産関係の残差から事後的に1次の自己相関係数を推定する関数 |
|
calculate F/Ftarget based on F_%SPR multiplier |
|
将来予測の結果オブジェクトから生物パラメータを取り出して、その生物パラメータをベースに、Fvectorで与えたF at ageに対応するSPRを計算して返す |
|
MSYを達成するときの%SPRを計算する(calc_future_perSPRのwrapper) |
|
Steepness (h) と関連するパラメータ (SB0,R0,B0)を計算する関数 |
|
res_futureから考えられるほぼすべてのパフォーマンス指標をとりだし、2行のtibbleで返す |
|
再生産関係の残差の確率分布に関するチェック図を出力する関数 |
|
再生産関係推定が収束しているかや最適解を得られているかを診断する関数 |
|
MSEの計算が正しいかを確認するための関数 |
|
fit.SRとfit.SRregimeの重み付け(どのデータを使うか)の設定方法は、引数wで与える場合とSRdata$weightで与える場合の2パターンある。両者が矛盾なく設定されているかを確認するための関数 |
|
漁獲量の上限設定をしたときの設定がちゃんと生きているかどうかを確かめる |
|
値の大小に応じて表の背景にグラデーションをつける |
|
複数の管理基準値の推定結果を重ね書きする |
|
複数の再生産関係を比較する関数 |
|
F一定の場合で平衡状態になったときの統計量をx軸、y軸にプロットして比較する |
|
短期的将来予測における複数の管理方策のパフォーマンスを比較する表を出力する |
|
複数のKobe II tableの結果を重ね書きする |
|
F at ageをVPAの結果から%SPRで変換したりするための関数 |
|
fit.SRの結果をtibble形式に治す |
|
vpaデータ, 選択率の参照年1, 漁獲圧の参照年2を与えて、参照年2の漁獲圧のもとで参照年1の選択率となるF at ageを作成する関数。漁獲圧の変換は%SPR換算 |
|
future_vpaの結果オブジェクトのリストをtibble形式に変換する関数 |
|
future_vpaの結果オブジェクトをtibble形式に変換する関数 |
|
半期VPAの結果オブジェクトをtibble形式に変換する関数 |
|
VPAの結果オブジェクトをtibble形式に変換する関数 |
|
再生産関係の推定パラメータの相関を出力する関数 |
|
VPAデータが1年分追加されたダミーデータを生成する |
|
csvデータを読み込んでvpa用のデータを作成する関数 |
|
csvデータを読み込んでvpa.hay用のデータを作成する関数 |
|
管理基準値表から目的の管理基準値を取り出す関数 |
|
将来予測やVPAの結果から生物パラメータをとりだす |
|
future_vpaの返り値に要約統計量を加えるための内部関数 |
|
与えられた個体群動態でプラスグループが考慮されているかどうか |
|
年齢別漁獲尾数の不確実性評価の関数 |
|
vpa結果の収束診断のための関数 |
|
ジャックナイフ法を行う関数 |
|
VPAのレトロスペクティブ解析結果を自動生成する関数 |
|
VPAの感度分析のための関数 |
|
Draw fitted SR line(s) over observed stock-recruitment plots |
|
future_vpaを使ってMSY管理基準値などを計算するwrapper関数 |
|
表を画像として保存 |
|
kobeIItableから任意の表を取得し画像として保存 |
|
β調整による管理効果を比較する表を画像として一括保存 |
|
Export tables to csv |
|
Extract value at age from vpadata |
|
再生産関係の推定 |
|
fit.SRのwrapper関数でhに制約を置いたもの |
|
再生産関係のフィット(頑健版) |
|
ある系群から別の系群への移入が想定される場合の再生産関係の推定(マアジ用) |
|
レジーム分けを考慮した再生産関係の推定 |
|
Fit and draw SR line(s) over observed stock-recruitment plots |
|
kobeIItableから取り出した表を整形 |
|
将来予測の実施関数 |
|
毎年のFの%SPRやターゲットした%SPRに相当するFの相対的な大きさを計算する |
|
VPAの結果から再生産関係推定用のデータを作成する |
|
beverton-holtのh,R0とbioparsを与えるとa,bを返す関数 |
|
将来予測の結果のリストを入れると、代表的なパフォーマンス指標をピックアップする |
|
会議用の図の出力関数(大きさ・サイズの指定済):通常サイズ |
|
会議用の図の出力関数(大きさ・サイズの指定済):大きいサイズ |
|
Give age names starts from zero and ends with plus group |
|
再生産関係のレジームシフトを隠れマルコフモデルで推定する |
|
半期VPAのレトロスペクティブ解析用の関数 |
|
半期VPA専用のレトロスペクティブ解析の作図関数 |
|
半期VPAの作図関数 |
|
半期VPAの作図関数 |
|
再生産関係のジャックナイフ解析 |
|
使うフォルダ名を与えると一連の結果の関数を読み込む関数 |
|
管理基準値の表を作成する |
|
Get fitted SR data with model names |
|
Make summary table of ABC estimation |
|
将来予測用の三次元行列(年齢×年×シミュレーション)を与えられたら, pars.yearで指定された期間のパラメータを平均するか、parで指定されたパラメータを、year_replace_future以降の年で置き換える |
|
future_vpaにインプットとして入れる将来予測の空のarrayを生成する関数 |
|
beta.simluationの結果などを読んで、kobeII talbeに整形する関数 |
|
Make Kobe ratio (?) for frasyr_tool |
|
Make summary table of MSY estimation |
|
Make table for stock assessment result 2 |
|
Make table depending on the class of given object |
|
再生産関係の推定結果をtxtファイルに出力する関数 |
|
VPA結果をcsvファイルに出力する |
|
半期VPA専用のジャックナイフ法の作図関数 |
|
F currentをプロットする |
|
ref.Fの出力をプロットするための関数 |
|
HCRを書く |
|
縦軸が漁獲量のHCRを書く(traceの結果が必要) |
|
再生産関係をプロットする関数 |
|
ARの効果を差し引いて再生産関係をプロットする |
|
モデル平均の場合の再生産関係のプロット |
|
SRdataをプロットする |
|
fit.SRregimeの結果で得られた再生産関係をプロットするための関数 |
|
一つの将来予測の結果をプロットする(ggplotは使わず) |
|
将来予測の複数の結果をggplotで重ね書きする |
|
複数の将来予測の結果をプロットする(ggplotは使わず) |
|
Kobe plotを書く |
|
ブートストラップ法による信頼区間の算出と作図を自動で行う関数 |
|
残差プロットの関数 |
|
Plot SPR figure |
|
複数のVPAの結果を重ね書きする |
|
漁獲量曲線(yield curve)を書く |
|
profile likelihood |
|
プロファイル尤度を計算する関数 |
|
VPAの推定値についてprofile likelihood (one parameter)を実施する |
|
Util functions for select SRmodels Pull data frame from list |
|
kobeIItable から任意の表を指名して取り出す |
|
リスクテーブルを作る |
|
csvファイルとしてまとめられた資源計算結果を読み込んでRのオブジェクトにする |
|
re-calculate projection with different arguments |
|
再生産関係を仮定しない管理基準値計算(SPR,YPR,F0.1,Fmax)のための関数 |
|
半期VPA専用の残差プロット作図の関数 |
|
Retrieve function argument settings to reuse |
|
レトロスペクティブ解析の実施 |
|
do.callのsafe版 |
|
過去にさかのぼってブロックサンプリングをおこなう |
|
将来予測用の再生産関係の設定を行う関数 |
|
Source specific lines in an R file |
|
再生産関係推定でSR=HSにおいて同一尤度を持つ複数のパラメータセットがえられたとき、SR=Mesnilとしてcheck.SRfitを呼び出し、gammaを変えながら尤度が一意に定まる値を探す関数 |
|
kobe.tableをさらにsummaryする |
|
Make table of latest SB and F |
|
Make tables of %SPRs |
|
Make tables of latest SB and F relative to MSY values |
|
@export |
|
将来予測においてFが非常に小さい場合には決定論的予測と確率論的予測の平均がほぼ一致するかを確認するための関数 |
|
会議用の図のフォーマット |
|
tidy形式のVPAへのインプットデータをdata_dandlerへ渡す形式に変換する(暫定版) |
|
再生産関係を網羅的にフィットする |
|
Ridge VPAの計算をTMBで実行するためにcppファイルのコンパイル等をする関数 |
|
再生産関係のパラメタが想定した形になっているかを確認 |
|
VPAによる資源計算を実施する |
|
半期VPAで資源量推定を行う関数 |