問題タブ [mixture-model]

パッケージ mclust の coordProj 関数を使用しています。オプション「errors」を使用してプロットすると、シンボルと色を変更したいと考えています。

私のコード

プロットはこちら( http://i.stack.imgur.com/Nmgbq.png )

黒い点を文字「F」で、その他を「T」で変更しようとしています。それが不可能なのか、それとも間違っているのかはわかりません。

MATLAB でガウス混合モデルを適合させたい。具体的には、さまざまな開始点から数百のモデルを適合させ、それらを数回の EM 反復 (たとえば 10 回) に適合させたいと考えています。最終的な目的は、遺伝的アルゴリズムを使用して適切なモデルを検索することです。

混合モデルの適合には時間がかかるため、MATLAB の並列処理ツールボックスを使用して gpu で実行される高速な実装があるかどうか疑問に思っていましたか? gpuArray を gmcluster 関数に渡し、gmcluster と関連ファイルの適切な変数を gpuArray に変更することで、単純な GPU バージョンを実装してみました。これにより、実際にランタイムが悪化しました (特定のデータセットで CPU 時間は 0.1 で、GPU は 0.2 でした)。

単純な GPU の実装にできる改善点はありますか? (私は GPU プログラミングの初心者です)。

ここにあるようないくつかの実装があることは承知していますが、平均値と共分散の開始パラメータを受け入れるほど柔軟ではありません。

ヘルプ/ポインタをいただければ幸いです。

ありがとう

編集:これは、かなり時間がかかるコードの一部であり、GPU で実行したいと考えています。このコードは、各 og k ガウスの確率メンバーシップ行列ポストが与えられた場合に、mu で平均を計算し、Sigma で共分散を計算します。

どこ：

私の目的では、n は 1000 から 10,000 の間、d は 10 前後、k は 2 から 15 の間です。

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

（fitdistrRのようなもの）

各観測値に重み (0 から 1) が割り当てられている単変量分布にガウス混合モデル (GMM) を当てはめようとしています。

重みの割り当ては、基本的な「密度」関数 (「統計」から) に対して簡単に実行できますが、「mixtools」または「mclust」などの GMM パッケージに対してはこれを行うことができませんでした。

後者にはそのための機能もあります-「me.weighted」-リンク「http://cran.open-source-solution.org/web/packages/mclust/vignettes/mclust-weights.pdf」に例示されていますが、私はそれを働かせることができません。彼らの例と私のデータの両方で、エラーが発生しました：

scikit Learn での混合モデリングで発生している問題を、次の 1D の例で抽出しました。混合モデルを 1D 密度に単純に適合させる。単純な GMM を使用すると良好なパフォーマンスが得られますが、Variational または Dirichlet Process GMM を使用すると奇妙な動作が発生します。第一に、適合が非常に悪く、第二に、推定された分布が正規化されていません (以下のコードと画像を参照)。

これと同じ動作は、例http://scikit-learn.org/stable/auto_examples/mixture/plot_gmm_pdf.htmlで明らかです。GMM クラスによって返される確率はほぼ正規化されますが、DPGMM から返される確率は正規化されません。正規化は、次の行を追加することによって計算されました。

これはバグですか、それとも何か間違っていますか?

コード：

混合分布モデルを値のベクトルに当てはめようとしています。混合は、2 つのガウス分布と 1 つの一様分布で構成される必要があります。これをWinbugsに実装しようとしています。ガウスの混合物を使用した例をたくさん見つけましたが、ユニフォームを追加する方法がわかりません。以下のコード ペーストは現在、0 から 1 の間でスケーリングされた値のベクトルに適合するようにパラメーター化されていますが、「ノード NSD[1] の複数の定義」が得られるため、私の構造はまだ間違っているようです。助言がありますか？

の実装Rを使用して、ハードル モデルの有限混合でコーディングしようとしています。私の関心は、カウント データの負の二項式のハードルに特にあります。psclhurdle

残念ながら、のflexmix組み込みドライバーを提供していないパッケージに出くわしました。何度も試みたにもかかわらず、ドライバーを正しくコーディングすることができませんでした。psclhurdle

でを使用したことがある人はいますか?ドライバーの動作する実装を共有できますかpscl?hurdleflexmix

前もって感謝します。