問題タブ [kernel-density]
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.
r - Rでの密度推定
以下を使用して、古い信者の密度を推定しました。
忠実ではない点の密度を見つけるにはどうすればよいですか? ありがとう!
python - Python での多変量カーネル密度推定
SciPy のgaussian_kde
関数を使用して、多変量データの密度を推定しようとしています。以下のコードでは、3D 多変量法線をサンプリングしてカーネル密度に適合させていますが、適合を評価する方法がわかりません。
これを見ましたが、3D に拡張する方法がわかりません。
また、適合密度の評価をどのように開始すればよいかわかりませんか? これを視覚化するにはどうすればよいですか?
r - 正規性に対する密度分布の比率のテスト
正規分布と一様分布があります。均一の密度に対する正規分布の密度の比率を計算したいと思います。次に、この比率の正規性をテストします。
上記の分布はhm
、ht
私の実験データが示す例です。実際に使用するベクトルは、R でランダムに生成されたものではありません。
QQ プロットの相関係数から、正規性の良いアイデアを得ることができることを知っています。
正規分布する の場合hm
、これは 1 に近い値になります。
密度ベクトルの正規性を判断する方法はありますか? 例えばhmd
とratio
.
(追加情報:hm
およびht
は、長さのゲノム全体でホモ接合およびヘテロ接合の SNP をモデル化しています18585056
)
r - R でループを高速化する
次の関数を使用して、一部のデータのカーネル密度を推定しています
これがループの高速化に関する最初の質問ではないことはわかっています。サイトをチェックしたところ、一部の機能を使用した方が高速な場合があることがわかりましたがapply
、ループを正しく設定できた場合は常にそうとは限りません。
上記のコードでは、「不要なもの」はすべてループから除外されています。これは、私の理解が正しければ、計算を高速化することが示唆されているためです。KDens
しかし、上記の関数とdensity
Rにデフォルトで実装されている関数を比較してみました。まあ、私のマシンでは〜30秒必要ですdensity
が、1〜2秒必要です。KDens
編集:私が提供した情報は完全ではありませんでした
KDens 関数を高速化したい場合、どうすればよいでしょうか?