問題タブ [gamma-distribution]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票する
0 に答える
823 参照

c++ - ランダム、ブースト、またはクロノなしで C++ のガンマ分布からサンプルを生成する

C++ でガンマ分布からサンプルを生成したいと思います。スーパーコンピューターでコードをコンパイルしていますが、これには g++ の最新バージョンがありません (また、更新することもできません)。クロノとランダムを使用してガンマ分布からサンプルを生成するコードを作成しました。残念ながら、g++ のバージョンではコードをコンパイルできません (フラグ -std=c++11 または -std=c++0x も使用できません)。また、コンパイラには Boost/TR1 がないため、使用できません。

基本 (指数の合計からガンマを生成するなど) に戻る以外に、この分布からサンプリングできる別の方法はありますか? MCMC を実装するのではなく、純粋な C++ を考えています。

0 投票する
2 に答える
2514 参照

java - Java でガンマ分布からランダム値を生成する

Javaでガンマ分布から数値を生成するにはどうすればよいですか?

0 投票する
0 に答える
161 参照

r - ガンマ分布の分散に antiD 関数を使用する

これは私の最初の投稿であり、コミュニティのすべてのルールに従うことを願っています.

モザイク パッケージの関数 antiD を使用して、R の形状パラメーター 2 とスケール パラメーター 3 でガンマ分布の分散を計算しようとしています。私が使用するRコードは次のとおりです

問題は、これらのパラメーターを使用したガンマ分布の分散が 2*3^2 = 18 (ガンマ分布に関する Wiki ページ)に等しくなければならないため、得られる数値が意味をなさないことです。さらに、variance() の上限 (デフォルトの下限は 0) として 10^4 を指定すると、次のように返されます。

10^4 から 10^5 までの積分は次のようになります。

variance(10^5)この場合、なぜ無意味な結果が生じるのか誰か知っていますか? また、投稿のスタイルに関する追加のコメントにも感謝します。

0 投票する
0 に答える
67 参照

r - SAS からの R での反復測定と分布のコード化

SAS に次のコードがあり、R で同じデータ セットに対してこのモデルを再作成する必要があります。モデルを再作成するための適切なコードを見つけるのに苦労しており、特にデータの分布をガンマ分布に変更しています。

SAS コード

どんな助けでも大歓迎です

0 投票する
0 に答える
344 参照

r - Matlab の負の二項分散パラメーター

関数 matlabは、負の二項式の値とをnbinfit返します。負の二項式の「生態学的」または「Polya」パラメーター化の代わりに、mu (平均) とシータ (分散パラメーター) の値を返す同等の MLE 関数が matlab にありますか?rp

fitdistrRのようなもの)

0 投票する
1 に答える
451 参照

r - トゥイーディーの有限混合

トゥイーディー (または複合ポアソンガンマ) 分布の有限混合を推定しようとしています。これを行う方法に関するリソースを見つけることなく、考えられるリソースを精査しました。

現在、R でflexmixパッケージを使用して、 12 ~ 14 ページのflexmixビネットで概説されているように、別の M ステップ ドライバーを作成しようとしています。cplmパッケージに依存する私のコードは次のとおりです。

ただし、これにより次のエラーが発生します。

dtweedie(y, xi = p, mu = mu, phi = phi) のエラー: 非適合配列でのバイナリ演算

tweedie ディストリビューションの有限混合を実行または見た人はいますか? flexmixなどを使用して、これを達成するための正しい方向を教えてもらえますか?