2変量確率分布に従った乱数の発生方法について

２列以上の行列データから多変量のカーネル分布を作成したのち、

その分布に従う乱数を任意の数、発生させるにはどのようにすればよいでしょうか？

Answer 1

[Hiro Yoshino]

変数間の関係を記述するコピュラを指定する方法があるようです：

https://jp.mathworks.com/help/stats/examples/simulating-dependent-random-variables-using-copulas.html

Answer 2

[Hiro Yoshino]

カーネル分布というのはカーネル密度推定で得られた確率密度分布の事でしょうか？

それでしたら：

https://jp.mathworks.com/matlabcentral/answers/464916-how-to-generate-random-sampling-from-a-probability-distribution-function-generated-by-kernel-density

が参考になります。

また、FileExchangeからはこんなモノもあります：

https://jp.mathworks.com/matlabcentral/fileexchange/27693-random-numbers-from-a-user-defined-piecewise-polynomial-distribution

多変量正規分布なら

https://jp.mathworks.com/help/stats/mvnrnd.html

を参考にしてみては？

[Hiro Yoshino]

補足：

すぐに思いつく普通のやり方としては

pd = fitdist(x,'Kernel','Kernel','epanechnikov');
x_sampled = random(pd,N,1);

のようにやるのは上手く行きませんでしたか？

Answer 3

[Kyotaro Nishimoto]

回答ありがとうございます。

1変量であるならばfitdistで実施できることは確認しております。

２変量以上の場合はksdensity,mvdensityを使用する必要があると認識しておりますが

そちらには”pd”を算出→randomに入れて乱数を発生させる機能はなかったように思います。

（正規分布など、特定の分布であれば乱数を発生させる関数が準備されているようですが。。。）

先に提案いただいたほうを詳しく調べてみます。