スキュー正規分布の乱数（NORMINV（）に相当）

Question

スキュー正規分布を使ってスキュー乱数ジェネレータを作成しようとしていますが、式（これはNORMDIST()を使用しています）私が望む数字を生成できるNORMINV()ベースの式に変換します。これは私の式である：

「x」は数ある

=2/s * NORMDIST((x-m)/s, 0, 1, 0) * NORMSDIST(k*(x-m)/s) 

、「S」は、スケールは「m」は位置であり、「k」は形状です。私は式onlineを見つけたと認めますが、間違った方法を除いて、私がやりたいことをするようです。

どうすればいいですか？

Answer 1

Adelchi Azzaliniは、スキューノーマル（SN）分布の疑似乱数の生成方法についてはpageです。ここでページから：

N（0,1）ステップ1 U0及びvはNORMINV(RAND(),0,1)としてExcelで計算することができる変量。形状パラメータαで表されるδはδ = α/SQRT(1+α*α)である。

注：PythonのSciPyは、SciPy 0.18.0rc2コードのhereを参照して、スキューノーマル確率変数（scipy.stats.skewnorm）を取得しています。 docsはAzzaliniとCapitanioの論文を参照しています。