Wolfram言語™

特異確率密度関数からのランダムサンプリング

数式分布を確率密度関数で定義する．密度関数は連続的ではなく，無限大の特異点を含む．

In[1]:=

dist = ProbabilityDistribution[ Piecewise[{{1/(4 Sqrt[Abs[x]]), -1 < x < 0}, {2/3, 1/4 < x <= 1}}, 0], {x, -1, 1}]

Out[1]=

分布からランダムなサンプルを生成し，そのヒストグラムを密度関数と比べる．

In[2]:=

sample = RandomVariate[dist, 10^5];

In[3]:=

Show[Histogram[sample, {-1, 1, 0.05}, "PDF", PlotTheme -> "Detailed"], Plot[PDF[dist, x], {x, -1, 1}, ImageSize -> Medium, PlotTheme -> "Marketing", PlotRange -> {0, 2.5}], ImageSize -> Medium]

Out[3]=