こちらからどうぞhttp://logics-of-blue.com/vpa%EF%BC%88%E5%8C%97%E5%A4%A7%E4%BB%95%E6%A7%98%EF%BC%89/資源量推定手法の王道、VPAのプログラムです。最終齢はプラスグループで漁業が定常な場合を想定しました。最近年のFは過去3年間のFの平均値だと仮定。最高齢のFは最高齢−1齢のFと等しいと仮定しています。これは管理人が適当に作ったものであり、実データではありません。このデータを使って計算する方法を説明します。Rに詳しい人は以下の説明は読まなくていいです。言う名前だったとします。ディレクトリをデスクトップに変更。したらCatch<-read.table("data.txt",header=T)を実行します。したら、仮想CatchデータがRに読み込まれたことになります。このプログラムをRに張り付けて実行します。あたらしスクリプトを開いてから張り付けると楽。実行は「F5」キーです。VPAにはほかにも色々の計算方法があります。このプログラムはその一例です。VPAの解説をしてあるPDF資料などはネット上にたくさん落ちています。このプログラムは、Rらしからぬforループの連続使用により作られています。