問題タブ [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.
datetime - Rで回転したxラベル
R には次のようなデータがあります。
class(event_ts) は「POSIXct」です。
異なる時間における各場所のバッグの密度を見つけたかったのです。そこで、関数密度を次のように使用しました。
最後のコメントは機能せず、関数「text」の「x」値をどのように認識すればよいかわかりません。
事前にコメントやご指導をよろしくお願いいたします。
さすがシマ。
r - ggplot2 の geom_density とベース R の密度の違い
R には次のようなデータがあります。
class(event_ts) はPOSIXctです。
異なる時間における各場所のバッグの密度を見つけたかったのです。
コマンドを使用したgeom_density(ggplot2)ところ、非常にうまくプロットできました。density(base)このコマンドとの間に違いがあるのだろうか。彼らが使用している方法や、使用しているデフォルトの帯域幅などに関する違いを意味します。
密度をデータ フレームに追加する必要があります。関数 を使用したことがあれば、関数を使用してこれらの値をデータ フレームに追加density(base)する方法を知っていましたが、 を使用する場合も同じでしょうか。approxfungeom_density(ggplot2)
python - Python 用の単純な固定帯域幅 1D Gaussian KDE 実装
固定帯域幅とガウス カーネルを使用した単純なカーネル密度推定が必要です。
残念ながら、scipy.stats.gaussian_kde私は帯域幅の自動選択しか見ませんでした。set_bandwidth私が見る限り、この方法は、自動選択された値にいくつかの修正比率を掛けるだけです。そして、設定したい帯域幅を設定するだけです。
手書きの Python コードを KDE に使用したくありません。動作が少し遅すぎます。
の代替品を知りませんscipy.stats.gaussian_kdeか?
r - カスタム密度で値をシミュレート
密度と値のシミュレーションに関係する理論的およびコーディング上の質問があります。
density(x)コマンドを使用してカスタム密度を構築しています。ただし、この密度から 1000 ～ 10000 のシミュレーション値を生成したいと考えています。全体的な目標は、密度 (x$y) 形式で構築された 2 つの密度を取得し、シミュレーションを実行して、この密度 A が密度 B x% の時間よりも大きいと言うことです。シミュレートされた各値を取得して、どちらが高いかを確認し、A が B よりも高い回数をカウントするコードを作成します。
これを達成する方法はありますか？または、これらの密度で同様のことを達成する方法はありますか? ありがとう！
python - Python で 2D の FFT ベースのカーネル密度推定器を実装し、それを SciPy の実装と比較する
2D カーネル密度推定 (KDE) を実行するコードが必要ですが、SciPy の実装が遅すぎることがわかりました。それで、私は FFT ベースの実装を書きましたが、いくつかのことが私を混乱させます。(FFT の実装では、周期的な境界条件も適用されます。これは私が望んでいるものです。)
実装は、サンプルから単純なヒストグラムを作成し、これをガウスで畳み込むことに基づいています。これを実行して SciPy の結果と比較するコードを次に示します。
上記の 2 つのサンプル セットがあります。1000 のランダム ポイントはベンチマーク用であり、コメント アウトされています。3 つのポイントはデバッグ用です。
後者の場合の結果のプロットは、この投稿の最後にあります。
ここに私の質問があります：
- ヒストグラムの .T と KDE1 の fftshift を避けることはできますか? なぜ必要なのかはわかりませんが、ガウス分布がないと間違った場所に表示されます。
- SciPy のスカラー帯域幅はどのように定義されていますか? 2 つの実装では、ガウス分布の幅が大きく異なります。
- 同様に、gaussian_kde にスカラー帯域幅を指定したにもかかわらず、SciPy 実装のガウス分布が放射状に対称でないのはなぜですか?
- FFT コード用に SciPy で利用可能な他の帯域幅メソッドを実装するにはどうすればよいですか?
(1000 ランダム ポイントの場合、FFT コードは SciPy コードよりも ~390 倍高速であることに注意してください。)
matlab - ksdensity は各データ ポイント xi をどのように計算しますか?
[f,xi] = ksdensity(x)Matlab で使用すると、確率密度推定値fと、 ksdensity がfを計算するxi評価点が得られます。
私の質問は次のとおりです。各xiポイントはどのように計算/決定されますか? 公式はありますか？
ドキュメンテーション センターには次のように記載されています。デフォルトは、 x のデータ範囲をカバーする 100 個の等間隔のポイントです。したがって、それらは範囲をカバーしていますが、これは がどのように計算されるかを説明していません。
どうもありがとうございました！
ファン