特異確率密度関数からのランダムサンプリング
数式分布を確率密度関数で定義する.密度関数は連続的ではなく,無限大の特異点を含む.
In[1]:=
![Click for copyable input](assets.ja/random-sampling-from-a-singular-pdf/In_44.png)
dist = ProbabilityDistribution[
Piecewise[{{1/(4 Sqrt[Abs[x]]), -1 < x < 0}, {2/3, 1/4 < x <= 1}},
0], {x, -1, 1}]
Out[1]=
![](assets.ja/random-sampling-from-a-singular-pdf/O_38.png)
分布からランダムなサンプルを生成し,そのヒストグラムを密度関数と比べる.
In[2]:=
![Click for copyable input](assets.ja/random-sampling-from-a-singular-pdf/In_45.png)
sample = RandomVariate[dist, 10^5];
完全なWolfram言語入力を表示する
Out[3]=
![](assets.ja/random-sampling-from-a-singular-pdf/O_39.png)