scipy.stats.確率密度関数を使って，skew distributionのランダムサンプリングはできないのでしょうか？

ここではvon mises分布を例にとります．
scipy.stats.vonmises
以下のように各種モーメントの計算例が載っているにも関わらず，.rvsメソッドで平均，集中度のパラメータしか指定しかできないことにもやもやします．

kappa = 3.99
mean, var, skew, kurt = vonmises.stats(kappa, moments='mvsk')

そもそも.statsメソッドの使い方もいまいちつかめません...
どなたかご教授いただければ幸いです．