|
|
◆ RanfArrGam_sub()
| Sub RanfArrGam_sub |
( |
Ret As |
Double, |
|
|
Optional Alpha As |
Double = 1, |
|
|
Optional Beta As |
Double = 1 |
|
) |
| |
53 ビット実数ガンマ乱数 (Marsaglia and Tsang の圧搾法) (ラグ付きフィボナッチ法) (サブルーチン版)
- 目的
- ガンマ分布 Γ(α, β) に従うガンマ乱数 (53 ビット実数) を生成する.
ラグ付きフィボナッチ法の乱数生成ルーチン RanfArrNext() により一様乱数を生成し, Marsaglia and Tsang の圧搾法によりガンマ乱数に変換する.
- 引数
-
| [out] | Ret | ガンマ乱数 (53 ビット実数).
Alpha あるいは Beta が正しくないときは0を返す. |
| [in] | Alpha | (省略可)
ガンマ分布 Γ(α, β) のパラメータ α(形状母数). (α > 0) (省略時 = 1) |
| [in] | Beta | (省略可)
ガンマ分布 Γ(α, β) のパラメータ β(尺度母数). (β > 0) (省略時 = 1) |
- 参考文献
- D. E. Knuth, “The Art of Computer Programming Third Edition, Vol.2 / Seminumerical Algorithms”, Addison Wesley Longman, 1998.
- 四辻哲章「確率分布乱数生成法」第2版, プレアデス出版, 2019.
|
1.9.7