ページ概要
このページは，確率日雨量の予測を行うプログラムを紹介しています．
- プログラムでは，以下の確率分布モデルが使用されています．
| No | モデル短縮表示 | 確率分布モデル | 母数の数 |
|---|
| 1 | Gumbel | グンベル分布 | 2 |
| 2 | GEV | 一般化極値分布 | 3 |
| 3 | Weibull (Goda) | ワイブル分布 (合田のL-モーメント法) | 3 |
| 4 | Weibull (LSM) | ワイブル分布 (最小二乗法) | 3 |
| 5 | Weibull (MLM) | ワイブル分布 (最尤法) | 3 |
| 6 | SQRT-ET | 平方根指数型最大値分布 | 2 |
| 7 | LP3 | 対数ピアソン III 型分布 | 3 |
| 8 | LN3 (Iwai) | 対数正規分布 (岩井法) | 3 |
| 9 | LN3 (Moment) | 対数正規分布 (積率法) | 3 |
| 10 | LN3 (Trial) | 対数正規分布 (試行法) | 3 |
| 11 | Exponential | 指数分布 | 2 |
| 12 | GPD | 一般化パレート分布 | 3 |
- プロッティングポジション公式としては，Cunnane 公式が用いられています．
- プログラム作成に当たり，以下の文献を参考にしました．
- http://thesis.ceri.go.jp/center/info/geppou/ceri/0005005050.html
(星清：水文統計解析，開発土木研究所月報，No.540，1998年5月)
- http://iss.ndl.go.jp/books/R000000004-I4497918-00
(星清：現場のための水文統計(1)，開発土木研究所月報，No.540，1998年5月)
- http://iss.ndl.go.jp/books/R000000004-I4522177-00
(星清・新目竜一・宮原雅幸：現場のための水文統計(2)，開発土木研究所月報，No.541，1998年6月)
- 入力サンプルデータの出典は，以下に示すとおりです．
- 札幌の日最大雨量データ : 文献：星清・新目竜一・宮原雅幸，現場のための水文統計(2)，開発土木研究所月報，No.541，1998年6月 (http://iss.ndl.go.jp/books/R000000004-I4522177-00)
- 前橋の日最大雨量データ : 気象庁ホームページ (Webpage of Japan Meteorological Agency)
- 京都の日最大雨量データ : 文献：椎葉充晴・立川康人・市川温，例題で学ぶ水文学，森北出版株式会社，20120年5月21日 (ISBN978-4-627-49631-6)
- 指数分布および一般化パレート分布は，閾値超過データに適用すべきものであるが，計算用の例題として，計算では前出の各都市の年最大降水量データを利用しています．
関連文書
入力サンプルデータ
プログラム
ソースコード
実行用バッチファイル
f90_QQ 用
gfortran -o f90_QQ.exe f90_QQ.f90
f90_QQ 1 inp_RF_M.txt out_QQ_M_01.csv
f90_QQ 2 inp_RF_M.txt out_QQ_M_02.csv
f90_QQ 3 inp_RF_M.txt out_QQ_M_03.csv
f90_QQ 4 inp_RF_M.txt out_QQ_M_04.csv
f90_QQ 5 inp_RF_M.txt out_QQ_M_05.csv
f90_QQ 6 inp_RF_M.txt out_QQ_M_06.csv
f90_QQ 7 inp_RF_M.txt out_QQ_M_07.csv
f90_QQ 8 inp_RF_M.txt out_QQ_M_08.csv
f90_QQ 9 inp_RF_M.txt out_QQ_M_09.csv
f90_QQ 10 inp_RF_M.txt out_QQ_M_10.csv
f90_QQ 11 inp_RF_M.txt out_QQ_M_11.csv
f90_QQ 12 inp_RF_M.txt out_QQ_M_12.csv
Format
| Arg. | Description |
|---|
| f90_QQ | 実行用プログラム |
| No | 確率分布モデルを指定する番号 |
| fnameR | 入力データファイル |
| fnameW | 出力データファイル |
上表中 'No' については，ページ上段「ページ概要」の表を参照願います．
f90_BTS 用
gfortran -o f90_BTS.exe f90_BTS.f90
f90_BTS 1 inp_RF_M.txt out_BTS_M_01.csv inp_BHI_M_01.txt
f90_BTS 2 inp_RF_M.txt out_BTS_M_02.csv inp_BHI_M_02.txt
f90_BTS 3 inp_RF_M.txt out_BTS_M_03.csv inp_BHI_M_03.txt
f90_BTS 4 inp_RF_M.txt out_BTS_M_04.csv inp_BHI_M_04.txt
f90_BTS 5 inp_RF_M.txt out_BTS_M_05.csv inp_BHI_M_05.txt
f90_BTS 6 inp_RF_M.txt out_BTS_M_06.csv inp_BHI_M_06.txt
f90_BTS 7 inp_RF_M.txt out_BTS_M_07.csv inp_BHI_M_07.txt
f90_BTS 8 inp_RF_M.txt out_BTS_M_08.csv inp_BHI_M_08.txt
f90_BTS 9 inp_RF_M.txt out_BTS_M_09.csv inp_BHI_M_09.txt
f90_BTS 10 inp_RF_M.txt out_BTS_M_10.csv inp_BHI_M_10.txt
f90_BTS 11 inp_RF_M.txt out_BTS_M_11.csv inp_BHI_M_11.txt
f90_BTS 12 inp_RF_M.txt out_BTS_M_12.csv inp_BHI_M_12.txt
Format
f90_BTS No fnameR fnameW1 fnameW2
| Arg. | Description |
|---|
| f90_BTS | 実行用プログラム |
| No | 確率分布モデルを指定する番号 |
| fnameR | 入力データファイル |
| fnameW1 | 出力データファイル |
| fnameW2 | ヒストグラム作成用データ出力ファイル |
上表中 'No' については，ページ上段「ページ概要」の表を参照願います．
GMT 実行用バッチファイル (前橋データ用)
- GMT での作図用に3種類のバッチファイルを用意しています．
- ひとつは制御用バッチファイル，その他はデータ作成用とプロットコマンド用バッチファイルです．
- データ作成用バッチファイルとプロットコマンド用バッチファイルは，制御用バッチファイルより呼び出されます．
- 同じような多くのグラフを作成する場合，この方法は大変便利です．
Q-Q プロット
Bootstrap 信頼区間
Bootstrap 点推定値ヒストグラム
